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)
