Arpack.jl

Arpack.jl copied to clipboard

complex-Hermitian support?

1

Looking at the source code, it seems that it explicitly checks for the real-symmetric case, but doesn't do anything for the analogous complex-Hermitian case. Is this a limitation of the...

stevengj

“eigs” function doesn’t work for a sparse matrix with one or more zero pivots

1

Hi, I am trying to obtain the lowest eigenvalue of a sparse hermitian matrix with one or more zero pivots. For example, my matrix has the following form: ``` ham...

k-tamuraphys

Generalized eigenvalue problem with symmetric matrices: complex solution!?

2

Arpack v0.5.3: The generalized eigenvalue problem with symmetric matrices should return real eigenvalues and real eigenvectors. The `eigs` function fails to do that: ``` using LinearAlgebra using SparseArrays using Arpack...

PetrKryslUCSD

Please refer to report https://discourse.julialang.org/t/is-arpack-really-slow-on-windows/73467 The issue in brief: Arpack call takes ~25% more time on Windows than in WSL2 running under Windows. Different Julia executables, but both 1.7.0.

PetrKryslUCSD

Arpack-ng bug preventing us upgrading from 3.5 to 3.8

23

If I have done my analysis right - the failure is introduced in https://github.com/opencollab/arpack-ng/compare/0e7d01d..7fc42e5 I'm not sure if we are doing something wrong in our wrappers or Arpack-ng introduced a...

ViralBShah

svds giving different singular vectors for same matrix

5

I am trying to use svds but svds for same matrix is giving different singular vectors after every run. Due to this reason my codes are non-reproducible. I also checked...

shashins

explicittransform does an unintuitive thing

2

I have two symmetric PD sparse matrices, and I call ``` d, v, nconv = eigs(Symmetric(K+OmegaShift*M), Symmetric(M); nev=neigvs, which=:SM) ``` Clearly it is reasonable to expect the results to be...

PetrKryslUCSD

Instead of switching to complex matrices if the input is nonsymmetric (in the case of standard eigendecomposition), one could have a straightforward version that retains the same eltype

matbesancon

svds StackOverfowError for Float16 matrices

1

`svds` throws a `StackOverflowError` on Julia 1.2-DEV.172 when applied to `Float16` matrices, sparse or not (could not test on Julia 1.0): ```julia julia> svds(Float16[1 2 3; 4 5 6;7 8...

zgornel

Using Arpack with MKL

14

```julia julia> versioninfo() Julia Version 0.7.0-beta.65 Commit e2de9b3d2e* (2018-06-28 03:45 UTC) Platform Info: OS: Linux (x86_64-linux-gnu) CPU: Intel(R) Core(TM) i5-3570 CPU @ 3.40GHz WORD_SIZE: 64 LIBM: libopenlibm LLVM: libLLVM-6.0.0 (ORCJIT,...

dmbates