def test_dB(self): bs = BSpline() x = np.linspace(bs.xmin, bs.xmax, 1000) N = len(x) M = bs.deg + bs.nknots + 1 for nu in xrange(1, 4): B = np.zeros((N,M)) for i in xrange(-bs.deg, bs.nknots + 1): for n in xrange(N): B[n,ix(i)] = bs.dB(nu, i, bs.deg + 1, x[n]) pass