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