checker('A2 = J', A2, J, x)
checker('exterior Proj.', Ce2, Ce, x)
checker('interior Proj.', Ci2, Ci, x)
checker('error-Calderon with random [no-sense]', A, J, x)
def dir_data(x, normal, dom_ind, result):
result[0] = -np.exp(1j * kRef * x[1])
def neu_data(x, normal, dom_ind, result):
result[0] = -1j * normal[1] * kRef * np.exp(1j * kRef * x[1])
b = mtf.rhs(dir_data, neu_data)
M = A - X
print('')
print(mtf.shape, flush=True)
print('')
#################################################
#################################################
#################################################
from krylov import gmres, bicgstab
#################################################
#################################################
#################################################
checker('interior Proj.', Ci2, Ci, x)
checker('error-Calderon with random [no-sense]', A, J, x)
#################################################
#################################################
def dir_data(x, normal, dom_ind, result):
result[0] = -np.exp( 1j * kRef * x[1])
def neu_data(x, normal, dom_ind, result):
result[0] = -1j * normal[1] * kRef * np.exp( 1j * kRef * x[1])
#################################################
#################################################
b = mtf.rhs(dir_data, neu_data)
M = A - X
print('')
print(mtf.shape, flush=True)
print('')
#################################################
#################################################
#################################################
iA = iJ * A * iJ
#################################################
Pjac = iA
