The Mossolov problem by the augmented Lagrangian method – solver class body
: Bi(0),
n(1), r(1), tol(1e-10), max_iter(1000000),
Xh(), Th(),
lh(),
a(),
b(), inv_mt(), pa()
{}
string grad_approx =
"P" +
itos(
Xh.degree()-1) +
"d";
Th =
space (omega, grad_approx,
"vector");
integrate_option iopt;
iopt.invert = true;
}
void
uh ["boundary"] = 0;
}
{
out << catchmark(
"Bi") <<
Bi << endl
<< catchmark(
"n") <<
n << endl
<< catchmark(
"r") <<
r << endl
<< catchmark("sigma") << sigma_h
<< catchmark("u") << uh;
}