def test_mul_inverse(self): ga = GeometricAlgebra(metric=dual_metric) # a = 2 a = ga.fill([], fill_value=2.0, kind="scalar") # b = 3 + 3e0 b = ga.fill([], fill_value=3.0, kind="mv") # a * b = 2 * (3 + 3e0) = 6 + 6e0 c = a * b self.assertEqual(c, 6.0 + 6.0 * ga.basis_mvs[0]) # a^-1 = 1 / 2 a_inv = a.inverse() self.assertEqual(a_inv, 0.5) # c = a * b # => a_inv * c = b self.assertEqual(a_inv * c, b) # Since a is scalar, should commute too. # => c * a_inv = b self.assertEqual(c * a_inv, b) # b is not invertible and will throw an exception self.assertRaises(Exception, b.inverse)
def test_mul_tf_mv(self): ga = GeometricAlgebra(metric=dual_metric) zero = ga.zeros([], kind="scalar") zero_tf = ga.as_mv(tf.convert_to_tensor(0.0, dtype=tf.float32)) one = ga.ones([], kind="scalar") one_tf = ga.as_mv(tf.convert_to_tensor(1.0, dtype=tf.float32)) eps = ga.ones([], kind="pseudoscalar") ten = ga.fill([], fill_value=10.0, kind="scalar") ten_tf = ga.as_mv(tf.convert_to_tensor([10.0], dtype=tf.float32)) self.assertEqual(one * one_tf, one) self.assertEqual(one_tf * one, one) self.assertEqual(zero * one_tf, zero) self.assertEqual(one_tf * zero, zero) self.assertEqual(zero_tf * one, zero) self.assertEqual(one * zero_tf, zero) self.assertEqual(one_tf * eps, eps) self.assertEqual(eps * one_tf, eps) self.assertEqual(zero_tf * zero, zero) self.assertEqual(zero * zero_tf, zero) self.assertEqual(ten_tf * zero, zero) self.assertEqual(zero * ten_tf, zero) self.assertEqual(ten * zero_tf, zero) self.assertEqual(zero_tf * ten, zero) self.assertEqual((ten_tf * eps) * eps, zero) self.assertEqual(ten_tf * one, ten) self.assertEqual(one * ten_tf, ten) self.assertEqual(ten * one_tf, ten) self.assertEqual(one_tf * ten, ten)
def test_mul_py_mv(self): ga = GeometricAlgebra(metric=dual_metric) zero = ga.zeros([], kind="scalar") zero_py = 0.0 one = ga.ones([], kind="scalar") one_py = 1.0 eps = ga.ones([], kind="pseudoscalar") ten = ga.fill([], fill_value=10.0, kind="scalar") ten_py = 10.0 self.assertEqual(one * one_py, one) self.assertEqual(one_py * one, one) self.assertEqual(zero * one_py, zero) self.assertEqual(one_py * zero, zero) self.assertEqual(zero_py * one, zero) self.assertEqual(one * zero_py, zero) self.assertEqual(one_py * eps, eps) self.assertEqual(eps * one_py, eps) self.assertEqual(zero_py * zero, zero) self.assertEqual(zero * zero_py, zero) self.assertEqual(ten_py * zero, zero) self.assertEqual(zero * ten_py, zero) self.assertEqual(ten * zero_py, zero) self.assertEqual(zero_py * ten, zero) self.assertEqual((ten_py * eps) * eps, zero) self.assertEqual(ten_py * one, ten) self.assertEqual(one * ten_py, ten) self.assertEqual(ten * one_py, ten) self.assertEqual(one_py * ten, ten)