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] <javascript:>> > : > >> 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 >
