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

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

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

revision 1973 by ksteube, Thu Sep 25 23:11:13 2008 UTC revision 1974 by jfenwick, Thu Nov 6 02:40:10 2008 UTC
# Line 169  err_t Paso_Solver_BiCGStab( Line 169  err_t Paso_Solver_BiCGStab(
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_Solver_solvePreconditioner(A,&phat[0], &p[0]);
172      Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(ONE, A, &phat[0],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)
175      for (i0 = 0; i0 < n; i0++) sum_2 += rtld[i0] * v[i0];      for (i0 = 0; i0 < n; i0++) sum_2 += rtld[i0] * v[i0];
# Line 201  err_t Paso_Solver_BiCGStab( Line 201  err_t Paso_Solver_BiCGStab(
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_Solver_solvePreconditioner(A,&shat[0], &s[0]);
204           Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(ONE, A, &shat[0],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)
207           for (i0 = 0; i0 < n; i0++) {           for (i0 = 0; i0 < n; i0++) {

Legend:
Removed from v.1973  
changed lines
  Added in v.1974

  ViewVC Help
Powered by ViewVC 1.1.26