3
Resolvendo um sistema esparso não simétrico e não diagonalmente dominante da melhor maneira
Lembro-me fracamente das minhas primeiras palestras "numéricas" de que os solucionadores lineares iterativos para geralmente exigem que quando seja decomposto comoAx=bAx=bAx=bAAA A=D+MA=D+MA=D + M onde D é uma matriz diagonal e tem diagonal zero, os elementos de devem ser dominantes sobre as entradas em para que os solucionadores iterativos tenham …