A Roughness Penalty Based on the Squared Second Derivative
A roughness penalty function based on squared second derivatives evaluated numerically. This is a possible template function for other types of penalties.
vector of length m, must be increasing
a function of one numeric variable
a vector of second derivatives evaluated at the points of xgrid