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

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

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

revision 1388 by trankine, Fri Jan 11 07:45:58 2008 UTC revision 1733 by gross, Thu Aug 28 01:46:07 2008 UTC
# Line 114  Finley_Mesh* Finley_RectangularMesh_Rec4 Line 114  Finley_Mesh* Finley_RectangularMesh_Rec4
114    }    }
115    offset0=e_offset0*N_PER_E;    offset0=e_offset0*N_PER_E;
116    offset1=e_offset1*N_PER_E;    offset1=e_offset1*N_PER_E;
117    local_N0=local_NE0*N_PER_E+1;    local_N0=local_NE0>0 ? local_NE0*N_PER_E+1 : 0;
118    local_N1=local_NE1*N_PER_E+1;    local_N1=local_NE1>0 ? local_NE1*N_PER_E+1 : 0;
119    
120    /* get the number of surface elements */    /* get the number of surface elements */
121    
122    NFaceElements=0;    NFaceElements=0;
123    if (!periodic[0]) {    if (!periodic[0] && (local_NE0>0)) {
124       NDOF0=N0;       NDOF0=N0;
125       if (e_offset0 == 0) NFaceElements+=local_NE1;       if (e_offset0 == 0) NFaceElements+=local_NE1;
126       if (local_NE0+e_offset0 == NE0) NFaceElements+=local_NE1;       if (local_NE0+e_offset0 == NE0) NFaceElements+=local_NE1;
127    } else {    } else {
128        NDOF0=N0-1;        NDOF0=N0-1;
129    }    }
130    if (!periodic[1]) {    if (!periodic[1]  && (local_NE1>0)) {
131       NDOF1=N1;       NDOF1=N1;
132       if (e_offset1 == 0) NFaceElements+=local_NE0;       if (e_offset1 == 0) NFaceElements+=local_NE0;
133       if (local_NE1+e_offset1 == NE1) NFaceElements+=local_NE0;       if (local_NE1+e_offset1 == NE1) NFaceElements+=local_NE0;
# Line 180  Finley_Mesh* Finley_RectangularMesh_Rec4 Line 180  Finley_Mesh* Finley_RectangularMesh_Rec4
180       NN=out->FaceElements->numNodes;       NN=out->FaceElements->numNodes;
181       totalNECount=NE0*NE1;       totalNECount=NE0*NE1;
182       faceNECount=0;       faceNECount=0;
183       if (!periodic[0]) {       if (!periodic[0]  && (local_NE0>0)) {
184          /* **  elements on boundary 001 (x1=0): */          /* **  elements on boundary 001 (x1=0): */
185            
186          if (e_offset0 == 0) {          if (e_offset0 == 0) {
# Line 231  Finley_Mesh* Finley_RectangularMesh_Rec4 Line 231  Finley_Mesh* Finley_RectangularMesh_Rec4
231           }           }
232           totalNECount+=NE1;           totalNECount+=NE1;
233       }       }
234       if (!periodic[1]) {       if (!periodic[1]  && (local_NE1>0)) {
235          /* **  elements on boundary 010 (x2=0): */          /* **  elements on boundary 010 (x2=0): */
236          if (e_offset1 == 0) {          if (e_offset1 == 0) {
237             #pragma omp parallel for private(i0,k,node0)             #pragma omp parallel for private(i0,k,node0)

Legend:
Removed from v.1388  
changed lines
  Added in v.1733

  ViewVC Help
Powered by ViewVC 1.1.26