On Tue, Jul 23, 2013 at 10:47 AM, Manav Bhatia <[email protected]>wrote:
> > > On Tue, Jul 23, 2013 at 12:42 PM, John Peterson <[email protected]>wrote: > >> >> Manav, >> >> There is some good stuff in here, particularly >> the _nonhermitian_eig_lapack stuff. >> >> I'd like to get it into libmesh, do you have any test codes that verify >> it is working? >> >> -- >> John >> > > I am not sure if I have verified this code. I had started to add this a > few months ago, but then it became unnecessary and I never got around to > veryfying the code. > I can test it at my end, and send you an updated patch if needed? > That would be good, but it would be even better if we had a very short libmesh code that also demonstrated it working (not a full application code). It looks like it's designed to work when Number=Real, but I wonder if you've tested that it compiles with Number=std::complex? Also, what do you envision the user-facing interface for _nonhermitian_eig_lapack should be? For example, lu_solve, cholesky_solve, and svd call some underlying Lapack routines. DenseMatrix::eig() maybe? -- John ------------------------------------------------------------------------------ See everything from the browser to the database with AppDynamics Get end-to-end visibility with application monitoring from AppDynamics Isolate bottlenecks and diagnose root cause in seconds. Start your free trial of AppDynamics Pro today! http://pubads.g.doubleclick.net/gampad/clk?id=48808831&iu=/4140/ostg.clktrk _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
