Skip to content

Commit

Permalink
patch
Browse files Browse the repository at this point in the history
  • Loading branch information
sbryngelson authored Nov 20, 2023
1 parent 3a5c9b8 commit 6a38aa4
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions src/simulation/m_qbmm.fpp
Original file line number Diff line number Diff line change
Expand Up @@ -687,14 +687,14 @@ contains
if(j == 1 .or. j == 2) then
drdt2 = -1d0 / (2d0 * dsqrt(moms(4) - moms(2)**2d0))
else
drdt2 = 1d0 / (2d0 * dqrt(moms(4) - moms(2)**2d0))
drdt2 = 1d0 / (2d0 * dsqrt(moms(4) - moms(2)**2d0))
end if
else
! Edge case where variance < 0
if(j == 1 .or. j == 2) then
drdt2 = -1d0 / (2d0 * dsqrt(verysmall))
else
drdt2 = 1d0 / (2d0 * dsqrt(verysmall)))
drdt2 = 1d0 / (2d0 * dsqrt(verysmall))
end if
end if

Expand Down

0 comments on commit 6a38aa4

Please sign in to comment.