diff --git a/pybamm/expression_tree/concatenations.py b/pybamm/expression_tree/concatenations.py
index 52549b374..1bef73f75 100644
--- a/pybamm/expression_tree/concatenations.py
+++ b/pybamm/expression_tree/concatenations.py
@@ -269,22 +269,6 @@ class SparseStack(Concatenation):
         children = list(children)
         super().__init__(*children, name="sparse stack", check_domain=False)
 
-    def evaluate(self, t=None, y=None, known_evals=None):
-        """ See :meth:`pybamm.Symbol.evaluate()`. """
-        children = self.cached_children
-        if known_evals is not None:
-            if self.id not in known_evals:
-                children_eval = [None] * len(children)
-                for idx, child in enumerate(children):
-                    children_eval[idx], known_evals = child.evaluate(t, y, known_evals)
-                known_evals[self.id] = self._concatenation_evaluate(children_eval)
-            return known_evals[self.id], known_evals
-        else:
-            children_eval = [None] * len(children)
-            for idx, child in enumerate(children):
-                children_eval[idx] = child.evaluate(t, y)
-            return self._concatenation_evaluate(children_eval)
-
     def _concatenation_evaluate(self, children_eval):
         """ See :meth:`Concatenation.evaluate()`. """
         return vstack(children_eval)
diff --git a/pybamm/solvers/base_solver.py b/pybamm/solvers/base_solver.py
index 1cadf7db6..18b2ad446 100644
--- a/pybamm/solvers/base_solver.py
+++ b/pybamm/solvers/base_solver.py
@@ -12,13 +12,22 @@ class BaseSolver(object):
         The tolerance for the solver (default is 1e-8).
     """
 
-    def __init__(self, tol=1e-8):
+    def __init__(self, method=None, tol=1e-8):
+        self._method = method
         self._tol = tol
 
         # Solutions
         self._t = None
         self._y = None
 
+    @property
+    def method(self):
+        return self._method
+
+    @method.setter
+    def method(self, value):
+        self._method = value
+
     @property
     def tol(self):
         return self._tol
diff --git a/pybamm/solvers/dae_solver.py b/pybamm/solvers/dae_solver.py
index 1d17ff7e1..95918aaef 100644
--- a/pybamm/solvers/dae_solver.py
+++ b/pybamm/solvers/dae_solver.py
@@ -19,8 +19,8 @@ class DaeSolver(pybamm.BaseSolver):
         The tolerance for the initial-condition solver (default is 1e-8).
     """
 
-    def __init__(self, tol=1e-8, root_method="lm", root_tol=1e-6):
-        super().__init__(tol)
+    def __init__(self, method=None, tol=1e-8, root_method="lm", root_tol=1e-6):
+        super().__init__(method, tol)
         self.root_method = root_method
         self.root_tol = root_tol
 
@@ -53,17 +53,7 @@ class DaeSolver(pybamm.BaseSolver):
 
         """
         pybamm.logger.info("Start solving {}".format(model.name))
-
-        # create simplified rhs algebraic and event expressions
-        concatenated_rhs = model.concatenated_rhs
-        concatenated_algebraic = model.concatenated_algebraic
-        events = model.events
-        if model.use_simplify:
-            # set up simplification object, for re-use of dict
-            simp = pybamm.Simplification()
-            concatenated_rhs = simp.simplify(concatenated_rhs)
-            concatenated_algebraic = simp.simplify(concatenated_algebraic)
-            events = [simp.simplify(event) for event in events]
+        concatenated_rhs, concatenated_algebraic, y0, events, jac = self.set_up(model)
 
         def residuals(t, y, ydot):
             pybamm.logger.debug(
@@ -78,12 +68,6 @@ class DaeSolver(pybamm.BaseSolver):
             alg_eval = alg_eval[:, 0]
             return np.concatenate((rhs_eval - ydot[: rhs_eval.shape[0]], alg_eval))
 
-        def rhs(t, y):
-            return concatenated_rhs.evaluate(t, y, known_evals={})[0][:, 0]
-
-        def algebraic(t, y):
-            return concatenated_algebraic.evaluate(t, y, known_evals={})[0][:, 0]
-
         # Create event-dependent function to evaluate events
         def event_fun(event):
             def eval_event(t, y):
@@ -93,20 +77,8 @@ class DaeSolver(pybamm.BaseSolver):
 
         events = [event_fun(event) for event in events]
 
-        y0 = self.calculate_consistent_initial_conditions(
-            rhs, algebraic, model.concatenated_initial_conditions[:, 0]
-        )
-
-        if model.use_jacobian:
-            # Create Jacobian from simplified rhs
-            y = pybamm.StateVector(slice(0, np.size(y0)))
-            jac_rhs = concatenated_rhs.jac(y)
-            jac_algebraic = concatenated_algebraic.jac(y)
-            if model.use_simplify:
-                jac_rhs = simp.simplify(jac_rhs)
-                jac_algebraic = simp.simplify(jac_algebraic)
-
-            jac = pybamm.SparseStack(jac_rhs, jac_algebraic)
+        # Create function to evaluate jacobian
+        if jac is not None:
 
             def jacobian(t, y):
                 return jac.evaluate(t, y, known_evals={})[0]
@@ -125,6 +97,77 @@ class DaeSolver(pybamm.BaseSolver):
 
         pybamm.logger.info("Finish solving {}".format(model.name))
 
+    def set_up(self, model):
+        """Unpack model, perform checks, simplify and calculate jacobian.
+
+        Parameters
+        ----------
+        model : :class:`pybamm.BaseModel`
+            The model whose solution to calculate. Must have attributes rhs and
+            initial_conditions
+
+        Returns
+        -------
+        concatenated_rhs : :class:`pybamm.Concatenation`
+            Right-hand side of differential equations
+        concatenated_algebraic : :class:`pybamm.Concatenation`
+            Algebraic equations, which should evaluate to zero
+        y0 : :class:`numpy.array`
+            Vector of initial conditions
+        events : list of :class:`pybamm.Symbol`
+            List of events at which the model should terminate
+        jac : :class:`pybamm.SparseStack`
+            Jacobian matrix for the differential and algebraic equations
+
+        Raises
+        ------
+        :class:`pybamm.SolverError`
+            If the model contains any algebraic equations (in which case a DAE solver
+            should be used instead)
+        """
+        # create simplified rhs algebraic and event expressions
+        concatenated_rhs = model.concatenated_rhs
+        concatenated_algebraic = model.concatenated_algebraic
+        events = model.events
+        if model.use_simplify:
+            # set up simplification object, for re-use of dict
+            simp = pybamm.Simplification()
+            concatenated_rhs = simp.simplify(concatenated_rhs)
+            concatenated_algebraic = simp.simplify(concatenated_algebraic)
+            events = [simp.simplify(event) for event in events]
+
+        # Calculate consistent initial conditions for the algebraic equations
+        def rhs(t, y):
+            return concatenated_rhs.evaluate(t, y, known_evals={})[0][:, 0]
+
+        def algebraic(t, y):
+            return concatenated_algebraic.evaluate(t, y, known_evals={})[0][:, 0]
+
+        if len(model.algebraic) > 0:
+            y0 = self.calculate_consistent_initial_conditions(
+                rhs, algebraic, model.concatenated_initial_conditions[:, 0]
+            )
+        else:
+            # can use DAE solver to solve ODE model
+            y0 = model.concatenated_initial_conditions[:, 0]
+
+        # Calculate jacobian
+        if model.use_jacobian:
+            # Create Jacobian from simplified rhs
+            y = pybamm.StateVector(slice(0, np.size(y0)))
+            jac_rhs = concatenated_rhs.jac(y)
+            jac_algebraic = concatenated_algebraic.jac(y)
+            if model.use_simplify:
+                jac_rhs = simp.simplify(jac_rhs)
+                jac_algebraic = simp.simplify(jac_algebraic)
+
+            jac = pybamm.SparseStack(jac_rhs, jac_algebraic)
+
+        else:
+            jac = None
+
+        return concatenated_rhs, concatenated_algebraic, y0, events, jac
+
     def calculate_consistent_initial_conditions(self, rhs, algebraic, y0_guess):
         """
         Calculate consistent initial conditions for the algebraic equations through
diff --git a/pybamm/solvers/ode_solver.py b/pybamm/solvers/ode_solver.py
index ee11d3971..cdc36dbc2 100644
--- a/pybamm/solvers/ode_solver.py
+++ b/pybamm/solvers/ode_solver.py
@@ -14,8 +14,8 @@ class OdeSolver(pybamm.BaseSolver):
         The tolerance for the solver (default is 1e-8).
     """
 
-    def __init__(self, tol=1e-8):
-        super().__init__(tol)
+    def __init__(self, method=None, tol=1e-8):
+        super().__init__(method, tol)
 
     def solve(self, model, t_eval):
         """Calculate the solution of the model at specified times.
@@ -30,16 +30,9 @@ class OdeSolver(pybamm.BaseSolver):
 
         """
         pybamm.logger.info("Start solving {}".format(model.name))
+        concatenated_rhs, y0, events, jac_rhs = self.set_up(model)
 
-        concatenated_rhs = model.concatenated_rhs
-        events = model.events
-        if model.use_simplify:
-            # set up simplification object, for re-use of dict
-            simp = pybamm.Simplification()
-            # create simplified rhs and event expressions
-            concatenated_rhs = simp.simplify(concatenated_rhs)
-            events = [simp.simplify(event) for event in events]
-
+        # Create function to evaluate rhs
         def dydt(t, y):
             pybamm.logger.debug("Evaluating RHS for {} at t={}".format(model.name, t))
             y = y[:, np.newaxis]
@@ -55,14 +48,8 @@ class OdeSolver(pybamm.BaseSolver):
 
         events = [event_fun(event) for event in events]
 
-        y0 = model.concatenated_initial_conditions[:, 0]
-
-        if model.use_jacobian:
-            # Create Jacobian from simplified rhs
-            y = pybamm.StateVector(slice(0, np.size(y0)))
-            jac_rhs = concatenated_rhs.jac(y)
-            if model.use_simplify:
-                jac_rhs = simp.simplify(jac_rhs)
+        # Create function to evaluate jacobian
+        if jac_rhs is not None:
 
             def jacobian(t, y):
                 return jac_rhs.evaluate(t, y, known_evals={})[0]
@@ -81,6 +68,61 @@ class OdeSolver(pybamm.BaseSolver):
 
         pybamm.logger.info("Finish solving {}".format(model.name))
 
+    def set_up(self, model):
+        """Unpack model, perform checks, simplify and calculate jacobian.
+
+        Parameters
+        ----------
+        model : :class:`pybamm.BaseModel`
+            The model whose solution to calculate. Must have attributes rhs and
+            initial_conditions
+
+        Returns
+        -------
+        concatenated_rhs : :class:`pybamm.Concatenation`
+            Right-hand side of differential equations
+        y0 : :class:`numpy.array`
+            Vector of initial conditions
+        events : list of :class:`pybamm.Symbol`
+            List of events at which the model should terminate
+        jac_rhs : :class:`pybamm.SparseStack`
+            Jacobian matrix for the differential equations
+
+        Raises
+        ------
+        :class:`pybamm.SolverError`
+            If the model contains any algebraic equations (in which case a DAE solver
+            should be used instead)
+
+        """
+        if len(model.algebraic) > 0:
+            raise pybamm.SolverError(
+                """Cannot use ODE solver to solve model with DAEs"""
+            )
+
+        concatenated_rhs = model.concatenated_rhs
+        events = model.events
+        if model.use_simplify:
+            # set up simplification object, for re-use of dict
+            simp = pybamm.Simplification()
+            # create simplified rhs and event expressions
+            concatenated_rhs = simp.simplify(concatenated_rhs)
+            events = [simp.simplify(event) for event in events]
+
+        y0 = model.concatenated_initial_conditions[:, 0]
+
+        if model.use_jacobian:
+            # Create Jacobian from simplified rhs
+            y = pybamm.StateVector(slice(0, np.size(y0)))
+            jac_rhs = concatenated_rhs.jac(y)
+            if model.use_simplify:
+                jac_rhs = simp.simplify(jac_rhs)
+
+        else:
+            jac_rhs = None
+
+        return concatenated_rhs, y0, events, jac_rhs
+
     def integrate(
         self, derivs, y0, t_eval, events=None, mass_matrix=None, jacobian=None
     ):
diff --git a/pybamm/solvers/scikits_dae_solver.py b/pybamm/solvers/scikits_dae_solver.py
index c92f8c52a..3aac5ef96 100644
--- a/pybamm/solvers/scikits_dae_solver.py
+++ b/pybamm/solvers/scikits_dae_solver.py
@@ -35,16 +35,7 @@ class ScikitsDaeSolver(pybamm.DaeSolver):
         if scikits_odes_spec is None:
             raise ImportError("scikits.odes is not installed")
 
-        super().__init__(tol, root_method, root_tol)
-        self._method = method
-
-    @property
-    def method(self):
-        return self._method
-
-    @method.setter
-    def method(self, value):
-        self._method = value
+        super().__init__(method, tol, root_method, root_tol)
 
     def integrate(
         self, residuals, y0, t_eval, events=None, mass_matrix=None, jacobian=None
diff --git a/pybamm/solvers/scikits_ode_solver.py b/pybamm/solvers/scikits_ode_solver.py
index 407f41b6a..c21b6cea3 100644
--- a/pybamm/solvers/scikits_ode_solver.py
+++ b/pybamm/solvers/scikits_ode_solver.py
@@ -33,18 +33,9 @@ class ScikitsOdeSolver(pybamm.OdeSolver):
         if scikits_odes_spec is None:
             raise ImportError("scikits.odes is not installed")
 
-        super().__init__(tol)
-        self._method = method
+        super().__init__(method, tol)
         self.linsolver = linsolver
 
-    @property
-    def method(self):
-        return self._method
-
-    @method.setter
-    def method(self, value):
-        self._method = value
-
     def integrate(
         self, derivs, y0, t_eval, events=None, mass_matrix=None, jacobian=None
     ):
diff --git a/pybamm/solvers/scipy_solver.py b/pybamm/solvers/scipy_solver.py
index 78312a567..3b1493076 100644
--- a/pybamm/solvers/scipy_solver.py
+++ b/pybamm/solvers/scipy_solver.py
@@ -20,16 +20,7 @@ class ScipySolver(pybamm.OdeSolver):
     """
 
     def __init__(self, method="BDF", tol=1e-8):
-        super().__init__(tol)
-        self.method = method
-
-    @property
-    def method(self):
-        return self._method
-
-    @method.setter
-    def method(self, value):
-        self._method = value
+        super().__init__(method, tol)
 
     def integrate(
         self, derivs, y0, t_eval, events=None, mass_matrix=None, jacobian=None

Solver check
Add checks to solver to make sure we don't try to solve an ODE model with a DAE solver, and vice versa
diff --git a/tests/unit/test_solvers/test_ode_solver.py b/tests/unit/test_solvers/test_ode_solver.py
new file mode 100644
index 000000000..5e2638958
--- /dev/null
+++ b/tests/unit/test_solvers/test_ode_solver.py
@@ -0,0 +1,30 @@
+#
+# Tests for the ODE Solver class
+#
+import pybamm
+import unittest
+
+
+class TestOdeSolver(unittest.TestCase):
+    def test_wrong_solver(self):
+        # Create model
+        model = pybamm.BaseModel()
+        var = pybamm.Variable("var")
+        model.rhs = {var: var}
+        model.algebraic = {var: var - 1}
+
+        # test errors
+        solver = pybamm.OdeSolver()
+        with self.assertRaisesRegex(
+            pybamm.SolverError, "Cannot use ODE solver to solve model with DAEs"
+        ):
+            solver.solve(model, None)
+
+
+if __name__ == "__main__":
+    print("Add -v for more debug output")
+    import sys
+
+    if "-v" in sys.argv:
+        debug = True
+    unittest.main()
diff --git a/tests/unit/test_solvers/test_scikits_solvers.py b/tests/unit/test_solvers/test_scikits_solvers.py
index bfbf65432..67f1ce417 100644
--- a/tests/unit/test_solvers/test_scikits_solvers.py
+++ b/tests/unit/test_solvers/test_scikits_solvers.py
@@ -457,6 +457,7 @@ class TestScikitsSolvers(unittest.TestCase):
         model.rhs = {var1: 0.1 * var1}
         model.algebraic = {var2: 2 * var1 - var2}
         model.initial_conditions = {var1: 1, var2: 2}
+        model.use_jacobian = False
         disc = get_discretisation_for_testing()
         disc.process_model(model)
 
@@ -551,6 +552,22 @@ class TestScikitsSolvers(unittest.TestCase):
         np.testing.assert_allclose(solver.y[0], np.exp(0.1 * solver.t))
         np.testing.assert_allclose(solver.y[-1], 2 * np.exp(0.1 * solver.t))
 
+    def test_solve_ode_model_with_dae_solver(self):
+        # Create model
+        model = pybamm.BaseModel()
+        var = pybamm.Variable("var")
+        model.rhs = {var: 0.1 * var}
+        model.initial_conditions = {var: 1}
+        disc = get_discretisation_for_testing()
+        disc.process_model(model)
+
+        # Solve
+        solver = pybamm.ScikitsDaeSolver(tol=1e-8)
+        t_eval = np.linspace(0, 1, 100)
+        solver.solve(model, t_eval)
+        np.testing.assert_array_equal(solver.t, t_eval)
+        np.testing.assert_allclose(solver.y[0], np.exp(0.1 * solver.t))
+
 
 if __name__ == "__main__":
     print("Add -v for more debug output")

pip install -e .[dev,docs]
parse_log_pytest
pytest
pytest
apt-get update apt-get install -y gfortran gcc libopenblas-dev liblapack-dev
3.6
pytest --no-header -rA --tb=line --color=no -p no:cacheprovider -W ignore::DeprecationWarning
alabaster==0.7.13
anyio==3.6.2
anytree==2.8.0
argon2-cffi==21.3.0
argon2-cffi-bindings==21.2.0
async-generator==1.10
attrs @ file:///opt/conda/conda-bld/attrs_1642510447205/work
autograd==1.6.2
Babel==2.11.0
backcall==0.2.0
bleach==4.1.0
certifi==2021.5.30
cffi==1.15.1
charset-normalizer==2.0.12
comm==0.1.4
contextvars==2.4
cycler==0.11.0
dataclasses==0.8
decorator==5.1.1
defusedxml==0.7.1
docutils==0.17.1
entrypoints==0.4
flake8==5.0.4
future==1.0.0
guzzle-sphinx-theme==0.7.11
idna==3.10
imagesize==1.4.1
immutables==0.19
importlib-metadata==4.2.0
iniconfig @ file:///home/linux1/recipes/ci/iniconfig_1610983019677/work
ipykernel==5.5.6
ipython==7.16.3
ipython-genutils==0.2.0
ipywidgets==7.8.5
jedi==0.17.2
Jinja2==3.0.3
json5==0.9.16
jsonschema==3.2.0
jupyter==1.1.1
jupyter-client==7.1.2
jupyter-console==6.4.3
jupyter-core==4.9.2
jupyter-server==1.13.1
jupyterlab==3.2.9
jupyterlab-pygments==0.1.2
jupyterlab-server==2.10.3
jupyterlab_widgets==1.1.11
kiwisolver==1.3.1
MarkupSafe==2.0.1
matplotlib==3.3.4
mccabe==0.7.0
mistune==0.8.4
more-itertools @ file:///tmp/build/80754af9/more-itertools_1637733554872/work
nbclassic==0.3.5
nbclient==0.5.9
nbconvert==6.0.7
nbformat==5.1.3
nest-asyncio==1.6.0
notebook==6.4.10
numpy==1.19.5
packaging @ file:///tmp/build/80754af9/packaging_1637314298585/work
pandas==1.1.5
pandocfilters==1.5.1
parso==0.7.1
pexpect==4.9.0
pickleshare==0.7.5
Pillow==8.4.0
pluggy @ file:///tmp/build/80754af9/pluggy_1615976315926/work
prometheus-client==0.17.1
prompt-toolkit==3.0.36
ptyprocess==0.7.0
py @ file:///opt/conda/conda-bld/py_1644396412707/work
-e git+https://github.com/pybamm-team/PyBaMM.git@8ef96429e9798d9225c062b380b883d36f998ab7#egg=pybamm
pycodestyle==2.9.1
pycparser==2.21
pyflakes==2.5.0
Pygments==2.14.0
pyparsing @ file:///tmp/build/80754af9/pyparsing_1635766073266/work
pyrsistent==0.18.0
pytest==6.2.4
python-dateutil==2.9.0.post0
pytz==2025.2
pyzmq==25.1.2
requests==2.27.1
scipy==1.5.4
Send2Trash==1.8.3
six==1.17.0
sniffio==1.2.0
snowballstemmer==2.2.0
Sphinx==4.3.2
sphinxcontrib-applehelp==1.0.2
sphinxcontrib-devhelp==1.0.2
sphinxcontrib-htmlhelp==2.0.0
sphinxcontrib-jsmath==1.0.1
sphinxcontrib-qthelp==1.0.3
sphinxcontrib-serializinghtml==1.1.5
terminado==0.12.1
testpath==0.6.0
toml @ file:///tmp/build/80754af9/toml_1616166611790/work
tornado==6.1
traitlets==4.3.3
typing_extensions @ file:///opt/conda/conda-bld/typing_extensions_1647553014482/work
urllib3==1.26.20
wcwidth==0.2.13
webencodings==0.5.1
websocket-client==1.3.1
widgetsnbextension==3.6.10
zipp @ file:///tmp/build/80754af9/zipp_1633618647012/work

name: PyBaMM
channels:
  - defaults
  - https://repo.anaconda.com/pkgs/main
  - https://repo.anaconda.com/pkgs/r
  - conda-forge
dependencies:
  - _libgcc_mutex=0.1=main
  - _openmp_mutex=5.1=1_gnu
  - attrs=21.4.0=pyhd3eb1b0_0
  - ca-certificates=2025.2.25=h06a4308_0
  - certifi=2021.5.30=py36h06a4308_0
  - importlib_metadata=4.8.1=hd3eb1b0_0
  - iniconfig=1.1.1=pyhd3eb1b0_0
  - ld_impl_linux-64=2.40=h12ee557_0
  - libffi=3.3=he6710b0_2
  - libgcc-ng=11.2.0=h1234567_1
  - libgomp=11.2.0=h1234567_1
  - libstdcxx-ng=11.2.0=h1234567_1
  - more-itertools=8.12.0=pyhd3eb1b0_0
  - ncurses=6.4=h6a678d5_0
  - openssl=1.1.1w=h7f8727e_0
  - packaging=21.3=pyhd3eb1b0_0
  - pip=21.2.2=py36h06a4308_0
  - pluggy=0.13.1=py36h06a4308_0
  - py=1.11.0=pyhd3eb1b0_0
  - pyparsing=3.0.4=pyhd3eb1b0_0
  - pytest=6.2.4=py36h06a4308_2
  - python=3.6.13=h12debd9_1
  - readline=8.2=h5eee18b_0
  - setuptools=58.0.4=py36h06a4308_0
  - sqlite=3.45.3=h5eee18b_0
  - tk=8.6.14=h39e8969_0
  - toml=0.10.2=pyhd3eb1b0_0
  - typing_extensions=4.1.1=pyh06a4308_0
  - wheel=0.37.1=pyhd3eb1b0_0
  - xz=5.6.4=h5eee18b_1
  - zipp=3.6.0=pyhd3eb1b0_0
  - zlib=1.2.13=h5eee18b_1
  - pip:
      - alabaster==0.7.13
      - anyio==3.6.2
      - anytree==2.8.0
      - argon2-cffi==21.3.0
      - argon2-cffi-bindings==21.2.0
      - async-generator==1.10
      - autograd==1.6.2
      - babel==2.11.0
      - backcall==0.2.0
      - bleach==4.1.0
      - cffi==1.15.1
      - charset-normalizer==2.0.12
      - comm==0.1.4
      - contextvars==2.4
      - cycler==0.11.0
      - dataclasses==0.8
      - decorator==5.1.1
      - defusedxml==0.7.1
      - docutils==0.17.1
      - entrypoints==0.4
      - flake8==5.0.4
      - future==1.0.0
      - guzzle-sphinx-theme==0.7.11
      - idna==3.10
      - imagesize==1.4.1
      - immutables==0.19
      - importlib-metadata==4.2.0
      - ipykernel==5.5.6
      - ipython==7.16.3
      - ipython-genutils==0.2.0
      - ipywidgets==7.8.5
      - jedi==0.17.2
      - jinja2==3.0.3
      - json5==0.9.16
      - jsonschema==3.2.0
      - jupyter==1.1.1
      - jupyter-client==7.1.2
      - jupyter-console==6.4.3
      - jupyter-core==4.9.2
      - jupyter-server==1.13.1
      - jupyterlab==3.2.9
      - jupyterlab-pygments==0.1.2
      - jupyterlab-server==2.10.3
      - jupyterlab-widgets==1.1.11
      - kiwisolver==1.3.1
      - markupsafe==2.0.1
      - matplotlib==3.3.4
      - mccabe==0.7.0
      - mistune==0.8.4
      - nbclassic==0.3.5
      - nbclient==0.5.9
      - nbconvert==6.0.7
      - nbformat==5.1.3
      - nest-asyncio==1.6.0
      - notebook==6.4.10
      - numpy==1.19.5
      - pandas==1.1.5
      - pandocfilters==1.5.1
      - parso==0.7.1
      - pexpect==4.9.0
      - pickleshare==0.7.5
      - pillow==8.4.0
      - prometheus-client==0.17.1
      - prompt-toolkit==3.0.36
      - ptyprocess==0.7.0
      - pycodestyle==2.9.1
      - pycparser==2.21
      - pyflakes==2.5.0
      - pygments==2.14.0
      - pyrsistent==0.18.0
      - python-dateutil==2.9.0.post0
      - pytz==2025.2
      - pyzmq==25.1.2
      - requests==2.27.1
      - scipy==1.5.4
      - send2trash==1.8.3
      - six==1.17.0
      - sniffio==1.2.0
      - snowballstemmer==2.2.0
      - sphinx==4.3.2
      - sphinxcontrib-applehelp==1.0.2
      - sphinxcontrib-devhelp==1.0.2
      - sphinxcontrib-htmlhelp==2.0.0
      - sphinxcontrib-jsmath==1.0.1
      - sphinxcontrib-qthelp==1.0.3
      - sphinxcontrib-serializinghtml==1.1.5
      - terminado==0.12.1
      - testpath==0.6.0
      - tornado==6.1
      - traitlets==4.3.3
      - urllib3==1.26.20
      - wcwidth==0.2.13
      - webencodings==0.5.1
      - websocket-client==1.3.1
      - widgetsnbextension==3.6.10
prefix: /opt/conda/envs/PyBaMM

tests/unit/test_solvers/test_ode_solver.py::TestOdeSolver::test_wrong_solver
pybamm/solvers/scikits_ode_solver.py pybamm/solvers/scipy_solver.py pybamm/expression_tree/concatenations.py pybamm/solvers/dae_solver.py pybamm/solvers/base_solver.py pybamm/solvers/ode_solver.py pybamm/solvers/scikits_dae_solver.py
pybamm/solvers/scikits_ode_solver.py pybamm/solvers/scipy_solver.py pybamm/expression_tree/concatenations.py pybamm/solvers/dae_solver.py pybamm/solvers/base_solver.py pybamm/solvers/ode_solver.py pybamm/solvers/scikits_dae_solver.py