Example #1
0
	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