function V=restriction(Eh,U)
  assert(fc_simesh.is_siMeshElt(Eh))
  assert(all(size(U)==[Eh.nqGlobal,1]))
  V=U(Eh.toGlobal);
end
