Tolerance of np.allclose() in check() function.

Tolerance of np.allclose(atol=..., rtol=...) in L.check() must depend on the dtype of array we check against L. What do you think ?

Would it be better to use np.linalg.norm(y - yref) / np.linalg.norm(yref) < rtol ?

Reference