src/GeomFill/GeomFill_Tensor.cxx

   1 // Created on: 1996-12-05
   2 // Created by: Philippe MANGIN
   3 // Copyright (c) 1996-1999 Matra Datavision
   4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
   5 //
   6 // This file is part of Open CASCADE Technology software library.
   7 //
   8 // This library is free software; you can redistribute it and/or modify it under
   9 // the terms of the GNU Lesser General Public License version 2.1 as published
  10 // by the Free Software Foundation, with special exception defined in the file
  11 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
  12 // distribution for complete text of the license and disclaimer of any warranty.
  13 //
  14 // Alternatively, this file may be used under the terms of Open CASCADE
  15 // commercial license or contractual agreement.
  16
  17
  18 #include <GeomFill_Tensor.hxx>
  19 #include <math_Matrix.hxx>
  20
  21 GeomFill_Tensor::GeomFill_Tensor(const Standard_Integer NbRow,
  22                                  const Standard_Integer NbCol,
  23                                  const Standard_Integer NbMat) :
  24                                    Tab(1,NbRow*NbMat*NbCol),
  25                                    nbrow( NbRow ),
  26                                    nbcol( NbCol ),
  27                                    nbmat( NbMat ),
  28                                    nbmtcl( NbMat*NbCol )
  29 {
  30 }
  31
  32 void GeomFill_Tensor::Init(const Standard_Real InitialValue)
  33 {
  34 // Standard_Integer I, T = nbrow * nbcol *  nbmat;
  35 // for (I=1; I<=T; I++) {Tab(I) = InitialValue;}
  36   Tab.Init(InitialValue);
  37 }
  38
  39 void GeomFill_Tensor::Multiply(const math_Vector& Right,
  40                                 math_Matrix& M) const
  41 {
  42   Standard_Integer i, j, k;
  43   Standard_Real Somme;
  44   for ( i=1; i<=nbrow; i++)
  45     {
  46       for ( j=1; j<=nbcol; j++)
  47         {
  48           Somme = 0;
  49           for ( k=1; k<=nbmat; k++)
  50             {
  51               Somme += Value(i,j,k)*Right(k);
  52             }
  53           M(i,j) = Somme;
  54         }
  55     }
  56 }
  57
  58