form(const space V, const space& M, "div");
Assembly the form associated to the divergence operator on a finite element space V:
/ | b(u,q) = | div(u) q dx | / OmegaThe V space may be a either P1 or P2 finite element space, while the M space may be P0 or P1d respectively. See also form(2) and space(2).
The following piece of code build the divergence form associated to the P1 approximation:
geo omega("square"); space V(omega, "P1", "vector"); space M(omega, "P0"); form b(V, M, "div");