Higher order to proper capture

KratosMultiphysics · Dec 9, 2022 · 46fdc94 · 46fdc94
1 parent 5690b3f
commit 46fdc94
Showing 1 changed file with 5 additions and 4 deletions.
diff --git a/kratos/geometries/line_3d_3.h b/kratos/geometries/line_3d_3.h
@@ -340,11 +340,12 @@ class Line3D3 : public Geometry<TPointType>
     double Length() const override
     {
         Vector temp;
-        this->DeterminantOfJacobian( temp, msGeometryData.DefaultIntegrationMethod() );
-        const IntegrationPointsArrayType& integration_points = this->IntegrationPoints( msGeometryData.DefaultIntegrationMethod() );
+        const auto& r_integration_order_required_analytical_solution = GeometryData::IntegrationMethod::GI_GAUSS_3;
+        this->DeterminantOfJacobian( temp, r_integration_order_required_analytical_solution );
+        const IntegrationPointsArrayType& r_integration_points = this->IntegrationPoints( r_integration_order_required_analytical_solution );
         double length = 0.0;
-        for ( unsigned int i = 0; i < integration_points.size(); i++ ) {
-            length += temp[i] * integration_points[i].Weight();
+        for ( unsigned int i = 0; i < r_integration_points.size(); i++ ) {
+            length += temp[i] * r_integration_points[i].Weight();
         }
 
         return length;