[eigen] [WIP] Power iteration algorithm

[ Thread Index | Date Index | More lists.tuxfamily.org/eigen Archives ]


Hello,

We have tried to implement the power iteration algorithm, and we came up with this version:
https://bitbucket.org/ensigll/eigen/src/eec44b78d44fbb57b0a4f2718ebccee5d6498b3e/Eigen/src/Eigenvalues/PowerIteration.h?at=ensigll

We have tested it using the existing eigensolver to make sure that the computed eigenvalues were correct.
Precision can be controlled with the parameter "threshold"; the smaller the threshold, the higher the precision and the longer the execution time.

Then we compared its effectiveness to eigensolver's. For the time being, program runs faster when the number of eigenvalues requested is lower than about 10% of the size of the input matrix.
However, we are still looking to improve the algorithm and make it usable in more general cases.

It's still a work in progress: we haven't finalized the algorithm (we'll try and make it more effective), neither the documentation.

Feel free to make any suggestions!
Best regards,
Imène, Bediss and Tristan.


Mail converted by MHonArc 2.6.19+ http://listengine.tuxfamily.org/