Add NinjaRS reduction and bare-loop UV/IR pole selection to MadLoop

.gitattributes

@@ -4,3 +4,6 @@
 .gitattributes export-ignore
 .pre-commit-config.yaml export-ignore
 .bzrignore export-ignore
+
+# Preserve the frozen report bytes and keep compressed PDF streams out of diffs.
+doc/benchmarkfolder/*.pdf binary

.github/workflows/unittest.yml

@@ -55,6 +55,15 @@ jobs:
       - uses: actions/checkout@v5
       - uses: ./.github/actions/clean_tmp
 
+      - name: "MadLoop NinjaRS interfaces, defaults, and reporting"
+        if: ${{ !cancelled() }}
+        run: |
+            cd $GITHUB_WORKSPACE
+            command -v gfortran
+            python3 -m unittest discover -s tests/unit_tests/loop -p 'test_ninjars*.py'
+            python3 -m unittest discover -s tests/unit_tests/various -p 'test_ninjars*.py'
+            python3 -m unittest tests.unit_tests.loop.test_loop_poles tests.unit_tests.interface.test_madevent_loop_warnings
+
       # ---- unittest_0 ----
       - name: "unittest_0"
         if: ${{ !cancelled() }}

.gitignore

@@ -1,4 +1,10 @@
 *~
+.benchmark-work/
+.ninjars-build/
+vendor/ninjars/target/
+vendor/ninjars/smoke
+vendor/ninjars/ninjars-native-libs.txt
+vendor/ninjars/ninjars-link-before.txt
 __pycache__
 *.o
 *.pyc
@@ -81,6 +87,12 @@ tests/input_files/IOTestsComparison/**/*.BackUp
 # output of ./tests/test_manager.py
 UNITTEST_proc/
 
+# Local development references and report plotting caches.
+/references/
+/doc/benchmarkfolder/final-report-v1/.mpl-cache/
+/doc/benchmarkfolder/final-report-v1/cli-report.json.gz
+/doc/benchmarkfolder/final-report-v1/fonts/*.otf
+
 # models/ only holds UFO models: downloaded ones are ignored, shipped ones kept
 models/*
 !models/sm/

Template/LO/Source/dsample.f

@@ -432,6 +432,9 @@ subroutine output_run_statistics(outUnit)
 c
 c     Global
 c
+      INTEGER NINJARS_DP_USAGE,NINJARS_QP_USAGE
+      COMMON/NINJARSSTATS/NINJARS_DP_USAGE,NINJARS_QP_USAGE
+      DATA NINJARS_DP_USAGE,NINJARS_QP_USAGE/0,0/
       INTEGER U_RETURN_CODES(0:9)
       INTEGER T_RETURN_CODES(0:9)
       INTEGER H_RETURN_CODES(0:9)
@@ -468,6 +471,10 @@ subroutine output_run_statistics(outUnit)
 
       write(outUnit,*) '<run_statistics> '
       write(outUnit,33) '<u_return_code>',U_RETURN_CODES,'</u_return_code>'
+      write(outUnit,'(A,I0,A)') '<ninjars_dp_usage>',NINJARS_DP_USAGE,
+     &     '</ninjars_dp_usage>'
+      write(outUnit,'(A,I0,A)') '<ninjars_qp_usage>',NINJARS_QP_USAGE,
+     &     '</ninjars_qp_usage>'
       write(outUnit,33) '<t_return_code>',T_RETURN_CODES,'</t_return_code>'
       write(outUnit,33) '<h_return_code>',H_RETURN_CODES,'</h_return_code>'
       write(outUnit,*) '<average_time>'//trim(toStr_real(AVG_TIMING))

Template/NLO/SubProcesses/BinothLHA.f

@@ -43,6 +43,8 @@ subroutine BinothLHA(p,born_wgt,virt_wgt)
       integer nsqso, MLResArrayDim
 c statistics for MadLoop
       double precision avgPoleRes(2),PoleDiff(2)
+      INTEGER NINJARS_DP_USAGE,NINJARS_QP_USAGE
+      COMMON/NINJARSSTATS/NINJARS_DP_USAGE,NINJARS_QP_USAGE
       integer ntot,nsun,nsps,nups,neps,n100,nddp,nqdp,nini,n10,n1(0:9)
       common/ups_stats/ntot,nsun,nsps,nups,neps,n100,nddp,nqdp,nini,n10,n1
       parameter (nbadmax = 5)
@@ -495,6 +497,10 @@ subroutine BinothLHA(p,born_wgt,virt_wgt)
          n10=n10+1              ! no known ret_code (10)
       endif
       n1(mod(ret_code,10))=n1(mod(ret_code,10))+1 ! unit ret code distribution
+      if (ret_code/100.eq.2.and.mod(ret_code,10).eq.8)
+     $     NINJARS_DP_USAGE=NINJARS_DP_USAGE+1
+      if (ret_code/100.eq.3.and.mod(ret_code,10).eq.2)
+     $     NINJARS_QP_USAGE=NINJARS_QP_USAGE+1
 
 c Write out the unstable, non-rescued phase-space points (MadLoop return
 c code is in the four hundreds) or the ones that are found by the pole

Template/NLO/SubProcesses/driver_mintFO.f

@@ -56,6 +56,8 @@ Program DRIVER
       common/SHevents/Hevents
       character*10 dum
 c statistics for MadLoop      
+      INTEGER NINJARS_DP_USAGE,NINJARS_QP_USAGE
+      COMMON/NINJARSSTATS/NINJARS_DP_USAGE,NINJARS_QP_USAGE
       integer ntot,nsun,nsps,nups,neps,n100,nddp,nqdp,nini,n10,n1(0:9)
       common/ups_stats/ntot,nsun,nsps,nups,neps,n100,nddp,nqdp,nini,n10,n1
 
@@ -127,6 +129,8 @@ Program DRIVER
       do i=0,9
         n1(i)=0
       enddo
+      NINJARS_DP_USAGE=0
+      NINJARS_QP_USAGE=0
 
       call setrun                !Sets up run parameters
       call setpara('param_card.dat')   !Sets up couplings and masses
@@ -248,6 +252,8 @@ Program DRIVER
               write(*,*) "#Unit ",j," = ",n1(j)
            endif
          enddo
+         write(*,*) '#NinjaRS DP = ',NINJARS_DP_USAGE
+         write(*,*) '#NinjaRS QP = ',NINJARS_QP_USAGE
       endif
 
       write (*,*) 'counters for the granny resonances'

Template/NLO/SubProcesses/driver_mintMC.f

@@ -72,6 +72,8 @@ Program DRIVER
       character*10 dum
       integer iFKS_picked
 c statistics for MadLoop      
+      INTEGER NINJARS_DP_USAGE,NINJARS_QP_USAGE
+      COMMON/NINJARSSTATS/NINJARS_DP_USAGE,NINJARS_QP_USAGE
       integer ntot,nsun,nsps,nups,neps,n100,nddp,nqdp,nini,n10,n1(0:9)
       common/ups_stats/ntot,nsun,nsps,nups,neps,n100,nddp,nqdp,nini,n10,n1
 
@@ -133,6 +135,8 @@ Program DRIVER
       do i=0,9
         n1(i)=0
       enddo
+      NINJARS_DP_USAGE=0
+      NINJARS_QP_USAGE=0
 
       call setrun                !Sets up run parameters
       call setpara('param_card.dat')   !Sets up couplings and masses
@@ -340,6 +344,8 @@ Program DRIVER
               write(*,*) "#Unit ",j," = ",n1(j)
            endif
          enddo
+         write(*,*) '#NinjaRS DP = ',NINJARS_DP_USAGE
+         write(*,*) '#NinjaRS QP = ',NINJARS_QP_USAGE
       endif
 
       write (*,*) 'counters for the granny resonances'

Template/loop_material/StandAlone/Cards/MadLoopParams.dat

@@ -14,13 +14,26 @@
 !     5 | Samurai
 !     6 | Ninja
 !     7 | COLLIER
+!     8 | NinjaRS (oneloopmaster scalar integrals)
 !     One can use the combinations to reduce integral,e.g.
 !     1|2|3 means first use CutTools, if it is not stable, use PJFry++, 
 !     if it is still unstable, use IREGI. If it failed, use QP of CutTools.
 !     Notice that any reduction tool not available on the system will be automatically
 !     skipped. 
 
-!     When using quadruple precision with Ninja or CutTools, the reduction will 
+#LoopPoleMode
+!-1
+! Default :: -1
+!     Select poles of the bare loop integral, leaving finite terms and
+!     counterterms/R2 unchanged: -1 = legacy backend behavior, 0 = UV+IR,
+!     1 = UV only, 2 = IR only, 3 = no poles.
+!     Modes 1, 2 and 3 require MLReductionLib containing only NinjaRS (8)
+!     and/or COLLIER (7), including all rescue backends.
+!     An explicit mode overrides COLLIERComputeUVpoles/IRpoles; -1 keeps
+!     those legacy switches active for COLLIER and both poles for NinjaRS.
+!     For example, use MLReductionLib=8|7 and LoopPoleMode=1 for UV poles.
+
+!     When using quadruple precision with Ninja, NinjaRS or CutTools, the reduction will
 !     always be done in quadruple precision, but the parameters below allow you to
 !     chose if you want to also recompute the *integrand* in quadruple precision.
 !     Doing so is slow but might improve the accuracy in some situation.
@@ -30,6 +43,9 @@
 #UseQPIntegrandForNinja
 !.TRUE.
 ! Default :: .TRUE.
+#UseQPIntegrandForNinjaRS
+!.TRUE.
+! Default :: .TRUE.
 !
 
 !     =================================================================================

Template/loop_material/StandAlone/SubProcesses/MadLoopCommons.inc

@@ -270,7 +270,7 @@ C         Not stable.
 C     U == 1
 C         Stable with CutTools in double precision.
 C     U == 2
-C         Stable with PJFry++.
+C         Stable with PJFry++ (H=2), or NinjaRS in quadruple precision (H=3).
 C     U == 3
 C         Stable with IREGI.
 C     U == 4
@@ -282,7 +282,7 @@ C         Stable with Ninja in double precision.
 C     U == 7
 C         Stable with COLLIER.
 C     U == 8
-C         Stable with Ninja in quadruple precision.
+C         Stable with Ninja in quadruple precision (H=3), or NinjaRS in double precision (H=2).
 C     U == 9
 C         Stable with CutTools in quadruple precision.
 C
@@ -315,14 +315,17 @@ C
          ELSEIF(MLRed.EQ.6)THEN
             SET_RET_CODE_U=8
             RETURN
-		 ELSE
-            STOP 'Only CutTools and Ninja can use quardruple precision'
+         ELSEIF(MLRed.EQ.8)THEN
+            SET_RET_CODE_U=2
+            RETURN
+         ELSE
+            STOP 'Only CutTools, Ninja and NinjaRS can use quadruple precision'
          ENDIF
       ENDIF
-      IF(MLRed.GE.1.AND.MLRed.LE.7)THEN
+      IF(MLRed.GE.1.AND.MLRed.LE.8)THEN
          SET_RET_CODE_U=MLRed
       ELSE
-         STOP 'Only CutTools, PJFry++, IREGI, Golem95, Samurai, Ninja and COLLIER are available'
+         STOP 'Unknown MadLoop reduction library; supported IDs are 1 through 8'
       ENDIF
       END
 
@@ -370,8 +373,11 @@ C Ninja
       ELSEIF(LIBNUM.EQ.7)THEN
 C Collier 
          CALL DETECT_COLLIER(NLOOPLINE,RANK,complex_mass,LPASS)
+      ELSEIF(LIBNUM.EQ.8)THEN
+C NinjaRS
+         CALL DETECT_NINJARS(NLOOPLINE,RANK,complex_mass,LPASS)
       ELSE
-         STOP 'Only CutTools, PJFry++, IREGI, Golem95, Samurai, Ninja and COLLIER are available'
+         STOP 'Unknown MadLoop reduction library; supported IDs are 1 through 8'
       ENDIF
       RETURN
       END
@@ -479,6 +485,15 @@ C     It can easily be increased if necessary.
       RETURN
       END
 
+      SUBROUTINE DETECT_NINJARS(NLOOPLINE,RANK,complex_mass,LPASS)
+C NinjaRS tensor reduction supports rank at most n+2 in both precisions.
+      IMPLICIT NONE
+      INTEGER NLOOPLINE,RANK
+      LOGICAL complex_mass,LPASS
+      LPASS=NLOOPLINE.GE.1.AND.RANK.GE.0.AND.
+     &      RANK.LE.NLOOPLINE+2
+      END
+
 ## if(collier_available){
       SUBROUTINE DETECT_COLLIER(NLOOPLINE,RANK,complex_mass,LPASS)
 C