diff --git a/qutip/solver/heom/bofin_baths.py b/qutip/solver/heom/bofin_baths.py
index ceba267232..7a0bd4c08f 100644
--- a/qutip/solver/heom/bofin_baths.py
+++ b/qutip/solver/heom/bofin_baths.py
@@ -1424,8 +1424,9 @@ def __init__(self, T, Q, alpha, wc, s, Nk=4, method="spectral", rmse=7e-5):
         self.Nk = Nk
         self.beta = 1 / T
         if _mpmath_available is False:
-            raise Exception("The mpmath module is needed for the description"
-                            "of Ohmic baths")
+            print(
+                "The mpmath module is needed for the description"
+                " of Ohmic baths")
         if method == "correlation":
             self.fit = FitCorr(self.Q)
             t = np.linspace(0, 15 / self.wc, 1000)