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

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

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

revision 3258 by jfenwick, Thu Jan 28 02:03:15 2010 UTC revision 3259 by jfenwick, Mon Oct 11 01:48:14 2010 UTC
# Line 139  void Paso_Pattern_reduceBandwidth(Paso_P Line 139  void Paso_Pattern_reduceBandwidth(Paso_P
139     index_t root, *AssignedLevel=NULL, *VerticesInTree=NULL, *firstVertexInLevel=NULL,k, *oldLabel=NULL;     index_t root, *AssignedLevel=NULL, *VerticesInTree=NULL, *firstVertexInLevel=NULL,k, *oldLabel=NULL;
140     /* check input */     /* check input */
141     if (N != pattern->numInput) {     if (N != pattern->numInput) {
142        Paso_setError(VALUE_ERROR,"Paso_Pattern_reduceBandwidth: pattern needs to be for a square matrix.");        Esys_setError(VALUE_ERROR,"Paso_Pattern_reduceBandwidth: pattern needs to be for a square matrix.");
143     } else if (N > 0) {     } else if (N > 0) {
144  /* printf("relabeling of %d DOFs started.\n",N); */  /* printf("relabeling of %d DOFs started.\n",N); */
145        degAndIdx=TMPMEMALLOC(N,Paso_DegreeAndIdx);        degAndIdx=TMPMEMALLOC(N,Paso_DegreeAndIdx);
# Line 147  void Paso_Pattern_reduceBandwidth(Paso_P Line 147  void Paso_Pattern_reduceBandwidth(Paso_P
147        AssignedLevel=TMPMEMALLOC(N,index_t);        AssignedLevel=TMPMEMALLOC(N,index_t);
148        VerticesInTree=TMPMEMALLOC(N,index_t);        VerticesInTree=TMPMEMALLOC(N,index_t);
149        firstVertexInLevel=TMPMEMALLOC(N+1,index_t);        firstVertexInLevel=TMPMEMALLOC(N+1,index_t);
150        if (! ( Paso_checkPtr(degAndIdx) || Paso_checkPtr(oldLabel) || Paso_checkPtr(AssignedLevel) || Paso_checkPtr(VerticesInTree) || Paso_checkPtr(firstVertexInLevel) ) ) {        if (! ( Esys_checkPtr(degAndIdx) || Esys_checkPtr(oldLabel) || Esys_checkPtr(AssignedLevel) || Esys_checkPtr(VerticesInTree) || Esys_checkPtr(firstVertexInLevel) ) ) {
151           /* get the initial bandwidth */           /* get the initial bandwidth */
152           #pragma omp parallel for private(i)           #pragma omp parallel for private(i)
153           for (i=0;i<N;++i) oldToNew[i]=i;           for (i=0;i<N;++i) oldToNew[i]=i;

Legend:
Removed from v.3258  
changed lines
  Added in v.3259

  ViewVC Help
Powered by ViewVC 1.1.26