Qu'est-ce qui peut mal tourner lors de l'utilisation de méthodes Krylov précondonifiées du KSP ( package de solveur linéaire de PETSc ) pour résoudre un système linéaire clairsemé tel que ceux obtenus en discrétisant et en linéarisant des équations différentielles partielles? Quelles mesures puis-je prendre pour déterminer ce qui ne …
Si je comprends bien, il existe deux grandes catégories de méthodes itératives pour résoudre des systèmes linéaires d'équations: Méthodes stationnaires (Jacobi, Gauss-Seidel, SOR, Multigrid) Méthodes Krylov Subspace (Gradient Conjugué, GMRES, etc.) Je comprends que la plupart des méthodes stationnaires fonctionnent en relaxant de manière itérative (lissage) les modes de Fourier …
J'ai des matrices et . est clairsemé et est avec très grand (peut être de l'ordre de plusieurs millions.) est une matrice haut avec plutôt petit ( ) et chaque colonne peut seulement un seul entrée avec le reste étant « s, de telle sorte que . est énorme, il …
Question: Supposons que vous ayez deux préconditionneurs différents (factorisés) pour une matrice définie positive symétrique : et où les inverses des facteurs sont facile à appliquer.A ≈ B T BUNEUNEAA ≈ BTBUNE≈BTBA \approx B^TBB , B T , C , C TA ≈ CTC,UNE≈CTC,A \approx C^TC,B , BT, C, CTB,BT,C,CTB, …
Je me demande ce qui est arrivé aux préconditionneurs polynomiaux. Je m'intéresse à eux, car ils semblent relativement élégants d'un point de vue mathématique, mais d'après ce que j'ai lu dans les enquêtes sur les méthodes de krylov, ils fonctionnent généralement très mal en tant que préconditionneurs. Pour reprendre les …
Pour autant que je sache, les solveurs multigrilles utilisent des lisseurs itératifs tels que Jacobi, Gauss-Seidel et SOR pour atténuer l'erreur à différentes fréquences. Une méthode de sous-espace de Krylov (comme le gradient conjugué, GMRES, etc.) pourrait-elle être utilisée à la place? Je ne pense pas qu'ils soient classés comme …
Disons que j'ai un système linéaire , qui converge rapidement en utilisant une méthode de Krylov appropriée (comme CG ou GMRES) pour tout . Si est une matrice de faible rang , la même méthode de Krylov sur le système convergera-t-elle également rapidement (idéalement avec un nombre supplémentaire d'itérations qui …
Multigrid (MG) peut être utilisé pour résoudre un système linéaire en construisant une supposition initiale et en répétant ce qui suit pour jusqu'à la convergence:A x = bUNEX=bAx=bX0X0x_0i = 0 , 1 ..je=0,1..i=0,1.. Calculer lerje= b - A xjerje=b-UNEXjer_i = b-Ax_i Appliquer un cycle multigrille pour obtenir une approximation , …
Dans une méthode comme le gmres ou le bicgstab, il pourrait être intéressant d'utiliser une autre méthode krylov comme préconditionneur. Après tout, ils sont faciles à mettre en œuvre sans matrice et dans un environnement parallèle. Par exemple, on peut utiliser quelques (disons ~ 5) itérations de bigcstab non conditionné …
Une équation de Poisson avec toutes les conditions aux limites de Neumann a un seul espace nul dimensionnel constant. Lors de la résolution via une méthode de Krylov, l'espace nul peut être supprimé soit en soustrayant la moyenne de la solution à chaque itération, soit en épinglant la valeur d'un …
Mon système est un problème FE symétrique avec des multiplicateurs de lagrange (par exemple le flux de Stokes incompressible): ( ABBTC)(ABTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} où est le cas typique (je me suis même assuré que les équations sont numérotées pour que les multiplicateurs de Lagrange apparaissent en …
Je résout pour une énorme matrice définie positive clairsemée utilisant la méthode du gradient conjugué (CG). Est-il possible de calculer le déterminant de utilisant les informations produites lors de la résolution?A AA x = bAx=bAx=bUNEAAUNEAA
En science informatique, nous rencontrons souvent de grands systèmes linéaires que nous devons résoudre par des moyens (efficaces), par exemple par des méthodes directes ou itératives. Si nous nous concentrons sur ces derniers, comment pouvons-nous établir qu'une méthode itérative pour résoudre un grand système linéaire est convergente dans la pratique? …
Considérez la situation où vous souhaitez résoudre un système linéaire à l'aide d'une méthode Krylov préconditionnée, mais l'application du préconditionneur lui-même implique la résolution d'un système auxiliaire, ce qui se fait avec une autre méthode Krylov préconditionnée. À un extrême, vous pouvez exécuter la résolution interne pour qu'elle converge à …
Je m'intéresse au calcul de la solution d'un système de lage d'OD en utilisant une méthode krylov comme dans [1]. Une telle méthode implique des fonctions liées à l'exponentielle (les fonctions dites ). Elle consiste essentiellement à calculer l'action de la fonction matricielle en construisant un sous-espace de Krylov en …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.