docs: def smeft_toarray(wc_name, wc_dict):
      """Construct a numpy array with Wilson coefficient values from a
      dictionary of label-value pairs corresponding to the non-redundant
      elements."""
      shape = smeftutil.C_keys_shape[wc_name]
      C = np.zeros(shape, dtype=complex)
      for k, v in wc_dict.items():
          if k.split('_')[0] != wc_name:
              continue
          indices = k.split('_')[-1]  # e.g. '1213'
          indices = tuple(int(s) - 1 for s in indices)  # e.g. (1, 2, 1, 3)
          C[indices] = v
      C = smeftutil.symmetrize({wc_name: C})[wc_name]
      return C