This code can read in a freeze out surface from 2+1D or 3+1D viscous hydro codes (cpu-vh, gpu-vh) or viscous anisotropic hydro codes (cpu-vah) and calculate 3D particle spectra. It can perform numerical integration to find the smooth particle spectra, or perform sampling for a particle list. It is accelerated via OpenMP and CUDA.
I used a boost invariant(2+1D) hypersurface file from MUSIC and tried to run the code. but it's showing an error of following type:
[.......]
Reading in freezeout surface in (new) public MUSIC boost invariant format
Loading particle info from urqmd/smash pdg file (check that only 1 blank line at end of file)
...using read_resonances_conventional()
There are 327 resonances: 77 mesons, 125 baryons and 125 antibaryons
Last particle: mcid = -12110, Anti-baryon-n-2250, m = 2.250000 GeV (please check this)
Reading in 14 moment and Chapman-Enskog coefficient tables...done
gsl: interp.c:83: ERROR: x values must be strictly increasing
Default GSL error handler invoked.
Aborted (core dumped)
When running with more than one chosen particle, something is wrong with resulting spectra for each particle (GPU version). Check indexing of particle species, momenta ...