[eigen] Implementation of eigs?
• To: eigen@xxxxxxxxxxxxxxxxxxx
• Subject: [eigen] Implementation of eigs?
• From: Mathieu Dutour <mathieu.dutour@xxxxxxxxx>
• Date: Thu, 18 Jun 2020 13:03:39 +0200
• Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20161025; h=mime-version:from:date:message-id:subject:to; bh=jg3Iob/sF705DnnYnql2zgOmdmHvhRwOxYGxaOFbt1E=; b=NPwfOG3dArYQHhFUjjxJEUNs1/m2OCwMUgu/xEt1/SfKcutU3wMt37Ockmf+XQHPrH kD4/c3QzSvInASQvs/upNWP51f3aVfupWU0G/q6ATUnrEOc1vfEhtvsq93CY7LtXsf0w 9CyCvzKM54yw1zTFP0Ffs+kMK+njtSoosamErsmPjIWmtC9btAOTtksozazbwWIuhXxp DPTkZtQHVVP88GONgPIxA6BfOWv7ebd5Pg3ZSW/azJ62Lq4KLT49qsRfsKGV+5++9K0F kw8OG1PTqLM7KL4aQj+/kHKrbdmIyk2sBBfCGqIs+RSx8ny8t6wT6wcmr5fySs71yyqY 2UhA==

For some computation I need to compute the say 10 highest eigenvalues
of a symmetric positive definite matrix.

The numerical technique is typically the power-method but this gives only
the highest eigenvector and more work is required for getting more than
one eigenvector.

In Matlab the command of choice is eigs(A, 10) and it is implemented as
well in the "octave" open source software.

Could we have such a functionality in Eigen? That would be fairly natural
for Eigen.

Mathieu

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