ReduceSpectralMesh.h

Go to the documentation of this file.

#ifndef __smoab_detail_ReduceSpectralMesh_h
#define __smoab_detail_ReduceSpectralMesh_h

#include "SimpleMoab.h"
#include "vtkCellType.h"

#include <algorithm>
#include <vector>

namespace smoab{ namespace detail{


class ReduceSpectralMesh
{
public:
  ReduceSpectralMeshReduceHighOrderCell(smoab::Range const& cells, moab::Interface* moab):
    LinearCellConnectivity(),
    HighOrderCellConnectivity()
    {
    int count = 0;
    const std::size_t cellSize=cells.size();
    while(count != cellSize)
      {
      EntityHandle* connectivity;
      int numVerts=0, iterationCount=0;
      //use the highly efficent calls, since we know that are of the same dimension
      moab->connect_iterate(cells.begin()+count,
                            cells.end(),
                            connectivity,
                            numVerts,
                            iterationCount);
      //if we didn't read anything, break!
      if(iterationCount == 0)
        {
        break;
        }

      //identify the cell type that we currently have,
      //store that along with the connectivity in a temp storage vector
      const moab::EntityType type = moab->type_from_handle(*cells.begin()+count);

      //instead of storing the connectivity in a single array,
      //what my goal is to have multiple vectors. The first vector
      //will hold all the linear cells, while the second vector
      //will hold the high order cells.

      //Since connect_iterate returns the number of cells of a given
      //type that we have read we can also store the RunLegthInfo
      //for all the cells allowing us to apply transformations to the
      //high order cells efficently in bulk.

      //while all these cells are contiously of the same type,
      //quadric hexs in vtk have 20 points, but moab has 21 so we
      //need to store this difference
      }
    }

private:
  std::vector<EntityHandle> LinearCellConnectivity;
  std::vector<EntityHandle> HighOrderCellConnectivity;

  std::vector<detail::ContinousCellInfo> LinearCellInfo;
  std::vector<detail::ContinousCellInfo> HighOrderCellInfo;


} } //namespace smoab::detail

#endif