Fix the issue with the adjoint in JOSS

JuliaSmoothOptimizers · Sep 22, 2023 · acefe34 · acefe34
1 parent 5945d70
commit acefe34
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/paper/paper.md b/paper/paper.md
@@ -210,8 +210,10 @@ if CUDA.functional()
 
   # Solve Py = x
   function ldiv_ic0!(P, x, y, z)
-    ldiv!(z, LowerTriangular(P), x)   # Forward substitution with L
-    ldiv!(y, LowerTriangular(P)', z)  # Backward substitution with Lᴴ
+    L = LowerTriangular(P)
+    Lᴴ = adjoint(L)
+    ldiv!(z, L, x)   # Forward substitution with L
+    ldiv!(y, Lᴴ, z)  # Backward substitution with Lᴴ
     return y
   end