diff --git a/Project.toml b/Project.toml
index e4c86156..fc5f0e1d 100644
--- a/Project.toml
+++ b/Project.toml
@@ -12,7 +12,7 @@ SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 [compat]
 FastClosures = "0.3.0"
 LinearAlgebra = "1.10"
-LinearOperators = "2"
+LinearOperators = "2.13"
 Printf = "1.10"
 SparseArrays = "1.10"
 Test = "1.10"
diff --git a/src/nlp/api.jl b/src/nlp/api.jl
index ca15e375..fb4e6a07 100644
--- a/src/nlp/api.jl
+++ b/src/nlp/api.jl
@@ -732,7 +732,7 @@ function jac_op!(
     end
     return res
   end
-  return LinearOperator{T}(nlp.meta.ncon, nlp.meta.nvar, false, false, prod!, ctprod!, ctprod!)
+  return LinearOperator{T,S}(nlp.meta.ncon, nlp.meta.nvar, false, false, prod!, ctprod!, ctprod!)
 end
 
 """
@@ -772,7 +772,7 @@ function jac_op!(
     end
     return res
   end
-  return LinearOperator{T}(nlp.meta.ncon, nlp.meta.nvar, false, false, prod!, ctprod!, ctprod!)
+  return LinearOperator{T,S}(nlp.meta.ncon, nlp.meta.nvar, false, false, prod!, ctprod!, ctprod!)
 end
 
 """
@@ -823,7 +823,7 @@ function jac_lin_op!(
     end
     return res
   end
-  return LinearOperator{T}(nlp.meta.nlin, nlp.meta.nvar, false, false, prod!, ctprod!, ctprod!)
+  return LinearOperator{T,S}(nlp.meta.nlin, nlp.meta.nvar, false, false, prod!, ctprod!, ctprod!)
 end
 
 """
@@ -863,7 +863,7 @@ function jac_lin_op!(
     end
     return res
   end
-  return LinearOperator{T}(nlp.meta.nlin, nlp.meta.nvar, false, false, prod!, ctprod!, ctprod!)
+  return LinearOperator{T,S}(nlp.meta.nlin, nlp.meta.nvar, false, false, prod!, ctprod!, ctprod!)
 end
 
 """
@@ -914,7 +914,7 @@ function jac_nln_op!(
     end
     return res
   end
-  return LinearOperator{T}(nlp.meta.nnln, nlp.meta.nvar, false, false, prod!, ctprod!, ctprod!)
+  return LinearOperator{T,S}(nlp.meta.nnln, nlp.meta.nvar, false, false, prod!, ctprod!, ctprod!)
 end
 
 """
@@ -954,7 +954,7 @@ function jac_nln_op!(
     end
     return res
   end
-  return LinearOperator{T}(nlp.meta.nnln, nlp.meta.nvar, false, false, prod!, ctprod!, ctprod!)
+  return LinearOperator{T,S}(nlp.meta.nnln, nlp.meta.nvar, false, false, prod!, ctprod!, ctprod!)
 end
 
 """
@@ -1347,7 +1347,7 @@ function hess_op!(
     end
     return res
   end
-  return LinearOperator{T}(nlp.meta.nvar, nlp.meta.nvar, true, true, prod!, prod!, prod!)
+  return LinearOperator{T,S}(nlp.meta.nvar, nlp.meta.nvar, true, true, prod!, prod!, prod!)
 end
 
 """
@@ -1379,7 +1379,7 @@ function hess_op!(
     end
     return res
   end
-  return LinearOperator{T}(nlp.meta.nvar, nlp.meta.nvar, true, true, prod!, prod!, prod!)
+  return LinearOperator{T,S}(nlp.meta.nvar, nlp.meta.nvar, true, true, prod!, prod!, prod!)
 end
 
 """
@@ -1411,7 +1411,7 @@ function hess_op!(
     end
     return res
   end
-  return LinearOperator{T}(nlp.meta.nvar, nlp.meta.nvar, true, true, prod!, prod!, prod!)
+  return LinearOperator{T,S}(nlp.meta.nvar, nlp.meta.nvar, true, true, prod!, prod!, prod!)
 end
 
 """
diff --git a/src/nls/api.jl b/src/nls/api.jl
index a21f2f93..4fc7b912 100644
--- a/src/nls/api.jl
+++ b/src/nls/api.jl
@@ -217,7 +217,7 @@ function jac_op_residual!(
     end
     return res
   end
-  return LinearOperator{T}(
+  return LinearOperator{T,S}(
     nls_meta(nls).nequ,
     nls_meta(nls).nvar,
     false,
@@ -264,7 +264,7 @@ function jac_op_residual!(
     end
     return res
   end
-  return LinearOperator{T}(
+  return LinearOperator{T,S}(
     nls_meta(nls).nequ,
     nls_meta(nls).nvar,
     false,
@@ -449,7 +449,7 @@ function hess_op_residual!(
     end
     return res
   end
-  return LinearOperator{T}(nls_meta(nls).nvar, nls_meta(nls).nvar, true, true, prod!, prod!, prod!)
+  return LinearOperator{T,S}(nls_meta(nls).nvar, nls_meta(nls).nvar, true, true, prod!, prod!, prod!)
 end
 
 """