/[escript]/trunk/dudley/src/IndexList.cpp
ViewVC logotype

Diff of /trunk/dudley/src/IndexList.cpp

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

revision 1721 by phornby, Fri Jul 11 13:12:46 2008 UTC revision 1722 by gross, Fri Aug 22 04:20:30 2008 UTC
# Line 99  void Finley_IndexList_insertElementsWith Line 99  void Finley_IndexList_insertElementsWith
99      }      }
100    }    }
101  }  }
102    void Finley_IndexList_insertElementsWithRowRangeNoMainDiagonal(Finley_IndexList* index_list, index_t firstRow, index_t lastRow,
103                                                                  Finley_ElementFile* elements, index_t* row_map, index_t* col_map)
104    {
105      index_t color;
106      dim_t e,kr,kc,icol,irow, NN,irow_loc;
107      if (elements!=NULL) {
108        NN=elements->numNodes;
109        for (color=elements->minColor;color<=elements->maxColor;color++) {
110               #pragma omp for private(e,irow,kr,kc,icol,irow_loc) schedule(static)
111               for (e=0;e<elements->numElements;e++) {
112                   if (elements->Color[e]==color) {
113                       for (kr=0;kr<NN;kr++) {
114                         irow=row_map[elements->Nodes[INDEX2(kr,e,NN)]];
115                         if ((firstRow<=irow) && (irow < lastRow)) {
116                              irow_loc=irow-firstRow;
117                              for (kc=0;kc<NN;kc++) {
118                                  icol=col_map[elements->Nodes[INDEX2(kc,e,NN)]];
119                                  if (icol != irow) Finley_IndexList_insertIndex(&(index_list[irow_loc]),icol);
120                              }
121                          }
122                      }
123                   }
124               }
125        }
126      }
127    }
128    
129  /* inserts row index row into the Finley_IndexList in if it does not exist */  /* inserts row index row into the Finley_IndexList in if it does not exist */
130    

Legend:
Removed from v.1721  
changed lines
  Added in v.1722

  ViewVC Help
Powered by ViewVC 1.1.26