This function computes regularization penalty term gradient based on the betas and difference order.
regfun_grad(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