I want to release conjugate gradient method to solve Ax = f by using the scipy function
from scipy.sparse.linalg import cg x3, _ = cg(A, f) print(np.mean(np.s