Julia 1.7.0 OpenBLAS 可能有问题..

johnnychen94 · 2021 年12 月 3 日 15:56

解决方案是将BLAS线程设为 1… BLAS.set_num_threads(1)

julia> n = 1000;

julia> A = randn(Float32, n, n);

julia> inv(A)
ERROR: StackOverflowError:

julia> using LinearAlgebra

julia> BLAS.set_num_threads(1)

julia> inv(A);

github.com/JuliaLang/julia

1.7.0-rc2: StackOverflowError with complex-valued matrix exp

opened 10:54PM - 08 Nov 21 UTC

closed 09:54PM - 23 Dec 21 UTC

daviehh

bug regression linear algebra external dependencies

On macos, using version `1.7.0-rc2`, julia just shows `StackOverflowError` with …no other info when taking the matrix exp with a ~ 300x300 complex-valued matrix. Minimum example: ``` using LinearAlgebra n = 300 m = rand(ComplexF64, n, n); mex = exp(m); ``` Alro ran with `--startup-file=no` to make sure it's not some clash with other packages. <img width="1187" alt="Screen Shot 2021-11-08 at 5 52 35 PM" src="https://user-images.githubusercontent.com/25255906/140830950-e321af22-8fa9-4175-8cac-b1412da98143.png"> my `versioninfo()`: ``` Julia Version 1.7.0-rc2 Commit f23fc0d27a (2021-10-20 12:45 UTC) Platform Info: OS: macOS (x86_64-apple-darwin19.5.0) CPU: Intel(R) Core(TM) i5-8259U CPU @ 2.30GHz WORD_SIZE: 64 LIBM: libopenlibm LLVM: libLLVM-12.0.1 (ORCJIT, skylake) Environment: JULIA_MATHLINK = /Applications/Mathematica.app/Contents/Frameworks/mathlink.framework JULIA_MATHKERNEL = /Applications/Mathematica.app/Contents/MacOS/MathKernel ``` Same code runs fine with `Version 1.6.3 (2021-09-23)`. Hardware: Macbook pro (intel) with 16 GB ram, activity monitor shows low memory usage. In addition, the same code sometimes gives `ERROR: LoadError: ReadOnlyMemoryError()`, not sure how to reproduce that one... Thanks!

github.com/JuliaLang/julia

1.7.0: `inv` Causes Stack Overflow

opened 03:34PM - 02 Dec 21 UTC

closed 09:54PM - 23 Dec 21 UTC

SebastianAment

linear algebra

Following up on [this Discourse post](https://discourse.julialang.org/t/inv-caus…es-stack-overflow-on-julia-1-7-0-and-mac-os/72411). According to @stevengj, this issue might be solved by the resolution [this Yggdrasil issue](https://github.com/JuliaPackaging/Yggdrasil/issues/3996). I am running into an odd error when using Julia 1.7.0 on my 2017 Mac Book Pro (Intel i7) on Mac OS 12.01. In particular, executing ```Julia n = 1000; A = randn(Float32, n, n); inv(A) ERROR: StackOverflowError: ``` With no further stack trace printed afterwards. - The same error occurs with `ComplexF32` and `ComplexF64` types, but not `Float64`. - According to @carstenbauer, the issue is also present in 1.6.4, 1.7.0.rc-2, 1.7.0.rc-3, 1.7.0. - 1.6.3 and the current master branch do not exhibit the issue. More detailed info: On my machine, inv works for all element types up to a certain size. For `Float32`, the failure size is 514 by 514 and for the complex types, it’s 258 by 258. Below these cutoffs, everything works as usual.

johnnychen94 · 2021 年12 月 4 日 14:23

[OpenBLAS_jll] Update to `v0.3.10+10`

JuliaLang:backports-release-1.6 ← giordano:mg/openblas-32-threads-16

opened 02:13PM - 04 Dec 21 UTC

giordano

+93 -93

This version has been rebuilt to have 32 threads by default, instead of 512 as …it accidentally happened before. The large number of threads caused problems on some platforms, including `StackOverflowError`s. Similar to #43300, but for Julia v1.6. I wasn't able to test this one because I'm not able to reproduce the error on machines I have access to. Usual memo to self: to update the library, change the version number in `stdlib/OpenBLAS_jll/Project.toml`, run `make -f contrib/refresh_checksums.mk -j openblas` and remove the old checksums from the updated file `deps/checksums/openblas`.

woclass · 2021 年12 月 6 日 09:16

PR 昨天合并了。