On Jan 17, 2014, at 10:10 AM, subramanya sadasiva
wrote:
> Thanks.I wanted to try and keep it simple to alternate between lumped and
> full mass matrices. I had a doubt about using them with adaptivity. will
> solving , followed by applying constraints exactly give me the correct
> results?
Not through the dof map, but you can manually assemble into a 'diagonal'
matrix you create - e.g. a vector.
-Ben
> On Jan 17, 2014, at 8:54 AM, "subramanya sadasiva"
> wrote:
>
> I am trying to compute the L2 projection of some gauss point quantities to
> the nodes, and I'd like to use a t
I am trying to compute the L2 projection of some gauss point quantities to the
nodes, and I'd like to use a trapezoid rule to integrate the mass matrix. Is
there anyway to optimize the allocation of the matrix, so that it is a strictly
diagonal matrix? Thanks, Subramanya
