This function computes regularization penalty term Hessian based on the betas and differencing order.
regfun_hess(beta, regularization_order = 2)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.
scalar regularization value