/[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 1424 by trankine, Fri Jan 11 07:45:58 2008 UTC revision 1425 by artak, Wed Feb 27 06:12:04 2008 UTC
# Line 209  err_t Paso_Solver_BiCGStab( Line 209  err_t Paso_Solver_BiCGStab(
209                  #pragma omp for private(i0) reduction(+:sum_4) schedule(static)                  #pragma omp for private(i0) reduction(+:sum_4) schedule(static)
210              for (i0 = 0; i0 < n; i0++) {              for (i0 = 0; i0 < n; i0++) {
211                x[i0] += alpha * phat[i0] + omega * shat[i0];                x[i0] += alpha * phat[i0] + omega * shat[i0];
212                r[i0] -= omega * t[i0];                r[i0] = s[i0]-omega * t[i0];
213                sum_4 += r[i0] * r[i0];                sum_4 += r[i0] * r[i0];
214              }              }
215              norm_of_residual = sqrt(sum_4);              norm_of_residual = sqrt(sum_4);
Legend:
Removed from v.1424  
changed lines
  Added in v.1425
  ViewVC Help
Powered by ViewVC 1.1.26