def METHOD_NAME(self):
    m = 5000
    n = 1000
    A = rand(m, n, random_state=0)
    b = self.rnd.randn(m)
    res = lsq_linear(A, b)
    assert_allclose(res.optimality, 0, atol=1e-6)
    A = aslinearoperator(A)
    res = lsq_linear(A, b)
    assert_allclose(res.optimality, 0, atol=1e-6)