106 |
} |
} |
107 |
} |
} |
108 |
|
|
109 |
void PenalizedMatrixFD::Discretize (ConstReferenceCounting<Problem> Pb) |
void PenalizedMatrixFD::Discretize (ConstReferenceCounting<Problem> givenProblem) |
110 |
{ |
{ |
111 |
__problem = Pb; |
ffout(2) << "Adding Characteristic function of the domain to PDEs ... " << std::flush; |
112 |
|
ConstReferenceCounting<UserFunction> u = new OneDomainUserFunction((*givenProblem).domain()); |
113 |
|
__problem = (*givenProblem) * u; |
114 |
|
ffout(2) << "done\n"; |
115 |
|
|
116 |
switch (mesh().type()) { |
switch (mesh().type()) { |
117 |
case Mesh::cartesianHexahedraMesh: { |
case Mesh::cartesianHexahedraMesh: { |