Gauss-Seidel"/>

Por fin el multigrid??

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 […]

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 […]