00001 00005 class MassMatrix : public HandleId 00006 00007 { 00008 friend class FEM; 00009 Handle(GridFE) mesh; 00010 int ndof_per_node; 00011 bool lumped; 00012 Handle(MatDiag(NUMT)) mass_matrix; 00013 00014 public: 00015 MassMatrix (); 00016 ~MassMatrix () {} 00017 00018 bool isLumped () const { return lumped; } 00019 const MatDiag(NUMT)& M () const; 00020 MatDiag(NUMT)& M (); 00021 bool make (FEM& simulator, GridFE& grid); 00022 void solve (const Vector(NUMT)& rhs, Vector(NUMT)& x); 00023 }; 00024 00025