Hi, thanks for taking the time to reply. You made me go back and find a reference/copy error deep in my code... because of badly designed tests and bad luck in the circumstances triggering an error I had managed to convince myself the issue came from this linear system. Thanks again. Ben
On Thursday, February 11, 2016 at 4:26:44 PM UTC-5, Christoph Ortner wrote: > > > How can a sparse direct solver fail? That seems strange to me. With 1000 x > 1000 even a full direct solver should take no time: > > julia> A = rand(1000,1000); b = rand(1000); @time A \ b; > 0.023178 seconds (15 allocations: 7.645 MB) >
