we should create a geomap context comptuting the point coordinates (POINT context) and use it in an element/face loops. The issues with dofpoints() are : - one call will build all dofpoints - when the mesh move, the developper should call the rebuild of all dofpoint else we get the node without the moving see #1879 and replace the use of dofPoints in `opselt()`