TAPsAdvSimSupport_DATA.hpp

Go to the documentation of this file.

/******************************************************************************
TAPsAdvSimSupport_DATA.hpp
******************************************************************************/
/******************************************************************************
SUKITTI PUNAK   (09/30/2010)
UPDATE          (02/28/2011)
******************************************************************************/
#ifndef TAPs_ADV_SIM_SUPPORT_DATA_HPP
#define TAPs_ADV_SIM_SUPPORT_DATA_HPP

#include <TAPs/Core/TAPsLib.hpp>

#include "TAPsListOfAdvSimConsts.hpp"
#include "TAPs/GeometricDataStructure/HalfEdge/TAPsHalfEdgeDataStructure.hpp"

//#include "TAPs/Support/TAPsAdvSimConstraint.hpp"
#include "AdvSimConstraints/TAPsAdvSimConstraint_HEVertexVsHEVertex.hpp"
#include "AdvSimConstraints/TAPsAdvSimConstraint_VertexIdVsVertexId.hpp"
#include "AdvSimConstraints/TAPsAdvSimConstraint_VertexIdVsHEVertex.hpp"

//#include <TAPs/OpenGLModel/TAPsModelSuture.hpp>
//#include <TAPs/ModelSpecific/TAPsModelForSurgery.hpp>
#include <TAPs/Model/TAPsModelElasticRod.hpp>
#include <TAPs/Model/TAPsModelDefBasedOnFEMTet.hpp>
#include <TAPs/Model/TAPsModelDefBasedOnFEMHex.hpp>
#include <TAPs/Model/RigidBodies/TAPsRigidBodyDynamics.hpp>


BEGIN_NAMESPACE_TAPs
//=============================================================================
template <typename T, typename DATA = unsigned char>
class AdvSimSupport_DATA {
//=============================================================================
public:
    //=========================================================================
    class IDForModel {
    public:
        IDForModel ( Enum::ModelType type, unsigned int slot ) : Type(type), Slot(slot) {}
        Enum::ModelType Type;   
        unsigned int    Slot;   
    };

    typedef std::vector< IDForModel >   ModelIDs;

    typedef std::list< PointForce<T> * >    ListOfPointForces;

    // Data for controlling forward and backward movements of an Elastic Rod
    class ERNodeProps {
    public:
        ERNodeProps () : DiffLength(0), IsConstrained(false), MoveStatus(0), GroupID(0) {}
        void Reset () { DiffLength=0; IsConstrained=false; MoveStatus=0; GroupID=0; }
        T               DiffLength;         
        bool            IsConstrained;      
        unsigned char   MoveStatus;         

        unsigned short  GroupID;    
    };

    typedef std::vector< ERNodeProps >  ERNodePropsList;

    class ModelLists {
    public:
        //std::vector< OpenGL::ModelSuture<T> * >   m_ListOfSutureModelsBasedOnMSS;     //!< list of suture models based on mass-spring system (currently is not added to this simulation)
        //std::vector< ModelForSurgery<T> * >       m_ListOfModelsForSurgeryBasedOnMSS; //!< list of models for surgery based on mass-spring system  (currently is not added to this simulation)
        //std::vector< OpenGL::HETriMeshOneModelMultiParts<T> * >   m_ListOfHETriMeshOneModelMultiParts;    //!< list of sub models for surgery based on mass-spring system  (currently is not added to this simulation)

        void ResetAllERNodePropsList ()
        {
            for ( unsigned int i = 0; i < AllERNodePropsList.size(); ++i ) {
                for ( unsigned int j = 0; j < AllERNodePropsList[i].size(); ++j ) {
                    AllERNodePropsList[i][j].Reset();
                }
            }
        }

        void ResetERNodePropsListOfERModel ( unsigned int ERloc )
        {
            if ( ERloc < AllERNodePropsList.size() ) {
                for ( unsigned int j = 0; j < AllERNodePropsList[ERloc].size(); ++j ) {
                    AllERNodePropsList[ERloc][j].Reset();
                }
            }
        }

        std::vector< ERNodePropsList >  AllERNodePropsList; 



        std::vector< ModelElasticRod<T> * >         BasedOnER;  
        std::vector< ModelDefBasedOnFEM<T,DATA> * > BasedOnFEM; 
        std::vector< RigidBodyDynamics<T> * >       BasedOnRBD; 
    };

    //-------------------------------------------------------------------------
    // (3D-array) list of model-model interaction data
    // SEMI-EXCLUSIVE DATA ARRAYS FOR INTERACTION

    // 1) ERModel vs HEModel
    //===============================================================
    typedef AdvSimConstraint_VertexIdVsHEVertex<T>          Constraint_ERModelVsHEModel;
    typedef std::vector< Constraint_ERModelVsHEModel * >    InteractionConstraints_ERModelVsHEModel;
    //---------------------------------------------------------------
    class InteractionConstants_ERModelVsHEModel {
    public:
        InteractionConstants_ERModelVsHEModel ()
            : DiffLengthForTriggeringForwardMove( 0.05 )
            , DiffLengthForTriggeringBackwardMove( 0.05 )
        {}

        T   DiffLengthForTriggeringForwardMove;     
        T   DiffLengthForTriggeringBackwardMove;    
    };
    class Interaction_ERModelVsHEModel {
    public:
        Interaction_ERModelVsHEModel (
            InteractionConstants_ERModelVsHEModel &     _Constants  = InteractionConstants_ERModelVsHEModel(),
            InteractionConstraints_ERModelVsHEModel &   _Constraints= InteractionConstraints_ERModelVsHEModel()
        ) : Constants( _Constants ), Constraints( _Constraints )
        {}

        InteractionConstants_ERModelVsHEModel   Constants;
        InteractionConstraints_ERModelVsHEModel Constraints;
    };
    typedef std::vector< std::vector< Interaction_ERModelVsHEModel > >  ListOfInteractions_ERModelVsHEModel;
    //---------------------------------------------------------------

    // 2) ERModel vs IPGModel
    //===============================================================
    typedef AdvSimConstraint_VertexIdVsVertexId<T>          Constraint_ERModelVsIPGModel;
    typedef std::vector< Constraint_ERModelVsIPGModel * >   InteractionConstraints_ERModelVsIPGModel;
    //---------------------------------------------------------------
    class InteractionConstants_ERModelVsIPGModel {
    public:
        InteractionConstants_ERModelVsIPGModel ()   {}
    };
    class Interaction_ERModelVsIPGModel {
    public:
        Interaction_ERModelVsIPGModel (
            InteractionConstants_ERModelVsIPGModel &    _Constants  = InteractionConstants_ERModelVsIPGModel(),
            InteractionConstraints_ERModelVsIPGModel &  _Constraints= InteractionConstraints_ERModelVsIPGModel()
        ) : Constants( _Constants ), Constraints( _Constraints )
        {}

        InteractionConstants_ERModelVsIPGModel      Constants;
        InteractionConstraints_ERModelVsIPGModel    Constraints;
    };
    typedef std::vector< std::vector< Interaction_ERModelVsIPGModel > > ListOfInteractions_ERModelVsIPGModel;
    //---------------------------------------------------------------

    // 3) IPGModel vs HEModel
    //===============================================================
    typedef AdvSimConstraint_VertexIdVsHEVertex<T>          Constraint_IPGModelVsHEModel;
    typedef std::vector< Constraint_IPGModelVsHEModel * >   InteractionConstraints_IPGModelVsHEModel;
    //---------------------------------------------------------------
    class InteractionConstants_IPGModelVsHEModel {
    public:
        InteractionConstants_IPGModelVsHEModel ()   {}
    };
    class Interaction_IPGModelVsHEModel {
    public:
        Interaction_IPGModelVsHEModel (
            InteractionConstants_IPGModelVsHEModel &    _Constants  = InteractionConstants_IPGModelVsHEModel(),
            InteractionConstraints_IPGModelVsHEModel &  _Constraints= InteractionConstraints_IPGModelVsHEModel()
        ) : Constants( _Constants ), Constraints( _Constraints )
        {}

        InteractionConstants_IPGModelVsHEModel      Constants;
        InteractionConstraints_IPGModelVsHEModel    Constraints;
    };
    typedef std::vector< std::vector< Interaction_IPGModelVsHEModel > > ListOfInteractions_IPGModelVsHEModel;
    //---------------------------------------------------------------

    // 4) ERModel vs ERModel
    //===============================================================
    typedef AdvSimConstraint_VertexIdVsVertexId<T>          Constraint_ERModelVsERModel;
    typedef std::vector< Constraint_ERModelVsERModel * >    InteractionConstraints_ERModelVsERModel;
    //---------------------------------------------------------------
    class InteractionConstants_ERModelVsERModel {
    public:
        InteractionConstants_ERModelVsERModel ()    {}
    };
    class Interaction_ERModelVsERModel {
    public:
        Interaction_ERModelVsERModel (
            InteractionConstants_ERModelVsERModel &     _Constants  = InteractionConstants_ERModelVsERModel(),
            InteractionConstraints_ERModelVsERModel &   _Constraints= InteractionConstraints_ERModelVsERModel()
        ) : Constants( _Constants ), Constraints( _Constraints )
        {}

        InteractionConstants_ERModelVsERModel   Constants;
        InteractionConstraints_ERModelVsERModel Constraints;
    };
    typedef std::vector< std::vector< Interaction_ERModelVsERModel > >  ListOfInteractions_ERModelVsERModel;
    //---------------------------------------------------------------

    // 5) HEModel vs HEModel
    //===============================================================
    typedef AdvSimConstraint_HEVertexVsHEVertex<T>          Constraint_HEModelVsHEModel;
    typedef std::vector< Constraint_HEModelVsHEModel * >    InteractionConstraints_HEModelVsHEModel;
    //---------------------------------------------------------------
    class InteractionConstants_HEModelVsHEModel {
    public:
        InteractionConstants_HEModelVsHEModel ()    {}
    };
    class Interaction_HEModelVsHEModel {
    public:
        Interaction_HEModelVsHEModel (
            InteractionConstants_HEModelVsHEModel &     _Constants  = InteractionConstants_HEModelVsHEModel(),
            InteractionConstraints_HEModelVsHEModel &   _Constraints= InteractionConstraints_HEModelVsHEModel()
        ) : Constants( _Constants ), Constraints( _Constraints )
        {}

        InteractionConstants_HEModelVsHEModel   Constants;
        InteractionConstraints_HEModelVsHEModel Constraints;
    };
    typedef std::vector< std::vector< Interaction_HEModelVsHEModel > >  ListOfInteractions_HEModelVsHEModel;
    //---------------------------------------------------------------

    // 6) IPGModel vs IPGModel
    //===============================================================
    typedef AdvSimConstraint_VertexIdVsVertexId<T>          Constraint_IPGModelVsIPGModel;
    typedef std::vector< Constraint_IPGModelVsIPGModel * >  InteractionConstraints_IPGModelVsIPGModel;
    //---------------------------------------------------------------
    class InteractionConstants_IPGModelVsIPGModel {
    public:
        InteractionConstants_IPGModelVsIPGModel ()  {}
    };
    class Interaction_IPGModelVsIPGModel {
    public:
        Interaction_IPGModelVsIPGModel (
            InteractionConstants_IPGModelVsIPGModel &   _Constants  = InteractionConstants_IPGModelVsIPGModel(),
            InteractionConstraints_IPGModelVsIPGModel & _Constraints= InteractionConstraints_IPGModelVsIPGModel();
        ) : Constants( _Constants ), Constraints( _Constraints )
        {}

        InteractionConstants_IPGModelVsIPGModel Constants;
        InteractionConstraints_IPGModelVsIPGModel   Constraints;
    };
    typedef std::vector< std::vector< Interaction_IPGModelVsIPGModel > >    ListOfInteractions_IPGModelVsIPGModel;
    //---------------------------------------------------------------

//#ifdef    TAPs_MODEL_SURGICAL_SUTURE_WITH_HEAD_NEEDLE_HPP
//  typedef std::vector< AdvSimConstraint_ForModelSurgicalSutureWithHeadNeedle<T> * >   ListOfInteraction_SwHDModelVsHEModel;
//#endif//TAPs_MODEL_SURGICAL_SUTURE_WITH_HEAD_NEEDLE_HPP

    class InteractionLists {
    public:
        ListOfInteractions_ERModelVsHEModel     ERModelVsHEModel;
        ListOfInteractions_ERModelVsIPGModel    ERModelVsIPGModel;
        ListOfInteractions_IPGModelVsHEModel    IPGModelVsHEModel;
        ListOfInteractions_ERModelVsERModel     ERModelVsERModel;
        ListOfInteractions_HEModelVsHEModel     HEModelVsHEModel;
        ListOfInteractions_IPGModelVsIPGModel   IPGModelVsIPGModel;

    //#ifdef    TAPs_MODEL_SURGICAL_SUTURE_WITH_HEAD_NEEDLE_HPP
    //  ListOfInteraction_SwHDModelVsHEModel    SwHDModelVsHEModel;
    //#endif//TAPs_MODEL_SURGICAL_SUTURE_WITH_HEAD_NEEDLE_HPP
    };
    //=========================================================================

    //=========================================================================
    // State variables
    //-------------------------------------------------------------------------

    // Used for interaction control of the suture's head needle for class ModelSurgicalSutureWithHeadNeedle<T>
    static bool STATE_EnforcePositionA; 
    static bool STATE_EnforcePositionB; 
    static T    STATE_ForceRatio;       
    static T    STATE_ForceScaleA;      
    static T    STATE_ForceScaleB;      
    static T    STATE_ForceThresholdA;  
    static T    STATE_ForceThresholdB;  

    // Used for interaction control of the suture's thread for class ModelSurgicalSutureWithHeadNeedle<T>
    static bool STATE_EnforcePositionA_2;   
    static bool STATE_EnforcePositionB_2;   
    static T    STATE_ForceRatio_2;         
    static T    STATE_ForceScaleA_2;        
    static T    STATE_ForceScaleB_2;        
    static T    STATE_ForceThresholdA_2;    
    static T    STATE_ForceThresholdB_2;    

    //-------------------------------------------------------------------------
    // State variables
    //=========================================================================


    //=========================================================================
    // Temporary data for helping processing interactions
    //-------------------------------------------------------------------------
    static std::vector< Vector3<T> >    TEMP_ListOfPoints;              
    static std::vector< HEFace<T>* >    TEMP_ListOfHEFaces;             
    static std::vector< HEVertex<T>* >  TEMP_ListOfHEVertices;          
    static std::vector< T >             TEMP_ListOfRatios;              
    static std::vector< T >             TEMP_ListOfIntersectionAngles;  
    //-------------------------------------------------------------------------
    // Temporary data for helping processing interactions
    //=========================================================================


//=============================================================================
}; // END CLASS AdvSimSupport_DATA
//=============================================================================


//=========================================================================
// State variables
//-------------------------------------------------------------------------
template <typename T, typename DATA> bool   AdvSimSupport_DATA<T,DATA>::STATE_EnforcePositionA = false;
template <typename T, typename DATA> bool   AdvSimSupport_DATA<T,DATA>::STATE_EnforcePositionB = false;
template <typename T, typename DATA> T      AdvSimSupport_DATA<T,DATA>::STATE_ForceRatio = 0.5;
template <typename T, typename DATA> T      AdvSimSupport_DATA<T,DATA>::STATE_ForceScaleA = 1;
template <typename T, typename DATA> T      AdvSimSupport_DATA<T,DATA>::STATE_ForceScaleB = 1;
template <typename T, typename DATA> T      AdvSimSupport_DATA<T,DATA>::STATE_ForceThresholdA = 1;
template <typename T, typename DATA> T      AdvSimSupport_DATA<T,DATA>::STATE_ForceThresholdB = 1;

template <typename T, typename DATA> bool   AdvSimSupport_DATA<T,DATA>::STATE_EnforcePositionA_2 = false;
template <typename T, typename DATA> bool   AdvSimSupport_DATA<T,DATA>::STATE_EnforcePositionB_2 = false;
template <typename T, typename DATA> T      AdvSimSupport_DATA<T,DATA>::STATE_ForceRatio_2 = 0.5;
template <typename T, typename DATA> T      AdvSimSupport_DATA<T,DATA>::STATE_ForceScaleA_2 = 1;
template <typename T, typename DATA> T      AdvSimSupport_DATA<T,DATA>::STATE_ForceScaleB_2 = 1;
template <typename T, typename DATA> T      AdvSimSupport_DATA<T,DATA>::STATE_ForceThresholdA_2 = 1;
template <typename T, typename DATA> T      AdvSimSupport_DATA<T,DATA>::STATE_ForceThresholdB_2 = 1;
//-------------------------------------------------------------------------
// State variables
//=========================================================================


//=========================================================================
// Static temporary data for helping processing interactions
//-------------------------------------------------------------------------
template <typename T, typename DATA> std::vector< Vector3<T> >      AdvSimSupport_DATA<T,DATA>::TEMP_ListOfPoints;
template <typename T, typename DATA> std::vector< HEFace<T>* >      AdvSimSupport_DATA<T,DATA>::TEMP_ListOfHEFaces;
template <typename T, typename DATA> std::vector< HEVertex<T>* >    AdvSimSupport_DATA<T,DATA>::TEMP_ListOfHEVertices;
template <typename T, typename DATA> std::vector< T >               AdvSimSupport_DATA<T,DATA>::TEMP_ListOfRatios;
template <typename T, typename DATA> std::vector< T >               AdvSimSupport_DATA<T,DATA>::TEMP_ListOfIntersectionAngles;
//-------------------------------------------------------------------------
// Static temporary data for helping processing interactions
//=========================================================================


END_NAMESPACE_TAPs
//-----------------------------------------------------------------------------
// Include definition if TAPs_USE_EXPORT is not defined
//#if !defined( TAPs_USE_EXPORT )
//  #include "TAPsAdvSimSupport_DATA.cpp"
//#endif
//-----------------------------------------------------------------------------
#endif//TAPs_ADV_SIM_SUPPORT_HPP_DS
//34567890123456789012345678901234567890123456789012345678901234567890123456789
//--+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----