Hi Roy, Thanks. Yes, my goal is to have a non-diagonal mass matrix. I want to take the easiest path possible for doing that using LibMesh. By mass lumping with a Lagrange basis, do you mean having cell centered finite volume using a 0th order Lagrange basis in finite element?
Thanks, Harshad On Wed, Dec 16, 2015 at 5:06 PM, Roy Stogner <royst...@ices.utexas.edu> wrote: > > On Wed, 16 Dec 2015, Harshad Sahasrabudhe wrote: > > Sorry I meant to say that the control volume will go beyond the >> element if the element has obtuse angles. I want to solve an >> elliptic second order PDE with vertex centered finite volume method. >> > > Ah! Vertex centered makes your question make sense, but also changes > things. libMesh doesn't have any built in features for dealing with > dual meshes, so you'd have to handle that all on your own. > > Are you sure you're not creating an XY problem here? Your original > concern was the non-diagonal mass matrix that FEM gives you, right? > Wouldn't simply mass lumping with a Lagrange basis be an easy fix? > --- > Roy > ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
