Thank you for the suggestions, Wolfgang. K, U, and A are correctly computed
(linear system is solved correctly) and they are verified by an analytical
example, it is S that I am not getting right when hanging nodes are
involved. I will do some hand computations by hand to see where I can get
f
On 1/8/23 08:40, Basca Jadamba wrote:
I am trying to assemble a rectangular matrix S with a property that K*U=S*A
where A/U is the coefficient/solution pair in an elliptic PDE (think -div(a
grad u) = f) and where K is the system matrix. S is a FullMatrix (no sparsity
pattern) and it works fine
I am trying to assemble a rectangular matrix S with a property that K*U=S*A
where A/U is the coefficient/solution pair in an elliptic PDE (think -div(a
grad u) = f) and where K is the system matrix. S is a FullMatrix (no
sparsity pattern) and it works fine on uniform meshes (correctness is
c
