What To Do When Julia's Lufact Runs Out Of Memory
Reid Atcheson, Accelerator Software Engineer at NAG has blogged about 'What to do when Julia's Lufact runs out of memory'. In his post he shows a sparse linear system which he was unable to solve with Julia’s
lufact because of its excessive memory requirements. Rather than resorting to an iterative algorithm he shows how to use nested dissection to get the same answer that
lufact would if it had enough memory. He starts by discussing the issues and why one may not want to use iterative solvers even though they solve the memory issue. Read his blog.