diff --git a/src/pyfmi/fmi_algorithm_drivers.py b/src/pyfmi/fmi_algorithm_drivers.py
index b67b4d3c..8a5e7c58 100644
--- a/src/pyfmi/fmi_algorithm_drivers.py
+++ b/src/pyfmi/fmi_algorithm_drivers.py
@@ -207,7 +207,7 @@ def __init__(self, *args, **kw):
             'sensitivities':None,
             'write_scaled_result':False,
             'result_file_name':'',
-            'with_jacobian':True,
+            'with_jacobian':"Default",
             'logging':False,
             'dynamic_diagnostics':False,
             'result_handling':"binary",
@@ -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: