grad(3) gradient operator


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
                / Omega

The 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");