code: def test_eigs_consistency(n, atol):
      vals = np.arange(1, n + 1, dtype=np.float64)
      A = spdiags(vals, 0, n, n)
      np.random.seed(345678)
      X = np.random.rand(n, 2)
      lvals, lvecs = lobpcg(A, X, largest=True, maxiter=100)
      vals, vecs = eigs(A, k=2)
  
      _check_eigen(A, lvals, lvecs, atol=atol, rtol=0)
      assert_allclose(np.sort(vals), np.sort(lvals), atol=1e-14)
  
imports: from __future__ import division, print_function, absolute_import
  import itertools
  import numpy
  from numpy.testing import assert_almost_equal, assert_equal, assert_allclose, assert_array_less
  from scipy import ones, rand, r_, diag, linalg, eye
  from scipy.linalg import eig, eigh, toeplitz
  import scipy.sparse
  from scipy.sparse.linalg.eigen.lobpcg import lobpcg
  from scipy.sparse.linalg import eigs
  from scipy.sparse import spdiags
  import pytest
label_desc: Use image node_id 15 for calling a global function with example usage: test_eigs_consistency(n, atol) without return types