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))