/[ff3d]/ff3d/solver/VectorialPDEOperator.hpp
ViewVC logotype

Diff of /ff3d/solver/VectorialPDEOperator.hpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1.1.1 by delpinux, Mon Feb 17 16:32:51 2003 UTC revision 1.2 by delpinux, Sat Sep 20 21:31:09 2003 UTC
# Line 21  Line 21 
21  // It means a PDE operator that can be applied to a vectorial function.  // It means a PDE operator that can be applied to a vectorial function.
22    
23    
24  #ifndef _VECTORIALPDEOPERATOR_HPP_  #ifndef VECTORIALPDEOPERATOR_HPP
25  #define _VECTORIALPDEOPERATOR_HPP_  #define VECTORIALPDEOPERATOR_HPP
26    
27  #include <vector>  #include <vector>
28  #include <PDEOperatorSum.hpp>  #include <PDEOperatorSum.hpp>
# Line 74  public: Line 74  public:
74      return os;      return os;
75    }    }
76    
77      ReferenceCounting<VectorialPDEOperator> operator * (ConstReferenceCounting<UserFunction> u) const
78      {
79        VectorialPDEOperator* newVectorialPDEOperator = new VectorialPDEOperator(this->size());
80    
81        for (size_t i=0; i<__vectorPDEOperator.size(); ++i)
82          (*newVectorialPDEOperator).__vectorPDEOperator[i] = (*__vectorPDEOperator[i]) * u;
83    
84        return newVectorialPDEOperator;
85      }
86    
87    const VectorialPDEOperator& operator += (const VectorialPDEOperator& VPDEOp)    const VectorialPDEOperator& operator += (const VectorialPDEOperator& VPDEOp)
88    {    {
89      assert(__vectorPDEOperator.size() == VPDEOp.__vectorPDEOperator.size());      assert(__vectorPDEOperator.size() == VPDEOp.__vectorPDEOperator.size());
# Line 85  public: Line 95  public:
95      return *this;      return *this;
96    }    }
97    
98    const VectorialPDEOperator& operator -= (VectorialPDEOperator& VPDEOp)    const VectorialPDEOperator& operator -= (const VectorialPDEOperator& VPDEOp)
99    {    {
100      assert(__vectorPDEOperator.size() == VPDEOp.__vectorPDEOperator.size());      assert(__vectorPDEOperator.size() == VPDEOp.__vectorPDEOperator.size());
101    
# Line 96  public: Line 106  public:
106      return *this;      return *this;
107    }    }
108    
109    VectorialPDEOperator(const size_t dimension)    VectorialPDEOperator(const size_t& dimension)
110        : __vectorPDEOperator(dimension)
111    {    {
     __vectorPDEOperator.resize(dimension);  
112      for(size_t i=0; i<dimension; ++i)      for(size_t i=0; i<dimension; ++i)
113        __vectorPDEOperator[i] = new PDEOperatorSum();        __vectorPDEOperator[i] = new PDEOperatorSum();
114    }    }
# Line 114  public: Line 124  public:
124    }    }
125  };  };
126    
127  #endif // _VECTORIALPDEOPERATOR_HPP_  #endif // VECTORIALPDEOPERATOR_HPP
   

Legend:
Removed from v.1.1.1.1  
changed lines
  Added in v.1.2

savannah-hackers-public@gnu.org
ViewVC Help
Powered by ViewVC 1.1.26