apply also updateEntitiesCoDimensionOnePermutation when nDim!=nRealDim
fix geotool cylinder (else connectivity problem between surface)
add doftable check in debug mode
add PointSet parameter with Pch,Pchv, Pdh (default is PointSetEquiSpaced ) :
    auto Xh = Pch<3,PointSetFekete>(mesh);
start update of expr eval based on loc+interp.  goal : use submeshrelation or conformal meshes
more tests+clean