animaGaussianMCMCost.h

Go to the documentation of this file.

#pragma once

#include <vnl/vnl_matrix.h>
#include <vnl/vnl_diag_matrix.h>

#include <animaBaseMCMCost.h>

#include <animaMultiCompartmentModel.h>
#include <AnimaMCMExport.h>

namespace anima
{    

class ANIMAMCM_EXPORT GaussianMCMCost : public anima::BaseMCMCost
{
public:
    typedef GaussianMCMCost Self;
    typedef anima::BaseMCMCost Superclass;
    typedef itk::SmartPointer<Self> Pointer;
    typedef itk::SmartPointer<const Self> ConstPointer;

    itkNewMacro(Self)

    
    itkTypeMacro(GaussianMCMCost, anima::BaseMCMCost)

    
    MeasureType GetValues(const ParametersType &parameters) ITK_OVERRIDE;

    double GetCurrentCostValue() ITK_OVERRIDE;

    void GetDerivativeMatrix(const ParametersType &parameters, DerivativeMatrixType &derivative) ITK_OVERRIDE;

    void GetCurrentDerivative(DerivativeMatrixType &derivativeMatrix, DerivativeType &derivative) ITK_OVERRIDE;

    void SetMarginalEstimation(bool val) {m_MarginalEstimation = val;}

protected:
    GaussianMCMCost()
    {
        m_MarginalEstimation = false;
    }

    virtual ~GaussianMCMCost() {}

private:
    GaussianMCMCost(const Self&); //purposely not implemented
    void operator=(const Self&); //purposely not implemented

    bool m_MarginalEstimation;

    MeasureType m_Residuals;
};

} // end namespace anima