/[ff3d]/ff3d/solver/PenalizedFictitousDomain.cpp
ViewVC logotype

Diff of /ff3d/solver/PenalizedFictitousDomain.cpp

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

revision 1.2 by delpinux, Sat Apr 12 16:12:55 2003 UTC revision 1.3 by delpinux, Sat Sep 20 21:31:09 2003 UTC
# Line 55  void PenalizedFictitousDomain::__discret Line 55  void PenalizedFictitousDomain::__discret
55  {  {
56    MemoryManager MM;    MemoryManager MM;
57    
58      
59    
60    bool performAssembling =MM.ReserveMatrix(__A,    bool performAssembling =MM.ReserveMatrix(__A,
61                                             problem().numberOfUnknown(),                                             problem().numberOfUnknown(),
62                                             __degreeOfFreedomSet.size());                                             __degreeOfFreedomSet.size());
# Line 114  void PenalizedFictitousDomain::__discret Line 116  void PenalizedFictitousDomain::__discret
116    }    }
117  }  }
118    
119  void PenalizedFictitousDomain::Discretize (ConstReferenceCounting<Problem> Pb)  void PenalizedFictitousDomain::Discretize (ConstReferenceCounting<Problem> givenProblem)
120  {  {
121    __problem = Pb;    ffout(2) << "Adding Characteristic function of the domain to PDEs ... " << std::flush;
122      ConstReferenceCounting<UserFunction> u = new ConstUserFunction(1.);
123        
124    //   __problem = (*givenProblem) * u;
125      __problem = givenProblem;
126      ffout(2) << "done\n";
127    
128    switch (mesh().type()) {    switch (mesh().type()) {
129    case Mesh::cartesianHexahedraMesh: {    case Mesh::cartesianHexahedraMesh: {
130      this->__discretize<Structured3DMesh>();      this->__discretize<Structured3DMesh>();

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

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