/[escript]/branches/domexper/paso/src/BiCGStab.c
ViewVC logotype

Diff of /branches/domexper/paso/src/BiCGStab.c

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

revision 3233 by jfenwick, Fri Oct 1 01:53:46 2010 UTC revision 3234 by jfenwick, Mon Oct 4 01:46:30 2010 UTC
# Line 168  err_t Paso_Solver_BiCGStab( Line 168  err_t Paso_Solver_BiCGStab(
168        
169      /*        Compute direction adjusting vector PHAT and scalar ALPHA. */      /*        Compute direction adjusting vector PHAT and scalar ALPHA. */
170        
171          Paso_Solver_solvePreconditioner(A,&phat[0], &p[0]);      Paso_SystemMatrix_solvePreconditioner(A,&phat[0], &p[0]);
172      Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(PASO_ONE, A, &phat[0],PASO_ZERO, &v[0]);      Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(PASO_ONE, A, &phat[0],PASO_ZERO, &v[0]);
173        
174          #pragma omp parallel for private(i0) reduction(+:sum_2) schedule(static)          #pragma omp parallel for private(i0) reduction(+:sum_2) schedule(static)
# Line 200  err_t Paso_Solver_BiCGStab( Line 200  err_t Paso_Solver_BiCGStab(
200           breakFlag = FALSE;           breakFlag = FALSE;
201         } else {         } else {
202           /*           Compute stabilizer vector SHAT and scalar OMEGA. */           /*           Compute stabilizer vector SHAT and scalar OMEGA. */
203               Paso_Solver_solvePreconditioner(A,&shat[0], &s[0]);           Paso_SystemMatrix_solvePreconditioner(A,&shat[0], &s[0]);
204           Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(PASO_ONE, A, &shat[0],PASO_ZERO,&t[0]);           Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(PASO_ONE, A, &shat[0],PASO_ZERO,&t[0]);
205        
206               #pragma omp parallel for private(i0) reduction(+:omegaNumtr,omegaDenumtr) schedule(static)               #pragma omp parallel for private(i0) reduction(+:omegaNumtr,omegaDenumtr) schedule(static)

Legend:
Removed from v.3233  
changed lines
  Added in v.3234

  ViewVC Help
Powered by ViewVC 1.1.26