Fill host_field with the right-hand side of the Poisson equation.
The routine uses the host-associated mesh, dims and n_pi.
Edit this routine to supply another right-hand side, then rebuild.
The coordinates are the cell centres from
mesh%get_coordinates(i, j, k, CELL).
The right-hand side must have zero mean.
In every supported combination of boundary conditions, the pressure
has periodic or homogeneous Neumann conditions.
The default is cos(2 pi x).
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(field_t), | intent(inout) | :: | host_field |