function U=relevement(Eh,V,varargin)
  assert(fc_simesh.is_siMeshElt(Eh))
  assert(all(size(V)==[Eh.nq,1]))
  p = inputParser;
  p.addParameter('funinit', @zeros, @(x) fc_tools.utils.isfunhandle(x) ); % can use rand, ones, ...
  p.parse(varargin{:});
  U=p.Results.funinit(Eh.nqGlobal,1);
  U(Eh.toGlobal)=V;
end
