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

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

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

revision 3488 by gross, Thu Mar 24 22:44:40 2011 UTC revision 3489 by caltinay, Wed Mar 30 00:46:04 2011 UTC
# Line 461  void Paso_Preconditioner_AMG_setStrongCo Line 461  void Paso_Preconditioner_AMG_setStrongCo
461       register double sum_row=0;       register double sum_row=0;
462       register double main_row=0;       register double main_row=0;
463       register dim_t kdeg=0;       register dim_t kdeg=0;
464           const register index_t koffset=A->mainBlock->pattern->ptr[i]+A->col_coupleBlock->pattern->ptr[i];           register const index_t koffset=A->mainBlock->pattern->ptr[i]+A->col_coupleBlock->pattern->ptr[i];
465                    
466       /* collect information for row i: */       /* collect information for row i: */
467       #pragma ivdep       #pragma ivdep
# Line 532  void Paso_Preconditioner_AMG_setStrongCo Line 532  void Paso_Preconditioner_AMG_setStrongCo
532                    
533            const double threshold = remote_threshold[2*i+1];            const double threshold = remote_threshold[2*i+1];
534            register dim_t kdeg=0;            register dim_t kdeg=0;
535                const register index_t koffset=koffset_0+A->row_coupleBlock->pattern->ptr[i];                register const index_t koffset=koffset_0+A->row_coupleBlock->pattern->ptr[i];
536                if (remote_threshold[2*i]>0) {                if (remote_threshold[2*i]>0) {
537               #pragma ivdep               #pragma ivdep
538               for (iptr=A->row_coupleBlock->pattern->ptr[i];iptr<A->row_coupleBlock->pattern->ptr[i+1]; ++iptr) {               for (iptr=A->row_coupleBlock->pattern->ptr[i];iptr<A->row_coupleBlock->pattern->ptr[i+1]; ++iptr) {
# Line 593  void Paso_Preconditioner_AMG_setStrongCo Line 593  void Paso_Preconditioner_AMG_setStrongCo
593       register double main_row=0;       register double main_row=0;
594       register index_t rtmp_offset=-A->mainBlock->pattern->ptr[i];       register index_t rtmp_offset=-A->mainBlock->pattern->ptr[i];
595       register dim_t kdeg=0;       register dim_t kdeg=0;
596       const register index_t koffset=A->mainBlock->pattern->ptr[i]+A->col_coupleBlock->pattern->ptr[i];       register const index_t koffset=A->mainBlock->pattern->ptr[i]+A->col_coupleBlock->pattern->ptr[i];
597            
598       /* collect information for row i: */       /* collect information for row i: */
599       for (iptr=A->mainBlock->pattern->ptr[i];iptr<A->mainBlock->pattern->ptr[i+1]; ++iptr) {       for (iptr=A->mainBlock->pattern->ptr[i];iptr<A->mainBlock->pattern->ptr[i+1]; ++iptr) {
# Line 685  void Paso_Preconditioner_AMG_setStrongCo Line 685  void Paso_Preconditioner_AMG_setStrongCo
685            
686       const double threshold2 = remote_threshold[2*i+1]*remote_threshold[2*i+1];       const double threshold2 = remote_threshold[2*i+1]*remote_threshold[2*i+1];
687       register dim_t kdeg=0;       register dim_t kdeg=0;
688       const register index_t koffset=koffset_0+A->row_coupleBlock->pattern->ptr[i];       register const index_t koffset=koffset_0+A->row_coupleBlock->pattern->ptr[i];
689       if (remote_threshold[2*i]>0) {       if (remote_threshold[2*i]>0) {
690          #pragma ivdep          #pragma ivdep
691          for (iptr=A->row_coupleBlock->pattern->ptr[i];iptr<A->row_coupleBlock->pattern->ptr[i+1]; ++iptr) {          for (iptr=A->row_coupleBlock->pattern->ptr[i];iptr<A->row_coupleBlock->pattern->ptr[i+1]; ++iptr) {

Legend:
Removed from v.3488  
changed lines
  Added in v.3489

  ViewVC Help
Powered by ViewVC 1.1.26