div(3) divergence operator


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
               / 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", "vector");
        space M(omega, "P0");
        form b(V, M, "div");