Implements the limiter-free third order logarithmic reconstruction (LDLR). In LDLR the values of the flux at the interface xi+1/2 are
reconstructed in the following way
where
where tol = 0.1h^q, q = 1.4 typically and the lateral derivatives d1 and d2 are obtained in the following way
Release Notes: