inv_dipole.h

Go to the documentation of this file.

//=============================================================================================================
 
#ifndef INV_DIPOLE_H
#define INV_DIPOLE_H
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
//=============================================================================================================
// EIGEN INCLUDES
//=============================================================================================================
 
#include <Eigen/Core>
 
#include <iostream>
 
//=============================================================================================================
// DEFINE NAMESPACE INVLIB
//=============================================================================================================
 
namespace INVLIB
{
 
//=============================================================================================================
// FORWARD DECLARATIONS
//=============================================================================================================
template<typename T>
class InvDipole;
 
//=============================================================================================================
template<typename T>
struct InvDipolePair
{
    int m_iIdx1;            
    InvDipole<T> m_Dipole1; 
 
    int m_iIdx2;            
    InvDipole<T> m_Dipole2; 
 
    T m_vCorrelation; 
};
 
//=============================================================================================================
template<class T>
class InvDipole
{
    //typedef Eigen::Matrix<T, 3, 1> Point3D;
 
public:
    //=========================================================================================================
    InvDipole();
 
    //=========================================================================================================
    virtual ~InvDipole();
 
    inline T& x()
    {
        return m_vecPosition[0];
    }
    inline T& y()
    {
        return m_vecPosition[1];
    }
    inline T& z()
    {
        return m_vecPosition[2];
    }
    inline T x() const
    {
        return m_vecPosition[0];
    }
    inline T y() const
    {
        return m_vecPosition[1];
    }
    inline T z() const
    {
        return m_vecPosition[2];
    }
 
    inline T& phi_x()
    {
        return m_vecDirection[0];
    }
    inline T& phi_y()
    {
        return m_vecDirection[1];
    }
    inline T& phi_z()
    {
        return m_vecDirection[2];
    }
    inline T phi_x() const
    {
        return m_vecDirection[0];
    }
    inline T phi_y() const
    {
        return m_vecDirection[1];
    }
    inline T phi_z() const
    {
        return m_vecDirection[2];
    }
 
    //=========================================================================================================
    void clean();
 
protected:
private:
    Eigen::Matrix<T, 3, 1> m_vecPosition;
    Eigen::Matrix<T, 3, 1> m_vecDirection;
 
    double m_dLength;
    double m_dFrequency;
 
    //TGreensFunction* green;
};
} // NAMESPACE
 
//TypeDefs
typedef INVLIB::InvDipole<int> InvDipole_INT;
typedef INVLIB::InvDipole<float> InvDipole_FLOAT;
typedef INVLIB::InvDipole<double> InvDipole_DOUBLE;
 
//Make the template definition visible to compiler in the first point of instantiation
#include "inv_dipole.cpp"
 
#endif // INV_DIPOLE_H