diff --git a/motile/solver.py b/motile/solver.py
index 6c9328a..f827f88 100644
--- a/motile/solver.py
+++ b/motile/solver.py
@@ -105,7 +105,9 @@ def add_constraints(self, constraints: Constraint) -> None:
         for constraint in constraints.instantiate(self):
             self.constraints.add(constraint)
 
-    def solve(self, timeout: float = 0.0, num_threads: int = 1) -> ilpy.Solution:
+    def solve(
+        self, timeout: float = 0.0, num_threads: int = 1, verbose: bool = False
+    ) -> ilpy.Solution:
         """Solve the global optimization problem.
 
         Args:
@@ -113,9 +115,10 @@ def solve(self, timeout: float = 0.0, num_threads: int = 1) -> ilpy.Solution:
                 The timeout for the ILP solver, in seconds. Default (0.0) is no
                 timeout. If the solver times out, the best solution encountered
                 so far is returned (if any has been found at all).
-
             num_threads:
                 The number of threads the ILP solver uses.
+            verbose:
+                If true, print more information from ILP solver. Defaults to False.
 
         Returns:
             :class:`ilpy.Solution`, a vector of variable values. Use
@@ -142,7 +145,7 @@ def solve(self, timeout: float = 0.0, num_threads: int = 1) -> ilpy.Solution:
         if timeout > 0:
             self.ilp_solver.set_timeout(timeout)
 
-        self.ilp_solver.set_verbose(False)
+        self.ilp_solver.set_verbose(verbose)
 
         self.solution = self.ilp_solver.solve()
         if message := self.solution.get_status():