Estrutura de classificação no complemento Schur


9

Estou pesquisando a estrutura dos complementos de Schur e encontro um fenômeno interessante:

Suponha que A seja de 5 - pt laplaciano. Se eu usar a ordenação de dissecações aninhadas e o método multifrontal para calcular a fatoração da LU e depois verificar o último bloco do complemento de schur, ele tem uma classificação baixa para os blocos fora da diagonal.

Mas, quando eu uso o mesmo método para fatorar , onde é algum valor positivo próximo aos valores próprios de A, então o último complemento schur não possui a propriedade low-rank.UMA-λEuλ

Não sei se o indefinido mudará a estrutura no complemento schur ou não. Alguém pode fornecer alguma referência para este tópico?

Respostas:


8

Bem-vindo ao maravilhoso mundo das equações de Helmholtz. Substitua por e você está descrevendo uma fatoração da equação de Helmholtz. Você pode estar interessado neste artigo , que aborda exatamente esse problema. Há também um bom artigo de revisão que explica por que as equações de Helmholtz são difíceis.λ0 0ω2


No artigo de Ying, ele mostrou que, para o problema 2D, o complemento schur deve ter a propriedade de baixo escalão. Ele apenas afirma que, para o problema 3D, a propriedade de baixa classificação não é significativa. Meu problema é um problema 2D, mas ele não tem uma classificação baixa.
precisa

@ Willowbrook: Eu acho que você deve dar uma leitura mais cuidadosa. Argumenta-se que a propriedade de baixo escalão é válida apenas para os subproblemas 1d do problema 2D, e somente no caso em que uma condição de contorno absorvente é usada. Se você introduzir um em sua formulação, acho que suas fileiras fora da diagonal diminuirão significativamente, embora ainda devam crescer significativamente com o tamanho do problema.
22612 Jack Poulson
Ao utilizar nosso site, você reconhece que leu e compreendeu nossa Política de Cookies e nossa Política de Privacidade.
Licensed under cc by-sa 3.0 with attribution required.