Comments (2)
Hi, I used shift invert to solve a generalized symmetric sparse problem, similar to https://github.com/m-reuter/arpackpp/blob/master/examples/superlu/sym/lsymgshf.cc
I had the first eigenvalue be 0, which lead to instabilities, so I moved it out of the way slightly (using sigma to shift the problem away from zero).
In your case, I don't know. You can play around with inverting (so that eigenvalues increase) and shift so that all are above zero.
About the missing method in ARdsNonSymMatrix, not everything is completely implemented it seems. I have never used non symmetric matrices, so not sure if this was omitted on purpose or if it is just missing. Sorry I can't be more helpful.
from arpackpp.
No worries! And thanks for a quick reply.
I understand why mathematically it should work, I'm just wondering if there is a ready-made solution (it seems not). I think they left out those implementations on purpose (which is good), but didn't provide a convenient interface to implement it yourself (not good).
The manual says the (A - sigma * I)^(-1) * v
solver is automatically provided by ARPACK++ for ARluSymStdEig
, but it uses LU-decomposition. This is too memory-hungry and inefficient for large matrices, I believe? I'd probably use some sort of an iterative solver.
I wonder if there's any theory for how to combine IRAM and an iterative solver that is also based on Krylov subspaces. Because it seems to me that something like MINRES could benefit greatly from reusing the already computed Krylov subspaces.
from arpackpp.
Related Issues (20)
- Wrong method signatures HOT 1
- Wrong eigenvalues with UMFpack HOT 5
- compiling on windows HOT 1
- Installing with dynamic gfortran library HOT 8
- Valgrind Error from ARrcStdEig::FindEigenvalues
- Access Violation exception HOT 1
- Error with multiple threads HOT 2
- install-superlu.sh download link fix
- Question about the arpackpp HOT 1
- Integer overflow HOT 4
- Is it possible to use `gemv` from BLAS directly on `ARdsNonSymMatrix`? HOT 2
- Generalised eigenvalue problem with shift-inverse mode using UMFPACK HOT 2
- Plans to make a release? HOT 2
- Update the ARPACK++ documentation HOT 3
- Error while using ARNonSymStdEig HOT 2
- Duplicate symbols when linking two objects including arsnsym.h HOT 2
- Error when compile with gcc 5.4 HOT 4
- superlu 5.0 should set $BLAS HOT 1
- Posible memory leak? HOT 6
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from arpackpp.