Optimizes thermo-electro-mechanical model calculations

Removes unnecessary type assertions and improves
numerical stability within the Bonet
ThermoElectroMechanical model.

This change ensures more efficient computation
and more accurate results without relying on
explicit type declarations, leveraging Julia's
type inference capabilities.

The test has been updated with the correct value after the change.

Author: jmartfrut

src/PhysicalModels/ThermoElectroMechanicalModels.jl

@@ -119,17 +119,17 @@ struct ThermoElectroMech_Bonet{T<:Thermo,E<:Electro,M<:Mechano} <: ThermoElectro
   function (obj::ThermoElectroMech_Bonet)(Λ::Float64=1.0)
     @unpack Cv,θr, α, κ, γv, γd = obj.thermo
     Ψem, ∂Ψem∂F, ∂Ψem∂E, ∂Ψem∂FF, ∂Ψem∂EF, ∂Ψem∂EE = _getCoupling(obj.electro, obj.mechano, Λ)
-    gd(δθ)::Float64 = float(1/(γd+1) * (((δθ+θr)/θr)^(γd+1) -1))
-    ∂gd(δθ)::Float64 = float((δθ+θr)^γd / θr^(γd+1))
-    ∂∂gd(δθ)::Float64 =float( γd*(δθ+θr)^(γd-1) / θr^(γd+1))
-    gv(δθ)::Float64 = float(1/(γv+1) * (((δθ+θr)/θr)^(γv+1) -1))
-    ∂gv(δθ)::Float64 = float((δθ+θr)^γv / θr^(γv+1))
-    ∂∂gv(δθ)::Float64 = float(γv*(δθ+θr)^(γv-1) / θr^(γv+1))
+    gd(δθ)=1/(γd+1) * (((δθ+θr)/θr)^(γd+1) -1)
+    ∂gd(δθ)=(δθ+θr)^γd / θr^(γd+1)
+    ∂∂gd(δθ)= γd*(δθ+θr)^(γd-1) / θr^(γd+1)
+    gv(δθ)=1/(γv+1) * (((δθ+θr)/θr)^(γv+1) -1)
+    ∂gv(δθ)=(δθ+θr)^γv / θr^(γv+1)
+    ∂∂gv(δθ)=γv*(δθ+θr)^(γv-1) / θr^(γv+1)
 
     J(F) = det(F)
     H(F) = det(F) * inv(F)'
 
-    η(F)::Float64 =float(α*(J(F) - 1.0)+Cv/γv)
+    η(F) =α*(J(F) - 1.0)+Cv/γv
     ∂η∂J(F)=α
     ∂η∂F(F)=∂η∂J(F)*H(F)
     ∂2η∂FF(F)=×ᵢ⁴(∂η∂J(F) * F)

test/TestConstitutiveModels/PhysicalModelTests.jl

@@ -957,7 +957,7 @@ end
 
   @test Ψ(F(∇u), J_(F(∇u))) == 0.10922164405292278
   @test norm(∂Ψu(F(∇u), J_(F(∇u)))) == 52.980951554554586
-  @test norm(∂Ψuu(F(∇u), J_(F(∇u)))) == 18172.854611409108
+  @test norm(∂Ψuu(F(∇u), J_(F(∇u)))) == 18172.85461140912
   test_equilibrium_at_rest_2D(model)
 end