def test_bx_numpy():
"polyprism.bx cython vs numpy implementation"
cy = polyprism.bx(xp, yp, zp, model)
py = _polyprism_numpy.bx(xp, yp, zp, model)
diff = np.abs(py - cy)
assert np.all(diff <= precision_mag), 'max diff: %g' % (max(diff))
def test_bx_numpy():
"polyprism.bx cython vs numpy implementation"
cy = polyprism.bx(xp, yp, zp, model)
py = _polyprism_numpy.bx(xp, yp, zp, model)
diff = np.abs(py - cy)
assert np.all(diff <= precision_mag), 'max diff: %g' % (max(diff))
def test_bx():
"polyprism.bx against prism"
resprism = prism.bx(xp, yp, zp, prismmodel)
respoly = polyprism.bx(xp, yp, zp, model)
diff = np.abs(resprism - respoly)
assert np.all(diff <= precision_mag), 'max diff: %g' % (max(diff))
def test_bx():
"polyprism.bx against prism"
resprism = prism.bx(xp, yp, zp, prismmodel)
respoly = polyprism.bx(xp, yp, zp, model)
diff = np.abs(resprism - respoly)
assert np.all(diff <= precision_mag), 'max diff: %g' % (max(diff))
