Cell< CellTrait > Class Template Reference

#include <Cell.h>

Public Types
enum	{ numFaces = CellTrait::numFaces }
enum	{ numNodes = CellTrait::numNodes }

Public Member Functions
	Cell ()
void	orderCellFacesAndNodes (const int c, CRConnectivity &cellFaces, CRConnectivity &cellNodes, const CRConnectivity &faceNodes, const CRConnectivity &faceCells, const Array< Vector< double, 3 > > &nodeCoordinates)

Private Attributes
unsigned int	_faceAllNodesSig [numFaces]
unsigned int	_nodeFirstFaceSig [numNodes]
unsigned int	_faceFirstFaceNodesSig [numNodes]
map< unsigned int, unsigned int >	_faceFirstFaceNodesSigMap

Detailed Description

template<class CellTrait>
class Cell< CellTrait >

Definition at line 16 of file Cell.h.

Member Enumeration Documentation

template<class CellTrait >

anonymous enum

Enumerator
numFaces

Definition at line 20 of file Cell.h.

{numFaces = CellTrait::numFaces};

template<class CellTrait >

anonymous enum

Enumerator
numNodes

Definition at line 21 of file Cell.h.

{numNodes = CellTrait::numNodes};

Constructor & Destructor Documentation

template<class CellTrait >

Cell< CellTrait >::Cell

(

)

Definition at line 52 of file Cell.cpp.

{
      for(int n=0; n<numNodes; n++) _nodeFirstFaceSig[n]  = 0;

    for(int f=0; f<numFaces; f++)
    {
        _faceAllNodesSig[f] = 0;
        for(unsigned int nf=0; nf<CellTrait::faceNodeCount[f]; nf++)
        {
            unsigned int n = CellTrait::faceNodes[f][nf];
            _faceAllNodesSig[f] |= (1<<n);
            if (f == 0)
              _nodeFirstFaceSig[n] = 1<<n;
        }
    }

    for(int f=0; f<numFaces; f++)
    {
        _faceFirstFaceNodesSig[f] = 0;
        for(unsigned int nf=0; nf<CellTrait::faceNodeCount[f]; nf++)
        {
            unsigned int n = CellTrait::faceNodes[f][nf];
            _faceFirstFaceNodesSig[f] |= _nodeFirstFaceSig[n];
        }
      _faceFirstFaceNodesSigMap[_faceFirstFaceNodesSig[f]] = f;
    }

#if 0
    for(int f=0; f<numFaces; f++)
      cout << _faceAllNodesSig[f] << " ";
    cout << endl;
    
    for(int n=0; n<numNodes; n++)
      cout << _nodeFirstFaceSig[n] <<  " " ;
    cout << endl;

    for(int f=0; f<numFaces; f++)
      cout << _faceFirstFaceNodesSig[f] <<  "   " ;
    cout << endl;
#endif
}

Member Function Documentation

template<class CellTrait >

void Cell< CellTrait >::orderCellFacesAndNodes	(	const int	c,
		CRConnectivity &	cellFaces,
		CRConnectivity &	cellNodes,
		const CRConnectivity &	faceNodes,
		const CRConnectivity &	faceCells,
		const Array< Vector< double, 3 > > &	nodeCoordinates
	)

Definition at line 97 of file Cell.cpp.

References CRConnectivity::getCount().

{

  int face0 = 0;
  int face0NodeCount=0;
  bool reverseFace0Nodes = false;

  for(int nf=0; nf<numFaces; nf++)
  {
      const int f = cellFaces(c,nf);
      if ((unsigned int)faceNodes.getCount(f) == CellTrait::faceNodeCount[0])
      {
          face0 = f;
          face0NodeCount = faceNodes.getCount(f);
          if (faceCells(f,0) != c)
          {
              if (faceCells(f,1) != c)
                cerr  << "malformed grid " << endl;
              reverseFace0Nodes = true;
          }
          break;
      }
  }

  map<int, unsigned int> thisNodesFirstFaceNodesSig;

  for(unsigned int nn=0; nn<(unsigned int)face0NodeCount; nn++)
  {
      int node = faceNodes(face0,nn);
      if (reverseFace0Nodes) node = faceNodes(face0,face0NodeCount-nn-1);
      thisNodesFirstFaceNodesSig[node] = 1<<nn;
  }

  int orderedFaces[numFaces];

  for(int nf=0; nf<numFaces; nf++)
  {
      unsigned int thisFaceFirstFaceNodesSig = 0;
      const int f = cellFaces(c,nf);

      for(int nn=0; nn<faceNodes.getCount(f); nn++)
      {
          const int n = faceNodes(f,nn);

          if (thisNodesFirstFaceNodesSig.find(n) !=
              thisNodesFirstFaceNodesSig.end())
            thisFaceFirstFaceNodesSig |= thisNodesFirstFaceNodesSig[n];
      }

      int faceOrder = _faceFirstFaceNodesSigMap[thisFaceFirstFaceNodesSig];
      orderedFaces[faceOrder] = f;
  }

  for(int nf=0; nf<numFaces; nf++)
    cellFaces(c,nf) = orderedFaces[nf];
  

  map<int, unsigned int> thisNodesAllFaceNodesSig;

  for(int nf=0; nf<numFaces; nf++)
  {
      //      unsigned int thisFaceFirstFaceNodeSig = 0;
      const int f = cellFaces(c,nf);

      //      bool reverseFaceNodes = (faceCells(f,0) != c);

      const int faceNodeCount = faceNodes.getCount(f);
      for(int nn=0; nn<faceNodeCount; nn++)
      {
          int n = faceNodes(f,nn);
          if (reverseFace0Nodes) n = faceNodes(f,faceNodeCount-nn-1);

          if (thisNodesAllFaceNodesSig.find(n) ==
              thisNodesAllFaceNodesSig.end())
          {
              thisNodesAllFaceNodesSig[n] =
                _faceAllNodesSig[nf];
          }
          else
          {
              thisNodesAllFaceNodesSig[n] &=
                _faceAllNodesSig[nf];
          }
      }
  }

  for( map<int,unsigned int>::const_iterator pos =
         thisNodesAllFaceNodesSig.begin();
       pos != thisNodesAllFaceNodesSig.end();
       ++pos)
  {
    const int node = pos->first;
    // log2 missing in msvc ?
          const int index = (int) log2(pos->second);
    // const int index = (int) (log( (double) pos->second)/log(2.));
      cellNodes(c,index) = node;
  }
   
}

Member Data Documentation

template<class CellTrait >

unsigned int Cell< CellTrait >::_faceAllNodesSig[numFaces]

private

Definition at line 35 of file Cell.h.

template<class CellTrait >

unsigned int Cell< CellTrait >::_faceFirstFaceNodesSig[numNodes]

private

Definition at line 37 of file Cell.h.

template<class CellTrait >

map<unsigned int, unsigned int> Cell< CellTrait >::_faceFirstFaceNodesSigMap

private

Definition at line 39 of file Cell.h.

template<class CellTrait >

unsigned int Cell< CellTrait >::_nodeFirstFaceSig[numNodes]

private

Definition at line 36 of file Cell.h.

---

The documentation for this class was generated from the following files:

• Cell.h
• Cell.cpp