.\"
.de Id
..
.de Sp
.if n .sp
.if t .sp 0.4
..
.TH div 3rheolef "rheolef-6.1" "rheolef-6.1" "rheolef-6.1"
.\" label: /*Class:div
.SH NAME
\fBdiv\fP -- divergence operator
.\" skip: @bfindex div
.\" skip: @cindex divergence
.\" skip: @apindex P1
.\" skip: @apindex P2
.\" skip: @apindex P1d
.\" skip: @apindex P0
.SH SYNOPSIS

.\" begin_example
.Sp
.nf
      form(const space V, const space& M, "div");
.Sp
.fi
.\" end_example
.SH DESCRIPTION

Assembly the form associated to the
divergence operator on a finite element space \fBV\fP:
.\" begin_example
.Sp
.nf
               /
               |
      b(u,q) = |  div(u) q dx
               |
               / Omega
.Sp
.fi
.\" end_example
.\" END IFINFO
The V space may be a either
\fBP1\fP or \fBP2\fP finite element space,
while the M space may be \fBP0\fP or \fBP1d\fP respectively.
See also form(2) and space(2).
.PP
.SH EXAMPLE

The following piece of code build the divergence form
associated to the \fBP1\fP approximation:
.\" begin_example
.Sp
.nf
        geo omega("square");
        space V(omega, "P1", "vector");
        space M(omega, "P0");
        form b(V, M, "div");
.Sp
.fi
.\" end_example
.\" skip start:AUTHOR:

.\" skip start:DATE:

.\" END

.\" LENGTH = 2
.SH SEE ALSO
form(2), space(2)