Say in the signatures that epochs data may arrive per epoch

BaseEpochs and EpochsArray both accept a list holding one array per
epoch, and tmin may be given per event, but the annotations still
described only the rectangular case, so passing what the docstrings
document was reported as a type error.

The fixed-duration branch of BaseEpochs.__init__ asserts the list is not
present before reshaping, which is the branch it already relies on.

Author: snesmaeili

mne/epochs.py

@@ -579,9 +579,11 @@ class BaseEpochs(
     Parameters
     ----------
     %(info_not_none)s
-    data : ndarray | None
+    data : ndarray | list of ndarray | None
         If ``None``, data will be read from the Raw object. If ndarray, must be
-        of shape (n_epochs, n_channels, n_times).
+        of shape (n_epochs, n_channels, n_times). A list of
+        (n_channels, n_times) arrays, one per epoch, gives epochs of differing
+        duration.
     %(events_epochs)s
     %(event_id)s
     %(epochs_tmin_tmax)s
@@ -639,7 +641,7 @@ class BaseEpochs(
     def __init__(
         self,
         info: Info,
-        data: np.ndarray | None,
+        data: np.ndarray | list[np.ndarray] | None,
         events: np.ndarray,
         event_id: int | list[int] | dict | str | list[str] | None = None,
         tmin: float | np.ndarray = -0.2,
@@ -807,6 +809,7 @@ def __init__(
             self._do_baseline = False
         else:
             assert decim == 1
+            assert not isinstance(data, list)  # only the branch above takes a list
             if (
                 data.ndim != 3
                 or data.shape[2] != round((tmax - tmin) * self.info["sfreq"]) + 1