MMSD_GaussianModel.h

Go to the documentation of this file.

#ifndef MMSD_GaussianModel_H
#define MMSD_GaussianModel_H


#include "MMSD_Model.h"
#include "MMSD_FluryGautschiGaussianLaw.h"
#include "MMSD_StiefelGaussianLaw.h"
#include "MMSD_GaussianCluster.h"

DEFINE_SPTR(MMSD_GaussianModel);
class MMSD_GaussianModel : public virtual MMSD_Model {
    SP_OBJECT(MMSD_GaussianModel);
    // ATTRIBUTES
    
public:
    static const tFlag FLURY_GAUTSCHI_LAW;
    static const tFlag STIEFEL_LAW;
    
private:
  
    tFlag mGaussianLaw;
    int mGaussianLawMaximumIterationsNumber;
    
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    MMSD_GaussianModel(void);
  
  
  
    // DESTRUCTORS 
  
  
    virtual ~MMSD_GaussianModel(void);   
  
  
  
public:
    //----------------
    // NEW constructors
    // ---------------

    inline static SP::MMSD_GaussianModel New() {
        return New(FLURY_GAUTSCHI_LAW,6);
    }
    inline static SP::MMSD_GaussianModel New(const tFlag& gaussianLaw,const int& nIterations) {
        SP::MMSD_GaussianModel p(new MMSD_GaussianModel(),
                                 MMSD_GaussianModel::Delete());
        p->setThis(p);
        p->setGaussianLawMaximumIterationsNumber(nIterations);
        p->setGaussianLaw(gaussianLaw);
        return p;
    }

    // SET & GET Methods
    // ==================
    
    inline void setGaussianLaw(const tFlag& f) {
        mGaussianLaw=f;
    }
    inline int getGaussianLaw() const {
        return mGaussianLaw;
    }
    inline void setGaussianLawMaximumIterationsNumber(const int& n) {
        mGaussianLawMaximumIterationsNumber=n;
    };
    inline int getGaussianLawMaximumIterationsNumber() const {
        return mGaussianLawMaximumIterationsNumber;
    };

    // SPECIALIZATION Methods
    // ========================
    
    virtual SP::MMSD_Law NewLawInstance() const {
        if (getGaussianLaw()==STIEFEL_LAW)
            return MMSD_StiefelGaussianLaw::New(getGaussianLawMaximumIterationsNumber(),getMinEigenValue());
        else 
            return MMSD_FluryGautschiGaussianLaw::New(getGaussianLawMaximumIterationsNumber(),getMinEigenValue());
    };
    virtual SP::MMSD_Cluster NewClusterInstance() const {
        return MMSD_GaussianCluster::New(); 
    };
   
};

#endif