def test_mod_c(): z1 = Bicomplex(np.linspace(0, np.pi, 5), 0) z2 = z1.mod_c() assert_array_equal(z2, np.sqrt(z1.z1**2 + z1.z2**2)) z3 = Bicomplex(0.1, np.linspace(0, np.pi, 5)) z4 = z3.mod_c() trueval = np.sqrt(z3*z3.conjugate()) assert_allclose(z4, np.sqrt(z3.z1**2 + z3.z2**2)) assert_allclose(z4, trueval.z1)