[SciPy-user] qz decomposition how to?

Sven Schreiber svetosch at gmx.net
Fri Nov 11 10:10:41 CST 2005


Nils Wagner schrieb:

> 
> Use
> 
> help (linalg.flapack)
> 
>       t,sdim,w,vs,work,info =
> cgees(cselect,a,compute_v=1,sort_t=0,lwork=3*n,cselect_extra_args=(),overwrite_a=0)
>       t,sdim,w,vs,work,info =
> zgees(zselect,a,compute_v=1,sort_t=0,lwork=3*n,zselect_extra_args=(),overwrite_a=0)
>       t,sdim,wr,wi,vs,work,info =
> dgees(dselect,a,compute_v=1,sort_t=0,lwork=3*n,dselect_extra_args=(),overwrite_a=0)
>       t,sdim,wr,wi,vs,work,info =
> sgees(sselect,a,compute_v=1,sort_t=0,lwork=3*n,sselect_extra_args=(),overwrite_a=0)
>       alphar,alphai,beta,vl,vr,work,info =
> sggev(a,b,compute_vl=1,compute_vr=1,lwork=8*n,overwrite_a=0,overwrite_b=0)
>       alphar,alphai,beta,vl,vr,work,info =
> dggev(a,b,compute_vl=1,compute_vr=1,lwork=8*n,overwrite_a=0,overwrite_b=0)
>       alpha,beta,vl,vr,work,info =
> cggev(a,b,compute_vl=1,compute_vr=1,lwork=8*n,overwrite_a=0,overwrite_b=0)
>       alpha,beta,vl,vr,work,info =
> zggev(a,b,compute_vl=1,compute_vr=1,lwork=8*n,overwrite_a=0,overwrite_b=0)
> 
> See also linalg.schur, linalg.rsf2csf and linalg.eig
> 
> Nils
>  
> 

That was a quick reply, and it seems to work. Thank you!
> _______________________________________________
> SciPy-user mailing list
> SciPy-user at scipy.net
> http://www.scipy.net/mailman/listinfo/scipy-user
> 
> 



More information about the SciPy-user mailing list