---------------------------------------------------------------------------
IndexError                                Traceback (most recent call last)
File ~/Desktop/04_github_repos/mne-opm/examples/reproduce_bug.py:59
     56 sim = SourceSimulator(src, tstep=1. / info['sfreq'],
     57                       duration=duration)
     58 sim.add_data(label_dipole, source_time_series, events)
---> 59 raw = mne.simulation.simulate_raw(
     60         info, sim, forward=None, src=src,
     61         head_pos=trans_matrices, bem=bem, trans=trans
     62 )

File <decorator-gen-408>:12, in simulate_raw(info, stc, trans, src, bem, head_pos, mindist, interp, n_jobs, use_cps, forward, first_samp, max_iter, verbose)

File ~/anaconda3/envs/mne_workshop/lib/python3.12/site-packages/mne/simulation/raw.py:364, in simulate_raw(info, stc, trans, src, bem, head_pos, mindist, interp, n_jobs, use_cps, forward, first_samp, max_iter, verbose)
    362 raw_datas.append(this_data)
    363 # Stim channel
--> 364 fwd, fi = interper.feed(this_stop - this_start)
    365 fi = fi[0]
    366 stc_data, stim_data, _ = _stc_data_event(
    367     stc_counted, fi, info["sfreq"], get_fwd.src, None if n == 0 else verts
    368 )

File ~/anaconda3/envs/mne_workshop/lib/python3.12/site-packages/mne/_ola.py:190, in _Interp2.feed(self, n_pts)
    188 # Convenience function for assembly
    189 out_arrays = None
--> 190 for o in self.feed_generator(n_pts):
    191     if out_arrays is None:
    192         out_arrays = [
    193             np.empty(v.shape + (n_pts,)) if v is not None else None
    194             for v in o[1]
    195         ]

File ~/anaconda3/envs/mne_workshop/lib/python3.12/site-packages/mne/_ola.py:131, in _Interp2.feed_generator(self, n_pts)
    129     eval_pt = self.control_points[self._left_idx + 1]
    130     logger.debug(f"  Eval @ {self._left_idx + 1} ({eval_pt})")
--> 131     self._right = self.values(eval_pt)
    132 assert self._right is not None
    133 left_point = self.control_points[self._left_idx]

File ~/anaconda3/envs/mne_workshop/lib/python3.12/site-packages/mne/simulation/raw.py:770, in _SimForwards.__call__(self, offset)
    768 assert self.offsets[self.idx] == offset
    769 self.idx += 1
--> 770 fwd = next(self.iter)
    771 self.src = fwd["src"]
    772 # XXX eventually we could speed this up by allowing the forward
    773 # solution code to only compute the normal direction

File ~/anaconda3/envs/mne_workshop/lib/python3.12/site-packages/mne/simulation/raw.py:863, in _iter_forward_solutions(***failed resolving arguments***)
    858     if not outside.all():
    859         raise RuntimeError(
    860             f"{np.sum(~outside)} MEG sensors collided with inner skull "
    861             f"surface for transform {ti}"
    862         )
--> 863     megfwd = _compute_forwards(
    864         rr, sensors=sensors, bem=bem, n_jobs=n_jobs, verbose=False
    865     )["meg"]
    866     megfwd = _to_forward_dict(megfwd, megnames)
    867 else:

File <decorator-gen-354>:10, in _compute_forwards(rr, bem, sensors, n_jobs, verbose)

File ~/anaconda3/envs/mne_workshop/lib/python3.12/site-packages/mne/forward/_compute_forward.py:838, in _compute_forwards(rr, bem, sensors, n_jobs, verbose)
    836 _check_option("solver", solver, ("mne", "openmeeg"))
    837 if bem["is_sphere"] or solver == "mne":
--> 838     fwd_data = _prep_field_computation(rr, sensors=sensors, bem=bem, n_jobs=n_jobs)
    839     Bs = _compute_forwards_meeg(
    840         rr, sensors=sensors, fwd_data=fwd_data, n_jobs=n_jobs
    841     )
    842 else:

File <decorator-gen-353>:12, in _prep_field_computation(rr, sensors, bem, n_jobs, verbose)

File ~/anaconda3/envs/mne_workshop/lib/python3.12/site-packages/mne/forward/_compute_forward.py:755, in _prep_field_computation(rr, sensors, bem, n_jobs, verbose)
    753     cf = FIFF.FIFFV_COORD_HEAD
    754     # multiply solution by "mults" here for simplicity
--> 755     solution = _bem_specify_coils(bem, coils, cf, mults, n_jobs)
    756 else:
    757     # Compute solution for EEG sensor
    758     logger.info("Setting up for EEG...")

File ~/anaconda3/envs/mne_workshop/lib/python3.12/site-packages/mne/forward/_compute_forward.py:199, in _bem_specify_coils(bem, coils, coord_frame, mults, n_jobs)
    179 """Set up for computing the solution at a set of MEG coils.
    180 
    181 Parameters
   (...)
    196     MEG solution
    197 """
    198 # Make sure MEG coils are in MRI coordinate frame to match BEM coords
--> 199 coils, coord_frame = _check_coil_frame(coils, coord_frame, bem)
    201 # leaving this in in case we want to easily add in the future
    202 # if method != 'simple':  # in ['ferguson', 'urankar']:
    203 #     raise NotImplementedError
   (...)
    207 
    208 # Process each of the surfaces
    209 rmags, cosmags, ws, bins = _triage_coils(coils)

File ~/anaconda3/envs/mne_workshop/lib/python3.12/site-packages/mne/forward/_compute_forward.py:53, in _check_coil_frame(coils, coord_frame, bem)
     50 if coord_frame != FIFF.FIFFV_COORD_MRI:
     51     if coord_frame == FIFF.FIFFV_COORD_HEAD:
     52         # Make a transformed duplicate
---> 53         coils, coord_frame = _dup_coil_set(coils, coord_frame, bem["head_mri_t"])
     54     else:
     55         raise RuntimeError(f"Bad coil coordinate frame {coord_frame}")

File ~/anaconda3/envs/mne_workshop/lib/python3.12/site-packages/mne/forward/_compute_forward.py:41, in _dup_coil_set(coils, coord_frame, t)
     39     if key in coil:
     40         coil[key] = apply_trans(t["trans"], coil[key], False)
---> 41 coil["r0"] = apply_trans(t["trans"], coil["r0"])
     42 coil["rmag"] = apply_trans(t["trans"], coil["rmag"])
     43 coil["cosmag"] = apply_trans(t["trans"], coil["cosmag"], False)

IndexError: only integers, slices (`:`), ellipsis (`...`), numpy.newaxis (`None`) and integer or boolean arrays are valid indices