CGS, un solucionador rápido tipo Lanczos para sistemas lineales no simétricos

Se presenta un método tipo Lanczos para sistemas lineales dispersos no simétricos que surgen de discretizaciones de ecuaciones diferenciales parciales elípticas. El método se basa en una variante polinómica del algoritmo de gradientes conjugados. Aunque está relacionado con el llamado algoritmo de gradientes bi-conjugados (Bi-CG), no implica multiplicaciones de matrices-vectores adjuntos, y la tasa de convergencia esperada es aproximadamente el doble que la del algoritmo Bi-CG. Se realiza una comparación numérica con otros solucionadores, probando el método en una familia de ecuaciones de convección-difusión, en varias rejillas, y con el uso de dos métodos de preacondicionamiento diferentes. Se utilizan diferencias hacia adelante, así como diferencias centrales, en los experimentos.