docs: def smeft_evolve(C_in, scale_in, scale_out, newphys=True, **kwargs):
      """Solve the SMEFT RGEs by numeric integration.
  
      Input C_in and output C_out are dictionaries of arrays."""
      sol = _smeft_evolve(C_in, scale_in, scale_out, newphys=newphys, **kwargs)
      return C_array2dict(sol.y[:, -1].view(complex))