/[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 3558 by artak, Mon Dec 13 23:15:01 2010 UTC revision 3559 by jfenwick, Fri Aug 26 07:47:27 2011 UTC
# Line 183  void Paso_Preconditioner_LocalSmoother_S Line 183  void Paso_Preconditioner_LocalSmoother_S
183     const dim_t n=A_p->numRows;     const dim_t n=A_p->numRows;
184     const dim_t n_block=A_p->row_block_size;     const dim_t n_block=A_p->row_block_size;
185     double *diag = smoother->diag;     double *diag = smoother->diag;
186     index_t* pivot = smoother->pivot;     index_t* pivot = smoother->pivot;
187     const dim_t block_len=A_p->block_size;     const dim_t block_len=A_p->block_size;
188    
189        
190     register dim_t i,k;     register dim_t i,k;
191     register index_t iptr_ik, mm;     register index_t iptr_ik, mm;
# Line 192  void Paso_Preconditioner_LocalSmoother_S Line 193  void Paso_Preconditioner_LocalSmoother_S
193     int failed = 0;     int failed = 0;
194        
195     const index_t* ptr_main = Paso_SparseMatrix_borrowMainDiagonalPointer(A_p);     const index_t* ptr_main = Paso_SparseMatrix_borrowMainDiagonalPointer(A_p);
196      
197       (void)pivot;         /* silence warning from var being unused by macros */
198       (void)block_len;    
199     /* forward substitution */     /* forward substitution */
200        
201     if (n_block==1) {     if (n_block==1) {
# Line 301  void Paso_Preconditioner_LocalSmoother_S Line 305  void Paso_Preconditioner_LocalSmoother_S
305     const dim_t n=A_p->numRows;     const dim_t n=A_p->numRows;
306     const dim_t n_block=A_p->row_block_size;     const dim_t n_block=A_p->row_block_size;
307     double *diag = smoother->diag;     double *diag = smoother->diag;
308     index_t* pivot = smoother->pivot;     index_t* pivot = smoother->pivot;  
309     const dim_t block_len=A_p->block_size;     const dim_t block_len=A_p->block_size;  
310     double *y;     double *y;
311        
312     register dim_t i,k;     register dim_t i,k;
# Line 313  void Paso_Preconditioner_LocalSmoother_S Line 317  void Paso_Preconditioner_LocalSmoother_S
317     const index_t* coloring = Paso_Pattern_borrowColoringPointer(A_p->pattern);     const index_t* coloring = Paso_Pattern_borrowColoringPointer(A_p->pattern);
318     const dim_t num_colors = Paso_Pattern_getNumColors(A_p->pattern);     const dim_t num_colors = Paso_Pattern_getNumColors(A_p->pattern);
319     const index_t* ptr_main = Paso_SparseMatrix_borrowMainDiagonalPointer(A_p);     const index_t* ptr_main = Paso_SparseMatrix_borrowMainDiagonalPointer(A_p);
320      
321       (void)pivot;         /* These vars are dropped by some macros*/
322       (void)block_len;
323    
324     #pragma omp parallel  private(mm, i,iptr_ik,k,rtmp, color, y)     #pragma omp parallel  private(mm, i,iptr_ik,k,rtmp, color, y)
325     {     {

Legend:
Removed from v.3558  
changed lines
  Added in v.3559

  ViewVC Help
Powered by ViewVC 1.1.26