[vecsizeFIX] madgraph5#629 in ggttmad, modify coupl.inc, auto_dsig1.f…

…, driver.f to make them easier for upstream codegen
valassi · Apr 11, 2023 · 04c085b · 04c085b
1 parent 8217c91
commit 04c085b
Show file tree

Hide file tree

Showing 3 changed files with 17 additions and 7 deletions.
diff --git a/epochX/cudacpp/gg_tt.mad/Source/MODEL/coupl.inc b/epochX/cudacpp/gg_tt.mad/Source/MODEL/coupl.inc
@@ -3,13 +3,16 @@ c      written by the UFO converter
 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 
 C     
-C     NB: VECSIZE_MEMMAX_COUPL must be set equal to VECSIZE_MEMMAX, which is defined in vector.inc
-C     NB: vector.inc is not included in coupl.inc because Fortran has no include guards
+C     NB: VECSIZE_MEMMAX is defined in vector.inc, but this is not
+C     included in coupl.inc because Fortran has no include guards
+C     
+C     NB: therefore coupl.inc introduces the definition of a separate
+C     variable VECSIZE_MEMMAX_COUPL, whose value must be the same as
+C     that of VECSIZE_MEMMAX (this is cross-checked in PROGRAM DRIVER)
 C     
 
       INTEGER VECSIZE_MEMMAX_COUPL
-      PARAMETER (VECSIZE_MEMMAX_COUPL=16384) ! NB: 16k events per GPU grid is the minimum required to fill a V100 GPU
-c     PARAMETER (VECSIZE_MEMMAX_COUPL=32) ! NB: workaround for out-of-memory on Juwels: 32 is enough for no-CUDA builds (issue #498)
+      PARAMETER (VECSIZE_MEMMAX_COUPL=16384)
 
       DOUBLE PRECISION G, ALL_G(VECSIZE_MEMMAX_COUPL)
       COMMON/STRONG/ G, ALL_G
@@ -26,10 +29,17 @@ c     PARAMETER (VECSIZE_MEMMAX_COUPL=32) ! NB: workaround for out-of-memory on
       PARAMETER(NL=2D0)
 
       DOUBLE PRECISION MDL_MB,MDL_MH,MDL_MT,MDL_MTA,MDL_MW,MDL_MZ
+
       COMMON/MASSES/ MDL_MB,MDL_MH,MDL_MT,MDL_MTA,MDL_MW,MDL_MZ
 
+
       DOUBLE PRECISION MDL_WH,MDL_WT,MDL_WW,MDL_WZ
+
       COMMON/WIDTHS/ MDL_WH,MDL_WT,MDL_WW,MDL_WZ
 
-      DOUBLE COMPLEX GC_10(VECSIZE_MEMMAX_COUPL), GC_11(VECSIZE_MEMMAX_COUPL)
+
+      DOUBLE COMPLEX GC_10(VECSIZE_MEMMAX_COUPL),
+     $  GC_11(VECSIZE_MEMMAX_COUPL)
+
       COMMON/COUPLINGS/ GC_10, GC_11
+
diff --git a/epochX/cudacpp/gg_tt.mad/SubProcesses/P1_gg_ttx/auto_dsig1.f b/epochX/cudacpp/gg_tt.mad/SubProcesses/P1_gg_ttx/auto_dsig1.f
@@ -80,7 +80,7 @@ DOUBLE PRECISION FUNCTION DSIG1(PP,WGT,IMODE)
 
       INTEGER SUBDIAG(MAXSPROC),IB(2)
       COMMON/TO_SUB_DIAG/SUBDIAG,IB
-      INCLUDE 'vector.inc' ! defines VECSIZE_MEMMAX
+      INCLUDE 'vector.inc'  ! defines VECSIZE_MEMMAX
       INCLUDE 'coupl.inc'  ! defines VECSIZE_MEMMAX_COUPL
       INCLUDE 'run.inc'
 C     Common blocks

diff --git a/epochX/cudacpp/gg_tt.mad/SubProcesses/P1_gg_ttx/driver.f b/epochX/cudacpp/gg_tt.mad/SubProcesses/P1_gg_ttx/driver.f
@@ -71,7 +71,7 @@ Program DRIVER
 c      integer ncols,ncolflow(maxamps),ncolalt(maxamps),ic
 c      common/to_colstats/ncols,ncolflow,ncolalt,ic
 
-      INCLUDE 'vector.inc' ! defines VECSIZE_MEMMAX
+      INCLUDE 'vector.inc'  ! defines VECSIZE_MEMMAX
       INCLUDE 'coupl.inc'  ! defines VECSIZE_MEMMAX_COUPL
       INTEGER VECSIZE_USED
       DATA VECSIZE_USED/VECSIZE_MEMMAX/ ! initial value