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

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

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

revision 2993 by artak, Tue Mar 23 04:29:30 2010 UTC revision 3000 by lgao, Tue Mar 30 03:10:13 2010 UTC
# Line 210  void Paso_Solver_solvePreconditioner(Pas Line 210  void Paso_Solver_solvePreconditioner(Pas
210                     #pragma omp parallel for private(i) schedule(static)                     #pragma omp parallel for private(i) schedule(static)
211                     for (i=0;i<prec->gs->n*prec->gs->n_block;++i) bnew[i]+=b[i];                     for (i=0;i<prec->gs->n*prec->gs->n_block;++i) bnew[i]+=b[i];
212                      /* Compute the residual b=b-Ax*/                      /* Compute the residual b=b-Ax*/
213                     Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(DBLE(-1), A, x, DBLE(1), bnew);                     Paso_SparseMatrix_MatrixVector_CSR_OFFSET0(DBLE(-1), A->mainBlock, x, DBLE(1), bnew);
214                     /* Go round again*/                     /* Go round again*/
215                     Paso_Solver_solveGS(prec->gs,x,bnew);                     Paso_Solver_solveGS(prec->gs,x,bnew);
216                     sweeps=sweeps-1;                     sweeps=sweeps-1;

Legend:
Removed from v.2993  
changed lines
  Added in v.3000

  ViewVC Help
Powered by ViewVC 1.1.26