TAPsCDFnsForInteractionPointGroup.cpp

Go to the documentation of this file.

/******************************************************************************
TAPsCDFnsForInteractionPointGroup.cpp
SUKITTI PUNAK   (08/27/2010)
UPDATE          (09/07/2010)
******************************************************************************/
#include "TAPsCDFnsForInteractionPointGroup.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.

#ifdef  TAPs_INTERACTION_POINT_GROUP_HPP
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

BEGIN_NAMESPACE_TAPs__CD__Fn
//=============================================================================
//-----------------------------------------------------------------------------




#ifdef  TAPs_BOUNDING_VOLUME_HPP
//=============================================================================
//-----------------------------------------------------------------------------
template <typename T>
bool CD_IPGroup_with_BoundingVolume (
    InteractionPointGroup<T> * const    pIPG,           
    BoundingVolume<T> * const   pBV,                    
    std::vector< PointForce<T> > & ListOfPointForces,   
    TransformationSupport<T> * const pTrx_BV            
)
{
    //---------------------------------------------------------------
    // For transforming the BV to world space
    TransformationSupport<T> BV_Trx;
    // pTrx_BV, if the extra transformation support is provided
    if ( pTrx_BV ) {
        BV_Trx.RefToMatrixTransform() *= pTrx_BV->GetMatrixTransform();
    }
    // BV's transformation support
    BV_Trx.RefToMatrixTransform() *= pBV->GetTransform().GetMatrixTransform();

    //---------------------------------------------------------------
    // CD with BV
    Vector3<T> Dist;
    std::vector< InteractionPoint<T> >::iterator pt = pIPG->GetTheSetOfInteractionPoints().begin();
    int countPF = 0;
    int i = 0;
    while ( pt != pIPG->GetTheSetOfInteractionPoints().end() ) {
        bool bCD = pBV->TestPointLocation( pt->GetPosition(), &BV_Trx, &Dist );
        if ( bCD ) {
            Vector3<T> Pos = pIPG->GetTheSetOfInteractionPoints()[i].GetPosition();
            Vector3<T> Force = Dist;
            ListOfPointForces.push_back( PointForce<T>( Pos, Force ) );
            ++countPF;
        }
        ++pt;
        ++i;
    }

    if ( countPF > 0 ) {
        //
        // CR can be implemented here!
        //
        return true;
    }
    return false;
}
// END: CD_IPGroup_with_BoundingVolume(...)
//-----------------------------------------------------------------------------
//=============================================================================
#endif//TAPs_BOUNDING_VOLUME_HPP




#ifdef  TAPs_MULTI_BOUNDING_VOLUME_HPP
//=============================================================================
//-----------------------------------------------------------------------------
template <typename T>
bool CD_IPGroup_with_MultiBoundingVolume (
    InteractionPointGroup<T> * const    pIPG,           
    MultiBoundingVolume<T> * const  pMBV,               
    std::vector< PointForce<T> > & ListOfPointForces,   
    TransformationSupport<T> * const pTrx_MBV           
)
{
    //ListOfPointForces.clear();

    //---------------------------------------------------------------
    // For transforming the MBV to world space
    TransformationSupport<T> MBV_Trx;
    // pTrx_MBV, if the extra transformation support is provided
    if ( pTrx_MBV ) {
        MBV_Trx.RefToMatrixTransform() *= pTrx_MBV->GetMatrixTransform();
    }
    // MBV's transformation support
    MBV_Trx.RefToMatrixTransform() *= pMBV->GetTransform().GetMatrixTransform();

    //---------------------------------------------------------------
    // Traverse the list of bounding volumes
    std::vector< BoundingVolume<T> * >::iterator pos = const_cast< MultiBoundingVolume<T> * >( pMBV )->GetBoundingVolumeList().begin();

    int countPF = 0;
    while ( pos != pMBV->GetBoundingVolumeList().end() ) {

        Vector3<T> Dist;
        std::vector< InteractionPoint<T> >::iterator pt = pIPG->GetTheSetOfInteractionPoints().begin();
        int i = 0;
        while ( pt != pIPG->GetTheSetOfInteractionPoints().end() ) {
            bool bCD = (*pos)->TestSphereLocation( pt->GetPosition(), pt->GetRadius(), &MBV_Trx, &Dist );
            if ( bCD ) {
                Vector3<T> Pos = pIPG->GetTheSetOfInteractionPoints()[i].GetPosition();
                Vector3<T> Force = Dist;
                ListOfPointForces.push_back( PointForce<T>( Pos, Force ) );
                ++countPF;
            }
            ++pt;
            ++i;
        }
        ++pos;  // next bounding volume
    }

    if ( countPF > 0 ) {
        //
        // CR can be implemented here!
        //
        return true;
    }
    return false;
}
// END: CD_IPGroup_with_MultiBoundingVolume(...)
//-----------------------------------------------------------------------------
//=============================================================================
#endif//TAPs_MULTI_BOUNDING_VOLUME_HPP




#ifdef  TAPs_MODEL_DEF_BASED_ON_FEM_HPP
//=============================================================================
//-----------------------------------------------------------------------------
template <typename T, typename DATA>
bool CD_IPGroup_with_ModelDefBasedOnFEM (
    InteractionPointGroup<T> * const    pIPG,           
    ModelDefBasedOnFEM<T,DATA> * const  pDMFEM,         
    std::vector< PointForce<T> > & ListOfPointForces    
)
{
    //if ( !pDMFEM->PtrToBVHTreeOfSurfaceMesh() )   return false;
    //return pDMFEM->PtrToBVHTreeOfSurfaceMesh()->TestOverlapWithTillLeafNodes( pIPG, ListOfPointForces );
    return pDMFEM->PtrToBVHTreeOfSurfaceMesh()->TestOverlapWith( pIPG, ListOfPointForces );
}
// END: CD_IPGroup_with_ModelDefBasedOnFEM(...)
//-----------------------------------------------------------------------------
//=============================================================================
#endif//TAPs_MODEL_DEF_BASED_ON_FEM_HPP




//-----------------------------------------------------------------------------
//=============================================================================
END_NAMESPACE_TAPs__CD__Fn
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#endif//TAPs_INTERACTION_POINT_GROUP_HPP
//34567890123456789012345678901234567890123456789012345678901234567890123456789
//--+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----