/[escript]/branches/diaplayground/ripley/src/Brick.cpp
ViewVC logotype

Diff of /branches/diaplayground/ripley/src/Brick.cpp

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

revision 3764 by caltinay, Tue Jan 10 06:31:15 2012 UTC revision 3766 by caltinay, Thu Jan 12 08:17:49 2012 UTC
# Line 684  IndexVector Brick::getNumFacesPerBoundar Line 684  IndexVector Brick::getNumFacesPerBoundar
684      return ret;      return ret;
685  }  }
686    
687    IndexVector Brick::getNumSubdivisionsPerDim() const
688    {
689        IndexVector ret;
690        ret.push_back(m_NX);
691        ret.push_back(m_NY);
692        ret.push_back(m_NZ);
693        return ret;
694    }
695    
696  pair<double,double> Brick::getFirstCoordAndSpacing(dim_t dim) const  pair<double,double> Brick::getFirstCoordAndSpacing(dim_t dim) const
697  {  {
698      if (dim==0)      if (dim==0)
# Line 1610  void Brick::createPattern() Line 1619  void Brick::createPattern()
1619      // The rest is assigned in the loop further down      // The rest is assigned in the loop further down
1620      m_dofMap.assign(getNumNodes(), 0);      m_dofMap.assign(getNumNodes(), 0);
1621  #pragma omp parallel for  #pragma omp parallel for
1622      for (index_t i=front; i<m_N2; i++) {      for (index_t i=front; i<front+nDOF2; i++) {
1623          for (index_t j=bottom; j<m_N1; j++) {          for (index_t j=bottom; j<bottom+nDOF1; j++) {
1624              for (index_t k=left; k<m_N0; k++) {              for (index_t k=left; k<left+nDOF0; k++) {
1625                  m_dofMap[i*m_N0*m_N1+j*m_N0+k]=(i-front)*nDOF0*nDOF1+(j-bottom)*nDOF0+k-left;                  m_dofMap[i*m_N0*m_N1+j*m_N0+k]=(i-front)*nDOF0*nDOF1+(j-bottom)*nDOF0+k-left;
1626              }              }
1627          }          }

Legend:
Removed from v.3764  
changed lines
  Added in v.3766

  ViewVC Help
Powered by ViewVC 1.1.26