Fix copy_triangle for Julia 1.11

JuliaSmoothOptimizers · May 22, 2024 · 7214572 · 7214572
1 parent 5605b36
commit 7214572
Showing 1 changed file with 7 additions and 1 deletion.
diff --git a/src/block_krylov_utils.jl b/src/block_krylov_utils.jl
@@ -164,7 +164,13 @@ function copy_triangle(Q::AbstractMatrix{FC}, R::AbstractMatrix{FC}, k::Int) whe
       end
     end
   else
-    copytrito!(R, Q, 'U')
+    mR, nR = size(R)
+    mQ, nQ = size(Q)
+    if (mR == mQ) && (nR == nQ)
+      copytrito!(R, Q, 'U')
+    else
+      copytrito!(R, view(Q, 1:k, 1:k), 'U')
+    end
   end
   return R
 end