119 |
void PenalizedFictitousDomain::Discretize (ConstReferenceCounting<Problem> givenProblem) |
void PenalizedFictitousDomain::Discretize (ConstReferenceCounting<Problem> givenProblem) |
120 |
{ |
{ |
121 |
ffout(2) << "Adding Characteristic function of the domain to PDEs ... " << std::flush; |
ffout(2) << "Adding Characteristic function of the domain to PDEs ... " << std::flush; |
122 |
ConstReferenceCounting<UserFunction> u = new ConstUserFunction(1.); |
ConstReferenceCounting<UserFunction> u = new OneDomainUserFunction((*givenProblem).domain()); |
123 |
|
__problem = (*givenProblem) * u; |
|
// __problem = (*givenProblem) * u; |
|
|
__problem = givenProblem; |
|
124 |
ffout(2) << "done\n"; |
ffout(2) << "done\n"; |
125 |
|
|
126 |
switch (mesh().type()) { |
switch (mesh().type()) { |