Minor updates to diff to improve type support

Author: KnutAM

src/differentiation.jl

@@ -128,11 +128,12 @@ is defined already.
 """
 function StressStateDerivatives(::AbstractStressState, m::AbstractMaterial)
     mderiv = MaterialDerivatives(m)
+    T = get_vector_eltype(m)
     np = get_vector_length(m)
     nt = get_num_tensorcomponents(m) # Should be changed to only save non-controlled entries
     # Filled with `NaN`s to ensure we don't rely on these being defined.
-    dϵdp = fill(NaN, nt, np)
-    dσdp = fill(NaN, nt, np)
+    dϵdp = fill(T(NaN), nt, np)
+    dσdp = fill(T(NaN), nt, np)
     # TODO: Should be changed to only save non-controlled entries
     vo = Tensors.DEFAULT_VOIGT_ORDER[3]
     if nt == 6

src/stressiterations.jl

@@ -414,7 +414,7 @@ end
 function get_full_tensor(state::GeneralStressState{Nσ}, ::TT, v::SVector{Nσ,T}) where {Nσ,T,TT}
     TB = Tensors.get_base(TT)
     shear_factor = if TB == SymmetricTensor{2,3}
-        T(1/√2)
+        1/sqrt(2 * one(T))
     elseif TB == Tensor{2,3}
         one(T)
     else
@@ -426,14 +426,14 @@ function get_full_tensor(state::GeneralStressState{Nσ}, ::TT, v::SVector{Nσ,T}
 end
 
 function get_unknowns(state::GeneralStressState{Nσ}, a::AbstractTensor{2,3,T}) where {Nσ, T}
-    shear_factor = a isa SymmetricTensor ? T(√2) : one(T)
+    shear_factor = a isa SymmetricTensor ? sqrt(2 * one(T)) : one(T)
     s(i,j) = i==j ? one(T) : shear_factor
     f(c) = ((i,j) = c; s(i,j)*(a[i,j]-state.σ[i,j]))
     return SVector{Nσ,T}((f(c) for c in state.σm_inds))
 end
 
 function get_unknowns(state::GeneralStressState{Nσ}, a::AbstractTensor{4,3,T}) where {Nσ,T}
-    shear_factor = a isa SymmetricTensor ? T(√2) : one(T)
+    shear_factor = a isa SymmetricTensor ? sqrt(2 * one(T)) : one(T)
     s(i,j) = i==j ? one(T) : shear_factor
     f(c1,c2) = ((i,j) = c1; (k,l) = c2; a[i,j,k,l]*s(i,j)*s(k,l))
     return SMatrix{Nσ,Nσ,T}((f(c1,c2) for c1 in state.σm_inds, c2 in state.σm_inds))