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
