Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

minres! crashes hard for a complex system #301

Closed
RaghuSivapuram opened this issue Jun 20, 2021 · 5 comments
Closed

minres! crashes hard for a complex system #301

RaghuSivapuram opened this issue Jun 20, 2021 · 5 comments

Comments

@RaghuSivapuram
Copy link

I am using minres! to obtain an improved solution than the given initial guess. The system matrices and the initial guess are complex.
I get this crash when I call: minres!(U, Kd, F; maxiter=20) where Kd, U are complex.

Please submit a bug report with steps to reproduce this fault, and any error messages that follow (in their entirety). Thanks.
Exception: EXCEPTION_ACCESS_VIOLATION at 0x1fec08e0 --
Please submit a bug report with steps to reproduce this fault, and any error messages that follow (in their entirety). Thanks.
Exception: EXCEPTION_ACCESS_VIOLATION at 0x1fec08e0 -- OLATION with steps to reproduce this fault, and any error messages that follow (in their entirety). Thanks.
Exception: EXCEPTION_ACCESS_VIOLATION with steps to reproduce this fault, and any error messages that follow (in their entirety). Thanks.
Exception: EXCEPTION_ACCESS_VDCABS164_ at C:\Users\sivap\AppData\Local\Programs\Julia-1.6.1\bin\libopenblas64_.DLL (unknown line)
in expression starting at C:\Users\sivap.julia\MyPackages\dev\FinEtoolsRapidHarmonicVA\batchr.jl:43
in expression starting at at 0x1fec08e0 -- at 0x1fec08e0 -- ing at at 0x1fec08e0 -- at 0x1fec08e0 -- rs\sivap\AppData\Local\Programs\Julia-1.6.1\bin\libopenblas64_.DLL (unknown line)

Please let me know if I'm misusing the function.

@haampie
Copy link
Member

haampie commented Jun 21, 2021

Can you provide something reproducible? I don't see this failing.

julia> using LinearAlgebra, IterativeSolvers

julia> A = rand(ComplexF64, 100, 100); A += A' + 100I; x = ones(ComplexF64, 100); b = A * x;

julia> A * minres(A, b) - b |> norm
1.7369022090081076e-5

@fredrikekre
Copy link
Member

JuliaLang/julia#40963

@RaghuSivapuram
Copy link
Author

using LinearAlgebra, IterativeSolvers, SparseArrays

A = sprand(19400, 19400, 0.1) + 1.0im * sprand(19400, 19400, 0.1)
F = rand(ComplexF64, 19400)
x = zeros(ComplexF64, 19400)

minres!(x, A, F; maxiter=20) blows up.

@haampie
Copy link
Member

haampie commented Jun 22, 2021

@RaghuSivapuram can you share more details (os, cpu, julia version...) in JuliaLang/julia#40963?

@RaghuSivapuram
Copy link
Author

This problem doesn't occur in Julia 1.7.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants