Viscoelastic: stabilized
(FENE-P works, though non-rigorously: c taken arbitrarily large)
Gavrilyuk: bad solver, bad choice of variables (entropy: full gas dynamics !)