def test_compute_N_04(self):
k = 10.0
p = np.array([0.5, 0.75, 1.0], dtype=np.float32)
normal_p = np.array([-np.sqrt(0.5), -np.sqrt(0.5), 0.0], dtype=np.float32)
a = np.array([0.0, 0.0, 0.0], dtype=np.float32)
b = np.array([0.1, 0.0, 0.0], dtype=np.float32)
c = np.array([0.0, 0.1, 0.0], dtype=np.float32)
pOnElement = True
zP = iops_pyx.n_3d(k, p, normal_p, a, b, c, pOnElement)
zC = iops_cpp.n_3d(k, p, normal_p, a, b, c, pOnElement)
self.assertAlmostEqual(zP, zC, 5)
def test_compute_N_02(self):
k = 10.0
p = np.array([0.5, 0.75, 1.0], dtype=np.float32)
normal_p = np.array([-np.sqrt(0.5), -np.sqrt(0.5), 0.0],
dtype=np.float32)
a = np.array([0.0, 0.0, 0.0], dtype=np.float32)
b = np.array([0.1, 0.0, 0.0], dtype=np.float32)
c = np.array([0.0, 0.1, 0.0], dtype=np.float32)
pOnElement = False
zP = iops_pyx.n_3d(k, self.p_off, self.n_p_off, self.a, self.b, c,
pOnElement)
zC = iops_cpp.n_3d(k, self.p_off, self.n_p_off, self.a, self.b, c,
pOnElement)
self.assertAlmostEqual(zP, zC)