################################################# ################################################# from krylov import gmres, bicgstab ################################################# ################################################# ################################################# iA = iJ * A * iJ ################################################# Pjac = iA E = mtf.upper() Pgs = iA + iA * E * iA Msigma = lambda sigma: (A - J) + sigma * (J - X) CCi = iJ * (0.5 * J + At) CCe = (0.5 * J - At) B = J - X BD = B * CCi + CCe F = B * CCi ################################################# ################################################# #################################################
#################################################
#################################################
b = mtf.rhs(dir_data, neu_data)
M = A - X
print('')
print(mtf.shape, flush=True)
print('')
#################################################
#################################################
#################################################
iA = iJ * A * iJ
#################################################
Pjac = iA
E = mtf.upper()
Pgs = iA + iA * E * iA
CCi = iJ * (0.5 * J + At)
CCe = (0.5 * J - At)
B = J - X
BD = B * CCi + CCe
F = B * CCi
