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