Numerical stochastic perturbation theory in the Schrödinger functional

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.