  Thanks to dimensional reduction, the infrared contributions to the QCD
pressure can be obtained from two different three-dimensional effective field
theories, called the Electrostatic QCD (Yang-Mills plus adjoint Higgs) and the
Magnetostatic QCD (pure Yang-Mills theory). Lattice measurements have been
carried out within these theories, but a proper interpretation of the results
requires renormalization, and in some cases also improvement, i.e. the removal
of terms of O(a) or O(a^2). We discuss how these computations can be
implemented and carried out up to 4-loop level with the help of Numerical
Stochastic Perturbation Theory.
