75 questions
- Bountied 0
- Unanswered
- Frequent
- Score
- Trending
- Week
- Month
- Unanswered (my tags)
2
votes
1
answer
133
views
PCA with arpack returns different values when the order of observations change, but why?
I have recently noticed that when I change the order of the observations in a sparse array, scikit-learn PCA with svd_solver="arpack" returns different floating point numbers. Is this an ...
5
votes
1
answer
103
views
Which differences between SciPy's and MATLAB's eigs function causes them to converge (or not) to different eigenvalues?
I am trying to solve a generalized eigenvalue problem in Python by translating a MATLAB script. Both A and B are not positive semi definite. In MATLAB, I'm guessing this makes the ARPACK library fall ...
0
votes
1
answer
69
views
Matrix free arnoldi
I want to do a "matrix-free eigenvalue evaluation" of a unitary matrix using ARPACK package. Unitary matrix comes from a sparse hermitian matrix, something that happens in many-body problem.
...
0
votes
0
answers
87
views
Can't link to C library despite explicit command in build.rs
I am trying to link a library against a local build of arpack. I have compiled arpack itself and I have it as a static binary. I then make this build.rs
use std::env;
use std::path::Path;
fn main() {
...
0
votes
1
answer
162
views
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 ...
1
vote
0
answers
127
views
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 ...
2
votes
0
answers
225
views
Why doesn't scipy.sparse.linalg.eigs produce a correct subset of the eigenvalues in this problem?
I have an eigenvalue problem which I would like to solve using SciPy. In this problem it is strongly preferable to use an iterative solver which will return a subset of the (largest/smallest) ...
0
votes
0
answers
208
views
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 ...
1
vote
1
answer
162
views
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 ...
2
votes
0
answers
193
views
Issue with finding eigenvalues with smallest magnitude using eigs
I have a large Hermitian sparse matrix, and I'd like to get a few of the eigenvalues with the smallest magnitude. I'm using the "eigs" function in Arpack however I get this error message:
...
0
votes
0
answers
84
views
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 ...
0
votes
1
answer
2k
views
I want to obtain eigenvalues of symmetric matrix in Julia in O(nmr) [closed]
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-...
2
votes
1
answer
608
views
Compute eigenvalues of complex-hermitian sparsematrix in Julia
Im working with some roughly 100000x100000 hermitian complex sparse-matrices, with roughly 5% of entries populated, and want to calculate the eigenvalues/eigenvectors.
Sofar ive been using Arpack.jl ...
2
votes
2
answers
830
views
Usage of arpack-ng/openblas/javacpp in the smile java library
I'm trying to set up the Smile machine learning library in Java but I'm having some issues getting some dependencies to work. Currently I am on a Fedora machine, but I'd like this to work on any ...
0
votes
1
answer
94
views
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 ...