Re: [julia-users] Schurfact ordered eigenvalues

2014-09-19 Thread Jiahao Chen
This is a tricky cfunction to write and I don't really know the right answer. (I hope others will correct me if I get anything wrong.) But as I understand it: 1. The LAPACK calls Julia uses are actually to the CLAPACK wrappers. BLAS.gges! calls the CLAPACK function ?gges and you want to pass it a

Re: [julia-users] Schurfact ordered eigenvalues

2014-09-19 Thread Chase Coleman
Thanks for the response. I'm glad that it isn't something trivial that I missed (although unfortunately that means that following this response will be the beginning of a potentiall long guessing game). If I make progress, I will provide an update. On Friday, September 19, 2014 1:44:26 PM

Re: [julia-users] Schurfact ordered eigenvalues

2014-09-19 Thread Jiahao Chen
You could also try specifying Ptr{Void} as the output type to see if it works.

Re: [julia-users] Schurfact ordered eigenvalues

2014-09-19 Thread Jameson Nash
http://lapackpp.sourceforge.net/html/f2c_8h.html#05fcfc78d8ede55ad78889ef2cfd6403 It looks like L_fp is a Ptr{Int} On Fri, Sep 19, 2014 at 6:00 PM, Jiahao Chen jia...@mit.edu wrote: You could also try specifying Ptr{Void} as the output type to see if it works. ​

[julia-users] Schurfact ordered eigenvalues

2014-09-18 Thread Chase Coleman
Hi, The function (`?gges`) being called in `schurfact` from LAPACK accepts an argument for sorting the eigenvalues according to a logical/boolean function. I am currently working on some code that would be valuable for me to sort eigenvalues into explosive and nonexplosive (modulus bounded by