diff --git a/meshwell/model.py b/meshwell/model.py
index 0bc3aa6..0547685 100644
--- a/meshwell/model.py
+++ b/meshwell/model.py
@@ -33,6 +33,7 @@ def __init__(
         gmsh.option.setNumber("Mesh.MaxNumThreads1D", n_threads)
         gmsh.option.setNumber("Mesh.MaxNumThreads2D", n_threads)
         gmsh.option.setNumber("Mesh.MaxNumThreads3D", n_threads)
+        gmsh.option.setNumber("Geometry.OCCParallel", n_threads>1)
 
         # Point snapping
         self.point_tolerance = point_tolerance