Just a minor thing: could you add more numerical accuracy here? Usually we do 1E-15