The Burgers equation – the Godonov flux
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