/[escript]/trunk/finley/src/Mesh_hex20.c
ViewVC logotype

Diff of /trunk/finley/src/Mesh_hex20.c

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

revision 782 by bcumming, Tue Jul 18 00:47:47 2006 UTC revision 787 by bcumming, Wed Jul 26 01:46:45 2006 UTC
# Line 829  Finley_Mesh* Finley_RectangularMesh_Hex2 Line 829  Finley_Mesh* Finley_RectangularMesh_Hex2
829        return NULL;        return NULL;
830    }    }
831    
832        k=0;
833    #pragma omp parallel for private(i0,i1,i2,k)    #pragma omp parallel for private(i0,i1,i2,k)
834    for (k=0,i2=0;i2<N2;i2++) {    for (i2=0;i2<N2;i2++) {
835      for (i1=0;i1<N1;i1++) {      for (i1=0;i1<N1;i1++) {
836        for (i0=0;i0<N0t;i0++,k++) {                for (i0=0;i0<N0t;i0++,k++) {        
837          out->Nodes->Coordinates[INDEX2(0,k,3)]=DBLE((i0+firstNodeConstruct) % N0)/DBLE(N0-1)*Length[0];          out->Nodes->Coordinates[INDEX2(0,k,3)]=DBLE((i0+firstNodeConstruct) % N0)/DBLE(N0-1)*Length[0];
# Line 960  Finley_Mesh* Finley_RectangularMesh_Hex2 Line 961  Finley_Mesh* Finley_RectangularMesh_Hex2
961        for (i0=0;i0<numElementsLocal;i0++,k++) {        for (i0=0;i0<numElementsLocal;i0++,k++) {
962                  node0 = (periodicLocal[0] && !i0) ? 2*(i1*N0t + i2*N1*N0t) : 2*(i1*N0t + i2*N1*N0t + i0) + periodicLocal[0];                  node0 = (periodicLocal[0] && !i0) ? 2*(i1*N0t + i2*N1*N0t) : 2*(i1*N0t + i2*N1*N0t + i0) + periodicLocal[0];
963    
964          out->Elements->Id[k]=k;                  out->Elements->Id[k]=((firstNodeConstruct/2+i0)%NE0)*NE1*NE2 + NE1*i2 + i1;
965          out->Elements->Tag[k]=0;          out->Elements->Tag[k]=0;
966          out->Elements->Color[k]=COLOR_MOD(i0)+3*COLOR_MOD(i1)+9*COLOR_MOD(i2);          out->Elements->Color[k]=COLOR_MOD(i0)+3*COLOR_MOD(i1)+9*COLOR_MOD(i2);
967          out->Elements->Dom[k]=ELEMENT_INTERNAL;          out->Elements->Dom[k]=ELEMENT_INTERNAL;

Legend:
Removed from v.782  
changed lines
  Added in v.787

  ViewVC Help
Powered by ViewVC 1.1.26