function A=Masse(Eh,varargin)
  assert(fc_simesh.is_siMeshElt(Eh))
  p = inputParser;
  p.addParameter('local', true, @islogical );
  p.parse(varargin{:});
  AL=FEM.Masse(Eh.nq,Eh.me,Eh.vols);
  if p.Results.local,A=AL;return;end
  A=sparse(Eh.nqGlobal,Eh.nqGlobal);
  I=Eh.toGlobal;
  A(I,I)=AL;
end
