diff --git a/NumLib/Fem/CoordinatesMapping/NaturalCoordinatesMapping.cpp b/NumLib/Fem/CoordinatesMapping/NaturalCoordinatesMapping.cpp
index 3b56c147fc947a6eb1cb39988e9ea49832e4933c..2a7adfb144f38af4c3a3e0d548094c61523217a2 100644
--- a/NumLib/Fem/CoordinatesMapping/NaturalCoordinatesMapping.cpp
+++ b/NumLib/Fem/CoordinatesMapping/NaturalCoordinatesMapping.cpp
@@ -143,8 +143,8 @@ inline void computeMappingMatrices(
             shapemat,
             FieldType<ShapeMatrixType::DNDR>());
 
-        auto const dim = ele.getDimension();
-        auto const nnodes = ele.getNumberOfNodes();  // TODO Taylor-Hood
+        auto const dim = T_SHAPE_FUNC::DIM;
+        auto const nnodes = T_SHAPE_FUNC::NPOINTS;
 
         // jacobian: J=[dx/dr dy/dr // dx/ds dy/ds]
         for (auto k = decltype(nnodes){0}; k < nnodes; k++)