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

SU(3) support #44

Open
wants to merge 129 commits into
base: master
Choose a base branch
from
Open

SU(3) support #44

wants to merge 129 commits into from

Conversation

simone-romiti
Copy link
Collaborator

@simone-romiti simone-romiti commented Mar 7, 2023

todo: write a program that can merge 2 .h5
…. 14 of https://arxiv.org/pdf/hep-lat/9804008.pdf. They are the same as before, no previous result needs to be changed
… format more easily readable by hadron.

To do:
- write the documentation on why this is done separately and not in the c++ code
- change the script such that the data in the old format can be removed
- generalize the script to do the same also for the correlators
script can be used for correlators and interpolators
HighFive dependence now optional
	modified:   hadron/glueball.py
verbose output + fixed wrong index in loop (some i_j pair were not analyzed)
Note: this is not an issue for our previous runs (even if we restarted the run from a different geometry) because:
- the geometry (and size) of the configuration is not determined by what is read from the binary file
- for SU(2): both the HMC and the Metropolis algorithms use the restoreSU() at the end of the trajectory (sweep). For U(1) this function acts trivially.
This implied that at the end of the trajectory everything
Still not working:

- sometimes gives segmentation fault (check weights size)
- acceptance rate almost zero
To do:
Test the behavior of the plaquette expectation value as a function of the volume. as $L \to \infty$ the obc and the pbc should give the same result
implementing support for SU(3) fields
implementing eq. 4.26 of Gattringer&Lang
1st version that compiles successfully. E and dH in the hmc give nan, still debugging.
- exponentiation of su(3) works. The formulas were found using sympy and I made a numerical test of the unitarity of the matrix.
- get_deriv() still implemented incorrectly
simone-romiti and others added 30 commits February 1, 2024 11:09
…nto su3_heatbath_overrelaxation on Christiane's Laptop
…nto su3_heatbath_overrelaxation on Christianes laptop
Prevent double use of seeds in metropolis with same mechanism as in heatbath.
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.

2 participants