Make the calc_force_density function dispatchable

Author: kaipartmann

src/auxiliary/nans.jl

@@ -5,8 +5,8 @@ function containsnan(K::T) where {T<:AbstractArray}
     return false
 end
 
-function nancheck(chunk::AbstractBodyChunk, t, Δt)
-    if containsnan(chunk.storage.b_int)
+function nancheck(storage::AbstractStorage, t, Δt)
+    if containsnan(storage.b_int)
         n = (t > 0 && Δt > 0) ? Int(t ÷ Δt) : 0
         msg = "NaN's found in field `b_int` at simulation time $(t), step $(n)!\n"
         error(msg)

src/discretization/bond_system.jl

@@ -257,15 +257,22 @@ end
 
 function calc_force_density!(chunk::AbstractBodyChunk{<:AbstractBondSystem}, t, Δt)
     (; system, mat, paramsetup, storage) = chunk
+    calc_force_density!(storage, system, mat, paramsetup, t, Δt)
+    return nothing
+end
+
+function calc_force_density!(storage::AbstractStorage, system::AbstractBondSystem,
+                             mat::AbstractBondSystemMaterial,
+                             paramsetup::AbstractParameterSetup, t, Δt)
     (; dmgmodel) = mat
-    storage.b_int .= 0
+    storage.b_int .= 0.0
     storage.n_active_bonds .= 0
-    for i in each_point_idx(chunk)
+    for i in each_point_idx(system)
         calc_failure!(storage, system, mat, dmgmodel, paramsetup, i)
         calc_damage!(storage, system, mat, dmgmodel, paramsetup, i)
         force_density_point!(storage, system, mat, paramsetup, t, Δt, i)
     end
-    nancheck(chunk, t, Δt)
+    nancheck(storage, t, Δt)
     return nothing
 end