LU Factorization

\ is left-associative in Julia, so this is parsed as

(LU.U\LU.L)\(LU.P*b)

whereas you want:

LU.U\(LU.L\(LU.P*b))

or equivalently LU.U\(LU.L\b[LU.p]), which gives the same answer as A \ b.

2 Likes