def test_kernelxx_numpy():
"polyprism.kernelxx cython vs numpy implementation"
for p in model:
cy = polyprism.kernelxx(xp, yp, zp, p)
py = _polyprism_numpy.kernelxx(xp, yp, zp, p)
diff = np.abs(py - cy)
assert np.all(diff <= precision), 'max diff: %g' % (max(diff))
def test_kernelxx_numpy():
"polyprism.kernelxx cython vs numpy implementation"
for p in model:
cy = polyprism.kernelxx(xp, yp, zp, p)
py = _polyprism_numpy.kernelxx(xp, yp, zp, p)
diff = np.abs(py - cy)
assert np.all(diff <= precision), 'max diff: %g' % (max(diff))
def test_kernelxx():
"polyprism.kernelxx against prism"
for pris, poly in zip(prismmodel, model):
resprism = prism.kernelxx(xp, yp, zp, pris)
respoly = polyprism.kernelxx(xp, yp, zp, poly)
diff = np.abs(resprism - respoly)
assert np.all(diff <= precision), \
'max diff: %g' % (max(diff))
def test_kernelxx():
"polyprism.kernelxx against prism"
for pris, poly in zip(prismmodel, model):
resprism = prism.kernelxx(xp, yp, zp, pris)
respoly = polyprism.kernelxx(xp, yp, zp, poly)
diff = np.abs(resprism - respoly)
assert np.all(diff <= precision), \
'max diff: %g' % (max(diff))
