pub fn minmod(a: f64, b: f64) -> f64
Minmod slope limiter: smallest slope among a and b with same sign.
a
b
Returns zero when a and b have opposite signs.