test fix

Author: CWillberg

src/Compute/compute_field_values.jl

@@ -125,8 +125,7 @@ function calculate_stresses(block_nodes::Dict{Int64,Vector{Int64}},
                                            material_parameter["Symmetry"],
                                            Data_Manager.get_dof())
             calculate_strain(active_nodes,
-                             invert(hookeMatrix,
-                                    "Hooke matrix not invertable"))
+                             inv(hookeMatrix))
         end
     end
 end

src/Models/Material/Material_Basis.jl

@@ -403,7 +403,7 @@ function get_2D_Hooke_matrix(aniso_matrix::MMatrix{T}, symmetry::String,
         matrix[3, 3] = aniso_matrix[6, 6]
         return matrix
     elseif occursin("plane stress", symmetry)
-        inv_aniso = invert(aniso_matrix, "Hooke matrix not invertable")
+        invert(inv_aniso, aniso_matrix, "Hooke matrix not invertable")
         matrix = get_MMatrix(36)
         matrix[1:2, 1:2] = inv_aniso[1:2, 1:2]
         matrix[3, 1:2] = inv_aniso[6, 1:2]

src/Support/Geometry.jl

@@ -296,7 +296,7 @@ Therefore,
 
 """
 
-function compute_left_stretch_tensor(deformation_gradient::Matrix{Float64})
+function compute_left_stretch_tensor(deformation_gradient::AbstractMatrix{Float64})
     return sqrt(deformation_gradient * deformation_gradient')
 end
 
@@ -328,12 +328,13 @@ function deformation_gradient_decomposition(nodes::Union{Base.OneTo{Int64},Vecto
                                             rot_tensor)
     for iID in nodes
         # EQ (15) in https://doi.org/10.1016/j.ijsolstr.2008.10.029
-        rot_tensor[iID, :,
-        :] = invert(compute_left_stretch_tensor(deformation_gradient[iID,
-                                                                                        :,
-                                                                                        :]),
-                                       "Inversion of left stretch tensor fails in function deformation_gradient_decomposition.") *
-                                deformation_gradient[iID, :, :]
+        @views rot_tensor[iID, :, :] = invert(rot_tensor[iID, :,
+                                                         :],
+                                              compute_left_stretch_tensor(deformation_gradient[iID,
+                                                                                               :,
+                                                                                               :]),
+                                              "Inversion of left stretch tensor fails in function deformation_gradient_decomposition.") *
+                                       deformation_gradient[iID, :, :]
     end
     return rot_tensor
 end

src/Support/Helpers.jl

@@ -943,6 +943,15 @@ function interpol_data(x::Union{Vector{Float64},Vector{Int64},Float64,Int64},
     return evaluate(values["spl"], x)
 end
 
+function invert(A::AbstractMatrix{Float64},
+                error_message::String = "Matrix is singular")
+    if abs(det(smat(A))) < 1e-15
+        @error "$error_message "
+        return
+    end
+    return inv(smat(A))
+end
+
 """
 	invert(A::Union{Matrix{Float64},Matrix{Int64}}, error_message::String="Matrix is singular")
 

test/unit_tests/Support/ut_helpers.jl

@@ -50,13 +50,16 @@ end
 end
 
 @testset "ut_invert" begin
+    test_inv = zeros(2, 2)
     PeriLab.Solver_Manager.Helpers.invert(zeros(Float64, 4, 4),
                                           "test Float is singular.")
 
-    @test PeriLab.Solver_Manager.Helpers.invert([1.0 0; 0 1]) == inv([1 0; 0 1])
-    @test PeriLab.Solver_Manager.Helpers.invert([1.0 0.0 0.0; 0.0 1.0 0.0; 0.0 0.0 1.0],
-                                                "test Int is singular.") ==
-          inv([1.0 0.0 0.0; 0.0 1.0 0.0; 0.0 0.0 1.0])
+    PeriLab.Solver_Manager.Helpers.invert(test_inv, [1.0 0; 0 1])
+    @test test_inv == inv([1 0; 0 1])
+    test_inv = zeros(3, 3)
+    PeriLab.Solver_Manager.Helpers.invert(test_inv, [1.0 0.0 0.0; 0.0 1.0 0.0; 0.0 0.0 1.0],
+                                          "test Int is singular.")
+    @test test_inv == inv([1.0 0.0 0.0; 0.0 1.0 0.0; 0.0 0.0 1.0])
 end
 
 @testset "ut_find_local_neighbors" begin
@@ -295,7 +298,7 @@ end
     @test PeriLab.Solver_Manager.Helpers.check_inf_or_nan(a, "a") == false
     a = NaN
     @test PeriLab.Solver_Manager.Helpers.check_inf_or_nan(a, "a")
-    a = 1/0
+    a = 1 / 0
     @test PeriLab.Solver_Manager.Helpers.check_inf_or_nan(a, "a") == true
 end
 @testset "get_matrix_style" begin