Improvements to shock indicator.

We can select between limiter, density, energy as shock indicator.
Dont use predictor for explicit scheme. Added missing call to compute_shock_indicator.
Minor changes.

Suppress warning about unsigned int due to trilinos.
Fixed bug in assert
Adding missing mapping from FEValues and some other places. By default it assumes MappingQ1.