  The Schr\"odinger functional (SF) is a powerful and widely used tool for the
treatment of a variety of problems in renormalization and related areas. Albeit
offering many conceptual advantages, one major downside of the SF scheme is the
fact that perturbative calculations quickly become cumbersome with the
inclusion of higher orders in the gauge coupling and hence the use of an
automated perturbation theory framework is desirable. We present the
implementation of the SF in numerical stochastic perturbation theory (NSPT) and
compare first results for the running coupling at two loops in pure SU(3)
Yang-Mills theory with the literature.
