/[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 3080 by jfenwick, Tue Aug 3 04:28:03 2010 UTC revision 3086 by jfenwick, Thu Aug 5 05:07:58 2010 UTC
# Line 14  Line 14 
14    
15  /**************************************************************/  /**************************************************************/
16    
17  /* Finley: Converting an element list into a matrix shape     */  /* Dudley: Converting an element list into a matrix shape     */
18    
19  /**************************************************************/  /**************************************************************/
20    
# Line 27  Line 27 
27     into the row index col. If symmetric is set, only the upper     into the row index col. If symmetric is set, only the upper
28     triangle of the matrix is stored. */     triangle of the matrix is stored. */
29    
30  void Finley_IndexList_insertElements(Finley_IndexList* index_list, Finley_ElementFile* elements,  void Dudley_IndexList_insertElements(Dudley_IndexList* index_list, Dudley_ElementFile* elements,
31                                         bool_t reduce_row_order, index_t* row_map,                                         bool_t reduce_row_order, index_t* row_map,
32                                         bool_t reduce_col_order, index_t* col_map) {                                         bool_t reduce_col_order, index_t* col_map) {
33    /* index_list is an array of linked lists. Each entry is a row (DOF) and contains the indices to the non-zero columns */    /* index_list is an array of linked lists. Each entry is a row (DOF) and contains the indices to the non-zero columns */
34    index_t color;    index_t color;
35    Finley_ReferenceElement*refElement;    Dudley_ReferenceElement*refElement;
36    dim_t e, kr, kc, NN_row, NN_col, icol, irow, NN, *row_node=NULL, *col_node=NULL, isub, numSub;    dim_t e, kr, kc, NN_row, NN_col, icol, irow, NN, *row_node=NULL, *col_node=NULL, isub, numSub;
37    if (elements!=NULL) {    if (elements!=NULL) {
38      NN=elements->numNodes;      NN=elements->numNodes;
39      refElement= Finley_ReferenceElementSet_borrowReferenceElement(elements->referenceElementSet, FALSE);      refElement= Dudley_ReferenceElementSet_borrowReferenceElement(elements->referenceElementSet, FALSE);
40                
41      if (reduce_col_order) {      if (reduce_col_order) {
42            numSub=1;            numSub=1;
# Line 67  void Finley_IndexList_insertElements(Fin Line 67  void Finley_IndexList_insertElements(Fin
67                             irow=row_map[elements->Nodes[INDEX2(row_node[INDEX2(kr,isub,NN_row)],e,NN)]];                             irow=row_map[elements->Nodes[INDEX2(row_node[INDEX2(kr,isub,NN_row)],e,NN)]];
68                             for (kc=0;kc<NN_col;kc++) {                             for (kc=0;kc<NN_col;kc++) {
69                                 icol=col_map[elements->Nodes[INDEX2(col_node[INDEX2(kc,isub,NN_col)],e,NN)]];                                 icol=col_map[elements->Nodes[INDEX2(col_node[INDEX2(kc,isub,NN_col)],e,NN)]];
70                                 Finley_IndexList_insertIndex(&(index_list[irow]),icol);                                 Dudley_IndexList_insertIndex(&(index_list[irow]),icol);
71                             }                             }
72                         }                         }
73                     }                     }
# Line 79  void Finley_IndexList_insertElements(Fin Line 79  void Finley_IndexList_insertElements(Fin
79  }  }
80    
81    
82  void Finley_IndexList_insertElementsWithRowRange(Finley_IndexList* index_list, index_t firstRow, index_t lastRow,  void Dudley_IndexList_insertElementsWithRowRange(Dudley_IndexList* index_list, index_t firstRow, index_t lastRow,
83                                                   Finley_ElementFile* elements, index_t* row_map, index_t* col_map)                                                   Dudley_ElementFile* elements, index_t* row_map, index_t* col_map)
84  {  {
85  /* this does not resolve macro elements */  /* this does not resolve macro elements */
86      index_t color;      index_t color;
# Line 97  void Finley_IndexList_insertElementsWith Line 97  void Finley_IndexList_insertElementsWith
97                            irow-=firstRow;                            irow-=firstRow;
98                            for (kc=0;kc<NN;kc++) {                            for (kc=0;kc<NN;kc++) {
99                                icol=col_map[elements->Nodes[INDEX2(kc,e,NN)]];                                icol=col_map[elements->Nodes[INDEX2(kc,e,NN)]];
100                                Finley_IndexList_insertIndex(&(index_list[irow]),icol);                                Dudley_IndexList_insertIndex(&(index_list[irow]),icol);
101                            }                            }
102                        }                        }
103                    }                    }
# Line 106  void Finley_IndexList_insertElementsWith Line 106  void Finley_IndexList_insertElementsWith
106      }      }
107    }    }
108  }  }
109  void Finley_IndexList_insertElementsWithRowRangeNoMainDiagonal(Finley_IndexList* index_list, index_t firstRow, index_t lastRow,  void Dudley_IndexList_insertElementsWithRowRangeNoMainDiagonal(Dudley_IndexList* index_list, index_t firstRow, index_t lastRow,
110                                                                Finley_ElementFile* elements, index_t* row_map, index_t* col_map)                                                                Dudley_ElementFile* elements, index_t* row_map, index_t* col_map)
111  {  {
112    /* this does not resolve macro elements */    /* this does not resolve macro elements */
113    index_t color;    index_t color;
# Line 124  void Finley_IndexList_insertElementsWith Line 124  void Finley_IndexList_insertElementsWith
124                            irow_loc=irow-firstRow;                            irow_loc=irow-firstRow;
125                            for (kc=0;kc<NN;kc++) {                            for (kc=0;kc<NN;kc++) {
126                                icol=col_map[elements->Nodes[INDEX2(kc,e,NN)]];                                icol=col_map[elements->Nodes[INDEX2(kc,e,NN)]];
127                                if (icol != irow) Finley_IndexList_insertIndex(&(index_list[irow_loc]),icol);                                if (icol != irow) Dudley_IndexList_insertIndex(&(index_list[irow_loc]),icol);
128                            }                            }
129                        }                        }
130                    }                    }
# Line 134  void Finley_IndexList_insertElementsWith Line 134  void Finley_IndexList_insertElementsWith
134    }    }
135  }  }
136    
137  /* inserts row index row into the Finley_IndexList in if it does not exist */  /* inserts row index row into the Dudley_IndexList in if it does not exist */
138    
139  void Finley_IndexList_insertIndex(Finley_IndexList* in, index_t index) {  void Dudley_IndexList_insertIndex(Dudley_IndexList* in, index_t index) {
140    dim_t i;    dim_t i;
141    /* is index in in? */    /* is index in in? */
142    for (i=0;i<in->n;i++) {    for (i=0;i<in->n;i++) {
# Line 146  void Finley_IndexList_insertIndex(Finley Line 146  void Finley_IndexList_insertIndex(Finley
146    if (in->n==INDEXLIST_LENGTH) {    if (in->n==INDEXLIST_LENGTH) {
147       /* if in->index is full check the extension */       /* if in->index is full check the extension */
148       if (in->extension==NULL) {       if (in->extension==NULL) {
149          in->extension=TMPMEMALLOC(1,Finley_IndexList);          in->extension=TMPMEMALLOC(1,Dudley_IndexList);
150          if (Finley_checkPtr(in->extension)) return;          if (Dudley_checkPtr(in->extension)) return;
151          in->extension->n=0;          in->extension->n=0;
152          in->extension->extension=NULL;          in->extension->extension=NULL;
153       }       }
154       Finley_IndexList_insertIndex(in->extension,index);       Dudley_IndexList_insertIndex(in->extension,index);
155    } else {    } else {
156       /* insert index into in->index*/       /* insert index into in->index*/
157       in->index[in->n]=index;       in->index[in->n]=index;
# Line 159  void Finley_IndexList_insertIndex(Finley Line 159  void Finley_IndexList_insertIndex(Finley
159    }    }
160  }  }
161    
162  /* counts the number of row indices in the Finley_IndexList in */  /* counts the number of row indices in the Dudley_IndexList in */
163    
164  dim_t Finley_IndexList_count(Finley_IndexList* in, index_t range_min,index_t range_max) {  dim_t Dudley_IndexList_count(Dudley_IndexList* in, index_t range_min,index_t range_max) {
165    dim_t i;    dim_t i;
166    dim_t out=0;    dim_t out=0;
167    register index_t itmp;    register index_t itmp;
# Line 172  dim_t Finley_IndexList_count(Finley_Inde Line 172  dim_t Finley_IndexList_count(Finley_Inde
172            itmp=in->index[i];            itmp=in->index[i];
173            if ((itmp>=range_min) && (range_max>itmp)) ++out;            if ((itmp>=range_min) && (range_max>itmp)) ++out;
174      }      }
175       return out+Finley_IndexList_count(in->extension, range_min,range_max);       return out+Dudley_IndexList_count(in->extension, range_min,range_max);
176    }    }
177  }  }
178    
179  /* count the number of row indices in the Finley_IndexList in */  /* count the number of row indices in the Dudley_IndexList in */
180    
181  void Finley_IndexList_toArray(Finley_IndexList* in, index_t* array, index_t range_min,index_t range_max, index_t index_offset) {  void Dudley_IndexList_toArray(Dudley_IndexList* in, index_t* array, index_t range_min,index_t range_max, index_t index_offset) {
182    dim_t i, ptr;    dim_t i, ptr;
183    register index_t itmp;    register index_t itmp;
184    if (in!=NULL) {    if (in!=NULL) {
# Line 191  void Finley_IndexList_toArray(Finley_Ind Line 191  void Finley_IndexList_toArray(Finley_Ind
191            }            }
192    
193      }      }
194      Finley_IndexList_toArray(in->extension,&(array[ptr]), range_min, range_max, index_offset);      Dudley_IndexList_toArray(in->extension,&(array[ptr]), range_min, range_max, index_offset);
195    }    }
196  }  }
197    
198  /* deallocates the Finley_IndexList in by recursive calls */  /* deallocates the Dudley_IndexList in by recursive calls */
199    
200  void Finley_IndexList_free(Finley_IndexList* in) {  void Dudley_IndexList_free(Dudley_IndexList* in) {
201    if (in!=NULL) {    if (in!=NULL) {
202      Finley_IndexList_free(in->extension);      Dudley_IndexList_free(in->extension);
203      TMPMEMFREE(in);      TMPMEMFREE(in);
204    }    }
205  }  }
206    
207  /* creates a Paso_pattern from a range of indices */  /* creates a Paso_pattern from a range of indices */
208  Paso_Pattern* Finley_IndexList_createPattern(dim_t n0, dim_t n,Finley_IndexList* index_list,index_t range_min,index_t range_max,index_t index_offset)  Paso_Pattern* Dudley_IndexList_createPattern(dim_t n0, dim_t n,Dudley_IndexList* index_list,index_t range_min,index_t range_max,index_t index_offset)
209  {  {
210     dim_t *ptr=NULL;     dim_t *ptr=NULL;
211     register dim_t s,i,itmp;     register dim_t s,i,itmp;
# Line 213  Paso_Pattern* Finley_IndexList_createPat Line 213  Paso_Pattern* Finley_IndexList_createPat
213     Paso_Pattern* out=NULL;     Paso_Pattern* out=NULL;
214    
215     ptr=MEMALLOC(n+1-n0,index_t);     ptr=MEMALLOC(n+1-n0,index_t);
216     if (! Finley_checkPtr(ptr) ) {     if (! Dudley_checkPtr(ptr) ) {
217         /* get the number of connections per row */         /* get the number of connections per row */
218         #pragma omp parallel for schedule(static) private(i)         #pragma omp parallel for schedule(static) private(i)
219         for(i=n0;i<n;++i) {         for(i=n0;i<n;++i) {
220                ptr[i-n0]=Finley_IndexList_count(&index_list[i],range_min,range_max);                ptr[i-n0]=Dudley_IndexList_count(&index_list[i],range_min,range_max);
221         }         }
222         /* accumulate ptr */         /* accumulate ptr */
223         s=0;         s=0;
# Line 229  Paso_Pattern* Finley_IndexList_createPat Line 229  Paso_Pattern* Finley_IndexList_createPat
229         ptr[n-n0]=s;         ptr[n-n0]=s;
230         /* fill index */         /* fill index */
231         index=MEMALLOC(ptr[n-n0],index_t);         index=MEMALLOC(ptr[n-n0],index_t);
232         if (! Finley_checkPtr(index)) {         if (! Dudley_checkPtr(index)) {
233                #pragma omp parallel for schedule(static)                #pragma omp parallel for schedule(static)
234                for(i=n0;i<n;++i) {                for(i=n0;i<n;++i) {
235                    Finley_IndexList_toArray(&index_list[i],&index[ptr[i-n0]],range_min,range_max,index_offset);                    Dudley_IndexList_toArray(&index_list[i],&index[ptr[i-n0]],range_min,range_max,index_offset);
236                }                }
237                out=Paso_Pattern_alloc(PATTERN_FORMAT_DEFAULT,n-n0,range_max+index_offset,ptr,index);                out=Paso_Pattern_alloc(PATTERN_FORMAT_DEFAULT,n-n0,range_max+index_offset,ptr,index);
238         }         }
239    }    }
240    if (! Finley_noError()) {    if (! Dudley_noError()) {
241          MEMFREE(ptr);          MEMFREE(ptr);
242          MEMFREE(index);          MEMFREE(index);
243          Paso_Pattern_free(out);          Paso_Pattern_free(out);

Legend:
Removed from v.3080  
changed lines
  Added in v.3086

  ViewVC Help
Powered by ViewVC 1.1.26