TAPsTetBarycentricCoordsForVertexRef.cpp

Go to the documentation of this file.

/******************************************************************************
TAPsTetBarycentricCoordsForVertexRef.cpp
******************************************************************************/
/******************************************************************************
SUKITTI PUNAK   (03/16/2010)
UPDATE          (09/28/2010)
******************************************************************************/
#include "TAPsTetBarycentricCoordsForVertexRef.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
//=============================================================================
//-----------------------------------------------------------------------------
//template <typename T>
//TetBarycentricCoordsForVertexRef<T>::TetBarycentricCoordsForVertexRef ()
//  :
//  m_pVertex( NULL ), m_pTetrahedron( NULL )
//{}
//-----------------------------------------------------------------------------
template <typename T>
TetBarycentricCoordsForVertexRef<T>::TetBarycentricCoordsForVertexRef ( 
        Vector3<T> & vertex,    
        FEM::Tetrahedron<T> & tetrahedron,  
        T u,    
        T v,    
        T w     
)   : BaseBarycentricCoordsForVertexRef( vertex, u, v, w), m_pTetrahedron( &tetrahedron )
{}
//-----------------------------------------------------------------------------
template <typename T>
TetBarycentricCoordsForVertexRef<T>::TetBarycentricCoordsForVertexRef ( 
        Vector3<T> & vertex,    
        FEM::Tetrahedron<T> * pTetrahedron, 
        T u,    
        T v,    
        T w     
)   : BaseBarycentricCoordsForVertexRef( vertex, u, v, w), m_pTetrahedron( pTetrahedron )
{}
//-----------------------------------------------------------------------------
template <typename T>
TetBarycentricCoordsForVertexRef<T>::TetBarycentricCoordsForVertexRef (
    TetBarycentricCoordsForVertexRef<T> const & orig
)
    : BaseBarycentricCoordsForVertexRef( orig ),
      m_pTetrahedron( orig.m_pTetrahedron )
{}
//-----------------------------------------------------------------------------
template <typename T>
TetBarycentricCoordsForVertexRef<T>::~TetBarycentricCoordsForVertexRef ()
{}  
//-----------------------------------------------------------------------------
template <typename T>
TetBarycentricCoordsForVertexRef<T> & TetBarycentricCoordsForVertexRef<T>::operator= (
    TetBarycentricCoordsForVertexRef<T> const & orig
)
{
    // Inherited properties
    m_pVertex = orig.m_pVertex;
    m_coords = orig.m_coords;
    m_uiElementID = orig.m_uiElementID;
    // Added properties
    m_pTetrahedron = orig.m_pTetrahedron;
    return *this;
}
//-----------------------------------------------------------------------------
template <typename T>
std::string TetBarycentricCoordsForVertexRef<T>::StrInfo () const
{
    std::ostringstream ss;
    ss << "TetBarycentricCoordsForVertexRef<" << typeid(T).name() << ">";
    ss << "ref vertex: " << m_pVertex << " " << *m_pVertex
        << ") with coordinates " << &m_coords << " " << m_coords 
        << " and elementID " << m_uiElementID
        << " associated with tetrahedron ";
    if ( m_pTetrahedron ) {
        ss << *m_pTetrahedron;
    }
    else { 
        ss << "NULL\n";
    }
    return ss.str();
}
//-----------------------------------------------------------------------------
template <typename T>
void TetBarycentricCoordsForVertexRef<T>::ComputeNewVertexLocation ()
{
    Vector3<T> &T0  = *m_pTetrahedron->DeformedNode(0);
    Vector3<T> &T10 = *m_pTetrahedron->DeformedNode(1) - T0;
    Vector3<T> &T20 = *m_pTetrahedron->DeformedNode(2) - T0;
    Vector3<T> &T30 = *m_pTetrahedron->DeformedNode(3) - T0;

    T x = T10[0]*m_coords[0] + T20[0]*m_coords[1] + T30[0]*m_coords[2] + T0[0];
    T y = T10[1]*m_coords[0] + T20[1]*m_coords[1] + T30[1]*m_coords[2] + T0[1];
    T z = T10[2]*m_coords[0] + T20[2]*m_coords[1] + T30[2]*m_coords[2] + T0[2];

    m_pVertex->SetXYZ( x, y, z );
}
//-----------------------------------------------------------------------------
template <typename T>
Vector3<T> TetBarycentricCoordsForVertexRef<T>::GetNewVertexLocation () const
{
    Vector3<T> &T0  = *m_pTetrahedron->DeformedNode(0);
    Vector3<T> &T10 = *m_pTetrahedron->DeformedNode(1) - T0;
    Vector3<T> &T20 = *m_pTetrahedron->DeformedNode(2) - T0;
    Vector3<T> &T30 = *m_pTetrahedron->DeformedNode(3) - T0;

    T x = T10[0]*m_coords[0] + T20[0]*m_coords[1] + T30[0]*m_coords[2] + T0[0];
    T y = T10[1]*m_coords[0] + T20[1]*m_coords[1] + T30[1]*m_coords[2] + T0[1];
    T z = T10[2]*m_coords[0] + T20[2]*m_coords[1] + T30[2]*m_coords[2] + T0[2];

    return Vector3<T>( x, y, z );
}
//-----------------------------------------------------------------------------
//=============================================================================
END_NAMESPACE_TAPs
//34567890123456789012345678901234567890123456789012345678901234567890123456789
//--+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----