Rheolef  7.1
an efficient C++ finite element environment
burgers_flux_godunov.icc

The Burgers equation – the Godonov flux

Float phi (const point& nu, Float a, Float b) {
if ((nu[0] >= 0 && a <= b) || (nu[0] <= 0 && a >= b))
return nu[0]*min(sqr(a),sqr(b))/2;
else
return nu[0]*max(sqr(a),sqr(b))/2;
}
Float phi(const point &nu, Float a, Float b)
see the Float page for the full documentation
see the point page for the full documentation
Definition: nu.h:26