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

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

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

revision 2806 by artak, Thu Dec 3 01:51:55 2009 UTC revision 2807 by artak, Mon Dec 7 00:02:55 2009 UTC
# Line 655  void Paso_Pattern_RS_MI(Paso_SparseMatri Line 655  void Paso_Pattern_RS_MI(Paso_SparseMatri
655    Paso_IndexList* index_list=NULL;    Paso_IndexList* index_list=NULL;
656        
657    index_t* counter;    index_t* counter;
658      /*dim_t lk;*/
659    
660    index_list=TMPMEMALLOC(A->pattern->numOutput,Paso_IndexList);    index_list=TMPMEMALLOC(A->pattern->numOutput,Paso_IndexList);
661     if (! Paso_checkPtr(index_list)) {     if (! Paso_checkPtr(index_list)) {
# Line 711  void Paso_Pattern_RS_MI(Paso_SparseMatri Line 712  void Paso_Pattern_RS_MI(Paso_SparseMatri
712    /*S_i={j \in N_i; i strongly coupled to j}*/    /*S_i={j \in N_i; i strongly coupled to j}*/
713    
714    /*    /*
715     #pragma omp parallel for private(i,k) schedule(static)    #pragma omp parallel for private(i,iptr,lk) schedule(static)
716    for (i=0;i<n;++i) {    for (i=0;i<n;++i) {
717          k=0;          lk=0;
718          for (iptr=A->pattern->ptr[i];iptr<A->pattern->ptr[i+1]; ++iptr) {          for (iptr=A->pattern->ptr[i];iptr<A->pattern->ptr[i+1]; ++iptr) {
719              if(ABS(A->val[iptr])>1.e-15 && A->pattern->index[iptr]!=i )              if(ABS(A->val[iptr])>1.e-15 && A->pattern->index[iptr]!=i )
720                 k++;                 lk++;
721          }          }
722            #pragma omp critical
723            k+=lk;
724          if(k==0) {          if(k==0) {
725              mis_marker[i]=IS_IN_F;              mis_marker[i]=IS_IN_F;
726          }          }

Legend:
Removed from v.2806  
changed lines
  Added in v.2807

  ViewVC Help
Powered by ViewVC 1.1.26