docs: def smeft_evolve_continuous(C_in, scale_in, scale_out, newphys=True, **kwargs):
      """Solve the SMEFT RGEs by numeric integration, returning a function that
      allows to compute an interpolated solution at arbitrary intermediate
      scales."""
      sol = _smeft_evolve(C_in, scale_in, scale_out, newphys=newphys,
                          dense_output=True, **kwargs)
      @np.vectorize
      def _rge_solution(scale):
          t = log(scale)
          y = sol.sol(t).view(complex)
          yd = C_array2dict(y)
          yw = arrays2wcxf_nonred(yd)
          return yw
      def rge_solution(scale):
          # this is to return a scalar if the input is scalar
          return _rge_solution(scale)[()]
      return rge_solution