inv_hpi_fit_data.h

Go to the documentation of this file.

//=============================================================================================================
 
#ifndef INV_HPI_FIT_DATA_H
#define INV_HPI_FIT_DATA_H
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "../inv_global.h"
#include "inv_hpi_fit.h"
#include "inv_sensor_set.h"
 
//=============================================================================================================
// EIGEN INCLUDES
//=============================================================================================================
 
#include <Eigen/Core>
 
//=============================================================================================================
// QT INCLUDES
//=============================================================================================================
 
#include <QSharedPointer>
 
//=============================================================================================================
// FORWARD DECLARATIONS
//=============================================================================================================
 
namespace FIFFLIB{
    class FiffInfo;
    class FiffCoordTrans;
    class FiffDigPointSet;
}
 
//=============================================================================================================
// DEFINE NAMESPACE HPILIB
//=============================================================================================================
 
namespace INVLIB
{
 
//=============================================================================================================
// 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 INVSHARED_EXPORT InvHpiFitData
{
 
public:
    typedef QSharedPointer<InvHpiFitData> SPtr;             
    typedef QSharedPointer<const InvHpiFitData> ConstSPtr;  
 
    //=========================================================================================================
    explicit InvHpiFitData();
 
    //=========================================================================================================
    void doDipfitConcurrent();
 
    Eigen::MatrixXd         m_coilPos;
    Eigen::RowVectorXd      m_sensorData;
    DipFitError             m_errorInfo;
    InvSensorSet               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 InvSensorSet& sensors);
 
    //=========================================================================================================
    DipFitError dipfitError(const Eigen::MatrixXd& matPos,
                            const Eigen::MatrixXd& matData,
                            const InvSensorSet& 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 InvSensorSet& sensors,
                               int &iSimplexNumitr);
};
 
//=============================================================================================================
// INLINE DEFINITIONS
//=============================================================================================================
} //NAMESPACE
 
#endif // INV_HPI_FIT_DATA_H