def test_InverseAdjointOperator(operator_with_arrays):
op, mu, U, V = operator_with_arrays
if not op.linear:
return
inv = InverseAdjointOperator(op)
rtol = atol = 1e-12
try:
assert np.all(almost_equal(inv.apply(U, mu=mu), op.apply_inverse_adjoint(U, mu=mu),
rtol=rtol, atol=atol))
except (InversionError, LinAlgError, NotImplementedError):
pass
try:
assert np.all(almost_equal(inv.apply_inverse(V, mu=mu), op.apply_adjoint(V, mu=mu),
rtol=rtol, atol=atol))
except (InversionError, LinAlgError, NotImplementedError):
pass
try:
assert np.all(almost_equal(inv.apply_adjoint(V, mu=mu), op.apply_inverse(V, mu=mu),
rtol=rtol, atol=atol))
except (InversionError, LinAlgError, NotImplementedError):
pass
try:
assert np.all(almost_equal(inv.apply_inverse_adjoint(U, mu=mu), op.apply(U, mu=mu),
rtol=rtol, atol=atol))
except (InversionError, LinAlgError, NotImplementedError):
pass
def test_InverseAdjointOperator(operator_with_arrays):
op, mu, U, V = operator_with_arrays
if not op.linear:
return
inv = InverseAdjointOperator(op)
rtol = atol = 1e-12
try:
assert np.all(almost_equal(inv.apply(U, mu=mu), op.apply_inverse_adjoint(U, mu=mu),
rtol=rtol, atol=atol))
except (InversionError, LinAlgError, NotImplementedError):
pass
try:
assert np.all(almost_equal(inv.apply_inverse(V, mu=mu), op.apply_adjoint(V, mu=mu),
rtol=rtol, atol=atol))
except (InversionError, LinAlgError, NotImplementedError):
pass
try:
assert np.all(almost_equal(inv.apply_adjoint(V, mu=mu), op.apply_inverse(V, mu=mu),
rtol=rtol, atol=atol))
except (InversionError, LinAlgError, NotImplementedError):
pass
try:
assert np.all(almost_equal(inv.apply_inverse_adjoint(U, mu=mu), op.apply(U, mu=mu),
rtol=rtol, atol=atol))
except (InversionError, LinAlgError, NotImplementedError):
pass