gps2GS function

Penalized B-splines estimation with automatic grid search of their smoothing parameter