On Jan 17, 2014, at 10:10 AM, subramanya sadasiva <[email protected]>
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?
> Subramanya
I'll see if I can find a reference, but I recall reading that since the mass
matrix is so nice it can be iteratively solved via Jacobi or Gauss-Siedel
iteration as
Mu = f
(D + OD)u = f
Du = f - ODu
Du^{l+1} = f - ODu^l
where D is diagonal and OD is the rest…
Now, I'd love a patch that basically implements that if you are feeling
ambitious…
For the lumped case of course OD=0 and everything converges in one iteration.
-Ben
------------------------------------------------------------------------------
CenturyLink Cloud: The Leader in Enterprise Cloud Services.
Learn Why More Businesses Are Choosing CenturyLink Cloud For
Critical Workloads, Development Environments & Everything In Between.
Get a Quote or Start a Free Trial Today.
http://pubads.g.doubleclick.net/gampad/clk?id=119420431&iu=/4140/ostg.clktrk
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
