TAPsBaseBarycentricCoordsForVertexRef.cpp

Go to the documentation of this file.

/******************************************************************************
TAPsBaseBarycentricCoordsForVertexRef.cpp
******************************************************************************/
/******************************************************************************
SUKITTI PUNAK   (04/16/2010)
UPDATE          (06/25/2010)
******************************************************************************/
#include "TAPsBaseBarycentricCoordsForVertexRef.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>
//BaseBarycentricCoordsForVertexRef<T>::BaseBarycentricCoordsForVertexRef ()
//  : m_pVertex( NULL )
//{}
//-----------------------------------------------------------------------------
template <typename T>
BaseBarycentricCoordsForVertexRef<T>::BaseBarycentricCoordsForVertexRef ( 
    Vector3<T> & vertex,    
    T u,    
    T v,    
    T w     
)   : m_pVertex( &vertex ), m_coords( u, v, w ), m_uiElementID( 0 )
{}
//-----------------------------------------------------------------------------
template <typename T>
BaseBarycentricCoordsForVertexRef<T>::BaseBarycentricCoordsForVertexRef (
    BaseBarycentricCoordsForVertexRef<T> const & orig
)   : m_pVertex( orig.m_pVertex ), m_coords( orig.m_coords ), m_uiElementID( orig.m_uiElementID )
{}
//-----------------------------------------------------------------------------
template <typename T>
BaseBarycentricCoordsForVertexRef<T>::~BaseBarycentricCoordsForVertexRef ()
{}  
//-----------------------------------------------------------------------------
template <typename T>
BaseBarycentricCoordsForVertexRef<T> & BaseBarycentricCoordsForVertexRef<T>::operator= (
    BaseBarycentricCoordsForVertexRef<T> const & orig
)
{
    m_pVertex = orig.m_pVertex;
    m_coords = orig.m_coords;
    m_uiElementID = orig.m_uiElementID;
    return *this;
}
//-----------------------------------------------------------------------------
template <typename T>
std::string BaseBarycentricCoordsForVertexRef<T>::StrInfo () const
{
    std::ostringstream ss;
    ss << "BaseBarycentricCoordsForVertexRef<" << typeid(T).name() << ">";
    ss << "ref vertex: " << *m_pVertex
        << ") with coordinates " << m_coords 
        << " and elementID " << m_uiElementID
        << "\n.";
    return ss.str();
}
//-----------------------------------------------------------------------------
template <typename T>
void BaseBarycentricCoordsForVertexRef<T>::ComputeNewVertexLocationBasedOnTetrahedralVertices (
    T x0, T y0, T z0,
    T x1, T y1, T z1,
    T x2, T y2, T z2,
    T x3, T y3, T z3
)
{
    T x = (x1-x0)*m_coords[0] + (x2-x0)*m_coords[1] + (x3-x0)*m_coords[2] + x0;
    T y = (y1-y0)*m_coords[0] + (y2-y0)*m_coords[1] + (y3-y0)*m_coords[2] + y0;
    T z = (z1-z0)*m_coords[0] + (z2-z0)*m_coords[1] + (z3-z0)*m_coords[2] + z0;
    m_pVertex->SetXYZ( x, y, z );
}
//-----------------------------------------------------------------------------
template <typename T>
void BaseBarycentricCoordsForVertexRef<T>::ComputeNewVertexLocationBasedOnHexahedralVertices (
    T x0, T y0, T z0,
    T x1, T y1, T z1,
    T x2, T y2, T z2,
    T x3, T y3, T z3,
    T x4, T y4, T z4,
    T x5, T y5, T z5,
    T x6, T y6, T z6,
    T x7, T y7, T z7
)
{
    T s0 = 1 - m_coords[0];
    T s1 = 1 + m_coords[0];
    T t0 = 1 - m_coords[1];
    T t1 = 1 + m_coords[1];
    T r0 = 1 - m_coords[2];
    T r1 = 1 + m_coords[2];

    T N0 = s0 * t0 * r0;
    T N1 = s1 * t0 * r0;
    T N2 = s1 * t1 * r0;
    T N3 = s0 * t1 * r0;
    T N4 = s0 * t0 * r1;
    T N5 = s1 * t0 * r1;
    T N6 = s1 * t1 * r1;
    T N7 = s0 * t1 * r1;

    T x = 0.125 * ( x0*N0 + x1*N1 + x2*N2 + x3*N3 + x4*N4 + x5*N5 + x6*N6 + x7*N7 );
    T y = 0.125 * ( y0*N0 + y1*N1 + y2*N2 + y3*N3 + y4*N4 + y5*N5 + y6*N6 + y7*N7 );
    T z = 0.125 * ( z0*N0 + z1*N1 + z2*N2 + z3*N3 + z4*N4 + z5*N5 + z6*N6 + z7*N7 );

    m_pVertex->SetXYZ( x, y, z );
}
//-----------------------------------------------------------------------------
//=============================================================================
END_NAMESPACE_TAPs
//34567890123456789012345678901234567890123456789012345678901234567890123456789
//--+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----