import numpy as np
from scipy.sparse import csc_matrix
M = csc_matrix(np.random.rand(5, 5))
B = np.random.rand(5, 4)
X, info = bicg_solve(M, B)
print(X.shape)
Итеративное решение разреженных систем линейных уравнений с (M, N) правой матрицей размера