def allclose_complex_numbers(a, b, atol=0.001):
    np.testing.assert_equal(np.shape(a), np.shape(b))
    np.testing.assert_allclose(np.abs(a), np.abs(b), rtol=1e-05, atol=atol)
    np.testing.assert_allclose(np.real(a), np.real(b), rtol=1e-05, atol=atol)
    np.testing.assert_allclose(np.imag(a), np.imag(b), rtol=1e-05, atol=atol)
def allclose_complex_numbers(a, b, atol=0.001):
    np.testing.assert_equal(np.shape(a), np.shape(b))
    np.testing.assert_allclose(np.abs(a), np.abs(b), rtol=1e-05, atol=atol)
    np.testing.assert_allclose(np.real(a), np.real(b), rtol=1e-05, atol=atol)
    np.testing.assert_allclose(np.imag(a), np.imag(b), rtol=1e-05, atol=atol)