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

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

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

revision 2881 by jfenwick, Thu Jan 28 02:03:15 2010 UTC revision 2998 by artak, Wed Mar 24 04:00:20 2010 UTC
# Line 203  void Paso_Solver_solveGS(Paso_Solver_GS Line 203  void Paso_Solver_solveGS(Paso_Solver_GS
203                     }                     }
204                }                }
205             } else if (n_block==2) {             } else if (n_block==2) {
206                #pragma omp parallel for schedule(static) private(i,iptr_ik,k,iptr_main,S1,S2,R1,R2)                #pragma omp parallel for schedule(static) private(i,iptr_ik,k,iptr_main,S1,S2,R1,R2,A11,A21,A12,A22,D,S11,S21,S12,S22)
207                for (i = 0; i < n; ++i) {                for (i = 0; i < n; ++i) {
208                     if (gs->colorOf[i]==color) {                     if (gs->colorOf[i]==color) {
209                       /* x_i=x_i-a_ik*x_k */                       /* x_i=x_i-a_ik*x_k */
# Line 239  void Paso_Solver_solveGS(Paso_Solver_GS Line 239  void Paso_Solver_solveGS(Paso_Solver_GS
239    
240                }                }
241             } else if (n_block==3) {             } else if (n_block==3) {
242                #pragma omp parallel for schedule(static) private(i,iptr_ik,iptr_main,k,S1,S2,S3,R1,R2,R3)                #pragma omp parallel for schedule(static) private(i,iptr_ik,iptr_main,k,S1,S2,S3,R1,R2,R3,A11,A21,A31,A12,A22,A32,A13,A23,A33,D,S11,S21,S31,S12,S22,S32,S13,S23,S33)
243                for (i = 0; i < n; ++i) {                for (i = 0; i < n; ++i) {
244                     if (gs->colorOf[i]==color) {                     if (gs->colorOf[i]==color) {
245                       /* x_i=x_i-a_ik*x_k */                       /* x_i=x_i-a_ik*x_k */
# Line 297  void Paso_Solver_solveGS(Paso_Solver_GS Line 297  void Paso_Solver_solveGS(Paso_Solver_GS
297       /* backward substitution */       /* backward substitution */
298       for (color=(gs->num_colors)-1;color>-1;--color) {       for (color=(gs->num_colors)-1;color>-1;--color) {
299             if (n_block==1) {             if (n_block==1) {
300                #pragma omp parallel for schedule(static) private(i,iptr_ik,k,S1,R1)                /*#pragma omp parallel for schedule(static) private(i,iptr_ik,k,S1,R1,iptr_main)*/
301                for (i = 0; i < n; ++i) {                for (i = 0; i < n; ++i) {
302                     if (gs->colorOf[i]==color) {                     if (gs->colorOf[i]==color) {
303                       /* x_i=x_i-a_ik*x_k */                       /* x_i=x_i-a_ik*x_k */
# Line 315  void Paso_Solver_solveGS(Paso_Solver_GS Line 315  void Paso_Solver_solveGS(Paso_Solver_GS
315                     }                     }
316                }                }
317             } else if (n_block==2) {             } else if (n_block==2) {
318                #pragma omp parallel for schedule(static) private(i,iptr_ik,k,S1,S2,R1,R2)                #pragma omp parallel for schedule(static) private(i,iptr_ik,k,S1,S2,R1,R2,iptr_main,D,A11,A21,A12,A22,S11,S21,S12,S22)
319                for (i = 0; i < n; ++i) {                for (i = 0; i < n; ++i) {
320                     if (gs->colorOf[i]==color) {                     if (gs->colorOf[i]==color) {
321                       /* x_i=x_i-a_ik*x_k */                       /* x_i=x_i-a_ik*x_k */
# Line 353  void Paso_Solver_solveGS(Paso_Solver_GS Line 353  void Paso_Solver_solveGS(Paso_Solver_GS
353                      }                      }
354                }                }
355             } else if (n_block==3) {             } else if (n_block==3) {
356                #pragma omp parallel for schedule(static) private(i,iptr_ik,k,S1,S2,S3,R1,R2,R3)                #pragma omp parallel for schedule(static) private(i,iptr_ik,k,S1,S2,S3,R1,R2,R3,iptr_main,D,A11,A21,A31,A12,A22,A32,A13,A23,A33,S11,S21,S31,S12,S22,S32,S13,S23,S33)
357                for (i = 0; i < n; ++i) {                for (i = 0; i < n; ++i) {
358                     if (gs->colorOf[i]==color) {                     if (gs->colorOf[i]==color) {
359                       /* x_i=x_i-a_ik*x_k */                       /* x_i=x_i-a_ik*x_k */

Legend:
Removed from v.2881  
changed lines
  Added in v.2998

  ViewVC Help
Powered by ViewVC 1.1.26