sensorset.h

Go to the documentation of this file.

//=============================================================================================================
#ifndef SENSORSET_H
#define SENSORSET_H
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "../inverse_global.h"
 
//=============================================================================================================
// QT INCLUDES
//=============================================================================================================
 
#include <QObject>
#include <QSharedPointer>
 
//=============================================================================================================
// EIGEN INCLUDES
//=============================================================================================================
 
#include <Eigen/Core>
 
//=============================================================================================================
// FORWARD DECLARATIONS
//=============================================================================================================
 
namespace FWDLIB{
    class FwdCoil;
    class FwdCoilSet;
}
 
namespace FIFFLIB{
    class FiffCoordTransOld;
    class FiffDigPointSet;
    class FiffChInfo;
}
 
//=============================================================================================================
// DEFINE NAMESPACE INVERSELIB
//=============================================================================================================
 
namespace INVERSELIB
{
enum class Accuracy : int{high = 2, medium = 1, low = 0};
 
class INVERSESHARED_EXPORT SensorSet {
 
public:
    typedef QSharedPointer<SensorSet> SPtr;            
    typedef QSharedPointer<const SensorSet> ConstSPtr; 
    //=========================================================================================================
    explicit SensorSet() = default;
 
    //=========================================================================================================
    explicit SensorSet(const QSharedPointer<FWDLIB::FwdCoilSet> pFwdCoilSet);
 
    inline int np() const;
    inline int ncoils() const;
    inline Eigen::VectorXd ez(int iSensor) const;
    inline Eigen::MatrixXd ez() const;
 
    inline Eigen::VectorXd r0(int iSensor) const;
    inline Eigen::MatrixXd r0() const;
 
    inline Eigen::MatrixXd rmag(int iSensor) const;
    inline Eigen::MatrixXd rmag() const;
 
    inline Eigen::MatrixXd cosmag(int iSensor) const;
    inline Eigen::MatrixXd cosmag() const;
 
    inline Eigen::MatrixXd tra(int iSensor) const;
    inline Eigen::MatrixXd tra() const;
 
    inline Eigen::RowVectorXd w(int iSensor) const;
    inline Eigen::RowVectorXd w() const;
 
    inline bool operator== (const SensorSet &b) const;
    inline bool operator!= (const SensorSet &b) const;
 
private:
    //=========================================================================================================
    void initFromFwdCoilSet(const QSharedPointer<FWDLIB::FwdCoilSet> pFwdCoilSet);
 
    //=========================================================================================================
    void initMatrices(int iNchan, int iNp);
 
    Eigen::MatrixXd m_ez{Eigen::MatrixXd(0,0)};
    Eigen::MatrixXd m_r0{Eigen::MatrixXd(0,0)};
    Eigen::MatrixXd m_rmag{Eigen::MatrixXd(0,0)};
    Eigen::MatrixXd m_cosmag{Eigen::MatrixXd(0,0)};
    Eigen::MatrixXd m_tra{Eigen::MatrixXd(0,0)};
    Eigen::RowVectorXd m_w{Eigen::RowVectorXd(0)};
    int m_ncoils{0};
    int m_np{0};
};
 
//=============================================================================================================
// INLINE DEFINITIONS
//=============================================================================================================
 
inline int SensorSet::np() const
{
    return m_np;
}
 
inline int SensorSet::ncoils() const
{
    return m_ncoils;
}
 
inline Eigen::VectorXd SensorSet::ez(int iSensor) const
{
    return m_ez.row(iSensor);
}
 
inline Eigen::MatrixXd SensorSet::ez() const
{
    return m_ez;
}
 
inline Eigen::VectorXd SensorSet::r0(int iSensor) const
{
    return m_r0.row(iSensor);
}
 
inline Eigen::MatrixXd SensorSet::r0() const
{
    return m_r0;
}
 
inline Eigen::RowVectorXd SensorSet::w(int iSensor) const
{
    return m_w.segment(iSensor*m_np,m_np);
}
 
inline Eigen::RowVectorXd SensorSet::w() const
{
    return m_w;
}
 
inline Eigen::MatrixXd SensorSet::rmag(int iSensor) const
{
    return m_rmag.block(iSensor*m_np,0,m_np,3);
}
 
inline Eigen::MatrixXd SensorSet::rmag() const
{
    return m_rmag;
}
 
inline Eigen::MatrixXd SensorSet::cosmag(int iSensor) const
{
    return m_cosmag.block(iSensor*m_np,0,m_np,3);
}
 
inline Eigen::MatrixXd SensorSet::cosmag() const
{
    return m_cosmag;
}
 
inline Eigen::MatrixXd SensorSet::tra() const
{
    return m_tra;
}
//=============================================================================================================
 
inline bool SensorSet::operator== (const SensorSet &b) const
{
    return (this->ez() == b.ez() &&
            this->r0() == b.r0() &&
            this->rmag() == b.rmag() &&
            this->cosmag() == b.cosmag() &&
            this->tra() == b.tra() &&
            this->w() == b.w() &&
            this->np() == b.np() &&
            this->ncoils() == b.ncoils());
}
 
//=============================================================================================================
 
inline bool SensorSet::operator!= (const SensorSet &b) const
{
    bool equal = this==&b;
    return !(equal);
}
 
//=============================================================================================================
class INVERSESHARED_EXPORT SensorSetCreator
{
 
public:
    typedef QSharedPointer<SensorSetCreator> SPtr;            
    typedef QSharedPointer<const SensorSetCreator> ConstSPtr; 
    //=========================================================================================================
    explicit SensorSetCreator();
 
    //=========================================================================================================
    SensorSet updateSensorSet(const QList<FIFFLIB::FiffChInfo>& channelList,
                              const Accuracy& accuracy);
 
private:
    QSharedPointer<FWDLIB::FwdCoilSet>  m_pCoilDefinitions{nullptr};    // the coil definitions as template
};
 
//=============================================================================================================
// INLINE DEFINITIONS
//=============================================================================================================
 
} // namespace INVERSELIB
 
#endif // SENSORSET_H