TAPsCGAlgo.cpp

Go to the documentation of this file.

/******************************************************************************
TAPsCGAlgo.hpp
******************************************************************************/
/******************************************************************************
SUKITTI PUNAK   (2/03/2006)
UPDATE          (2/15/2010)
******************************************************************************/
#include "TAPsCGAlgo.hpp"
// Using Inclusion Model (i.e. definitions are included in declarations)
//                       (this name.cpp is included in name.hpp)
// Each friend is defined directly inside its declaration.

BEGIN_NAMESPACE_TAPs
//=============================================================================
// Subdivision Fns
//=============================================================================
//-----------------------------------------------------------------------------
template <typename T>
bool CGAlgo<T>::SubdivisionLoop ( OpenGL::OpenGLModel<T> * prModel )
{
    //-------------------------------------------------------------------
    // Loop Subdivision for HalfEdge Model
    {
        OpenGL::HalfEdgeModel<T> * const  model = 
            dynamic_cast< OpenGL::HalfEdgeModel<T> * const >( prModel );
        if ( model ) {
            std::cout << "Do Loop Subdivision to A HalfEdge Model\n";
            return SubdivisionLoop( model );
        }
    }
    //-------------------------------------------------------------------
    // Loop Subdivision for XPolygonal Model
    {
        OpenGL::XPolygonalModel<T> * const  model = 
            dynamic_cast< OpenGL::XPolygonalModel<T> * const >( prModel );
        if ( model ) {
            std::cout << "Do Loop Subdivision to An XPolygonal Model\n";
            return SubdivisionLoop( model );
        }
    }
    //-------------------------------------------------------------------
    // Loop Subdivision for Polygonal Model
    {
        OpenGL::PolygonalModel<T> * const  model = 
            dynamic_cast< OpenGL::PolygonalModel<T> * const >( prModel );
        if ( model ) {
            std::cout << "Do Loop Subdivision to A Polygonal Model\n";
            return SubdivisionLoop( model );
        }
    }
    //-------------------------------------------------------------------
    std::cout << "WARNING: Failed Doing Loop Subdivision --> the model is unrecognized!"
        << std::endl;
    return false;
}
//-----------------------------------------------------------------------------
template <typename T>
bool CGAlgo<T>::SubdivisionLoop ( OpenGL::PolygonalModel<T> * prModel )
{
    return false;
}
//-----------------------------------------------------------------------------
template <typename T>
bool CGAlgo<T>::SubdivisionLoop ( OpenGL::XPolygonalModel<T> * prModel )
{
    return false;
}
//-----------------------------------------------------------------------------
template <typename T>
bool CGAlgo<T>::SubdivisionLoop ( OpenGL::HalfEdgeModel<T> * prModel )
{
    return false;
}
//-----------------------------------------------------------------------------
//=============================================================================




//=============================================================================
//-----------------------------------------------------------------------------
template <typename T>
void CGAlgo<T>::Physics::CalVolumeAndCentroid_BasedOnTetrahedron (
    OpenGL::OpenGLModel<T> *    prModel,    
    T &                         Volume,     
    Vector3<T> &                Centroid    
)
{
    std::cout << "CalVolumeAndCentroid_BasedOnTetrahedron\n";
    OpenGL::HalfEdgeModel<T> * pHEModel = dynamic_cast< OpenGL::HalfEdgeModel<T> * >( prModel );
    if ( pHEModel ) {
        CalVolumeAndCentroid_BasedOnTetrahedron( pHEModel, Volume, Centroid );
        return;
    }
    OpenGL::XPolygonalModel<T> * pXPolyModel = dynamic_cast< OpenGL::XPolygonalModel<T> * >( prModel );
    if ( pXPolyModel ) {
        std::cout << "CGAlgo<T>::Physics::CalVolumeAndCentroid_BasedOnTetrahedron(...) -- ";
        std::cout << "(" << __FILE__ << ", " << __LINE__ << ") -- ";
        std::cout << "NOT IMPLEMENTED YET\n";
        return;
    }
    OpenGL::PolygonalModel<T> * pPolyModel = dynamic_cast< OpenGL::PolygonalModel<T> * >( prModel );
    if ( pPolyModel ) {
        std::cout << "CGAlgo<T>::Physics::CalVolumeAndCentroid_BasedOnTetrahedron(...) -- ";
        std::cout << "(" << __FILE__ << ", " << __LINE__ << ") -- ";
        std::cout << "NOT IMPLEMENTED YET\n";
        return;
    }
}
//-----------------------------------------------------------------------------
template <typename T>
void CGAlgo<T>::Physics::CalVolumeAndCentroid_BasedOnTetrahedron (
    OpenGL::HalfEdgeModel<T> *  prModel,    
    T &                         Volume,     
    Vector3<T> &                Centroid    
)
{
    T TetVol = 0;
    Vector3<T> sumOfTetCentroidMultWithTetVol;
    Volume = 0;

    //std::vector< HEVertex<T> const * > vertices;
    std::vector< HEVertex<T> * const > vertices;
    TAPs::HEFace<T> * pFace = prModel->GetFaceList()->Head();
    while ( pFace ) {

        // volume of tetrahedron = |(a-d)*((b-d)^(c-d))| / 6
        // d := (0,0,0)
        vertices = pFace->GetPtrsToVertices();
        TetVol = vertices[0]->GetPosition() * ( vertices[1]->GetPosition() ^ vertices[2]->GetPosition() ) / T(6.0);
        Volume += TetVol;

        // centroid of tetrahedron = (a + b + c + d ) / 4
        sumOfTetCentroidMultWithTetVol += ( vertices[0]->GetPosition() + vertices[1]->GetPosition() + vertices[2]->GetPosition() ) * T(0.25) * TetVol;

        pFace = pFace->Next();
    }
    Centroid = sumOfTetCentroidMultWithTetVol / Volume;
}
//-----------------------------------------------------------------------------
//=============================================================================




//=============================================================================
//=============================================================================
//-----------------------------------------------------------------------------
//-----------------------------------------------------------------------------
//=============================================================================
END_NAMESPACE_TAPs
//-----------------------------------------------------------------------------
//34567890123456789012345678901234567890123456789012345678901234567890123456789
//--+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----