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: