Hi Peter, good luck. I completely agree with Rasmus, it really depends on your problem, what the best/fastest/smallest solution is. I am sorry to hear that the Spectra thing did not work.
As an intermediate solution, you can also just adapt your current lapack wrapper to work with eigen types. It is pretty straightforward to convert from ublas. Cheers Jens | Jens Wehner, PhD | eScience Research Engineer | Email: [email protected]<https://outlook.office.com/mail/options/mail/[email protected]> | Tel: +31(0)6 438 666 87 | | Netherlands eScience Center<https://www.esciencecenter.nl/> | Science Park 140 | 1098 XG Amsterdam | The Netherlands | | Twitter<https://twitter.com/eScienceCenter> | LinkedIn<https://www.linkedin.com/company/netherlands-escience-center> | Facebook<https://www.facebook.com/NLeScienceCenter/> | YouTube<https://www.youtube.com/channel/UCYLAIMi62d8tx3Ru6DOSPeg> | Newsletter<https://esciencecenter.us8.list-manage.com/subscribe/post?u=a0a563ca342f1949246a9f92f&id=31bfc2303d> | ________________________________ From: Peter <[email protected]> Sent: Thursday, January 21, 2021 9:42 To: [email protected] <[email protected]> Subject: Re: [eigen] lowest eigen value only Dear Jens, >> if your system is reasonably large iterative methods could be what you are >> looking for >> >> Libraries like https://github.com/yixuan/spectra >> <https://github.com/yixuan/spectra> implement such things on top of Eigen. Sadly it gives me a segfault at eig.init(), except for tiny matrices, although I just followed the simple example based on a dense matrix. I'll try write a bug report for spectra later, if I can nail it down to a simple example. Best regards, Peter
