Stable direct estimation for GPLSIAMs using P-splines with dynamically updated boundaries | Steady Practice | SteadyPractice