TAPsAdvSimConstraints.hpp

Go to the documentation of this file.

/******************************************************************************
TAPsAdvSimConstraints.hpp
******************************************************************************/
/******************************************************************************
SUKITTI PUNAK   (03/16/2011)
UPDATE          (03/31/2011)
******************************************************************************/
#ifndef TAPs_ADV_SIM_CONSTRAINTS_HPP
#define TAPs_ADV_SIM_CONSTRAINTS_HPP

#include <TAPs/Model/TAPsModelElasticRod.hpp>

BEGIN_NAMESPACE_TAPs
//=============================================================================


//-----------------------------------------------------------------------------
template <typename T>
class AdvSimConstraint_ConstrainedPoint {
public:
    AdvSimConstraint_ConstrainedPoint ();

    AdvSimConstraint_ConstrainedPoint ( AdvSimConstraint_ConstrainedPoint<T> const &orig );

    AdvSimConstraint_ConstrainedPoint<T> & operator= ( AdvSimConstraint_ConstrainedPoint<T> const &orig );

    // Destructor
    ~AdvSimConstraint_ConstrainedPoint ();

    std::string StrInfo () const;

    int         PointID;            
    Vector3<T>  HomePosition;       
    Vector3<T>  TargetPosition;     
};
//-----------------------------------------------------------------------------


//-----------------------------------------------------------------------------
template <typename T>
class AdvSimConstraint_GroupOfConstrainedPoints {
public:
    AdvSimConstraint_GroupOfConstrainedPoints ();

    AdvSimConstraint_GroupOfConstrainedPoints ( AdvSimConstraint_GroupOfConstrainedPoints<T> const &orig );

    AdvSimConstraint_GroupOfConstrainedPoints<T> & operator= ( AdvSimConstraint_GroupOfConstrainedPoints<T> const &orig );

    // Destructor
    ~AdvSimConstraint_GroupOfConstrainedPoints ();

    std::string StrInfo () const;

    void AddConstrainedPoint ( AdvSimConstraint_ConstrainedPoint<T> & CPt = AdvSimConstraint_ConstrainedPoint<T>() )
    {
        ListOfConstrainedPoints.push_back( CPt );
    }

    AdvSimConstraint_ConstrainedPoint<T> & ConstrainedPoint ( unsigned int i ) const
    {
        assert( i < ListOfConstrainedPoints.size() );
        return ListOfConstrainedPoints[i];
    }

    AdvSimConstraint_ConstrainedPoint<T> & ConstrainedPoint ( unsigned int i )
    {
        assert( i < ListOfConstrainedPoints.size() );
        return ListOfConstrainedPoints[i];
    }

    std::vector< AdvSimConstraint_ConstrainedPoint<T> > ListOfConstrainedPoints;    
    int     ModelID;    
    T       Ratio;      


#if defined(__gl_h_) || defined(__GL_H__)
public:
    virtual void Draw ( Vector3<T> & color = Vector3<T>(1,0,0), T PtSize = 20 ) const;
#endif
};
//-----------------------------------------------------------------------------


//-----------------------------------------------------------------------------
template <typename T>
class AdvSimConstraint_ConstrainedER {
public:
    AdvSimConstraint_ConstrainedER () : pModelER( NULL ) {}

    ~AdvSimConstraint_ConstrainedER () {}

    ModelElasticRod<T> * pModelER;                  
    std::vector< unsigned int > ListOfCPtsGrpIDs;   
    //std::vector< T >          ListOfRatios;       //!< the ratio of the enforcement; 0 means the suture points strictly follow the constrained points
};
//-----------------------------------------------------------------------------


//=============================================================================
template <typename T>
class AdvSimConstraints {
//=============================================================================
public:

    typedef std::vector< AdvSimConstraint_ConstrainedPoint<T> >         VectorOfCPts;
    typedef std::vector< AdvSimConstraint_GroupOfConstrainedPoints<T> > VectorOfCPtGroups;

    typedef std::vector< AdvSimConstraint_ConstrainedER<T> >    VectorOfConstrainedERs;

    // Member Functions -------------------------------------------------------
    friend std::ostream & operator<< ( std::ostream &output, AdvSimConstraints<T> const &obj )
    {
        output << obj.StrInfo();
        return output;
    }

    //-------------------------------------------------------------------------
    AdvSimConstraints ();

    //AdvSimConstraints ( AdvSimConstraints<T> const &orig );

    ~AdvSimConstraints ();
    //-------------------------------------------------------------------------
    std::string StrInfo () const;
    //-------------------------------------------------------------------------
    //AdvSimConstraints<T> & operator= ( AdvSimConstraints<T> const &orig );

    //-------------------------------------------------------------------------
    // Get/Set Functions

    const VectorOfCPtGroups & RefToListOfConstrainedPtGrps () const { return m_ListOfConstrainedPtGrps; }
    VectorOfCPtGroups   & RefToListOfConstrainedPtGrps ()           { return m_ListOfConstrainedPtGrps; }

    AdvSimConstraint_GroupOfConstrainedPoints<T> & RefToConstrainedPtGrp ( unsigned int i ) const
    {
        assert( i < m_ListOfConstrainedPtGrps.size() );
        return m_ListOfConstrainedPtGrps[i];
    }

    AdvSimConstraint_GroupOfConstrainedPoints<T> & RefToConstrainedPtGrp ( unsigned int i )
    {
        assert( i < m_ListOfConstrainedPtGrps.size() );
        return m_ListOfConstrainedPtGrps[i];
    }

    unsigned int AddConstrainedPointGroup ( AdvSimConstraint_GroupOfConstrainedPoints<T> & CPtsGrp = AdvSimConstraint_GroupOfConstrainedPoints<T>() )
    {
        m_ListOfConstrainedPtGrps.push_back( CPtsGrp );
        return m_ListOfConstrainedPtGrps.size() - 1;
    }

    //-------------------------------------------------------------------------
    // ModelElasticRod

    unsigned int AddConstrainedModelER (
        ModelElasticRod<T> & ModelER    
    );

    void AssociateCPtsGrpWithModelER (
        unsigned int CPtsGrp_ID,    
        unsigned int ModelER_ID     
    );

    //-------------------------------------------------------------------------
    // Functions for Advanced Simulation

    inline void EnforceAllConstrainedModels_A ();
    inline void EnforceAllConstrainedModels_B ();
    inline void EnforceAllConstrainedModels_B__ByForce ( Vector3<T> & Force = Vector3<T>(1,1,1) );

    inline void EnforceAllConstrainedModelERs_A ();
    inline void EnforceAllConstrainedModelERs_B ();
    inline void EnforceAllConstrainedModelERs_B__ByForce ( Vector3<T> & Force = Vector3<T>(1,1,1) );


    void ClearAllConstraints ();

    // Data Members -----------------------------------------------------------
//=============================================================================
protected:
    // Member Functions -------------------------------------------------------

    inline void EnforceConstrainedModelER_A (
        AdvSimConstraint_ConstrainedER<T> & ConstrainedER   
    );

    inline void EnforceConstraints_A (
        ModelElasticRod<T> * pER,   
        int groupID     
    );
    inline void EnforceConstraints_B (
        ModelElasticRod<T> * pER,   
        int groupID     
    );
    inline void EnforceConstraints_B__ByForce (
        ModelElasticRod<T> * pER,   
        int groupID,                
        Vector3<T> & Force          
    );

    inline void EnforceConstraints_A__NotForMultiGrp (
        ModelElasticRod<T> * pER,   
        int groupID     
    );
    inline void EnforceConstraints_B__NotForMultiGrp (
        ModelElasticRod<T> * pER,   
        int groupID     
    );

    void ClearAllConstriantsOnERatID ( unsigned int ERid );

    // Data Members -----------------------------------------------------------
    VectorOfCPtGroups       m_ListOfConstrainedPtGrps;  
    VectorOfConstrainedERs  m_ListOfConstrainedERs;     
//=============================================================================
private:
    // Member Functions -------------------------------------------------------
    // Data Members -----------------------------------------------------------
//=============================================================================

#if defined(__gl_h_) || defined(__GL_H__)
public:
    virtual void Draw ( Vector3<T> & color = Vector3<T>(1,0,0), T PtSize = 20 ) const;
#endif

//=============================================================================
}; // END CLASS AdvSimConstraints
//=============================================================================
END_NAMESPACE_TAPs
//-----------------------------------------------------------------------------
#include "TAPsAdvSimConstraints.cpp"

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