Exemplo n.º 1
0
def test_potential():
    "gravmag.prism.potential python vs cython implementation"
    py = _prism.potential(xp, yp, zp, model)
    cy = _cprism.potential(xp, yp, zp, model)
    diff = np.abs(py - cy)
    assert np.all(diff <= precision), 'max diff: %g' % (max(diff))
Exemplo n.º 2
0
def test_potential():
    "gravmag.prism.potential python vs cython implementation"
    py = _prism.potential(xp, yp, zp, model)
    cy = _cprism.potential(xp, yp, zp, model)
    diff = np.abs(py - cy)
    assert np.all(diff <= precision), 'max diff: %g' % (max(diff))
Exemplo n.º 3
0
def test_potential_ne():
    "gravmag.prism.potential python vs numexpr implementation"
    py = _prism.potential(xp, yp, zp, model)
    ne = _neprism.potential(xp, yp, zp, model)
    diff = np.abs(py - ne)
    assert np.all(diff <= precision), 'max diff: %g' % (max(diff))
Exemplo n.º 4
0
def test_potential_ne():
    "gravmag.prism.potential python vs numexpr implementation"
    py = _prism.potential(xp, yp, zp, model)
    ne = _neprism.potential(xp, yp, zp, model)
    diff = np.abs(py - ne)
    assert np.all(diff <= precision), 'max diff: %g' % (max(diff))