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

Event by event random choice of one helicity #570

Merged
merged 65 commits into from
Dec 15, 2022
Merged

Conversation

valassi
Copy link
Member

@valassi valassi commented Dec 14, 2022

Hi @oliviermattelaer @roiser I have finally ~completed the patch for the random choice of helicity #403

This MR does NOT include the random choice of color, which I have not yet done (and may be a bit more complex).

The MR is essentially complete with code regeneration of all processes. I am now rerunning the full tests to check performance and see functionalities.

…omenta' (prepare to add RndNum for color and helicity)
…mented out for now)

There are seven other types that are equivalent to this one and could be removed...
…atrixElements was used instead of PinnedHostBufferGs

HOWEVER this was harmless because the two types are the same!
This is an example why BufferOneFp might simplify the code...
…and device

Also improve the existing copies of gs and mes
./tmad/teeMadX.sh -ggtt +10x -makeclean
./tmad/teeMadX.sh -ggtt +10x -makeclean -fltonly
./tmad/teeMadX.sh -ggtt +10x -makeclean -mixonly

Note that the bridge throughputs are a bit lower as there are many more copies of rnd/sel hel/col
… to simplify helicity choice (madgraph5#403)

Disable OMP loops for the moment
./tmad/teeMadX.sh -ggtt +10x -makeclean -mix
… and new code madgraph5#568 - ok but cuda is 30% slower

I guess the assert is not helping?

Anyway, will revert and do it later, maybe.
…as "new" mulrichannel enabled madgraph5#568

Will revert the last changes anyway and move madgraph5#568 to a separate MR if any...
…st 4 commits

Revert "[lhe] rerun tput tee ggtt with multichannel enabled, all ok, as slow as "new" mulrichannel enabled madgraph5#568"
This reverts commit 445b321.

Revert "[lhe] rerun tput ggtt mad with "#undef MGONGPU_SUPPORTS_MULTICHANNEL" and new code madgraph5#568 - ok but cuda is 30% slower"
This reverts commit 3458706.

Revert "[lhe] in ggtt.mad prototype a simplification of multichannel ifdefs (madgraph5#568) - will revert it and do it separately"
This reverts commit 63e04b9.

Revert "[lhe] rebuild ggtt with #undef MGONGPU_SUPPORTS_MULTICHANNEL and rerun tput, looks good (a bit faster?)"
This reverts commit b9f534b.
./tmad/teeMadX.sh -ggtt +10x -makeclean -mix
./CODEGEN/generateAndCompare.sh gg_tt --mad --nopatch
git diff --no-ext-diff -R gg_tt.mad/Source/dsample.f gg_tt.mad/Source/genps.inc gg_tt.mad/Source/vector.inc gg_tt.mad/SubProcesses/makefile > CODEGEN/MG5aMC_patches/PROD/patch.common
git diff --no-ext-diff -R gg_tt.mad/SubProcesses/P1_gg_ttx/auto_dsig1.f gg_tt.mad/SubProcesses/P1_gg_ttx/driver.f gg_tt.mad/SubProcesses/P1_gg_ttx/matrix1.f > CODEGEN/MG5aMC_patches/PROD/patch.P1
git checkout gg_tt.mad
… gg_tt.mad to CODEGEN

The Fortran codegen was using allow_reverse=True and the cpp codegen allow_reverse=False.
Now moved to allow_reverse=True also in cudacpp.
@valassi
Copy link
Member Author

valassi commented Dec 14, 2022

A few comments

While writing this I see that some tests failed in the CI, so there is more work to be done...

@valassi
Copy link
Member Author

valassi commented Dec 14, 2022

Peculiar, it is some tests that fail, I can reproduce this interactively

  • eemumu fails
  • ggtt (where I developed the code) succeeds
  • ggttg also succeeds
  • ggttgg fails

The failures are that the events with given momenta give MEs which are not those in the refernce files?....

@valassi
Copy link
Member Author

valassi commented Dec 14, 2022

Ah voila, as I thought I had forgotten some hardcoded parameters in the codegen...

epochX/cudacpp/gg_ttgg.mad/SubProcesses/P1_gg_ttxgg/fcheck_sa.f /tmp/avalassi/jgigIx_fcheck_sa.f e43625dcf872b3faa292e8b5df57ecf103fdd8f0 100644 epochX/cudacpp/gg_ttgg.mad/SubProcesses/P1_gg_ttxgg/fcheck_sa.f c0bbf580efa6892ab11c9a340b82f15338a2ce7c 100644
7c7
<       PARAMETER(NEVTMAX=2048*256, NEXTERNAL=6, NP4=4)
---
>       PARAMETER(NEVTMAX=2048*256, NEXTERNAL=4, NP4=4)

@valassi
Copy link
Member Author

valassi commented Dec 14, 2022

OK I fixed the bug and regenerated all processes. It looks better in ggttgg interactively

… (slow)

Also reorder the arguments to keep makeclean and sa last (easier to copy/paste)
…es short (skip 6 ggttggg to gain time)

./tput/allTees.sh -short

This completes the random helicity choice madgraph5#403
@valassi
Copy link
Member Author

valassi commented Dec 15, 2022

This MR is now complete. The CI tests are succeeding, I will self merge

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