Perguntas com a marcação «preconditioning»

O que pode dar errado ao usar métodos Krylov pré-condicionados da KSP ( pacote de solucionadores lineares do PETSc ) para resolver um sistema linear esparso, como os obtidos pela discretização e linearização de equações diferenciais parciais? Que etapas posso tomar para determinar o que está...

Para a solução de grandes sistemas lineares usando métodos iterativos, é frequentemente interessante introduzir o pré-condicionamento, por exemplo, resolva M - 1 ( A x = b ) , onde M é usado aqui para o pré-condicionamento à esquerda do sistema. Normalmente, devemos ter esse M - 1 ≈ A - 1 e...

Isso é direcionado principalmente para PDEs elípticos sobre domínios convexos, para que eu possa ter uma boa visão geral dos dois

Estou muito impressionado com o desempenho serial de pré-condicionadores de ILU de nível inverso multinível , particularmente para Helmholtz heterogêneo , mas estou surpreso por não conseguir encontrar nenhuma implementação de código aberto. Em particular, o ILUPACK disponibiliza binários...

Eu me pergunto o que aconteceu com os pré-condicionadores polinomiais. Estou interessado neles, porque eles parecem ser comparativamente elegantes do ponto de vista matemático, mas, tanto quanto eu li em pesquisas sobre métodos de krylov, eles geralmente têm um desempenho muito ruim como...

Questão: Suponha que você tenha dois pré-condicionadores diferentes (fatorados) para uma matriz definida positiva simétrica : e onde os inversos dos fatores são fácil de aplicar.A ≈ B T B A ≈ C T C , B , B T , C , C TAAAA≈BTBA≈BTBA \approx B^TBA ≈ CTC,A≈CTC,A \approx C^TC,B , BT, C,...

Em métodos como gmres ou bicgstab, pode ser atraente usar outro método de krylov como pré-condicionador. Afinal, eles são fáceis de implementar de maneira livre de matriz e em um ambiente paralelo. Por exemplo, um exemplo pode usar algumas (digamos ~ 5) iterações da bigcstab não condicionada como...

Eu quero resolver um problema não linear com restrições de igualdade não lineares e estou usando um Lagrangiano aumentado com um termo de regularização de penalidade que, como se sabe, estraga o número de condições dos meus sistemas linearizados (a cada iteração de Newton, quero dizer) . Quanto...

Meu sistema é um problema simétrico de FE com multiplicadores de lagrange (por exemplo, fluxo incompressível de Stokes): ( ABBTC)(ABTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} onde é o caso típico (eu até me certifiquei de que as equações sejam numeradas para que os multiplicadores de...

Os métodos Newton-Krylov livre de Jacobiano (JFNK) e métodos de Krylov em geral podem ser muito úteis porque não requerem armazenamento explícito ou construção de uma matriz, apenas os resultados de produtos vetoriais de matriz. Se você realmente formar o sistema esparso, existem muitos...

Eu tenho um sistema linear com matriz cujos valores próprios são distribuídos uniformemente no círculo unitário como este: É possível resolver esse tipo de sistema de maneira eficaz pelo método iterativo, talvez com algum

Considere a situação em que você deseja resolver um sistema linear usando um método pré-condicionado de Krylov, mas a aplicação do pré-condicionador envolve a solução de um sistema auxiliar, o que é feito com outro método pré-condicionado de Krylov. Em um extremo, você pode executar a resolução...

Em esta discussão anterior da seguinte maneira multiplicativo para combinar precondicionadores simétricas e P 2 para o sistema simétrico Um x = b foi sugerido: P - 1 combinação : = P - 1 1 + P - 1 2 ( I - A P - 1 1 ) = P - 1 1 + P - 1 2 - P -P1P1P_1P2P2P_2A x = bAx=bAx=bP- 1combinação: ==P- 11+ P-...

Eu tenho um problema generalizado de matriz: de um método espectral em um problema de análise de estabilidade linear. Minha matriz B é diagonal e positiva semi-definida. A é não-eremita e complexo.Ax=λBxAx=λBxA x = \lambda B x Meu problema é essencialmente que, ao usar o solucionador de autovalor...

Qual é uma boa maneira de extrair a diagonal de uma matriz simétrica que já é quase diagonal, mas onde você não possui os elementos da matriz (apenas a capacidade de aplicá-la a vetores)? Outras restrições são: (1) aplicar a matriz n-por-n n-vezes para construir explicitamente a diagonal seria...

Implementei um Navier Stokes baseado na discretização FEM e PETSc para resolver o sistema linear de equações. Para criar um procedimento de solução eficiente, sigo o artigo "Pré-condicionamento eficiente das equações de Navier-Stokes linearizadas para fluxo incompressível" (Silvester et al.), Que...

Aprendi que para o solucionador BiCGStab para sistemas lineares esparsos, é quase sempre necessário usar um pré-condicionador. Eu já percebi que escolher um bom depende do problema. Navegando na web, descobri que existem muitos pré-condicionadores baseados em ILU (ILUT, MILU, etc.) Foi aqui que...

Quando penso em resolver equações não lineares, geralmente penso em linearizar primeiro e depois aplicar um pré-condicionador à matriz linear. Ocorreu-me que poderia ser possível condicionar as equações não lineares ANTES de serem linearizadas. Existem estratégias gerais para pré-condicionar...