General scalar renormalisation group equations at three-loop order

For arbitrary scalar QFTs in four dimensions, renormalisation group equations of quartic and cubic interactions, mass terms, as well as field anomalous dimensions are computed at three-loop order in the $overline{text{MS}}$ scheme. Utilising pre-existing literature expressions for a specific model, loop integrals are avoided and templates for general theories are obtained. We reiterate known four-loop expressions, and derive $beta$ functions for scalar masses and cubic interactions from it. As an example, the results are applied to compute all renormalisation group equations in $U(n) times U(n)$ scalar theories.