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