fix

Author: miguelmaso

src/WeakForms/ElectroMechanics.jl

@@ -27,15 +27,15 @@ end
 # Stagered jacobian
 # -----------------
 
-function jacobian(physicalmodel::ElectroMechano, ::Type{Mechano}, kine::NTuple{2,KinematicModel}, (u, φ), du, v, dΩ, Λ, vars...; kwargs...)
+function jacobian(physicalmodel::ElectroMechano, ::Type{Mechano}, kine::NTuple{2,KinematicModel}, (u, φ), du, v, dΩ, Λ=1.0, vars...; kwargs...)
     DΨ    = physicalmodel(; kwargs...)
     F,_,_ = get_Kinematics(kine[1]; Λ=Λ)
     E     = get_Kinematics(kine[2]; Λ=Λ)
     ∂Ψuu  = DΨ[4]
     ∫(∇(v)' ⊙ ((∂Ψuu ∘ (F∘∇(u)', E∘∇(φ), vars...)) ⊙ ∇(du)'))dΩ
 end
 
-function jacobian(physicalmodel::ElectroMechano, ::Type{Electro}, kine::NTuple{2,KinematicModel}, (u, φ), dφ, vφ, dΩ, Λ, vars...; kwargs...)
+function jacobian(physicalmodel::ElectroMechano, ::Type{Electro}, kine::NTuple{2,KinematicModel}, (u, φ), dφ, vφ, dΩ, Λ=1.0, vars...; kwargs...)
     DΨ    = physicalmodel(; kwargs...)
     F,_,_ = get_Kinematics(kine[1]; Λ=Λ)
     E     = get_Kinematics(kine[2]; Λ=Λ)

test/TestConstitutiveModels/ElectroMechanicalTests.jl

@@ -115,8 +115,8 @@ end
   λvn = 1e-3
   Avn = VectorValue(Uvn..., λvn)
   Ψ, ∂Ψu, ∂Ψφ, ∂Ψuu, ∂Ψφu, ∂Ψφφ = model(Δt=0.01)
-  @test norm(∂Ψu(F(∇u), F(∇un), E(∇φ), Avn)) ≈ 25.049301121178615
-  @test norm(∂Ψuu(F(∇u), F(∇un), E(∇φ), Avn)) ≈ 3110.7607787445168
+  @test norm(∂Ψu(F(∇u), E(∇φ), F(∇un), Avn)) ≈ 25.049301121178615
+  @test norm(∂Ψuu(F(∇u), E(∇φ), F(∇un), Avn)) ≈ 3110.7607787445168
 end
 
 
@@ -137,6 +137,6 @@ end
   λvn = 1e-3
   Avn = VectorValue(Uvn..., λvn)
   Ψ, ∂Ψu, ∂Ψφ, ∂Ψuu, ∂Ψφu, ∂Ψφφ = model(Δt=0.01)
-  @test norm(∂Ψu(F(∇u), F(∇un), E(∇φ), Avn, Avn)) ≈ 25.102080194257017
-  @test norm(∂Ψuu(F(∇u), F(∇un), E(∇φ), Avn, Avn)) ≈ 3110.9722775475557
+  @test norm(∂Ψu(F(∇u), E(∇φ), F(∇un), Avn, Avn)) ≈ 25.102080194257017
+  @test norm(∂Ψuu(F(∇u), E(∇φ), F(∇un), Avn, Avn)) ≈ 3110.9722775475557
 end