Spezialisierte Methoden für komplexe symmetrische tridiagonale verallgemeinerte Eigenwertprobleme

Ich muss verallgemeinerte Eigenwertprobleme lösen, wobei A und B beide tridiagonal sind, B symmetrisch positiv definit und real ist, aber A nur komplex symmetrisch ist (nicht definit oder hermitisch). Außerdem brauche ich die volle eigendecomposition. Ich rufe gerade Lapacks verallgemeinerten Eigensolver an, aber ich frage mich, ob es für dieses spezielle, stark strukturierte Problem bessere Methoden gibt. Insbesondere wäre frei verfügbarer Code (C ++) das Beste.$Ax = \lambda Bx$$A$$B$$B$$A$ZGGEV

Der Pole Expansions- und Selected Inversion ( PEXSI ) Methode könnte die Antwort sein. Ich habe diese Methode nicht verwendet, aber sie bietet eine Inversionsroutine für komplexe symmetrische Matrizen. Es ist nicht spezifisch für tridiagonale Matrizen, nutzt aber die Sparsamkeit.