pdSplineReg function

Cubic smoothing spline regression for HPD matrices