I prefer the present solution where we use dispatch instead of keywords. I think it is transparent what is happening: real matrix -> real schur form, complex matrix -> complex schur form. However, we could add a line in the documentation explaining it. If you are okay with that idea, please open a pr.
2014-04-17 6:58 GMT+02:00 Joey Huchette <[email protected]>: > Ahh simple enough, thanks. It might be nice to have a keyword argument > that just dispatches on this---it's not completely obvious that's the right > thing to do unless you dig around lapack.jl (or get a hint). I can PR. > > > On Thursday, April 17, 2014 12:18:00 AM UTC-4, Andreas Noack Jensen wrote: > >> The trick is to convert your matrix to complex before the calculation. >> >> schurfact(complex(A))[:T] >> >> gives the triangular part. >> >> >> 2014-04-17 1:18 GMT+02:00 Joey Huchette <[email protected]>: >> >> Is there an implementation in Base (or elsewhere) of a Schur >>> decomposition that returns a complex matrix matrix T that is triangular? >>> For reference, MATLAB has a optional switch between the two forms. I didn't >>> do enough digging to see if this option is exposed by Lapack, so maybe the >>> conversion could be done at the Julia level? Apologies if I overlooked it >>> in the docs/source. >>> >> >> >> >> -- >> Med venlig hilsen >> >> Andreas Noack Jensen >> > -- Med venlig hilsen Andreas Noack Jensen
