diff --git a/Project.toml b/Project.toml
index 2dc2807e4..f44d9025e 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "LinearSolve"
 uuid = "7ed4a6bd-45f5-4d41-b270-4a48e9bafcae"
 authors = ["SciML"]
-version = "1.41.0"
+version = "1.41.1"
 
 [deps]
 ArrayInterface = "4fba245c-0d91-5ea0-9b3e-6abc04ee57a9"
diff --git a/src/factorization.jl b/src/factorization.jl
index 8d0fa6ba1..421f42fa3 100644
--- a/src/factorization.jl
+++ b/src/factorization.jl
@@ -91,7 +91,7 @@ end
 
 function do_factorization(alg::QRFactorization, A, b, u)
     A = convert(AbstractMatrix, A)
-    if alg.inplace && !(A isa SparseMatrixCSC)
+    if alg.inplace && !(A isa SparseMatrixCSC) && !(A isa GPUArraysCore.AbstractGPUArray)
         fact = qr!(A, alg.pivot)
     else
         fact = qr(A) # CUDA.jl does not allow other args!