Stack Overflow
Stack Overflow

Questions tagged [arpack]

ARPACK is a package for solving large scale eigen value problems. Question tagged arpack relate to usage of this library, linking with it and compiling code using its subroutines.

ARPACK is a collection of Fortran77 subroutines designed to solve large scale eigenvalue problems.

For more information, see ARPACK's website

72 questions
1
vote
1 answer

scipy generalized eigenproblem with positive semidefinite

Hi, guys!!! I want to compute generalized eigendecomposition of the form: Lf = lambda Af by using scipy.sparse.linalg.eigs function, but get this error: /usr/local/lib/python2.7/dist-packages/scipy/linalg/decomp_lu.py:61: RuntimeWarning: Diagonal…
maksmara
  • 27
  • 6
0
votes
1 answer

Error in integers when building Arpack-ng with mkl

I'm trying to write a program that is going to use both Arpack-ng and MKL routines. I was testing the behavior by running a simple code which use a single call to the function dsaupd, which should return ido=-1 and info = 0, but it actually…
VINCENZO BISOGNO
  • 1
  • 3
0
votes
0 answers

Best generalised eigenvalue solver symmetric real matrices Julia

I need to solve generalised eigenvalue problems that stem from solid mechanics (mass/stiffness matrices, (K-l*M)x=0). The systems dimension range from 500k to 15M degrees of freedom. Matrices are real and symmetric. I am currenty using Arpack.jl to…
Smoop94
  • 73
  • 3
0
votes
0 answers

Detected problem when computing eigenvalues using octave eigs function

I need to perform Principal Components Analysis in a dataset, I am using octave 6.4.0, so it should not be a problem related to outdated software. Anyway, the function I am using relies on a very well tested so I proceeded as usual: Compute the…
user4518773
  • 1
0
votes
0 answers

Is nx.eigenvector_centrality_numpy() using the Arnoldi iteration instead of the basic power method?

Since nx.eigenvector_centrality_numpy() using ARPACK, is it mean that nx.eigenvector_centrality_numpy() using Arnoldi iteration instead of the basic power method? because when I try to compute manually using the basic power method, the result of my…
Nurul Jannah
  • 1
  • 1
0
votes
1 answer

I want to obtain eigenvalues of symmetric matrix in Julia in O(nmr)

I am a beginner at Julia. I want to obtain r eigenvalues and eigenvectors of input symmetric n times n matrix X in increasing order. I heard the computational complexity is O(n^2 r). n is around 1000-20000, r is around 100-1000. How can I obtain the…
Sakurai.JJ
  • 553
  • 3
  • 10
0
votes
1 answer

How to perturb datas for quadratic optimization wiht CVXPY

I have the same problem described here CVXPY Quadratic Programming; ArpackNoConvergence error. I would like to try the pertubation/trasformation solution, but i have no idea how to perform it. I have a quadratic problem like this…
Wodes
  • 1
0
votes
0 answers

PARPACK implementation runs into memory errors

I am making a module in Fortran 90 to run PARPACK on a given matrix. I have an existing ARPACK code which functions normally as expected. I tried converting it into PARPACK and it runs into memory clear errors. I am fairly new to coding and fortran,…
Srivatsank Pullabhotla
  • 11
  • 4
0
votes
0 answers

ARPACK error with eigenvector on igraph - python

I am recieving this error when I try to get the communities from somes Erdos graphs: InternalError: Error at d:\bld\python-igraph_1588168334861\work\vendor\build\igraph\igraph-0.9.0-pre+-msvc\src\arpack.c:1001: ARPACK error, Maximum number of…
Felipe
  • 1
0
votes
2 answers

Strange behaviour using Arpack EigenLibrary Wrapper for small matrices

I want to solve the following generalized EVP using the ArpackWrapper of the Eigen library: K_e is SPD. Normally K_g is indefinite and singular but for this MVP it is just indefinite. Furthermore, I'm interested in the smallest eigenvalues. For…
rath3t
  • 101
  • 4
0
votes
1 answer

counterintuitive speed difference between LM and shift-invert modes in scipy.sparse.linalg.eigsh?

I'm trying to find the smallest (as in most negative, not lowest magnitude) several eigenvalues of a list of sparse Hermitian matrices in Python using scipy.sparse.linalg.eigsh. The matrices are ~1000x1000, and the list length is ~500-2000. In…
lhk
  • 1
0
votes
1 answer

Problem linking my code with ARPACK on OSX (using MacPorts for ARPACK)

I am trying to compile a C++ program which invokes the ARPACK library. My problem is that when everything is linked, some of the symbols in the ARPACK library do not get resolved. They are __gfortran_transfer_integer __gfortran_transfer_character…
Marcus P S
  • 871
  • 10
  • 16
0
votes
0 answers

Trying to build ARPACK, but many source files are not being compiled

Let me start by saying I am not a superuser and somewhat of a noob here, but I have a basic understanding of what I'm doing. I have installed MSYS2 and MinGW-w64 following the directions here. I am working on 64-bit Windows 10, and I want to…
Kai
  • 213
  • 1
  • 12
0
votes
1 answer

Strange behaviour of ARPACK for hermitian matrix

I want to obtain numerically the energy of the ground state of some hermitian matrix (see the definition of this matrix in the following code) and plot it in terms of the matrix-parametter "phase". import scipy.sparse as sparse import scipy import…
Joe
  • 189
  • 10
0
votes
1 answer

How to compute the first eigenvalue and eigenvector in Fortran

I tried to use ARPACK using some examples (here), but I could not even figure out how to input my matrix. From this question, It seems that the implementations in Python and Matlab are the only solution to avoid the complexity of ARPACK. Is there…
Garini
  • 1,088
  • 16
  • 29
1 2 3
4
5