diff --git a/src/pyfmi/fmi_algorithm_drivers.py b/src/pyfmi/fmi_algorithm_drivers.py
index 923c8a01..8a5e7c58 100644
--- a/src/pyfmi/fmi_algorithm_drivers.py
+++ b/src/pyfmi/fmi_algorithm_drivers.py
@@ -632,7 +632,7 @@ def _set_options(self):
                 self.with_jacobian = True
             else:
                 fnbr, gnbr = self.model.get_ode_sizes()
-                if fnbr >= PYFMI_JACOBIAN_LIMIT and solver == "CVode":
+                if fnbr >= PYFMI_JACOBIAN_LIMIT and (solver == "CVode" or solver == "Radau5ODE"):
                     self.with_jacobian = True
                     if fnbr >= PYFMI_JACOBIAN_SPARSE_SIZE_LIMIT:
                         try: