/[escript]/trunk/finley/src/NodeFile.cpp
ViewVC logotype

Diff of /trunk/finley/src/NodeFile.cpp

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

revision 4934 by jfenwick, Tue May 13 00:28:11 2014 UTC revision 5037 by caltinay, Thu Jun 12 04:39:23 2014 UTC
# Line 1090  void NodeFile::createNodeMappings(const Line 1090  void NodeFile::createNodeMappings(const
1090    
1091      // get the Ids for DOFs and reduced nodes      // get the Ids for DOFs and reduced nodes
1092      if (noError()) {      if (noError()) {
1093            const int rnTargets = reducedNodesMapping.getNumTargets();
1094            const int dofTargets = degreesOfFreedomMapping.getNumTargets();
1095            const int rdofTargets = reducedDegreesOfFreedomMapping.getNumTargets();
1096  #pragma omp parallel  #pragma omp parallel
1097          {          {
1098  #pragma omp for  #pragma omp for
1099           for (int i=0; i<reducedNodesMapping.getNumTargets(); ++i)           for (int i=0; i<rnTargets; ++i)
1100               reducedNodesId[i]=Id[reducedNodesMapping.map[i]];               reducedNodesId[i]=Id[reducedNodesMapping.map[i]];
1101  #pragma omp for  #pragma omp for
1102           for (int i=0; i<degreesOfFreedomMapping.getNumTargets(); ++i)           for (int i=0; i<dofTargets; ++i)
1103               degreesOfFreedomId[i]=Id[degreesOfFreedomMapping.map[i]];               degreesOfFreedomId[i]=Id[degreesOfFreedomMapping.map[i]];
1104  #pragma omp for  #pragma omp for
1105           for (int i=0; i<reducedDegreesOfFreedomMapping.getNumTargets(); ++i)           for (int i=0; i<rdofTargets; ++i)
1106               reducedDegreesOfFreedomId[i]=Id[reducedDegreesOfFreedomMapping.map[i]];               reducedDegreesOfFreedomId[i]=Id[reducedDegreesOfFreedomMapping.map[i]];
1107          }          }
1108      } else {      } else {

Legend:
Removed from v.4934  
changed lines
  Added in v.5037

  ViewVC Help
Powered by ViewVC 1.1.26