This is a toolbox to model color glass condensates in the dense-dilute approximation.
You can access the documentation for the project here.
The simplest usage example would be to create nucleus and proton objects to calculate the momenta of particles produced:
First, import the toolbox and whatever other packages you will use for post-processing
import cgc
import matplotlib.pyplot as plt
Define constants for the dilute and dense objects.
N = 128 # Number of lattice sites
delta = .1 # Inter-lattice spacing
mu = 2 # Gaussian distribution width for color field generation
radius = 1 # Radius of dilute proton
Create the dilute and dense objects for (SU(3)), as well as a collision between them. Note that no actual calculations are done in this step; rather, calculations are only performed as they are deemed necessary (ie. when a particular quantity is called)
nucelus = cgc.Nucleus(3, N, delta, mu)
proton = cgc.Proton(3, N, delta, mu, radius)
col = cgc.Collision(nucleus, proton)
From here, quantities can be generated from calling any of the three objects above.
plt.plot(col.momentaBins(), col.particlesProduced())
Computationally expensive calculations (nucleus gauge field, wilson line, and adjoint wilson line) are optimized using the numba library, which transpiles the functions into C code.
(not complete list yet)
Higham, N. J. (2008). Functions of matrices: Theory and computation. Society for Industrial and Applied Mathematics.
Bertlmann, R. A., & Krammer, P. (2008). Bloch vectors for qudits. Journal of Physics A: Mathematical and Theoretical, 41(23), 235303. 10.1088/1751-8113/41/23/235303