def test_kernelxz_numpy():
    "polyprism.kernelxz cython vs numpy implementation"
    for p in model:
        cy = polyprism.kernelxz(xp, yp, zp, p)
        py = _polyprism_numpy.kernelxz(xp, yp, zp, p)
        diff = np.abs(py - cy)
        assert np.all(diff <= precision), 'max diff: %g' % (max(diff))
Example #2
0
def test_kernelxz_numpy():
    "polyprism.kernelxz cython vs numpy implementation"
    for p in model:
        cy = polyprism.kernelxz(xp, yp, zp, p)
        py = _polyprism_numpy.kernelxz(xp, yp, zp, p)
        diff = np.abs(py - cy)
        assert np.all(diff <= precision), 'max diff: %g' % (max(diff))
def test_kernelxz():
    "polyprism.kernelxz against prism"
    for pris, poly in zip(prismmodel, model):
        resprism = prism.kernelxz(xp, yp, zp, pris)
        respoly = polyprism.kernelxz(xp, yp, zp, poly)
        diff = np.abs(resprism - respoly)
        assert np.all(diff <= precision), \
            'max diff: %g' % (max(diff))
Example #4
0
def test_kernelxz():
    "polyprism.kernelxz against prism"
    for pris, poly in zip(prismmodel, model):
        resprism = prism.kernelxz(xp, yp, zp, pris)
        respoly = polyprism.kernelxz(xp, yp, zp, poly)
        diff = np.abs(resprism - respoly)
        assert np.all(diff <= precision), \
            'max diff: %g' % (max(diff))