Benchmark-Probleme für Eigenwert-Neuordnungsalgorithmen gesucht

10

Jede reelle Matrix kann unter Verwendung einer orthogonalen Ähnlichkeitstransformation U auf die reelle Schurform T = U T A U reduziert werden . Hier ist die Matrix T quasi dreieckig mit 1 mal 1 oder 2 mal 2 Blöcken auf der Hauptdiagonale. Jeweils 1 x 1 Block entspricht einem realen Eigenwert von A und jeweils 2 x 2 Block entspricht einem Paar von konjugiert - komplexen Eigenwerte von A .AT=UTAUUTAA

Das Problem der Neuwertumordnung besteht darin, eine orthogonale Ähnlichkeitstransformation so dass die Auswahl der Eigenwerte von A durch den Benutzer entlang der Diagonale der oberen linken Ecke von S = V T T V erscheint .VAS=VTTV

In LAPACK heißt die relevante Routine mit doppelter Genauigkeit DTRSEN. Daniel Kressner hat eine blockierte Version mit dem Namen BDTRSEN geschrieben. Die ScaLAPACK-Routine ist PDTRSEN.

Ich suche nach Anwendungen und Algorithmen, bei denen Fortschritte bei der Lösung des Problems der Eigenwertumordnung echte Vorteile haben.

Wir können leicht Testmatrizen in quasi dreieckiger Form erzeugen, aber wir haben Schwierigkeiten, die Form einer realistischen Verteilung der Eigenwertauswahl des Benutzers zu bestimmen.

Aus meiner Sicht ist die Subraumiteration mit Ritz-Beschleunigung ein idealer Algorithmus zum Testen von Verbesserungen des Neuordnungsalgorithmus. Es benötigt eine (spärliche) Matrixvektormultiplikation, einen hohen QR-Algorithmus und einen Neuordnungsalgorithmus.

Es fällt mir jedoch schwer, Probleme im wirklichen Leben zu finden, bei denen klar ist, dass eine bestimmte Menge von Eigenpaaren physikalisch interessant ist.

Wir können eine Eigenwertumordnung für dichte Matrizen der Dimension 40.000 unter Verwendung einer gemeinsam genutzten Speichermaschine durchführen. Die beste Leistung wird erzielt, wenn der Benutzer etwa 50% aller Eigenwerte auswählt.

Carl Christian
quelle

Antworten:

0

Ich bin mir sicher, dass ich die Nützlichkeit des Eigenwert-Neuordnungsalgorithmus nicht richtig einschätze, aber für diesen Teil Ihrer Frage fallen mir viele Antworten ein:

Es fällt mir jedoch schwer, Probleme im wirklichen Leben zu finden, bei denen klar ist, dass eine bestimmte Menge von Eigenpaaren physikalisch interessant ist.

Beispielsweise haben Sie bei bestimmten hydrodynamischen Stabilitätsproblemen eindeutige Eigenwerte, die mit physikalischen Phänomenen wie Kelvin-Helmholtz-Modi oder Tollmien-Schlichting-Wellen verbunden sind. Bei Fluid-Struktur-Wechselwirkungsproblemen kann der Resonanzmodus mit einer Schlaginstabilität verbunden sein.

Entspricht dies dem, wonach Sie suchen? Wenn ja, bin ich sicher, dass andere mit Beispielen aus ihren Bereichen zwitschern werden; Wenn nicht, können Sie die Frage schärfen?

Spencer Bryngelson
quelle
Es ist sehr spät in meiner Zeitzone, ich werde antworten, wenn ich geschlafen habe.
Carl Christian
Verzeih mir, es erwies sich als weniger trivial, die Frage zu schärfen, und andere Dinge lenkten mich ab. Schließlich kann ich auf die Sache zurückkommen.
Carl Christian