derekeverett/is3d

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.

is3d's Issues

Generalize delta-f coefficient calculation

Make a script to automatically generate the delta-f coefficient tables for arbitrary pdg file and chosen resonance file.

Resonance Decays

Still under development. Need to get compilation working and testing.

GSL interpolation error

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)

particle indexing in CUDA version

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 ...