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

Diff of /ff3d/solver/FictitiousDomainMethod.cpp

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

revision 1.1 by delpinux, Mon Nov 3 18:56:05 2003 UTC revision 1.2 by delpinux, Tue Nov 4 19:24:35 2003 UTC
# Line 105  __computesDegreesOfFreedom(ConstReferenc Line 105  __computesDegreesOfFreedom(ConstReferenc
105    
106    for (size_t i=0; i<m.numberOfCells();++i) {    for (size_t i=0; i<m.numberOfCells();++i) {
107      typename MeshType::ElementsGeometry& h = m.cell(i);      typename MeshType::ElementsGeometry& h = m.cell(i);
108        h.isFictitious() = false; // reset fictitious state!
109    
110      unsigned numberIn = 0;      unsigned numberIn = 0;
111    
112      for (unsigned j=0; j<MeshType::ElementsGeometry::NumberOfVertices; ++j) {      for (unsigned j=0; j<MeshType::ElementsGeometry::NumberOfVertices; ++j) {
# Line 127  __computesDegreesOfFreedom(ConstReferenc Line 129  __computesDegreesOfFreedom(ConstReferenc
129    ConstReferenceCounting<UserFunction> u = new FEM0UserFunction(__mesh, kiValues);    ConstReferenceCounting<UserFunction> u = new FEM0UserFunction(__mesh, kiValues);
130    __problem = (*givenProblem) * u;    __problem = (*givenProblem) * u;
131    ffout(2) << "done\n";    ffout(2) << "done\n";
       
   for (size_t i=0; i<dofVertices.size(); ++i) {  
     dofVertices[i] = (correspondance[i] != -1);  
   }  
132    
133    __degreeOfFreedomSet = new DegreeOfFreedomFDMSet(__givenDegreeOfFreedomSet.numberOfVariables(),    __degreeOfFreedomSet = new DegreeOfFreedomFDMSet(__givenDegreeOfFreedomSet.numberOfVariables(),
134                                                     numberOfUsedVertices,                                                     numberOfUsedVertices,

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

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