docs: def beta_array(C, HIGHSCALE, *args, **kwargs):
      """Return the beta functions of all SM parameters and SMEFT Wilson
      coefficients as a 1D numpy array."""
      beta_odict = beta(C, HIGHSCALE, *args, **kwargs)
      return np.hstack([np.asarray(b).ravel() for b in beta_odict.values()])