docs: def _smeft_evolve(C_in, scale_in, scale_out, newphys=True, **kwargs):
      """Axuliary function used in `smeft_evolve` and `smeft_evolve_continuous`"""
      def fun(t0, y):
          return beta.beta_array(C=C_array2dict(y.view(complex)),
                                 newphys=newphys).view(float) / (16 * pi**2)
      y0 = C_dict2array(C_in).view(float)
      sol = solve_ivp(fun=fun,
                      t_span=(log(scale_in), log(scale_out)),
                      y0=y0, **kwargs)
      return sol