Added some comments to evalJacobian and evalHessian. If evaluated
outside the B-spline domain they now return a zero matrix.