Skip to content

Commit

Permalink
[gpucpp] regenerate all 15 processes after Olivier's latest upstream …
Browse files Browse the repository at this point in the history
…changes, merging to 3.5.2

Most changes are in the version comments (from 3.5.1 to 3.5.2)
There are also some minor changes in genps.f but they look like bug fixes (nincming instead of hardcoded 2)
  • Loading branch information
valassi committed Nov 10, 2023
1 parent aae8ef1 commit 56308e9
Show file tree
Hide file tree
Showing 307 changed files with 1,017 additions and 753 deletions.
31 changes: 16 additions & 15 deletions epochX/cudacpp/ee_mumu.mad/CODEGEN_mad_ee_mumu_log.txt
Original file line number Diff line number Diff line change
Expand Up @@ -14,7 +14,7 @@ Running MG5 in debug mode
* * * * * *
* * * *
* *
* VERSION 3.5.1_lo_vect 2023-08-08 *
* VERSION 3.5.2_lo_vect 2023-11-08 *
* *
* WARNING: UNKNOWN DEVELOPMENT VERSION. *
* WARNING: DO NOT USE FOR PRODUCTION *
Expand Down Expand Up @@ -62,7 +62,7 @@ generate e+ e- > mu+ mu-
No model currently active, so we import the Standard Model
INFO: load particles
INFO: load vertices
DEBUG: model prefixing takes 0.005647420883178711 
DEBUG: model prefixing takes 0.005372047424316406 
INFO: Restrict model sm with file models/sm/restrict_default.dat .
DEBUG: Simplifying conditional expressions 
DEBUG: remove interactions: u s w+ at order: QED=1 
Expand Down Expand Up @@ -161,10 +161,10 @@ Load PLUGIN.CUDACPP_OUTPUT
Addition matrix-element will be done with PLUGIN: CUDACPP_OUTPUT
Output will be done with PLUGIN: CUDACPP_OUTPUT
DEBUG: cformat =  standalone_cudacpp [export_cpp.py at line 3071] 
DEBUG: Entering PLUGIN_ProcessExporter.__init__ (initialise the exporter) [output.py at line 158] 
DEBUG: Entering PLUGIN_ProcessExporter.__init__ (initialise the exporter) [output.py at line 155] 
INFO: initialize a new directory: CODEGEN_mad_ee_mumu
INFO: remove old information in CODEGEN_mad_ee_mumu
DEBUG: Entering PLUGIN_ProcessExporter.copy_template (initialise the directory) [output.py at line 163] 
DEBUG: Entering PLUGIN_ProcessExporter.copy_template (initialise the directory) [output.py at line 160] 
WARNING: File exists /data/avalassi/GPU2023/madgraph4gpuX/MG5aMC/TMPOUT/CODEGEN_mad_ee_mumu 
INFO: Creating subdirectories in directory /data/avalassi/GPU2023/madgraph4gpuX/MG5aMC/TMPOUT/CODEGEN_mad_ee_mumu
WARNING: File exists /data/avalassi/GPU2023/madgraph4gpuX/MG5aMC/TMPOUT/CODEGEN_mad_ee_mumu/Cards 
Expand All @@ -174,7 +174,7 @@ INFO: Generating Helas calls for process: e+ e- > mu+ mu- WEIGHTED<=4 @1
INFO: Processing color information for process: e+ e- > mu+ mu- @1
INFO: Creating files in directory P1_epem_mupmum
DEBUG: kwargs[prefix] = 0 [model_handling.py at line 1058] 
DEBUG: process_exporter_cpp =  <PLUGIN.CUDACPP_OUTPUT.model_handling.PLUGIN_OneProcessExporter object at 0x7f8e1c521af0> [export_v4.py at line 6240] 
DEBUG: process_exporter_cpp =  <PLUGIN.CUDACPP_OUTPUT.model_handling.PLUGIN_OneProcessExporter object at 0x7f64657c44f0> [export_v4.py at line 6262] 
INFO: Creating files in directory .
FileWriter <class 'PLUGIN.CUDACPP_OUTPUT.model_handling.PLUGIN_CPPWriter'> for ././CPPProcess.h
FileWriter <class 'PLUGIN.CUDACPP_OUTPUT.model_handling.PLUGIN_CPPWriter'> for ././CPPProcess.cc
Expand All @@ -191,19 +191,19 @@ INFO: Created files CPPProcess.h and CPPProcess.cc in directory ./.
INFO: Generating Feynman diagrams for Process: e+ e- > mu+ mu- WEIGHTED<=4 @1
INFO: Finding symmetric diagrams for subprocess group epem_mupmum
Generated helas calls for 1 subprocesses (2 diagrams) in 0.004 s
Wrote files for 8 helas calls in 0.100 s
Wrote files for 8 helas calls in 0.098 s
ALOHA: aloha starts to compute helicity amplitudes
ALOHA: aloha creates FFV1 routines
ALOHA: aloha creates FFV2 routines
ALOHA: aloha creates FFV4 routines
ALOHA: aloha creates 3 routines in 0.203 s
DEBUG: Entering PLUGIN_ProcessExporter.convert_model (create the model) [output.py at line 200] 
ALOHA: aloha creates 3 routines in 0.200 s
DEBUG: Entering PLUGIN_ProcessExporter.convert_model (create the model) [output.py at line 197] 
ALOHA: aloha starts to compute helicity amplitudes
ALOHA: aloha creates FFV1 routines
ALOHA: aloha creates FFV2 routines
ALOHA: aloha creates FFV4 routines
ALOHA: aloha creates FFV2_4 routines
ALOHA: aloha creates 7 routines in 0.267 s
ALOHA: aloha creates 7 routines in 0.255 s
<class 'aloha.create_aloha.AbstractRoutine'> FFV1
<class 'aloha.create_aloha.AbstractRoutine'> FFV1
<class 'aloha.create_aloha.AbstractRoutine'> FFV2
Expand All @@ -226,12 +226,13 @@ save configuration file to /data/avalassi/GPU2023/madgraph4gpuX/MG5aMC/TMPOUT/CO
INFO: Use Fortran compiler gfortran
INFO: Use c++ compiler g++
INFO: Generate web pages
DEBUG: 'Entering PLUGIN_ProcessExporter.finalize', self.in_madevent_mode, type(self) =  Entering PLUGIN_ProcessExporter.finalize True <class 'PLUGIN.CUDACPP_OUTPUT.output.PLUGIN_ProcessExporter'> [output.py at line 209] 
DEBUG: 'Entering PLUGIN_ProcessExporter.finalize', self.in_madevent_mode, type(self) =  Entering PLUGIN_ProcessExporter.finalize True <class 'PLUGIN.CUDACPP_OUTPUT.output.PLUGIN_ProcessExporter'> [output.py at line 206] 
DEBUG: cd /data/avalassi/GPU2023/madgraph4gpuX/MG5aMC/TMPOUT/CODEGEN_mad_ee_mumu; patch -p4 -i /data/avalassi/GPU2023/madgraph4gpuX/epochX/cudacpp/CODEGEN/PLUGIN/CUDACPP_SA_OUTPUT/MG5aMC_patches/PROD/patch.common
patching file Source/genps.inc
patching file Source/makefile
patching file SubProcesses/makefile
patching file bin/internal/gen_ximprove.py
Hunk #1 succeeded at 391 (offset 6 lines).
patching file bin/internal/madevent_interface.py
DEBUG: cd /data/avalassi/GPU2023/madgraph4gpuX/MG5aMC/TMPOUT/CODEGEN_mad_ee_mumu/SubProcesses/P1_epem_mupmum; patch -p6 -i /data/avalassi/GPU2023/madgraph4gpuX/epochX/cudacpp/CODEGEN/PLUGIN/CUDACPP_SA_OUTPUT/MG5aMC_patches/PROD/patch.P1
patching file auto_dsig1.f
Expand All @@ -241,16 +242,16 @@ patching file matrix1.f
Hunk #3 succeeded at 230 (offset 9 lines).
Hunk #4 succeeded at 267 (offset 18 lines).
Hunk #5 succeeded at 312 (offset 18 lines).
DEBUG: p.returncode =  0 [output.py at line 235] 
DEBUG: p.returncode =  0 [output.py at line 232] 
Output to directory /data/avalassi/GPU2023/madgraph4gpuX/MG5aMC/TMPOUT/CODEGEN_mad_ee_mumu done.
Type "launch" to generate events from this process, or see
/data/avalassi/GPU2023/madgraph4gpuX/MG5aMC/TMPOUT/CODEGEN_mad_ee_mumu/README
Run "open index.html" to see more information about this process.
quit

real 0m2.189s
real 0m4.853s
user 0m1.653s
sys 0m0.232s
sys 0m0.201s
************************************************************
* *
* W E L C O M E to *
Expand All @@ -263,7 +264,7 @@ sys 0m0.232s
* * * * * *
* * * *
* *
* VERSION 3.5.1_lo_vect *
* VERSION 3.5.2_lo_vect *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* https://server06.fynu.ucl.ac.be/projects/madgraph *
Expand Down Expand Up @@ -297,7 +298,7 @@ launch in debug mode
* * * * * *
* * * *
* *
* VERSION 3.5.1_lo_vect *
* VERSION 3.5.2_lo_vect *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* https://server06.fynu.ucl.ac.be/projects/madgraph *
Expand Down
2 changes: 1 addition & 1 deletion epochX/cudacpp/ee_mumu.mad/Cards/proc_card_mg5.dat
Original file line number Diff line number Diff line change
Expand Up @@ -8,7 +8,7 @@
#* * * *
#* *
#* *
#* VERSION 3.5.1_lo_vect 2023-08-08 *
#* VERSION 3.5.2_lo_vect 2023-11-08 *
#* *
#* WARNING: UNKNOWN DEVELOPMENT VERSION. *
#* WARNING: DO NOT USE FOR PRODUCTION *
Expand Down
2 changes: 1 addition & 1 deletion epochX/cudacpp/ee_mumu.mad/MGMEVersion.txt
Original file line number Diff line number Diff line change
@@ -1 +1 @@
3.5.1_lo_vect
3.5.2_lo_vect
2 changes: 1 addition & 1 deletion epochX/cudacpp/ee_mumu.mad/SubProcesses/MGVersion.txt
Original file line number Diff line number Diff line change
@@ -1 +1 @@
3.5.1_lo_vect
3.5.2_lo_vect
Original file line number Diff line number Diff line change
Expand Up @@ -7,7 +7,7 @@
// Further modified by: S. Hageboeck, O. Mattelaer, S. Roiser, A. Valassi, Z. Wettersten (2020-2023) for the MG5aMC CUDACPP plugin.
//==========================================================================
// This file has been automatically generated for CUDA/C++ standalone by
// MadGraph5_aMC@NLO v. 3.5.1_lo_vect, 2023-08-08
// MadGraph5_aMC@NLO v. 3.5.2_lo_vect, 2023-11-08
// By the MadGraph5_aMC@NLO Development Team
// Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
//==========================================================================
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -7,7 +7,7 @@
// Further modified by: O. Mattelaer, S. Roiser, A. Valassi (2020-2023) for the MG5aMC CUDACPP plugin.
//==========================================================================
// This file has been automatically generated for CUDA/C++ standalone by
// MadGraph5_aMC@NLO v. 3.5.1_lo_vect, 2023-08-08
// MadGraph5_aMC@NLO v. 3.5.2_lo_vect, 2023-11-08
// By the MadGraph5_aMC@NLO Development Team
// Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
//==========================================================================
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -359,7 +359,7 @@ SUBROUTINE DSIG_VEC(ALL_P,ALL_WGT,ALL_XBK,ALL_Q2FACT,ALL_CM_RAP
DOUBLE PRECISION FUNCTION DSIG(PP,WGT,IMODE)
C ****************************************************
C
C Generated by MadGraph5_aMC@NLO v. 3.5.1_lo_vect, 2023-08-08
C Generated by MadGraph5_aMC@NLO v. 3.5.2_lo_vect, 2023-11-08
C By the MadGraph5_aMC@NLO Development Team
C Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
C
Expand Down
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
DOUBLE PRECISION FUNCTION DSIG1(PP,WGT,IMODE)
C ****************************************************
C
C Generated by MadGraph5_aMC@NLO v. 3.5.1_lo_vect, 2023-08-08
C Generated by MadGraph5_aMC@NLO v. 3.5.2_lo_vect, 2023-11-08
C By the MadGraph5_aMC@NLO Development Team
C Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
C
Expand Down Expand Up @@ -225,7 +225,7 @@ DOUBLE PRECISION FUNCTION DSIG1_VEC(ALL_PP, ALL_XBK, ALL_Q2FACT,
$ ALL_CM_RAP, ALL_WGT, IMODE, ALL_OUT, VECSIZE_USED)
C ****************************************************
C
C Generated by MadGraph5_aMC@NLO v. 3.5.1_lo_vect, 2023-08-08
C Generated by MadGraph5_aMC@NLO v. 3.5.2_lo_vect, 2023-11-08
C By the MadGraph5_aMC@NLO Development Team
C Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
C
Expand Down
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
SUBROUTINE SMATRIX1(P, RHEL, RCOL, CHANNEL, IVEC, ANS, IHEL,
$ ICOL)
C
C Generated by MadGraph5_aMC@NLO v. 3.5.1_lo_vect, 2023-08-08
C Generated by MadGraph5_aMC@NLO v. 3.5.2_lo_vect, 2023-11-08
C By the MadGraph5_aMC@NLO Development Team
C Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
C
Expand Down Expand Up @@ -319,7 +319,7 @@ SUBROUTINE SMATRIX1(P, RHEL, RCOL, CHANNEL, IVEC, ANS, IHEL,
REAL*8 FUNCTION MATRIX1(P,NHEL,IC, IHEL,AMP2, JAMP2, IVEC)
C
C Generated by MadGraph5_aMC@NLO v. 3.5.1_lo_vect, 2023-08-08
C Generated by MadGraph5_aMC@NLO v. 3.5.2_lo_vect, 2023-11-08
C By the MadGraph5_aMC@NLO Development Team
C Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
C
Expand Down
4 changes: 2 additions & 2 deletions epochX/cudacpp/ee_mumu.mad/SubProcesses/genps.f
Original file line number Diff line number Diff line change
Expand Up @@ -1877,12 +1877,12 @@ double precision function get_channel_cut(p, config)
d1 = iforest(1, -i, config)
d2 = iforest(2, -i, config)
do j=0,3
if (d1.gt.0.and.d1.le.2) then
if (d1.gt.0.and.d1.le.nincoming) then
ptemp(j,-i) = ptemp(j,-i) - ptemp(j, d1)
else
ptemp(j,-i) = ptemp(j,-i)+ptemp(j, d1)
endif
if (d2.gt.0.and.d2.le.2) then
if (d2.gt.0.and.d2.le.nincoming) then
ptemp(j,-i) = ptemp(j,-i) - ptemp(j, d2)
else
ptemp(j,-i) = ptemp(j,-i)+ptemp(j, d2)
Expand Down
1 change: 1 addition & 0 deletions epochX/cudacpp/ee_mumu.mad/bin/internal/__init__.py
Original file line number Diff line number Diff line change
Expand Up @@ -26,6 +26,7 @@ class aMCatNLOError(MadGraph5Error):
import os
import logging
import time
pjoin = os.path.join

#Look for basic file position MG5DIR and MG4DIR
MG5DIR = os.path.realpath(os.path.join(os.path.dirname(__file__),
Expand Down
5 changes: 4 additions & 1 deletion epochX/cudacpp/ee_mumu.mad/bin/internal/banner.py
Original file line number Diff line number Diff line change
Expand Up @@ -4877,6 +4877,9 @@ def create_default_for_process(self, proc_characteristic, history, proc_def):
continue
break

if proc_characteristic['ninitial'] == 1:
self['SDE_strategy'] =1

if 'MLM' in proc_characteristic['limitations']:
if self['dynamical_scale_choice'] == -1:
self['dynamical_scale_choice'] = 3
Expand Down Expand Up @@ -5942,7 +5945,7 @@ def default_setup(self):
self.add_param("CheckCycle", 3)
self.add_param("MaxAttempts", 10)
self.add_param("ZeroThres", 1e-9)
self.add_param("OSThres", 1.0e-13)
self.add_param("OSThres", 1.0e-8)
self.add_param("DoubleCheckHelicityFilter", True)
self.add_param("WriteOutFilters", True)
self.add_param("UseLoopFilter", False)
Expand Down
17 changes: 14 additions & 3 deletions epochX/cudacpp/ee_mumu.mad/bin/internal/common_run_interface.py
Original file line number Diff line number Diff line change
Expand Up @@ -4906,6 +4906,7 @@ def __init__(self, question, cards=[], from_banner=None, banner=None, mode='auto
self.load_default()
self.define_paths(**opt)
self.last_editline_pos = 0
self.update_dependent_done = False

if 'allow_arg' not in opt or not opt['allow_arg']:
# add some mininal content for this:
Expand Down Expand Up @@ -6585,7 +6586,9 @@ def postcmd(self, stop, line):
self.check_card_consistency()
if self.param_consistency:
try:
self.do_update('dependent', timer=20)
if not self.update_dependent_done:
self.do_update('dependent', timer=20)
self.update_dependent_done = False
except MadGraph5Error as error:
if 'Missing block:' in str(error):
self.fail_due_to_format +=1
Expand Down Expand Up @@ -6638,6 +6641,8 @@ def do_update(self, line, timer=0):
self.update_dependent(self.mother_interface, self.me_dir, self.param_card,
self.paths['param'], timer, run_card=self.run_card,
lhapdfconfig=self.lhapdf)
self.update_dependent_done = True


elif args[0] == 'missing':
self.update_missing()
Expand Down Expand Up @@ -6717,12 +6722,13 @@ class TimeOutError(Exception):
def handle_alarm(signum, frame):
raise TimeOutError
signal.signal(signal.SIGALRM, handle_alarm)

if timer:
signal.alarm(timer)
log_level=30
else:
log_level=20


if run_card:
as_for_pdf = {'cteq6_m': 0.118,
'cteq6_d': 0.118,
Expand Down Expand Up @@ -6781,6 +6787,10 @@ def handle_alarm(signum, frame):
logger.log(log_level, "update the strong coupling value (alpha_s) to the value from the pdf selected: %s", as_for_pdf[pdlabel])
modify = True

if timer:
signal.alarm(timer)


# Try to load the model in the limited amount of time allowed
try:
model = mecmd.get_model()
Expand Down Expand Up @@ -6909,7 +6919,8 @@ def check_block(self, blockname):
def check_answer_consistency(self):
"""function called if the code reads a file"""
self.check_card_consistency()
self.do_update('dependent', timer=20)
if not self.update_dependent_done:
self.do_update('dependent', timer=20)

def help_set(self):
'''help message for set'''
Expand Down
18 changes: 12 additions & 6 deletions epochX/cudacpp/ee_mumu.mad/bin/internal/gen_ximprove.py
Original file line number Diff line number Diff line change
Expand Up @@ -157,12 +157,16 @@ def get_helicity(self, to_submit=True, clean=True):

(stdout, _) = p.communicate(''.encode())
stdout = stdout.decode('ascii',errors='ignore')
try:
if stdout:
nb_channel = max([math.floor(float(d)) for d in stdout.split()])
except Exception as error:
misc.sprint(stdout, 'no channel or error for %s' % Pdir)
continue

else:
for matrix_file in misc.glob('matrix*orig.f', Pdir):
files.cp(matrix_file, matrix_file.replace('orig','optim'))
P_zero_result.append(Pdir)
if os.path.exists(pjoin(self.me_dir, 'error')):
os.remove(pjoin(self.me_dir, 'error'))
continue # bypass bad process

self.cmd.compile(['madevent_forhel'], cwd=Pdir)
if not os.path.exists(pjoin(Pdir, 'madevent_forhel')):
raise Exception('Error make madevent_forhel not successful')
Expand All @@ -183,11 +187,13 @@ def get_helicity(self, to_submit=True, clean=True):
#sym_input = "%(points)d %(iterations)d %(accuracy)f \n" % self.opts
(stdout, _) = p.communicate(" ".encode())
stdout = stdout.decode('ascii',errors='ignore')
if os.path.exists(pjoin(self.me_dir,'error')):
if os.path.exists(pjoin(self.me_dir, 'error')):
raise Exception(pjoin(self.me_dir,'error'))
# note a continue is not enough here, we have in top to link
# the matrixX_optim.f to matrixX_orig.f to let the code to work
# after this error.
# for matrix_file in misc.glob('matrix*orig.f', Pdir):
# files.cp(matrix_file, matrix_file.replace('orig','optim'))

if 'no events passed cuts' in stdout:
raise Exception
Expand Down
4 changes: 2 additions & 2 deletions epochX/cudacpp/ee_mumu.mad/bin/internal/launch_plugin.py
Original file line number Diff line number Diff line change
Expand Up @@ -57,7 +57,7 @@ def reset_simd(self, old_value, new_value, name):
return
Sourcedir = pjoin(os.path.dirname(os.path.dirname(self.path)), 'Source')
subprocess.call(['make', 'cleanavx'], cwd=Sourcedir, stdout=subprocess.DEVNULL, stderr=subprocess.DEVNULL)

def plugin_input(self, finput):
return

Expand All @@ -79,7 +79,7 @@ def check_validity(self):
self['sde_strategy'] = 1
if self['hel_recycling']:
self['hel_recycling'] = False

class GPURunCard(CPPRunCard):
def default_setup(self):
super(CPPRunCard, self).default_setup()
Expand Down
Loading

0 comments on commit 56308e9

Please sign in to comment.