/[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

trunk/esys2/paso/src/Solvers/BiCGStab.c revision 150 by jgs, Thu Sep 15 03:44:45 2005 UTC trunk/paso/src/Solvers/BiCGStab.c revision 584 by gross, Thu Mar 9 23:03:38 2006 UTC
# Line 66  err_t Paso_Solver_BiCGStab( Line 66  err_t Paso_Solver_BiCGStab(
66      double * r,      double * r,
67      double * x,      double * x,
68      dim_t *iter,      dim_t *iter,
69      double * tolerance) {      double * tolerance,
70        Paso_Performance* pp) {
71    
72    
73    /* Local variables */    /* Local variables */
# Line 158  err_t Paso_Solver_BiCGStab( Line 159  err_t Paso_Solver_BiCGStab(
159      /*        Compute direction adjusting vector PHAT and scalar ALPHA. */      /*        Compute direction adjusting vector PHAT and scalar ALPHA. */
160        
161          Paso_Solver_solvePreconditioner(A,&phat[0], &p[0]);          Paso_Solver_solvePreconditioner(A,&phat[0], &p[0]);
162      Paso_SystemMatrix_MatrixVector(ONE, A, &phat[0],ZERO, &v[0]);      Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(ONE, A, &phat[0],ZERO, &v[0]);
163        
164          #pragma omp for private(i0) reduction(+:sum_2) schedule(static)          #pragma omp for private(i0) reduction(+:sum_2) schedule(static)
165      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 182  err_t Paso_Solver_BiCGStab( Line 183  err_t Paso_Solver_BiCGStab(
183         } else {         } else {
184           /*           Compute stabilizer vector SHAT and scalar OMEGA. */           /*           Compute stabilizer vector SHAT and scalar OMEGA. */
185               Paso_Solver_solvePreconditioner(A,&shat[0], &s[0]);               Paso_Solver_solvePreconditioner(A,&shat[0], &s[0]);
186           Paso_SystemMatrix_MatrixVector(ONE, A, &shat[0],ZERO,&t[0]);           Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(ONE, A, &shat[0],ZERO,&t[0]);
187        
188               #pragma omp for private(i0) reduction(+:omegaNumtr,omegaDenumtr) schedule(static)               #pragma omp for private(i0) reduction(+:omegaNumtr,omegaDenumtr) schedule(static)
189           for (i0 = 0; i0 < n; i0++) {           for (i0 = 0; i0 < n; i0++) {

Legend:
Removed from v.150  
changed lines
  Added in v.584

  ViewVC Help
Powered by ViewVC 1.1.26