Stability of pdsygst versus pdsyngst

Post here if you have a question about LAPACK or ScaLAPACK algorithm or data format

Stability of pdsygst versus pdsyngst

Postby ldm001 » Fri Feb 15, 2019 4:24 pm

I am working with some code which uses pdsygst to reduce general eigenproblems, with 'U'. I recently switched it to use 'L' with pdsyngst. Beyond the well established speed improvement of pdsyngst, I have provisional evidence [1] that pdsyngst is more stable. Is this something that is expected?



[1] The code I am using employs a simplex derivative method to find a fixed-point solution to a long sequence of iterations, each of which involves solving an eigenproblem. These problems are susceptible to numerical noise (instabilities) even with double precision. They appear to be converging better with pdsyngst. (I can add more details than I expect anyone will want to know, so won't unless asked.)
ldm001
 
Posts: 9
Joined: Tue Sep 16, 2008 10:21 pm

Return to Algorithm / Data

Who is online

Users browsing this forum: No registered users and 5 guests