MMSD_GaussianLaw.h

Go to the documentation of this file.

#ifndef MMSD_GaussianLaw_H
#define MMSD_GaussianLaw_H


#include "MMSD_Law.h"
#include "MMSD_IntegerVector.h"
#include "MMSD_DoubleVector.h"
#include "MMSD_DoubleSymmetricMatrix.h"
#include "MMSD_DoubleFullMatrix.h"
#include "MMSD_IntegerFullMatrix.h"
#include "MMSD_DataSet.h"


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


  
    // work vector
    MMSD_DoubleVector mWork;
    

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

   

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



    inline MMSD_DoubleVector& getWork() {
        return mWork;
    }

protected:
   
   
public:
  
    
    // -----------------------------
    // Parameters estimation methods
    // ------------------------------

   
protected:
    
    virtual void updateMean(const MMSD_DoubleFullMatrix& P,
                            const MMSD_DoubleFullMatrix& YP,
                            const MMSD_DoubleFullMatrix& W,
                            const MMSD_DoubleVector& probs,
                            MMSD_DoubleVector& mean) ;

      

    virtual void updateCovarianceDecompositionDiagonal( const MMSD_DoubleFullMatrix& P, // PxP
                                                        const MMSD_DoubleFullMatrix& YP, //NxP
                                                        const MMSD_DoubleVector& mean, // P
                                                        const MMSD_DoubleFullMatrix& W, // NxP
                                                        const MMSD_DoubleVector& probs, // N
                                                        MMSD_DoubleVector& D){
        updateCovarianceDecompositionDiagonal(P,YP,mean,W,probs,getMinEigenValue(),mWork,D);
    }

private:
    
    static void updateCovarianceDecompositionDiagonal(const MMSD_DoubleFullMatrix& P, // PxP
                                                      const MMSD_DoubleFullMatrix& YP, //NxP
                                                      const MMSD_DoubleVector& mean, // P
                                                      const MMSD_DoubleFullMatrix& W, // NxP
                                                      const MMSD_DoubleVector& probs, // N
                                                      const double& minValue,
                                                      MMSD_DoubleVector& work,
                                                      MMSD_DoubleVector& D);
 
    // ----------------------
    // String representation
    // ----------------------
public:
    virtual tString toString() const {
        tString ret=MMSD_Law::toString();
        return ret;
    }
    
};

#endif