Added density to model definition
Added a generic mass matrix routine