The 232nd R-CCS Cafe - part 2

The exploration of the structure of the phase diagram of strongly interacting matter is an active field of research both in the experimental and theoretical nuclear physics communities with important phenomenological implications for heavy-ion collisions, early universe physics. These calculations are done by discretized the space-time in four dimensional lattice and exploring different parameter regimes like mass, temperature and chemical potential and finally take the vanishing lattice spacing limit(continuum limit). Lattice calculations requires high-performance, highly parallelized code as for larger lattice sizes, smaller masses the computational cost increases drastically such that it pushes the limits of our available computational power. Currently, the algorithm that we use is known as Hamiltonian(hybrid) monte carlo for simulations and in various steps we need to invert a large higher dimensional sparse matrix. The inversion is the most time consuming or most computer resources demanding part of our calculations. In this talk, I will mainly focus on different solvers which we use to accelerate our simulations.