This function computes regularization penalty term based on the betas and a difference.
regfun(beta, regularization_order = 2)A spline parameter vector length num_params.
An integer (typically 0, 1, 2), indicating differencing order for L2 regularization of spline parameters. Default is 2 for second derivative penalty.
A scalar regularization value.