Ранкова структура в комплексі Шура


9

Я займаюся дослідженням структури в доповненнях Шура і знаходжу цікаве явище:

Припустимо, що A походить від 5 - pt лаплаціан. Якщо я використовую вкладений впорядкований розсічення та мультифронтальний метод для обчислення LU-факторизації, а потім перевіряю останній блок доповнення Щура, він має низький ранг для позадіагональних блоків.

Але, коли я використовую той же метод для факторизації , де - деяке позитивне значення біля власних значень A, то останній доповнення Щура не має властивості низького рангу.А-λЯλ

Я не знаю, змінить чи невизначеність структуру в доповненнях до Щура чи ні. Чи може хтось надати певну посилання на цю тему?

Відповіді:


8

Ласкаво просимо у чудовий світ рівнянь Гельмгольца. Замінітьλ0 з ω2і ви описуєте факторизацію рівняння Гельмгольца. Вас може зацікавити цей документ , який висвітлює саме цю проблему. Також є приємний оглядовий документ, який пояснює, чому рівняння Гельмгольца важкі.


У статті Ін він показав, що для 2D проблеми доповнення Щура повинно мати властивість низького рангу. Він стверджує лише, що для 3D-задачі властивість низького рангу не є істотною. Моя проблема - це 2D проблема, але вона не має низького рангу.
Willowbrook

@Willowbrook: Я думаю, що ви повинні уважніше прочитати. Властивість низького рангу утримується лише для підпроблем 1d задачі 2D і лише у випадку, коли використовується гранична умова поглинання. Якщо ви введете його у свою формулювання, я думаю, що ваші позадіагональні ранги значно зменшаться, хоча вони все одно значно зростуть із розміром проблеми.
Джек Поульсон
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.