The Poisson problem by the hybrid discontinuous Galerkin method – local averaging function.
The Poisson problem by the hybrid discontinuous Galerkin method – local averaging function
size_t k = uh.get_space().degree();
size_t d = uh.get_geo().dimension();
space Zh (uh.get_geo(),
"P0");
integrate_option iopt;
iopt.invert = true;
if (k >= 1) {
form inv_mz = integrate (zeta*xi, iopt);
iopt.set_order(2*k+2);
field lh = integrate (uh*xi, iopt);
}
const space& Mh = lambda_h.get_space();
form inv_ms = integrate(
"sides",
lambda*mu, iopt);
field lh = integrate (on_local_sides(inv_sh*lambda_h*xi));
}
field lh(Float epsilon, Float t, const test &v)
see the field page for the full documentation
see the space page for the full documentation
see the test page for the full documentation
see the test page for the full documentation
field dirichlet_hdg_average(field uh, field lambda_h)