The UpperTriangular and LinAlg.UnitLowerTriangular parts are not necessary. Julia will detect that but since we already know that they are triangular we can save the costs of the checks.
A_ldiv_B!(A,b) is indeed ugly. It's just A\b where the result is stored in b. That is useful here because it will save us a couple of temporary allocations. On Tue, Jul 19, 2016 at 11:12 AM, Kurolong <twaldm...@gmail.com> wrote: > Thank you, this was very helpful. > Just a small question: Are the UpperTriangular and LowerTriangular > commands actually necessary in some case? > The other question i have is what A_div_B! actually does in this case, but > i have a feeling that this is hard to explain. >
