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
2
votes
1 answer

Linking to ARPACK from VS2010

I am trying to get the ARPACK library to run on VS2010. I would like to use the C++ wrappers provided by ARPACK++ (some background - i need to get eigenvalues and eigenvectors of huge matrices). There is a very good tutorial on the subject…
Silentbit
  • 21
  • 3
1
vote
0 answers

Calculating full set of eigenvectors with arbitrary error tolerance

I'm trying to numerically compute the full set of eigenvectors of a large dense matrix, with the freedom to specify the error tolerance and maximum number of iterations. (The default error tolerance has proven too large for my purposes.) So far I…
Higglet
  • 11
  • 1
1
vote
1 answer

Custom C++ library dependent on ARPACK wrapped in pybind11 segfaults when NumPy is also imported

I am creating a custom library (written in C++) that does some numerical stuff with ARPACK-NG. The function is wrapped in pybind11 to provide access to the method from Python in a package. I observe strange behavior. Overview of the issue When NumPy…
Sock
  • 43
  • 4
1
vote
1 answer

Complex eigenvalues computation using scipy.sparse.linalg.eigs

Given the following input numpy 2d-array A that may be retrieved with the following link through the file hill_mat.npy, it would be great if I can compute only a subset of its eigenvalues using an iterative solver like…
Yacola
  • 2,873
  • 1
  • 10
  • 27
1
vote
1 answer

Has anyone had success installing Calculix 2.16 using the preCICE adapter?

I am trying to install the 2.16 version of Calculix. When installing according to that described at (https://github.com/precice/calculix-adapter/wiki/Installation-instructions-for-CalculiX)to my system (Ubuntu 16), I get an error…
FleaInTheSea11
  • 11
  • 3
1
vote
0 answers

Using an arpack tool from apt-get repository: undefined reference, segmentation fault

I installed Arpack using "sudo apt install libarpack2-dev". Next libraries are present in /usr/lib: libarpack.a, libarpack.so, libarpack.so.2, libarpack.so.2.0.0 I want to call the dsaupd subroutine in this library from a cpp file called…
Josja
  • 131
  • 7
1
vote
0 answers

ARPACK Eigenvalues with 16-Byte integer indexing

I have code that works fine to compute eigenvalues in my test case for ARPACK Shamelessly taken from here and adapted to a quick 4x4 matrix. (Comments at the top removed in my sample code for brevity). Okay, my problem. I have very large matrices,…
Jared Bland
  • 11
  • 3
1
vote
1 answer

Can I use parallel ARPACK in scipy?

I've been using scipy.sparse.linalg.eigs on some large matrices, and not surprisingly, it takes a while. So, I've been looking for ways to speed it up. My understand is that, under the hood, the scipy code use ARPACK, and there is a parallel version…
lnmaurer
  • 1,687
  • 4
  • 22
  • 32
1
vote
1 answer

how to run an MPI example in PARPACK

I am able to compile psndrv1.f, which is an example in ARPACK/PARPACK/EXAMPLE/MPI folder. However, when I run the program with the command mpirun -np 4 a.out, I got the following error. [.th:9951] *** An error occurred in MPI_Allreduce [.th:9951]…
wannik
  • 12,212
  • 11
  • 46
  • 58
1
vote
1 answer

Undefined reference to to `blacs_gridinfo_'

I am getting a undefined reference errors while building LSDMap (Locally Scaled Diffusion Map). What libraries contain the missing references? demo@argo-2 LSDMap]$ make mpifort -O3 -ffree-line-length-none src/p_wlsdmap.o -o p_wlsdmap…
AXP
  • 11
  • 2
1
vote
0 answers

stability of scipy.sparse.linalg.arpack.eigsh for eigen-decomposition of positive semi-definite sparse matrix?

I'm working on a variation of spectral clustering algorithm which performs standard spectral embedding on updated similarity matrix W for n_iters times, so I have to check the PSD of normalized laplacian matrix L_rw generated by updated W each…
LittleLittleQ
  • 440
  • 3
  • 7
  • 21
1
vote
2 answers

Using ARPACK solving eigenvalueproblem, but getting inconsistent results with Matlab

I'm new to ARPACK, I downloaded a script like the following import time import numpy as np from scipy.linalg import eigh from scipy.sparse.linalg import eigs np.set_printoptions(suppress=True) n=30 rstart=0 rend=n A=np.zeros(shape=(n,n)) # first…
jengmge
  • 77
  • 1
  • 5
1
vote
0 answers

Finding eigenvectors and eigenvalues of a sparse matrix with ARPACK ( called form PYTHON, MATLAB or as a FORTRAN subroutine)

Few days ago I asked a question how to find the eigenvalues of a large sparse matrix. I got no answers, so I decided to describe a potential solution. One question remains: Can I use the python implementation of ARPACK to compute the eigenvalues…
Alexander Cska
  • 738
  • 1
  • 7
  • 29
1
vote
2 answers

Inconsistent eigenvalues from eigs function in Scipy sparse

I am using eigs function from scipy.sparse.linalg module and found some inconsistent results. Running twice the same code gives different results, ie, the output of np.allclose is False. Can anyone explain why is that? from scipy.sparse.linalg…
poeticcapybara
  • 555
  • 1
  • 5
  • 15
1
vote
0 answers

Largescale sparse matrix inversion and using it in ARPACK OP

I'm trying rewrite matlab eigs function as C++ wrapper of ARPACK because arpack++ (C++ wrapper for arpack) is very slow in regular mode, shift-invert mode doesn't work properly. But I'm in stuck and I need a kick. :) LU Factorization Matlab…
user47779
  • 411
  • 1
  • 4
  • 9
1 2
3
4 5