[Getdp] lumping method
Christophe Geuzaine
geuzaine at acm.caltech.edu
Thu Apr 18 17:20:44 CEST 2002
Lin Ji wrote:
>
> Dear Christophe,
> Have you heared of the lumping method that simply replace the mass matrix
> by a diagonal matrix when solving the wave equation? It makes solving the
Sure.
> linear system at each time step explicit and therefore much faster. It seems
Yes, but you
1) loose some accuracy
2) are subjected to the classical CFL stability condition on your time
step
3) have constraints on the quality of your mesh ("no flat elements") to
have a stable scheme
It would be nice to implement such a lumping in getdp, though.
Unfortunately I do not have the time to do it myself. Would you be
interested?
Christophe
--
Christophe Geuzaine
Tel: (626) 395-4552 http://www.geuz.org
Fax: (626) 578-0124 mailto:geuzaine at acm.caltech.edu