form(const space V, const space& M, "grad");
Assembly the form associated to the gradient operator on a finite element space V:
/ | b(u, q) = | grad(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"); space M(omega, "P0", "vector"); form b(V, M, "grad");