Rangstruktur im Schur-Komplement

9

Ich recherchiere über die Struktur in den Schur-Ergänzungen und finde ein interessantes Phänomen:

Angenommen, A stammt von 5 - Punkt Laplace. Wenn ich zur Berechnung der LU-Faktorisierung eine verschachtelte Dissektionsreihenfolge und eine multifrontale Methode verwende und dann den letzten Schur-Komplementblock überprüfe, hat er einen niedrigen Rang für die nicht diagonalen Blöcke.

Wenn ich jedoch dieselbe Methode verwende, um zu faktorisieren , wobei λ ein positiver Wert nahe den Eigenwerten von A ist, dann hat das letzte Schur-Komplement nicht die Eigenschaft mit niedrigem Rang.EIN- -λichλ

Ich weiß nicht, ob das Unbestimmte die Struktur im Schur-Komplement verändern wird oder nicht. Kann jemand eine Referenz zu diesem Thema geben?

Willowbrook
quelle

Antworten:

8

λ0ω2

Jack Poulson
quelle
In Yings Arbeit zeigte er, dass für ein 2D-Problem das Schur-Komplement die Eigenschaft mit niedrigem Rang haben sollte. Er behauptet nur, dass für das 3D-Problem die niedrigrangige Eigenschaft nicht signifikant ist. Mein Problem ist ein 2D-Problem, aber es hat keinen niedrigen Rang.
Willowbrook
@ Willowbrook: Ich denke, dass Sie es genauer lesen sollten. Es wird argumentiert, dass die Eigenschaft mit niedrigem Rang nur für 1d- Teilprobleme des 2D-Problems gilt, und nur in dem Fall, in dem eine absorbierende Randbedingung verwendet wird. Wenn Sie eine in Ihre Formulierung einführen, werden Ihre nicht diagonalen Ränge meiner Meinung nach erheblich abnehmen, obwohl sie mit der Problemgröße immer noch erheblich wachsen sollten.
Jack Poulson