diff --git a/NumLib/Fem/FiniteElement/TemplateIsoparametric.h b/NumLib/Fem/FiniteElement/TemplateIsoparametric.h
index cf93e3af5381c8a302da175401a25126c948a88d..d986d22fc1ddea07398ef1717833bd0d1eeaeb4d 100644
--- a/NumLib/Fem/FiniteElement/TemplateIsoparametric.h
+++ b/NumLib/Fem/FiniteElement/TemplateIsoparametric.h
@@ -56,25 +56,25 @@ public:
     /**
      * Constructor without specifying a mesh element. resetMeshElement() must be called afterwards.
      *
-     * @param integration_points_level      the sampling level (default 2)
+     * @param integration_order      default 2.
      */
-    explicit TemplateIsoparametric(std::size_t integration_points_level=2)
+    explicit TemplateIsoparametric(std::size_t integration_order=2)
     : _ele(nullptr)
     {
-        this->_integration.setSamplingLevel(integration_points_level);
+        this->_integration.setIntegrationOrder(integration_order);
     }
 
     /**
      * Construct this object for the given mesh element.
      *
-     * @param e                             Mesh element object
-     * @param integration_points_level      the sampling level (default 2)
+     * @param e                      Mesh element object
+     * @param integration_order      default 2.
      */
-    TemplateIsoparametric(const MeshElementType &e, std::size_t integration_points_level=2)
+    TemplateIsoparametric(const MeshElementType &e, std::size_t integration_order=2)
     : _ele(&e)
     {
         this->resetMeshElement(e);
-        this->_integration.setSamplingLevel(integration_points_level);
+        this->_integration.setIntegrationOrder(integration_order);
     }
 
     ///
diff --git a/Tests/NumLib/TestFeQuad4.cpp b/Tests/NumLib/TestFeQuad4.cpp
index 71df5f4da780012531207a399cb0ec4698f35c8d..b906e6e8d3fc9c3258b6f51cba21e76ae7939d75 100644
--- a/Tests/NumLib/TestFeQuad4.cpp
+++ b/Tests/NumLib/TestFeQuad4.cpp
@@ -279,7 +279,7 @@ TYPED_TEST(NumLibFemIsoQuad4Test, CheckGaussIntegrationLevel)
     ASSERT_ARRAY_NEAR(this->expectedM.data(), M.data(), M.size(), this->eps);
 
     // Change gauss quadrature level to 3
-    q.setSamplingLevel(3);
+    q.setIntegrationOrder(3);
     M *= .0;
     for (std::size_t i=0; i<q.getNPoints(); i++) {
         shape.setZero();