MMSD_StiefelGaussianLaw.h

Go to the documentation of this file.

#ifndef MMSD_StiefelGaussianLaw_H
#define MMSD_StiefelGaussianLaw_H


#include "MMSD_GaussianLaw.h"
#include "MMSD_StiefelFunction.h"
#include "MATH_StiefelOptimizer.h"


DEFINE_SPTR(MMSD_StiefelGaussianLaw);
class MMSD_StiefelGaussianLaw : public virtual MMSD_GaussianLaw {
    SP_OBJECT(MMSD_StiefelGaussianLaw);
    // ATTRIBUTES
public:
  private:


    SP::MATH_StiefelOptimizer mOptimizer;
    SP::MMSD_StiefelFunction mFunction;
   

    
  
    
    
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    MMSD_StiefelGaussianLaw(void);
  
  
  
    // DESTRUCTORS 
  
  
    virtual ~MMSD_StiefelGaussianLaw(void);   
  
  
  
public:
    //----------------
    // NEW constructors
    // ---------------

    inline static SP::MMSD_StiefelGaussianLaw New() {
        return New(6,1.e-16);
    }
    inline static SP::MMSD_StiefelGaussianLaw New(const int& nIterations,const double& minEigenValue) {
        SP::MMSD_StiefelGaussianLaw p(new MMSD_StiefelGaussianLaw(),
                                      MMSD_StiefelGaussianLaw::Delete());
        p->setThis(p);
        p->setMaxIterationsNumber(nIterations);
        p->setMinEigenValue(minEigenValue);
        return p;
    }

    // -----------------
    // SET & GET Methods
    // -----------------


    inline void setMaxIterationsNumber(const int& n) {
        mOptimizer->setMaxIterationsNumber(n);
    };
    inline int getMaxIterationsNumber() const {
        return mOptimizer->getMaxIterationsNumber();
    };


protected:
   
   
public:
    virtual void restore(const MMSD_DoubleFullMatrix& properties) {
        MMSD_GaussianLaw::restore(properties);
    }
    
    virtual void initialize(const int& clusterIndex,
                            const int& samplesNumber,
                            const MMSD_IntegerVector& clustersIndex,
                            const MMSD_DoubleFullMatrix& Y) {
        MMSD_Law::initialize(clusterIndex,samplesNumber,clustersIndex,Y);

        
        
     
    }
    
    // -----------------------------
    // Parameters estimation methods
    // ------------------------------

 
    
protected:
  

    virtual void updateCovarianceDecompositionMatrix(const MMSD_DoubleFullMatrix& Y,
                                                     const MMSD_DoubleVector& Diag,
                                                     const MMSD_DoubleVector& mean,
                                                     const MMSD_DoubleFullMatrix& W,
                                                     const MMSD_DoubleVector& probs,
                                                     MMSD_DoubleFullMatrix& P,
                                                     MMSD_DoubleFullMatrix& YP);
    

  

    // ----------------------
    // String representation
    // ----------------------
public:
    virtual tString toString() const {
        tString ret=MMSD_GaussianLaw::toString();
        return ret;
    }
    
};

#endif