Build the UW stabilized stiffness matrix for an advection problem.
The equation taken into account is:
(beta u)’ = f
where beta is an element-wise constant.
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: bim1a_rhs, bim1a_reaction, bim1a_laplacian, bim2a_advection_diffusion.
Package: bim