And Staro Pickle writes: > I am doing some matrix factorizations like LU, LDL'. I find that the > julia's lufact() is much faster than my code. I wonder why my code is > slow and why lufact is fast(is this one question or two?). I use > vectorization(is that vectorization?).
Others have given some reasons... If you want to see contortions that were (perhaps still are) necessary to be competitive, see: https://groups.google.com/forum/#!msg/julia-users/DJRxApmpJmU/NqepNaX5ARIJ > n=5, a=10 > Matrix to be decomposed: > 1 0 0 0 a > 0 1 0 0 a > 0 0 1 0 a > 0 0 0 1 a > a a a a 1 However, if your matrices always have this form, I would highly recommend solving a few by hand to notice the pattern.
