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) { |
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, |