/[escript]/branches/domexper/dudley/src/ElementFile_markNodes.c
ViewVC logotype

Diff of /branches/domexper/dudley/src/ElementFile_markNodes.c

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

revision 1769 by gross, Mon Sep 8 02:47:55 2008 UTC revision 1770 by ksteube, Mon Sep 8 06:58:47 2008 UTC
# Line 25  Line 25 
25    
26  /**************************************************************/  /**************************************************************/
27    
28  void Finley_ElementFile_markNodes(index_t* mask,index_t offset,Finley_ElementFile* in,bool_t useLinear) {  void Finley_ElementFile_markNodes(index_t* mask,index_t offset,dim_t numNodes,Finley_ElementFile* in,bool_t useLinear) {
29     dim_t i,NN,NN2,e;     dim_t i,NN,NN2,e;
30     index_t *lin_node,*id=NULL;     index_t *lin_node,*id=NULL;
31    
# Line 44  void Finley_ElementFile_markNodes(index_ Line 44  void Finley_ElementFile_markNodes(index_
44          #pragma omp parallel for private(e,i) schedule(static)          #pragma omp parallel for private(e,i) schedule(static)
45          for (e=0;e<in->numElements;e++) {          for (e=0;e<in->numElements;e++) {
46              for (i=0;i<NN;i++) {              for (i=0;i<NN;i++) {
47    #ifdef BOUNDS_CHECK
48    if ((in->Nodes[INDEX2(lin_node[i],e,NN2)]-offset) >= numNodes) { printf("BOUNDS_CHECK %s %d i=%d e=%d NN2=%d offset=%d index=%d\n", __FILE__, __LINE__, i, e, NN2, offset, in->Nodes[INDEX2(lin_node[i],e,NN2)]-offset); exit(1); }
49    #endif
50                  mask[in->Nodes[INDEX2(lin_node[i],e,NN2)]-offset]=1;                  mask[in->Nodes[INDEX2(lin_node[i],e,NN2)]-offset]=1;
51             }             }
52          }          }

Legend:
Removed from v.1769  
changed lines
  Added in v.1770

  ViewVC Help
Powered by ViewVC 1.1.26