Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Next patch for cudacpp-only MEs: add dummy color choice to fix 'failed to reduce to color index' #489

Merged
merged 24 commits into from
Jun 15, 2022

Conversation

valassi
Copy link
Member

@valassi valassi commented Jun 15, 2022

This is a workaround for oliviermattelaer/mg5amc_test#14

It is the first step of #402: it implements a dummy color choice (as opposed to a random color choice)

Essentially I have just added JAMP2_MULTI(0,IVEC) = NCOLOR. I have not even initialised JAMP2_MULTI(1:MAXFLOW,IVEC), it seems that addmothers.f manages this by itself somehow

valassi added 24 commits June 15, 2022 17:59
… is not called (dummy random choice of color madgraph5#402)

This is a WORKAROUND for "Error: failed to reduce to color indices" oliviermattelaer/mg5amc_test#14

Note: in this patch I only set JAMP2(0), while all other JAMP2(1:NCOLOR) are zero.
It seems that addmothers.f already has a workaround to cover this case.
…iermattelaer/mg5amc_test#14

(En passant: ggttg addmothers is regenerated without debug printouts)
… cudacpp MEs) - no change as expected

Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/ee_mumu.mad/SubProcesses/P1_ll_ll
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/ee_mumu.mad/SubProcesses/P1_ll_ll
 [XSECTION] Cross section = 0.09017
 [XSECTION] Cross section = 0.09017 [9.0170633677521428E-002]
 [UNWEIGHT] Wrote 1009 events (found 1010 events)
 [XSECTION] Cross section = 0.09017 [9.0170633677521442E-002]
 [UNWEIGHT] Wrote 1009 events (found 1010 events)
 [XSECTION] Cross section = 0.09017 [9.0170633677521442E-002]
 [UNWEIGHT] Wrote 1009 events (found 1010 events)
 [XSECTION] Cross section = 0.2175 [0.21754196695805308]
 [UNWEIGHT] Wrote 966 events (found 967 events)
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttggg.mad/SubProcesses/P1_gg_ttxggg
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttggg.mad/SubProcesses/P1_gg_ttxggg
 [XSECTION] Cross section = 3.447e-07
 [XSECTION] Cross section = 3.447e-07 [3.4471868628442136E-007]
 [UNWEIGHT] Wrote 4 events (found 15 events)
 [XSECTION] Cross section = 3.447e-07 [3.4471908676045729E-007]
 [UNWEIGHT] Wrote 4 events (found 15 events)
 [XSECTION] Cross section = 3.447e-07 [3.4471908676045713E-007]
 [UNWEIGHT] Wrote 4 events (found 15 events)
 [XSECTION] Cross section = 2.008e-06 [2.0075279247605889E-006]
 [UNWEIGHT] Wrote 2 events (found 11 events)
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttgg.mad/SubProcesses/P1_gg_ttxgg
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttgg.mad/SubProcesses/P1_gg_ttxgg
 [XSECTION] Cross section = 0.0001289
 [XSECTION] Cross section = 0.0001289 [1.2885825323149218E-004]
 [UNWEIGHT] Wrote 4 events (found 74 events)
 [XSECTION] Cross section = 0.0001289 [1.2885954287258873E-004]
 [UNWEIGHT] Wrote 4 events (found 74 events)
 [XSECTION] Cross section = 0.0001289 [1.2885954287258873E-004]
 [UNWEIGHT] Wrote 4 events (found 74 events)
 [XSECTION] Cross section = 0.0002221 [2.2210365310367753E-004]
 [UNWEIGHT] Wrote 7 events (found 79 events)
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttg.mad/SubProcesses/P1_gg_ttxg
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttg.mad/SubProcesses/P1_gg_ttxg
 [XSECTION] Cross section = 0.08045
 [XSECTION] Cross section = 0.08045 [8.0445416635721884E-002]
 [UNWEIGHT] Wrote 56 events (found 449 events)
 [XSECTION] Cross section = 0.08045 [8.0446406133394599E-002]
 [UNWEIGHT] Wrote 56 events (found 449 events)
 [XSECTION] Cross section = 0.08045 [8.0446406133394599E-002]
 [UNWEIGHT] Wrote 56 events (found 449 events)
 [XSECTION] Cross section = 0.08402 [8.4018030529591323E-002]
 [UNWEIGHT] Wrote 16 events (found 397 events)
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_tt.mad/SubProcesses/P1_gg_ttx
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_tt.mad/SubProcesses/P1_gg_ttx
 [XSECTION] Cross section = 45.91
 [XSECTION] Cross section = 45.91 [45.914216281363188]
 [UNWEIGHT] Wrote 788 events (found 2238 events)
 [XSECTION] Cross section = 45.91 [45.914455838006880]
 [UNWEIGHT] Wrote 788 events (found 2238 events)
 [XSECTION] Cross section = 45.91 [45.914455838006880]
 [UNWEIGHT] Wrote 788 events (found 2238 events)
 [XSECTION] Cross section = 47.67 [47.668452211994605]
 [UNWEIGHT] Wrote 804 events (found 2280 events)
…ith dummy color choice - better in ggttg*, but eemumu crashes!

The eemumu issue, not captured in the log, is:

Executing ' ./cmadevent_cudacpp < /tmp/avalassi/input_eemumu_cpp > /tmp/avalassi/output_eemumu_cpp'
Program received signal SIGSEGV: Segmentation fault - invalid memory reference.
Backtrace for this error:
./madX.sh: line 201: 20065 Segmentation fault      $timecmd $cmd < ${tmpin} > ${tmp}
ERROR! ' ./cmadevent_cudacpp < /tmp/avalassi/input_eemumu_cpp > /tmp/avalassi/output_eemumu_cpp' failed
 Renormalization scale set on event-by-event basis
 Factorization   scale set on event-by-event basis
 getting user params
Enter number of events and max and min iterations:
 Number of events and iterations         2048           1           1
Enter desired fractional accuracy:
 Desired fractional accuracy:    9.9999999999999995E-007
Enter 0 for fixed, 2 for adjustable grid:

This is what is in the logs:

cat tmad/logs_*/* | egrep '(Cross|Working|Wrote|Error|ERROR)'
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/ee_mumu.mad/SubProcesses/P1_ll_ll
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/ee_mumu.mad/SubProcesses/P1_ll_ll
 [XSECTION] Cross section = 0.09017
 [XSECTION] Cross section = 0.09017 [9.0170633677521428E-002]
 [UNWEIGHT] Wrote 1009 events (found 1010 events)
ERROR! ' ./cmadevent_cudacpp < /tmp/avalassi/input_eemumu_cpp > /tmp/avalassi/output_eemumu_cpp' failed
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttggg.mad/SubProcesses/P1_gg_ttxggg
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttggg.mad/SubProcesses/P1_gg_ttxggg
 [XSECTION] Cross section = 3.447e-07
 [XSECTION] Cross section = 3.447e-07 [3.4471868628442136E-007]
 [UNWEIGHT] Wrote 4 events (found 15 events)
 [XSECTION] Cross section = 3.447e-07 [3.4471908676045729E-007]
 [UNWEIGHT] Wrote 4 events (found 15 events)
 [XSECTION] Cross section = 3.447e-07 [3.4471908676045713E-007]
 [UNWEIGHT] Wrote 4 events (found 15 events)
 [XSECTION] Cross section = 2.008e-06 [2.0075279247605889E-006]
 [UNWEIGHT] Wrote 2 events (found 11 events)
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttgg.mad/SubProcesses/P1_gg_ttxgg
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttgg.mad/SubProcesses/P1_gg_ttxgg
 [XSECTION] Cross section = 0.0001289
 [XSECTION] Cross section = 0.0001289 [1.2885825323149218E-004]
 [UNWEIGHT] Wrote 4 events (found 74 events)
 [XSECTION] Cross section = 0.0001289 [1.2885954287258873E-004]
 [UNWEIGHT] Wrote 4 events (found 74 events)
 [XSECTION] Cross section = 0.0001289 [1.2885954287258873E-004]
 [UNWEIGHT] Wrote 4 events (found 74 events)
 [XSECTION] Cross section = 0.0002221 [2.2210365310367753E-004]
 [UNWEIGHT] Wrote 7 events (found 79 events)
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttg.mad/SubProcesses/P1_gg_ttxg
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttg.mad/SubProcesses/P1_gg_ttxg
 [XSECTION] Cross section = 0.08045
 [XSECTION] Cross section = 0.08045 [8.0445416635721884E-002]
 [UNWEIGHT] Wrote 56 events (found 449 events)
 [XSECTION] Cross section = 0.08045 [8.0446406133394599E-002]
 [UNWEIGHT] Wrote 56 events (found 449 events)
 [XSECTION] Cross section = 0.08045 [8.0446406133394599E-002]
 [UNWEIGHT] Wrote 56 events (found 449 events)
 [XSECTION] Cross section = 0.08402 [8.4018030529591323E-002]
 [UNWEIGHT] Wrote 16 events (found 397 events)
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_tt.mad/SubProcesses/P1_gg_ttx
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_tt.mad/SubProcesses/P1_gg_ttx
 [XSECTION] Cross section = 45.91
 [XSECTION] Cross section = 45.91 [45.914216281363188]
 [UNWEIGHT] Wrote 788 events (found 2238 events)
 [XSECTION] Cross section = 45.91 [45.914455838006880]
 [UNWEIGHT] Wrote 788 events (found 2238 events)
 [XSECTION] Cross section = 45.91 [45.914455838006880]
 [UNWEIGHT] Wrote 788 events (found 2238 events)
 [XSECTION] Cross section = 47.67 [47.668452211994605]
 [UNWEIGHT] Wrote 804 events (found 2280 events)
… madX, revert to previous logs)

Revert "[fvsc] try another hack to collect stderr in madX and rerun, but not better, will revert"
This reverts commit 3614722.

Revert "[fvsc] rerun tmad alltees with new madX that should capture stderr... actually worse than before?"
This reverts commit a11c9ab.

Revert "[fvsc] improve tmad madX.sh script, keep stderr in the log"
This reverts commit e9307aa.
…has random values?!

This is probably because NCOLOR is taken by Fortran to be an IMPLICIT (and uninitialised!) variable, ouf...
Will revert and add IMPLICIT NONE...

 RESET CUMULATIVE VARIABLE
 DEBUG JAMP2_MULTI           1  -524355579
 DEBUG JAMP2_MULTI           2  -524355579
 DEBUG JAMP2_MULTI           3  -524355579
 DEBUG JAMP2_MULTI           4  -524355579
…ements

Indeed now NCOLOR no longer builds...

ccache /cvmfs/sft.cern.ch/lcg/releases/gcc/10.2.0-c44b3/x86_64-centos7/bin/gfortran -w -fPIC  -ffixed-line-length-132 -w -cpp  -c -DMG5AMC_MEEXPORTER_CUDACPP auto_dsig1.f -I../../Source/ -fopenmp -o auto_dsig1_cudacpp.o
auto_dsig1.f:584:38:

  584 |           JAMP2_MULTI(0,IVEC) = NCOLOR ! workaround for oliviermattelaer/mg5amc_test#14
      |                                      1
Error: Symbol ‘ncolor’ at (1) has no IMPLICIT type
…ded NCOLOR (from matrix1.f) into auto_dsig1.f
…sue is fixed for gggttg* generation

This shows that the workaround is ok for oliviermattelaer/mg5amc_test#14

cat tmad/logs_*/* | egrep '(Cross|Working|Wrote|Error|ERROR)'
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/ee_mumu.mad/SubProcesses/P1_ll_ll
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/ee_mumu.mad/SubProcesses/P1_ll_ll
 [XSECTION] Cross section = 0.09017
 [XSECTION] Cross section = 0.09017 [9.0170633677521428E-002]
 [UNWEIGHT] Wrote 1009 events (found 1010 events)
 [XSECTION] Cross section = 0.09017 [9.0170633677521442E-002]
 [UNWEIGHT] Wrote 1009 events (found 1010 events)
 [XSECTION] Cross section = 0.09017 [9.0170633677521442E-002]
 [UNWEIGHT] Wrote 1009 events (found 1010 events)
 [XSECTION] Cross section = 0.2175 [0.21754196695805308]
 [UNWEIGHT] Wrote 966 events (found 967 events)
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttggg.mad/SubProcesses/P1_gg_ttxggg
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttggg.mad/SubProcesses/P1_gg_ttxggg
 [XSECTION] Cross section = 3.447e-07
 [XSECTION] Cross section = 3.447e-07 [3.4471868628442136E-007]
 [UNWEIGHT] Wrote 4 events (found 15 events)
 [XSECTION] Cross section = 3.447e-07 [3.4471908676045729E-007]
 [UNWEIGHT] Wrote 4 events (found 15 events)
 [XSECTION] Cross section = 3.447e-07 [3.4471908676045713E-007]
 [UNWEIGHT] Wrote 4 events (found 15 events)
 [XSECTION] Cross section = 2.008e-06 [2.0075279247605889E-006]
 [UNWEIGHT] Wrote 2 events (found 11 events)
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttgg.mad/SubProcesses/P1_gg_ttxgg
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttgg.mad/SubProcesses/P1_gg_ttxgg
 [XSECTION] Cross section = 0.0001289
 [XSECTION] Cross section = 0.0001289 [1.2885825323149218E-004]
 [UNWEIGHT] Wrote 4 events (found 74 events)
 [XSECTION] Cross section = 0.0001289 [1.2885954287258873E-004]
 [UNWEIGHT] Wrote 4 events (found 74 events)
 [XSECTION] Cross section = 0.0001289 [1.2885954287258873E-004]
 [UNWEIGHT] Wrote 4 events (found 74 events)
 [XSECTION] Cross section = 0.0002221 [2.2210365310367753E-004]
 [UNWEIGHT] Wrote 7 events (found 79 events)
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttg.mad/SubProcesses/P1_gg_ttxg
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttg.mad/SubProcesses/P1_gg_ttxg
 [XSECTION] Cross section = 0.08045
 [XSECTION] Cross section = 0.08045 [8.0445416635721884E-002]
 [UNWEIGHT] Wrote 56 events (found 449 events)
 [XSECTION] Cross section = 0.08045 [8.0446406133394599E-002]
 [UNWEIGHT] Wrote 56 events (found 449 events)
 [XSECTION] Cross section = 0.08045 [8.0446406133394599E-002]
 [UNWEIGHT] Wrote 56 events (found 449 events)
 [XSECTION] Cross section = 0.08402 [8.4018030529591323E-002]
 [UNWEIGHT] Wrote 16 events (found 397 events)
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_tt.mad/SubProcesses/P1_gg_ttx
Working directory: /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_tt.mad/SubProcesses/P1_gg_ttx
 [XSECTION] Cross section = 45.91
 [XSECTION] Cross section = 45.91 [45.914216281363188]
 [UNWEIGHT] Wrote 788 events (found 2238 events)
 [XSECTION] Cross section = 45.91 [45.914455838006880]
 [UNWEIGHT] Wrote 788 events (found 2238 events)
 [XSECTION] Cross section = 45.91 [45.914455838006880]
 [UNWEIGHT] Wrote 788 events (found 2238 events)
 [XSECTION] Cross section = 47.67 [47.668452211994605]
 [UNWEIGHT] Wrote 804 events (found 2280 events)
…1, more compact printouts

cat tmad/logs_*/* | egrep '(run|Cross|Wrote|Error|ERROR)'
Working directory (run): /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/ee_mumu.mad/SubProcesses/P1_ll_ll
 [XSECTION] Cross section = 0.09017
 [XSECTION] Cross section = 0.09017 [9.0170633677521428E-002]
 [UNWEIGHT] Wrote 1009 events (found 1010 events)
 [XSECTION] Cross section = 0.09017 [9.0170633677521442E-002] fbridge_mode=1
 [UNWEIGHT] Wrote 1009 events (found 1010 events)
 [XSECTION] Cross section = 0.09017 [9.0170633677521442E-002] fbridge_mode=1
 [UNWEIGHT] Wrote 1009 events (found 1010 events)
 [XSECTION] Cross section = 0.2175 [0.21754196695805308] fbridge_mode=1
 [UNWEIGHT] Wrote 966 events (found 967 events)
Working directory (run): /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttggg.mad/SubProcesses/P1_gg_ttxggg
 [XSECTION] Cross section = 3.447e-07
 [XSECTION] Cross section = 3.447e-07 [3.4471868628442136E-007]
 [UNWEIGHT] Wrote 4 events (found 15 events)
 [XSECTION] Cross section = 3.447e-07 [3.4471908676045729E-007] fbridge_mode=1
 [UNWEIGHT] Wrote 4 events (found 15 events)
 [XSECTION] Cross section = 3.447e-07 [3.4471908676045713E-007] fbridge_mode=1
 [UNWEIGHT] Wrote 4 events (found 15 events)
 [XSECTION] Cross section = 2.008e-06 [2.0075279247605889E-006] fbridge_mode=1
 [UNWEIGHT] Wrote 2 events (found 11 events)
Working directory (run): /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttgg.mad/SubProcesses/P1_gg_ttxgg
 [XSECTION] Cross section = 0.0001289
 [XSECTION] Cross section = 0.0001289 [1.2885825323149218E-004]
 [UNWEIGHT] Wrote 4 events (found 74 events)
 [XSECTION] Cross section = 0.0001289 [1.2885954287258873E-004] fbridge_mode=1
 [UNWEIGHT] Wrote 4 events (found 74 events)
 [XSECTION] Cross section = 0.0001289 [1.2885954287258873E-004] fbridge_mode=1
 [UNWEIGHT] Wrote 4 events (found 74 events)
 [XSECTION] Cross section = 0.0002221 [2.2210365310367753E-004] fbridge_mode=1
 [UNWEIGHT] Wrote 7 events (found 79 events)
Working directory (run): /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_ttg.mad/SubProcesses/P1_gg_ttxg
 [XSECTION] Cross section = 0.08045
 [XSECTION] Cross section = 0.08045 [8.0445416635721884E-002]
 [UNWEIGHT] Wrote 56 events (found 449 events)
 [XSECTION] Cross section = 0.08045 [8.0446406133394599E-002] fbridge_mode=1
 [UNWEIGHT] Wrote 56 events (found 449 events)
 [XSECTION] Cross section = 0.08045 [8.0446406133394599E-002] fbridge_mode=1
 [UNWEIGHT] Wrote 56 events (found 449 events)
 [XSECTION] Cross section = 0.08402 [8.4018030529591323E-002] fbridge_mode=1
 [UNWEIGHT] Wrote 16 events (found 397 events)
Working directory (run): /data/avalassi/GPU2020/madgraph4gpuX/epochX/cudacpp/gg_tt.mad/SubProcesses/P1_gg_ttx
 [XSECTION] Cross section = 45.91
 [XSECTION] Cross section = 45.91 [45.914216281363188]
 [UNWEIGHT] Wrote 788 events (found 2238 events)
 [XSECTION] Cross section = 45.91 [45.914455838006880] fbridge_mode=1
 [UNWEIGHT] Wrote 788 events (found 2238 events)
 [XSECTION] Cross section = 45.91 [45.914455838006880] fbridge_mode=1
 [UNWEIGHT] Wrote 788 events (found 2238 events)
 [XSECTION] Cross section = 47.67 [47.668452211994605] fbridge_mode=1
 [UNWEIGHT] Wrote 804 events (found 2280 events)
@valassi
Copy link
Member Author

valassi commented Jun 15, 2022

This is a followup to th eprevious PR #486.

Amongst the things that remain to be done before a first milestone (xsections ok, and events ~okish except for random color and helicity) there are these two:

All tests are passing, I am self merging

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant