Fix shear heating bug (#86)

Fixes bug in the shear heating kernel for single phase materials.
PTsolvers · Jan 5, 2024 · 1de2830 · 1de2830
1 parent 973d613
commit 1de2830
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/thermal_diffusion/Shearheating.jl b/src/thermal_diffusion/Shearheating.jl
@@ -1,7 +1,7 @@
 @parallel_indices (I...) function compute_shear_heating!(
     shear_heating, τ::NTuple{N,T}, τ_old::NTuple{N,T}, ε::NTuple{N,T}, rheology, dt
 ) where {N,T}
-    _Gdt = inv(fn_ratio(get_G, rheology) * dt)
+    _Gdt = inv(get_G(rheology) * dt)
     τij, τij_o, εij = cache_tensors(τ, τ_old, ε, I...)
     εij_el = @. 0.5 * ((τij - τij_o) * _Gdt)
     shear_heating[I...] = compute_shearheating(rheology, τij, εij, εij_el)