/[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 2726 by artak, Wed Oct 21 23:50:05 2009 UTC revision 2728 by artak, Thu Oct 22 00:13:10 2009 UTC
# Line 391  void Paso_Pattern_greedy_color(Paso_Patt Line 391  void Paso_Pattern_greedy_color(Paso_Patt
391  /*For testing */  /*For testing */
392  void Paso_Pattern_greedy_diag(Paso_SparseMatrix* A, index_t* mis_marker, double threshold) {  void Paso_Pattern_greedy_diag(Paso_SparseMatrix* A, index_t* mis_marker, double threshold) {
393    
394    dim_t i,j=0,k;    dim_t i,j,k;
395    double *theta;    double *theta;
396    index_t iptr;    index_t iptr;
397    dim_t n=A->numRows;    dim_t n=A->numRows;
# Line 408  void Paso_Pattern_greedy_diag(Paso_Spars Line 408  void Paso_Pattern_greedy_diag(Paso_Spars
408      return;      return;
409    }    }
410        
411      j=0;
412    
     #pragma omp parallel for private(i,iptr,j,rsum) schedule(static)  
413      for (i=0;i<n;++i) {      for (i=0;i<n;++i) {
414          rsum=0;      rsum=0;
415          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) {
416              j=A->pattern->index[iptr];              j=A->pattern->index[iptr];
417              if(j!=i) {              if(j!=i) {
# Line 429  void Paso_Pattern_greedy_diag(Paso_Spars Line 429  void Paso_Pattern_greedy_diag(Paso_Spars
429            
430      while (Paso_Util_isAny(n,mis_marker,IS_AVAILABLE)) {      while (Paso_Util_isAny(n,mis_marker,IS_AVAILABLE)) {
431           k=0;           k=0;
          #pragma omp parallel for private(i,j,k) schedule(static)  
432           for (i=0;i<n;++i) {           for (i=0;i<n;++i) {
433             if(mis_marker[i]==IS_AVAILABLE) {             if(mis_marker[i]==IS_AVAILABLE) {
434                  if(k==0) {                  if(k==0) {
# Line 454  void Paso_Pattern_greedy_diag(Paso_Spars Line 453  void Paso_Pattern_greedy_diag(Paso_Spars
453                            
454           }           }
455                            
          #pragma omp parallel for private(i,iptr,j,rsum) schedule(static)  
456          for (i=0;i<n;++i) {          for (i=0;i<n;++i) {
457              if(AvADJ[i]) {              if(AvADJ[i]) {
458                  rsum=0;                  rsum=0;

Legend:
Removed from v.2726  
changed lines
  Added in v.2728

  ViewVC Help
Powered by ViewVC 1.1.26