(test_ptmass) bug fix with previous fix for zero-mass sinks; tighten …

…tolerances in disc tests
danieljprice · Dec 20, 2024 · c003cb6 · c003cb6
1 parent 7948ddc
commit c003cb6
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/src/main/ptmass.F90 b/src/main/ptmass.F90
@@ -1032,7 +1032,7 @@ subroutine update_ptmass(dptmass,xyzmh_ptmass,vxyz_ptmass,fxyz_ptmass,nptmass)
                                 - xyzmh_ptmass(2,1:nptmass)*vxyz_ptmass(1,1:nptmass))
  ! Calculate new masses
  newptmass(1:nptmass)           =xyzmh_ptmass(4,1:nptmass)+dptmass(idmsi,1:nptmass)
- where (newptmass1(1:nptmass) > 0)
+ where (newptmass(1:nptmass) > 0)
     newptmass1(1:nptmass)       = 1./newptmass(1:nptmass)
  elsewhere
     newptmass1(1:nptmass)       = 1.

diff --git a/src/tests/test_ptmass.f90 b/src/tests/test_ptmass.f90
@@ -416,12 +416,14 @@ subroutine test_binary(ntests,npass,string)
           tolang = 6.e-10
        endif
        tolen = 1.2e-2
+       if (use_fourthorder) tolen = 5.e-4
     case(2)
        tolen = 1.2e-3
        if (gravity) tolen = 3.1e-3
 
        if (use_fourthorder) then
           tolang = 2.e-11
+          tolen = 8.e-4
        endif
     case default
        if (calc_gravitwaves .and. itest==1) then