computeHermiteInterpolation.m
041b8992f8b7master

computeHermiteInterpolation.m
View Options

			function [ output ] = computeHermiteInterpolation(herm_c, herm_d, X, H, Npts)
			% computeHermiteInterpolation(herm_c, herm_d, X, H, Npts) evaluates the
			% cubic Hermite spline interpolation at points X. Interpolation
			% coefficients are given by herm_c and herm_d.

			output = zeros(size(X));
			precision = 12;

			leftIndices = floor(round(X * (Npts - 1), precision)) + 1;
			leftShifts = H * (leftIndices - 1);
			[~, upperBoundary] = find(leftIndices == Npts);

			output = output + herm_c(leftIndices) .* round(cubicHermite1(X, leftShifts, H), precision);
			output = output + herm_d(leftIndices) .* round(cubicHermite2(X, leftShifts, H), precision);

			% rightIndices = ceil(round(X *(Npts - 1), precision)) + 1;
			rightIndices = leftIndices + 1;
			rightShifts = H * (rightIndices -1);

			% Trick for x-values corresponding to upper boundary b. Has to be done
			% after defining rightShifts in order to get the correct solution.
			rightIndices(upperBoundary) = Npts;

			output = output + herm_c(rightIndices) .* round(cubicHermite1(X, rightShifts, H), precision);
			output = output + herm_d(rightIndices) .* round(cubicHermite2(X, rightShifts, H), precision);

			end