Estou tentando resolver uma equação de Poisson 2D por diferenças finitas. No processo, obtenho uma matriz esparsa com apenas variáveis em cada equação. Por exemplo, se as variáveis fossem U , a discretização
Estou interessado em resolver a equação de Poisson usando a abordagem de diferenças finitas. Eu gostaria de entender melhor como escrever a equação da matriz com as condições de contorno de Neumann. Alguém revisaria o seguinte, está correto? A matriz de diferenças finitas A equação de...
No meu projeto, tenho que resolver duas matrizes tridiagonais a cada passo do tempo, por isso é crucial ter um bom solucionador para elas. Fiz minha própria implementação, exatamente da maneira clássica de fazê-lo, descrita na Wikipedia. Então tentei usar o Lapack e, para minha surpresa, foi mais...