Build the UW stabilized stiffness matrix for an advection problem.
The equation taken into account is:
div (beta u) = f
where beta is an element-wise constant vector function.
Instead of passing the vector field beta directly one can pass a piecewise linear conforming scalar function phi as the last input. In such case beta = grad phi is assumed.
If phi is a single scalar value beta is assumed to be 0 in the whole domain.
See also: bim2a_rhs, bim2a_reaction, bim2c_mesh_properties.
Package: bim