/[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 3448 by gross, Thu Jan 20 06:47:07 2011 UTC revision 3449 by gross, Fri Jan 21 05:48:15 2011 UTC
# Line 159  Paso_Preconditioner_AMG* Paso_Preconditi Line 159  Paso_Preconditioner_AMG* Paso_Preconditi
159       } else {       } else {
160             Paso_Preconditioner_AMG_setStrongConnections(A_p, degree_S, offset_S, S, theta,tau);             Paso_Preconditioner_AMG_setStrongConnections(A_p, degree_S, offset_S, S, theta,tau);
161       }       }
162    {
163       dim_t p;
164       for (i=0; i< my_n; ++i) {
165             printf("%d: ",i);
166            for (p=0; p<degree_S[i];++p) printf("%d ",S[offset_S[i]+p]);
167            printf("\n");
168       }
169    }
170          
171    
172  /*MPI:  /*MPI:
173       Paso_Preconditioner_AMG_RungeStuebenSearch(n, A_p->pattern->ptr, degree_S, S, F_marker, options->usePanel);       Paso_Preconditioner_AMG_RungeStuebenSearch(n, A_p->pattern->ptr, degree_S, S, F_marker, options->usePanel);
# Line 786  void Paso_Preconditioner_AMG_CIJPCoarsen Line 795  void Paso_Preconditioner_AMG_CIJPCoarsen
795    const dim_t my_n;    const dim_t my_n;
796    const dim_t overlap_n;    const dim_t overlap_n;
797    const dim_t n= my_n + overlap_n;    const dim_t n= my_n + overlap_n;
798    
799    
800          /* initialize  split_marker and split_marker :*/
801          /* those unknows which are not influenced go into F, the rest is available for F or C */
802          #pragma omp parallel for private(i) schedule(static)
803          for (i=0;i<n;++i) {
804         degree_ST[i]=0;
805         if (degree_S[i]>0) {
806            lambda[i]=0;
807            split_marker[i]=PASO_AMG_UNDECIDED;
808         } else {
809            split_marker[i]=PASO_AMG_IN_F;
810            lambda[i]=-1;
811         }
812          }
813    
814     /* set local lambda + overlap */     /* set local lambda + overlap */
815     #pragma omp parallel for private(i)     #pragma omp parallel for private(i)
816     for (i=0; i<n ++i) {     for (i=0; i<n ++i) {

Legend:
Removed from v.3448  
changed lines
  Added in v.3449

  ViewVC Help
Powered by ViewVC 1.1.26