@yvesrevaz, grackle 3.0 is implemented with the different primordial chemistry modes. Unfortunately, the code fails at low density and high temperature with float.
In mode 3 with float:
In mode 3 with double:
