docs: def symmetrize_JMS_dict(C):
      """For a dictionary with JMS Wilson coefficients but keys that might not be
      in the non-redundant basis, return a dictionary with keys from the basis
      and values conjugated if necessary."""
      wc_keys = set(wcxf.Basis['WET', 'JMS'].all_wcs)
      Cs = {}
      for op, v in C.items():
          if '_' not in op or op in wc_keys:
              Cs[op] = v
              continue
          name, ind = op.split('_')
          if name in C_symm_keys[5]:
              i, j, k, l = ind
              indnew = ''.join([j, i, l, k])
              Cs['_'.join([name, indnew])] = v.conjugate()
          elif name in C_symm_keys[41]:
              i, j, k, l = ind
              indnew = ''.join([k, l, i, j])
              Cs['_'.join([name, indnew])] = v
          elif name in C_symm_keys[4]:
              i, j, k, l = ind
              indnew = ''.join([l, k, j, i])
              newname = '_'.join([name, indnew])
              if newname in wc_keys:
                  Cs[newname] = v.conjugate()
              else:
                  indnew = ''.join([j, i, l, k])
                  newname = '_'.join([name, indnew])
                  if newname in wc_keys:
                      Cs[newname] = v.conjugate()
                  else:
                      indnew = ''.join([k, l, i, j])
                      newname = '_'.join([name, indnew])
                      Cs[newname] = v
          elif name in C_symm_keys[9]:
              i, j, k, l = ind
              indnew = ''.join([j, i, k, l])
              Cs['_'.join([name, indnew])] = -v
      return Cs