More changes in use compile time matrix sizes in LIE

For these changes, a new alias type of

#ifdef OGS_EIGEN_DYNAMIC_SHAPE_MATRICES
...
template <typename ShapeFunction>
using ShapeMatrixPolicyTypeNew =
    EigenDynamicShapeMatrixPolicy<ShapeFunction, ShapeFunction::DIM>;
...
#else
template <typename ShapeFunction>
using ShapeMatrixPolicyTypeNew =
    EigenFixedShapeMatrixPolicy<ShapeFunction, ShapeFunction::DIM>;
...
#endif

is introduced, which is also used in MR. !3632. The plan is to replace ShapeMatrixPolicyType with ShapeMatrixPolicyTypeNew assembler by assembler, and at the end to remove the template argument of int GlobalDim from all template entities in ShapeMatrixPolicy.h. Once this is done, the type name can be changed back.

NumLib/Fem/CoordinatesMapping/NaturalCoordinatesMapping.cpp

@@ -214,11 +214,10 @@ inline void computeMappingMatrices(
         // J^-1, dshape/dx
         shapemat.invJ.noalias() = shapemat.J.inverse();
 
-        auto const nnodes(shapemat.dNdr.cols());
-        auto const ele_dim(shapemat.dNdr.rows());
         assert(shapemat.dNdr.rows() == ele.getDimension());
-        shapemat.dNdx.topLeftCorner(ele_dim, nnodes).noalias() =
-            shapemat.invJ * shapemat.dNdr;
+        shapemat.dNdx
+            .template topLeftCorner<T_SHAPE_FUNC::DIM, T_SHAPE_FUNC::NPOINTS>()
+            .noalias() = shapemat.invJ * shapemat.dNdr;
     }
 }