Questions marquées «banded-matrix»

J'essaie de résoudre une équation de Poisson 2D par des différences finies. Dans le processus, j'obtiens une matrice clairsemée avec seulement variables dans chaque équation. Par exemple, si les variables étaient U , alors la discrétisation donnerait:555UUU Ui−1,j+Ui+1,j−4Ui,j+Uje,j−1+Uje,j+1=fi ,...

Je suis intéressé à résoudre l'équation de Poisson en utilisant l'approche des différences finies. Je voudrais mieux comprendre comment écrire l'équation matricielle avec les conditions aux limites de Neumann. Quelqu'un pourrait-il revoir ce qui suit, est-ce exact? La matrice des différences finies...

Dans mon projet, je dois résoudre quelques matrices tridiagonales à chaque pas de temps, il est donc crucial d'avoir un bon solveur pour celles-ci. J'ai fait ma propre implémentation, juste la manière classique de le faire décrite sur Wikipedia. J'ai ensuite essayé d'utiliser Lapack à la place, et...