Update example2.jl

JuliaSmoothOptimizers · Sep 22, 2023 · 53ef564 · 53ef564
1 parent acefe34
commit 53ef564
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/paper/example2/example2.jl b/paper/example2/example2.jl
@@ -32,8 +32,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