Calculate force density for t=0 inside data handler initialization

Author: kaipartmann

src/core/data_handler.jl

@@ -1,4 +1,5 @@
 function initialize!(dh, solver)
+    calc_force_density!(dh, 0.0, solver.Δt)
     return nothing
 end
 

src/core/mpi_body_data_handler.jl

@@ -421,7 +421,8 @@ function export_reference_results(dh::MPIBodyDataHandler, options::AbstractJobOp
     return nothing
 end
 
-function initialize!(::AbstractMPIBodyDataHandler, ::AbstractTimeSolver)
+function initialize!(dh::AbstractMPIBodyDataHandler, solver::AbstractTimeSolver)
+    calc_force_density!(dh, 0.0, solver.Δt)
     return nothing
 end
 

src/core/threads_multibody_data_handler.jl

@@ -85,8 +85,9 @@ function update_caches!(dh::ThreadsMultibodyDataHandler)
     return nothing
 end
 
-function initialize!(dh::AbstractThreadsMultibodyDataHandler, ::AbstractTimeSolver)
+function initialize!(dh::AbstractThreadsMultibodyDataHandler, solver::AbstractTimeSolver)
     init_contact_nhs!(dh)
+    calc_force_density!(dh, 0.0, solver.Δt)
     return nothing
 end
 

src/discretization/bond_associated_system.jl

@@ -105,20 +105,22 @@ end
 @inline get_hood_volume(chunk::AbstractBodyChunk) = chunk.system.hood_volume
 
 function initialize!(dh::AbstractThreadsBodyDataHandler{<:BondAssociatedSystem},
-                     ::AbstractTimeSolver)
+                     solver::AbstractTimeSolver)
     @threads :static for chunk in dh.chunks
         calc_hood_volumes!(chunk)
     end
     @threads :static for chunk_id in eachindex(dh.chunks)
         exchange_loc_to_halo!(get_hood_volume, dh, chunk_id)
     end
+    calc_force_density!(dh, 0.0, solver.Δt)
     return nothing
 end
 
 function initialize!(dh::AbstractMPIBodyDataHandler{<:BondAssociatedSystem},
-                     ::AbstractTimeSolver)
+                     solver::AbstractTimeSolver)
     calc_hood_volumes!(dh.chunk)
     exchange_loc_to_halo!(get_hood_volume, dh)
+    calc_force_density!(dh, 0.0, solver.Δt)
     return nothing
 end
 

src/discretization/bond_system_corrections.jl

@@ -53,22 +53,24 @@ end
 end
 
 function initialize!(dh::AbstractThreadsBodyDataHandler{BondSystem{EnergySurfaceCorrection}},
-                     ::AbstractTimeSolver)
+                     solver::AbstractTimeSolver)
     @threads :static for chunk in dh.chunks
         calc_mfactor!(chunk)
     end
     @threads :static for chunk_id in eachindex(dh.chunks)
         exchange_loc_to_halo!(get_mfactor, dh, chunk_id)
         calc_scfactor!(dh.chunks[chunk_id])
     end
+    calc_force_density!(dh, 0.0, solver.Δt)
     return nothing
 end
 
 function initialize!(dh::AbstractMPIBodyDataHandler{BondSystem{EnergySurfaceCorrection}},
-                     ::AbstractTimeSolver)
+                     solver::AbstractTimeSolver)
     calc_mfactor!(dh.chunk)
     exchange_loc_to_halo!(get_mfactor, dh)
     calc_scfactor!(dh.chunk)
+    calc_force_density!(dh, 0.0, solver.Δt)
     return nothing
 end