A follow-up comment on this: We added lapack-based SVD in DenseMatrix before we included Eigen in contrib, and I don't think it was ever used very extensively (and it seems to be buggy, as shown by Jonas's example).
I find the lapack API to be quite opaque, so the less we rely on it the better IMO. So now that we include Eigen, I would be in favor of deprecating the lapack-based SVD and EVD in DenseMatrix and suggesting that users use Eigen instead, as Jonas did in his example. I'm not sure that we need a wrapper for Eigen's SVD or EVD in DenseMatrix --- it's just as easy (and more general) for users to call Eigen themselves however they need to. David On Mon, Nov 23, 2015 at 8:36 PM, John Peterson <jwpeter...@gmail.com> wrote: > > A quick thought is that we may be computing the SVD of A^T, since we have > to treat A as though it is column-major to interface with Fortran. In that > case, > > A^T = V S U^T > > that is, we need to treat the "U" handed back to us as "V" and the "V" as > "U". > > I'm out for the next week on vacation but will check this when I get back. > > ------------------------------------------------------------------------------ > Go from Idea to Many App Stores Faster with Intel(R) XDK > Give your users amazing mobile app experiences with Intel(R) XDK. > Use one codebase in this all-in-one HTML5 development environment. > Design, debug & build mobile apps & 2D/3D high-impact games for multiple > OSs. > http://pubads.g.doubleclick.net/gampad/clk?id=254741551&iu=/4140 > _______________________________________________ > Libmesh-users mailing list > Libmesh-users@lists.sourceforge.net > https://lists.sourceforge.net/lists/listinfo/libmesh-users > ------------------------------------------------------------------------------ Go from Idea to Many App Stores Faster with Intel(R) XDK Give your users amazing mobile app experiences with Intel(R) XDK. Use one codebase in this all-in-one HTML5 development environment. Design, debug & build mobile apps & 2D/3D high-impact games for multiple OSs. http://pubads.g.doubleclick.net/gampad/clk?id=254741551&iu=/4140 _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
