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();
 
    //=========================================================================================================
    /*  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