Yes, it looks like the same thing. If I do norm(zeros(129,129)) on my machine, I also get an Abort. Also eigfact(zeros(129,129)).
My machine is a Mac: versioninfo() Julia Version 0.4.2 Commit bb73f34 (2015-12-06 21:47 UTC) Platform Info: System: Darwin (x86_64-apple-darwin13.4.0) CPU: Intel(R) Core(TM) i7-2600 CPU @ 3.40GHz WORD_SIZE: 64 BLAS: libopenblas (USE64BITINT DYNAMIC_ARCH NO_AFFINITY Sandybridge) LAPACK: libopenblas64_ LIBM: libopenlibm LLVM: libLLVM-3.3 On Monday, January 18, 2016 at 3:13:52 PM UTC-8, Kevin Squire wrote: > > Hi Steven, > > This seems to be a problem with specific machines--I don't run into it, > for example. > > Are you on a Mac? If so, check out > https://github.com/JuliaLang/julia/issues/14507 and > https://github.com/staticfloat/homebrew-julia/issues/194 and see if > they're the same issue. If they are, can you comment there? (Additional > comments on open issues sometimes helps them get resolved faster) > > Cheers, > Kevin > > On Mon, Jan 18, 2016 at 2:57 PM, Steven White <[email protected] > <javascript:>> wrote: > >> I discovered this with matrices that were not random, but this simple test >> illustrates the problem: julia aborts when it tries to diagonalize a >> 705x705 matrix >> >> mat = rand(704,704) >> mat = mat' + mat >> e = eigfact(mat) >> @show "done 704" >> >> mat = rand(705,705) >> mat = mat' + mat >> e = eigfact(mat) >> @show "done 705" >> >> Results: >> "done 704" = "done 704" >> Abort >> >> I also get this when I try an svd: >> mat = rand(705,705) >> svd(mat) >> >> This gives >> Abort >> >> So, to find what the key problem sizes are: >> for i=1:705 >> mat = rand(i,i) >> svd(mat) >> @show i >> end >> >> This dies on i=129. Note that 129 is odd. >> But replacing the for loop to count by twos is fine: >> for i=2:2:1705 >> mat = rand(i,i) >> svd(mat) >> @show i >> end >> No problem here. >> >> >
