/[escript]/trunk/paso/src/Solver_GS.c
ViewVC logotype

Diff of /trunk/paso/src/Solver_GS.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1841 by artak, Wed Oct 1 05:56:05 2008 UTC revision 1842 by artak, Fri Oct 3 04:18:26 2008 UTC
# Line 405  void Paso_Solver_solveGS(Paso_Solver_GS Line 405  void Paso_Solver_solveGS(Paso_Solver_GS
405           }           }
406       }       }
407            
      if (gs->sweeps>1) {  
      /* Compute the residual b=b-Ax*/  
      Paso_SparseMatrix_MatrixVector_CSR_OFFSET0(DBLE(-1), gs->factors, x, DBLE(2), b);  
      /* Go round again*/  
      gs->sweeps=gs->sweeps-1;  
      Paso_Solver_solveGS(gs,x,b);  
      }  
       
408       return;       return;
409  }  }
410    

Legend:
Removed from v.1841  
changed lines
  Added in v.1842

  ViewVC Help
Powered by ViewVC 1.1.26