To construct a multivariate B-spline, we simply take the Kronecker tensor product between the univariate basis functions constructed for each individual dimension.
What I'd like to know is how do you construct a B-spline derivative for a multivariate function? Is it the same method of using the tensor product, but with the univariate B-spline derivatives?