Fixed a bug in the routine for scaline the bond force after a fusion …

…event.
PNNL-CompBio · Sep 19, 2024 · 7a465f4 · 7a465f4
1 parent 1271a81
commit 7a465f4
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/des/bmx_pc_fusion.cpp b/src/des/bmx_pc_fusion.cpp
@@ -159,6 +159,7 @@ bool BMXParticleContainer::EvaluateTipFusion (const Vector<MultiFab*> cost,
 
           }
           } // end of neighbor loop
+          incrementBondScale(&particle.rdata(0),&particle.idata(0),fpar);
           }); // end of loop over particles
 
       amrex::Gpu::Device::synchronize();
@@ -586,7 +587,6 @@ void BMXParticleContainer::CleanupFusion (const Vector<MultiFab*> cost,
           auto& particle = pstruct[i];
 
           // increment bond scaling parameter
-          incrementBondScale(&particle.rdata(0),&particle.idata(0),fpar);
           RealVect pos1(particle.pos());
 
           const auto neighbs = nbor_data.getNeighbors(i);