30 pA.solve(mryh, delta_yh);
36 field rhs =
b.trans_mult(rh);
38 pa.solve (rhs, delta_vh);
45 return sqrt (
mb(rh,rh));
50 return sqrt (dual (mrh,rh));
54 return (
Cf+
r)*uh[
"boundary"];
57 field rhs =
lh -
b.trans_mult(beta_h);
60 lambda_h = beta_h -
r*uh[
"boundary"];
field gh(Float epsilon, Float t, const field &uh, const test &v)
see the Float page for the full documentation
see the field page for the full documentation
field derivative_trans_mult(const field &mrh) const
Float dual_space_norm(const field &) const
Float space_norm(const field &) const
field derivative_solve(const field &mrh) const
void post(const field &beta_h, field &uh, field &lambda_h) const
The Poisson problem with Robin boundary condition – solver function.
field poisson_robin(Float Cf, const geo &boundary, const field &lh)