Patrick Alken writes: > I have run my code through Brian Gough's test program > on 12x12 and 200x200 random matrices, as well as some integer > matrix cases and find no significant disagreement with LAPACK. > ..... > I propose adding the double-shift method to gsl (assuming it > passes any further tests needed) which would be a perfectly > fine eigenvalue solver (EISPACK used it for a number of years) > until I (or someone else) can find the time to implement method 3.
Excellent work! Useful summary. I will add the code at the next opportunity. -- Brian Gough
