Fixed solver by changing order of temp b.c. and temp gather

ORNL · Nov 6, 2024 · 0be8643 · 0be8643
1 parent 0e4cc00
commit 0be8643
Showing 1 changed file with 4 additions and 3 deletions.
diff --git a/src/CabanaPD_Solver.hpp b/src/CabanaPD_Solver.hpp
@@ -281,14 +281,15 @@ class SolverElastic
                     computeHeatTransfer( *heat_transfer, *particles, *neighbors,
                                          neigh_iter_tag{}, dt );
             }
-            if constexpr ( is_temperature_dependent<
-                               typename force_model_type::thermal_type>::value )
-                comm->gatherTemperature();
 
             // Add non-force boundary condition.
             if ( !boundary_condition.forceUpdate() )
                 boundary_condition.apply( exec_space(), *particles, step * dt );
 
+            if constexpr ( is_temperature_dependent<
+                               typename force_model_type::thermal_type>::value )
+                comm->gatherTemperature();
+
             // Compute internal forces.
             updateForce();