fix type of temperature field

PetrKryslUCSD · Oct 10, 2023 · e6591a0 · e6591a0 · PetrKryslUCSD · Oct 10, 2023
1 parent 85cb7c6
commit e6591a0
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/FEMMBaseModule.jl b/src/FEMMBaseModule.jl
@@ -879,7 +879,7 @@ function elemfieldfromintegpoints(
     component::IT;
     context...,
 ) where {FEMM<:AbstractFEMM, GFT<:Number, UFT<:Number, IT<:Integer}
-    dT = NodalField(zeros(FT, nnodes(geom), 1)) # zero difference in temperature
+    dT = NodalField(zeros(GFT, nnodes(geom), 1)) # zero difference in temperature
     return elemfieldfromintegpoints(self, geom, u, dT, quantity, [component]; context...)
 end
 
@@ -891,7 +891,7 @@ function elemfieldfromintegpoints(
     component::AbstractVector{IT};
     context...,
 ) where {FEMM<:AbstractFEMM, GFT<:Number, UFT<:Number, IT<:Integer}
-    dT = NodalField(zeros(FT, nnodes(geom), 1)) # zero difference in temperature
+    dT = NodalField(zeros(GFT, nnodes(geom), 1)) # zero difference in temperature
     return elemfieldfromintegpoints(self, geom, u, dT, quantity, component; context...)
 end