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