Llego ya un tiempo implementando los diferentes esquemas multigrid en C++ y, tras un periodo de larga oscuridad, por fin aparece un poco de luz: ayer , por primera vez, pareció funcionar todo (en una dimensión y con fronteras Dirichlet, claro… Aunque todo lo demás ya está también implementado). En el primer nivel necesitamos un …
Archivos de la etiqueta: Gauss-Seidel
Métodos iterativos para la resolución de sistemas
Ya comentamos en el post anterior que formato tienen los sistemas $latex Au=f$ que surgen a partir de la discretización de un problema de Poisson y que tendremos que resolver. Esta resolución la podriamos hacer de manera directa pero son mucho mas eficaces los métodos iterativos. Se puede demostrar que la manera mas sencilla para …
Seguir leyendo «Métodos iterativos para la resolución de sistemas»