FORTRAN90: Matrizen-System loesen (LAPACK)

[zahor]

joa, ueber SVD koennte man es auch berechnen. Danke fuer die Idee! Ich hab mal gegoogelt, wie ich das mit LAPACK anstellen muesste, und bin auch fuendig geworden. Das wuerde dann aber ziemlich umstaendlich:
1. ZGEBRD aufrufen, um die Matrix in bidiagonale form zu bringen
2. ZUNGBR aufrufen (2 mal), um Q sowie P^H (P**H) zu kriegen
3. Dann das ganze durch ZBDSQR jagen um die Singulaerwerte zu erhalten
4. Letzendlich das Pseudinverse mithilfe von ZGELSS und den Singulaerwerten erhalten
5. Dieses Ergebnis mithilfe von MATMUL mit der anderen Matrix multiplizieren.
Da wird man doch doof davon, vor allem weil das ja ungemein viel Zeit beansprucht! Da wird ja alles moegliche berechnet, was ich weiter dann gar nicht brauche. Da ist ZGESV wohl doch schneller. Der Fakt, das MATLAB es fuer seine Berechnungen benutzt, will ja doch was heissen. Die werden ja alle Varianten durchdacht haben.
Ich wuerde schon lieber bei ZGESV oder notfalls ZGERTF / ZGETRI / ZGETRS bleiben... Danke trotzdem
Meine Frage war eher so gedacht: Wie mache ich ZGESV meine Matrizen schmackhaft?