def METHOD_NAME(self):
    m = 5000
    n = 1000
    A = rand(m, n, random_state=0)
    b = self.rnd.randn(m)
    lb = self.rnd.randn(n)
    ub = lb + 1
    res = lsq_linear(A, b, (lb, ub))
    assert_allclose(res.optimality, 0.0, atol=1e-6)
    res = lsq_linear(A, b, (lb, ub), lsmr_tol=1e-13)
    assert_allclose(res.optimality, 0.0, atol=1e-6)
    res = lsq_linear(A, b, (lb, ub), lsmr_tol='auto')
    assert_allclose(res.optimality, 0.0, atol=1e-6)