def test_gxy(): "gravmag.prism.gxy python vs cython implementation" py = _prism.gxy(xp, yp, zp, model) cy = _cprism.gxy(xp, yp, zp, model) diff = np.abs(py - cy) assert np.all(diff <= precision), 'max diff: %g' % (max(diff))
def test_gxy_ne(): "gravmag.prism.gxy python vs numexpr implementation" py = _prism.gxy(xp, yp, zp, model) ne = _neprism.gxy(xp, yp, zp, model) diff = np.abs(py - ne) assert np.all(diff <= precision), 'max diff: %g' % (max(diff))