Skip to content
This repository was archived by the owner on Aug 28, 2024. It is now read-only.

Commit

Permalink
Fixed problem that was causing fraggle to fail when the number of fra…
Browse files Browse the repository at this point in the history
…gments was reduced to find a valid solution
  • Loading branch information
daminton committed Feb 1, 2024
1 parent 16479d2 commit 2f23a61
Showing 1 changed file with 1 addition and 0 deletions.
1 change: 1 addition & 0 deletions src/fraggle/fraggle_generate.f90
Original file line number Diff line number Diff line change
Expand Up @@ -884,6 +884,7 @@ module subroutine fraggle_generate_vel_vec(collider, nbody_system, param, lfailu
L_mag_factor = .mag.(collider%L_total(:,1) + collider%L_total(:,2))
L_residual(:) = (collider%L_total(:,2) / L_mag_factor - collider%L_total(:,1)) / L_mag_factor
call collision_util_velocity_torque(-L_residual(:) * L_mag_factor, collider%fragments%mtot, impactors%rbcom, impactors%vbcom)
nfrag = collider%fragments%nbody

#ifdef DOCONLOC
do concurrent(i = 1:nfrag) shared(collider, impactors)
Expand Down

0 comments on commit 2f23a61

Please sign in to comment.