minor stuff

Author: NickNick9

iga/use_cases/external_boundary_circle_with_nurbs/plot_surrogate_boundaries.py

@@ -99,8 +99,6 @@
                     # plot gp, arrows, etc...
 
 
-
-    
     # === ACTIVE GAUSS POINTS: FROM ELEMENT CENTERS ===
     mp = model["IgaModelPart.StructuralAnalysisDomain"]
     active_gauss_points = []

iga/use_cases/external_boundary_circle_with_nurbs/run_and_post.py

@@ -18,11 +18,7 @@ def InitializeSolutionStep(self):
 
 # Define the analytical solution
 def analytical_temperature(x, y, t):
-    # return x**2*y
     return -np.cos(x)*np.sinh(y)
-    # return x**3+y**3
-    # return x**2+y**2
-    # return x+y
 
 # Main simulation function
 if __name__ == "__main__":
@@ -61,25 +57,6 @@ def analytical_temperature(x, y, t):
             curr_temperature += N[0, i] * node.GetSolutionStepValue(KratosMultiphysics.DISPLACEMENT_X, 0)
         computed_temperature.append(curr_temperature)
 
-
-    # # Loop over elements to gather computed solution
-    # for cond in mp.Conditions:
-    #     geom = cond.GetGeometry()
-    #     N = geom.ShapeFunctionsValues()
-    #     center = geom.Center()
-    #     weight = elem.GetValue(KratosMultiphysics.INTEGRATION_WEIGHT)
-    #     weights.append(weight)
-    #     # Extract Gauss point (center) coordinates
-    #     x_coord.append(center.X)
-    #     y_coord.append(center.Y)
-
-    #     # Initialize solution values at the center
-    #     curr_temperature = 0
-    #     # Compute nodal contributions using shape functions
-    #     for i, node in enumerate(geom):
-    #         curr_temperature += N[0, i] * node.GetSolutionStepValue(KratosMultiphysics.TEMPERATURE, 0)
-    #     computed_temperature.append(curr_temperature)
-
     # Get the current time after simulation run
     current_time = simulation._GetSolver().GetComputingModelPart().ProcessInfo[KratosMultiphysics.TIME]
     print("Current time after simulation:", current_time)