I am looking to solve 1D burgers equation with various random initial conditions. What is the best algorithm to find the exact solution?
One method that is covered in literature is the equal area principle (R. J. LeVeque, Numerical Methods for Conservation Laws). Is this the standard procedure employed in numerics papers on limiters (for ex: this paper) for the exact solutions?
Can anyone provide a pseudo-code or algorithm to implement this method to find the shock location and hence the exact solution?