diff --git a/src/model.jl b/src/model.jl
index fb4e420..185be4a 100644
--- a/src/model.jl
+++ b/src/model.jl
@@ -198,13 +198,15 @@ end
 ###### J product ######
 #######################
 
+import Base.Threads
+
 function add_jprod!(
     model::Model,
     V::AbstractMatrix,
     Jv::AbstractVector,
     i::MatrixIndex,
 )
-    for j in 1:model.meta.ncon
+    Threads.@threads for j in 1:model.meta.ncon
         Jv[j] += LinearAlgebra.dot(model.A[i.value, j], V)
     end
 end