LoadGeometry.h

Go to the documentation of this file.

#ifndef __smoab_detail_LoadGeometry_h
#define __smoab_detail_LoadGeometry_h

#include "SimpleMoab.h"
#include "MixedCellConnectivity.h"

#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
namespace smoab
{
namespace detail
{

    namespace
    {

        //empty definition, only allow specializations to compile
        template < typename VTKGridType >
        struct setCells;

        template <>
        struct setCells< vtkUnstructuredGrid >
        {
            typedef vtkUnstructuredGrid VTKGridType;
            void operator()( VTKGridType* grid,
                             int,
                             vtkUnsignedCharArray* types,
                             vtkIdTypeArray* locations,
                             vtkCellArray* cells ) const
            {
                grid->SetCells( types, locations, cells, NULL, NULL );
            }
        };

        template <>
        struct setCells< vtkPolyData >
        {
            typedef vtkPolyData VTKGridType;
            void operator()( VTKGridType* grid,
                             int dim,
                             vtkUnsignedCharArray*,
                             vtkIdTypeArray*,
                             vtkCellArray* cells ) const
            {
                if( dim == 0 )
                {
                    grid->SetVerts( cells );
                }
                else if( dim == 1 )
                {
                    grid->SetLines( cells );
                }
                else if( dim == 2 )
                {
                    grid->SetPolys( cells );
                }
            }
        };

    }  // namespace
    //basic class that given a range of moab cells will convert them to am unstructured grid.
    //holds only references to input cells, so they can't go out of scope

    //The Topology tag is used to describe the dimensoniality of the cells we are reading
    // 0 == verts, 1 == lines, 2 == tri/quads, 3 = volume elements
    class LoadGeometry
    {
        moab::Interface* Interface;
        int TopologyDim;
        smoab::detail::MixedCellConnectivity MixConn;
        smoab::Range Points;

      public:
        //----------------------------------------------------------------------------
        //warning cells input to constructor is held by reference
        LoadGeometry( const smoab::Range& cells, int topologyDim, const smoab::Interface& interface )
            : Interface( interface.Moab ), TopologyDim( topologyDim ), MixConn( cells, interface.Moab ), Points()
        {
        }

        const smoab::Range& moabPoints() const
        {
            return this->Points;
        }

        //----------------------------------------------------------------------------
        //todo: have support for using only a subsection of the input cells
        template < typename vtkDataSetType >
        void fill( vtkDataSetType* dataSet )
        {
            //now that mixConn has all the cells properly stored, lets fixup
            //the ids so that they start at zero and keep the same logical ordering
            //as before.
            vtkIdType numCells, connLen;
            this->MixConn.compactIds( numCells, connLen );
            this->addGridsTopology( dataSet, numCells, connLen );
            this->addCoordinates( dataSet );
        }

      private:
        //----------------------------------------------------------------------------
        void addCoordinates( vtkPointSet* grid )
        {
            //this is sorta of hackish as moabPoints is only valid
            //after compactIds has been called
            this->MixConn.moabPoints( this->Points );

            //since the smoab::range are always unique and sorted
            //we can use the more efficient coords_iterate
            //call in moab, which returns moab internal allocated memory
            vtkNew< vtkPoints > newPoints;
            newPoints->SetDataTypeToDouble();
            newPoints->SetNumberOfPoints( this->Points.size() );

            //need a pointer to the allocated vtkPoints memory so that we
            //don't need to use an extra copy and we can bypass all vtk's check
            //on out of bounds
            double* rawPoints = static_cast< double* >( newPoints->GetVoidPointer( 0 ) );
            this->Interface->get_coords( this->Points, rawPoints );

            grid->SetPoints( newPoints.GetPointer() );
        }

        //----------------------------------------------------------------------------
        template < typename vtkDataSetType >
        void addGridsTopology( vtkDataSetType* grid, vtkIdType numCells, vtkIdType numConnectivity ) const
        {
            //correct the connectivity size to account for the vtk padding
            const vtkIdType vtkConnectivity = numCells + numConnectivity;

            vtkNew< vtkIdTypeArray > cellArray;
            vtkNew< vtkIdTypeArray > cellLocations;
            vtkNew< vtkUnsignedCharArray > cellTypes;

            cellArray->SetNumberOfValues( vtkConnectivity );
            cellLocations->SetNumberOfValues( numCells );
            cellTypes->SetNumberOfValues( numCells );

            vtkIdType* rawArray     = static_cast< vtkIdType* >( cellArray->GetVoidPointer( 0 ) );
            vtkIdType* rawLocations = static_cast< vtkIdType* >( cellLocations->GetVoidPointer( 0 ) );
            unsigned char* rawTypes = static_cast< unsigned char* >( cellTypes->GetVoidPointer( 0 ) );

            this->MixConn.copyToVtkCellInfo( rawArray, rawLocations, rawTypes );

            vtkNew< vtkCellArray > cells;
            cells->SetCells( numCells, cellArray.GetPointer() );

            setCells< vtkDataSetType >()( grid, this->TopologyDim, cellTypes.GetPointer(), cellLocations.GetPointer(),
                                          cells.GetPointer() );
        }
    };

}  // namespace detail
}  // namespace smoab

#endif  // __smoab_detail_LoadGeometry_h