FMEwork/Release/Sources/Shared/DTICalculator/DTICalculator.h

Go to the documentation of this file.

// **InsertLicense** code
//----------------------------------------------------------------------------------

//----------------------------------------------------------------------------------

#ifndef __DTICalculator__H_
#define __DTICalculator__H_

#include "IDTICalculator.h"


#define DEFINE_DIM(DIM)                                        \
const MLuint32 DTICalculator<DIM>::ms_Dim         = DIM;            \
const MLuint32 DTICalculator<DIM>::ms_SymmTensDim = SYMM_MAT_SIZE(DIM)

#define DEFINE_EPSILONS(DIM)                                        \
const double DTICalculator<DIM>::m_sEpsilonRel  = 1.0e-00;  \
const double DTICalculator<DIM>::m_sEpsilonAbs  = 0.0e-08

template<MLuint32 DIM>
class DTICalculator : public IDTICalculator
{

public:

    DTICalculator();
    virtual ~DTICalculator();

    virtual void DeleteDTICalculator ();

    DTICalculator(const DTICalculator& calcObj);
    DTICalculator& operator=(const DTICalculator& calcObj);

    virtual MLuint32 getDimension() const;


    virtual Error_Type setGradients(const double **aGradients, const MLuint32 lNrOfGradients);
    virtual MLuint32 getNrOfGradients() const;

    virtual double getTensorError(const double *pLNNorm, const MLuint32 lNrOfGradients);

    virtual const double* calcDiffusionTensor(const double *pLNNorm);

    virtual Error_Type eigenAnaly();
    virtual Error_Type eigenAnaly(double regularizationParam);
    virtual const double* getEigenvalues() const;
    virtual const double* getEigenvectors() const;


private:

    static const MLuint32 ms_SymmTensDim; //(DIM*(DIM+1))/2;

    static const double m_sEpsilonRel;
    static const double m_sEpsilonAbs;

    float *m_pGradientMem;
    float *m_pGradPseudoInvMem;
    float *m_aGradPseudoInv[MLMISCMATH_SYMM_MAT_SIZE(DIM)];

    double m_aEigenvalue[DIM];
    double m_aEigenvector[DIM*DIM];
    double m_aDiffusionTensor[MLMISCMATH_SYMM_MAT_SIZE(DIM)];
    double m_aOriginalDiffusionTensor[MLMISCMATH_SYMM_MAT_SIZE(DIM)]; // remains equal after eigenAnaly() !!!!!!!!

    MLuint32 m_lNrOfGradients;

};




#endif