алгоритм — интерполяция с B-сплайнами различной степени в переполнении стека

Я пытаюсь интерполировать через набор точек данных движения, используя кусочные движения, как определено Вагнером, поэтому каждая точка данных содержит 8 чисел. Каждое из этих чисел формирует параметр движения, поэтому мне нужно сгенерировать функции Bspline для каждого параметра (всего 8 параметров, 8 функций Bpsline). 3 числа относятся к переводу, 4 к повороту и 1 к весу перевода. Каждая из этих трех групп определена, чтобы иметь определенную степень для соответствующих кривых, поэтому общая кубическая B-сплайн-интерполяция не работает.

Формула для кубической линии Bspline определена явно, но это не относится к более высоким степеням кривой, так как я могу интерполировать набор данных, используя степень, скажем, 4? Я знаю, что можно использовать декомпозицию LU, но я не уверен, как реализовать это в C ++ или как определить необходимые значения для декомпозиции (у меня есть вектор узла и время каждой точки данных, но это все).