hpifitdata.h

Go to the documentation of this file.

//=============================================================================================================
#ifndef HPIFITDATA_H
#define HPIFITDATA_H
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "../inverse_global.h"
#include "hpifit.h"
#include "sensorset.h"
 
//=============================================================================================================
// EIGEN INCLUDES
//=============================================================================================================
 
#include <Eigen/Core>
 
//=============================================================================================================
// QT INCLUDES
//=============================================================================================================
 
#include <QSharedPointer>
 
//=============================================================================================================
// FORWARD DECLARATIONS
//=============================================================================================================
 
namespace FIFFLIB{
    class FiffInfo;
    class FiffCoordTrans;
    class FiffDigPointSet;
}
 
//=============================================================================================================
// DEFINE NAMESPACE INVERSELIB
//=============================================================================================================
 
namespace INVERSELIB
{
 
//=============================================================================================================
// Declare all structures to be used
//=============================================================================================================
struct DipFitError {
    double error;
    Eigen::MatrixXd moment;
    int numIterations;
};
 
//=========================================================================================================
struct HPISortStruct {
    double base_arr;
    int idx;
};
 
//=============================================================================================================
// FORWARD DECLARATIONS
//=============================================================================================================
 
//=============================================================================================================
class INVERSESHARED_EXPORT HPIFitData
{
 
public:
    typedef QSharedPointer<HPIFitData> SPtr;             
    typedef QSharedPointer<const HPIFitData> ConstSPtr;  
    //=========================================================================================================
    explicit HPIFitData();
 
    //=========================================================================================================
    void doDipfitConcurrent();
 
    Eigen::MatrixXd         m_coilPos;
    Eigen::RowVectorXd      m_sensorData;
    DipFitError             m_errorInfo;
    SensorSet               m_sensors;
    Eigen::MatrixXd         m_matProjector;
 
    int                     m_iMaxIterations;
    float                   m_fAbortError;
 
protected:
    //=========================================================================================================
    Eigen::MatrixXd magnetic_dipole(Eigen::MatrixXd matPos,
                                    Eigen::MatrixXd matPnt,
                                    Eigen::MatrixXd matOri);
 
    //=========================================================================================================
    Eigen::MatrixXd compute_leadfield(const Eigen::MatrixXd& matPos,
                                      const SensorSet& sensors);
 
    //=========================================================================================================
    DipFitError dipfitError(const Eigen::MatrixXd& matPos,
                            const Eigen::MatrixXd& matData,
                            const SensorSet& sensors,
                            const Eigen::MatrixXd& matProjectors);
 
    //=========================================================================================================
    static bool compare(HPISortStruct a, HPISortStruct b);
 
    //=========================================================================================================
    Eigen::MatrixXd fminsearch(const Eigen::MatrixXd& matPos,
                               int iMaxiter,
                               int iMaxfun,
                               int iDisplay,
                               const Eigen::MatrixXd& matData,
                               const Eigen::MatrixXd& matProjectors,
                               const SensorSet& sensors,
                               int &iSimplexNumitr);
};
 
//=============================================================================================================
// INLINE DEFINITIONS
//=============================================================================================================
} //NAMESPACE
 
#endif // HPIFITDATA_H