MMSD_GaussianCluster.h

Go to the documentation of this file.

#ifndef MMSD_GaussianCluster_H
#define MMSD_GaussianCluster_H


#include "MMSD_Cluster.h"

#include "MMSD_GaussianFDFunction.h"

#include  "MATH_ZeroEquation.h"
#include  "MATH_NewtonEquation.h"

DEFINE_SPTR(MMSD_GaussianCluster);
DEFINE_SVPTR(MMSD_GaussianCluster);
class MMSD_GaussianCluster : public virtual MMSD_Cluster {
    SP_OBJECT(MMSD_GaussianCluster);
    // ATTRIBUTES
  
private:



    // function for estimating freedom degrees
    SP::MMSD_GaussianFDFunction mFunction;
    // solver for estimating freedom degrees
    SP::MATH_Equation mSolver;
    
  
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    MMSD_GaussianCluster(void);
  
  
  
    // DESTRUCTORS 
  
  
    virtual ~MMSD_GaussianCluster(void);   
  
  
  
public:
    //----------------
    // NEW constructors
    // ---------------

    static inline SP::MMSD_GaussianCluster New() {
        SP::MMSD_GaussianCluster p(new MMSD_GaussianCluster(),MMSD_GaussianCluster::Delete());
        p->setThis(p);
        return p;
    };
    

    // -----------------
    // SET & GET Methods
    // -----------------
    
protected:
    inline MATH_Equation& getSolver() {
        return *mSolver.get();
    }
    inline MMSD_GaussianFDFunction& getFunction() {
        return *mFunction.get();
    }

protected:
    virtual void computeMultivariateDensity(const MMSD_DoubleFullMatrix& YP, //size N x P
                                            const MMSD_DoubleFullMatrix& P, // size P x P
                                            const MMSD_DoubleVector& D, // size P
                                            const MMSD_DoubleVector& mean, // size P
                                            const MMSD_DoubleVector& nu, // size P
                                            double *Pk) const;
    
    virtual void updateWeights(const MMSD_DoubleFullMatrix& YP, // N x P
                               const MMSD_DoubleFullMatrix& P, // P x P
                               const MMSD_DoubleVector& Diag, // P
                               const MMSD_DoubleVector& mean, // P 
                               const MMSD_DoubleVector& nus, // P
                               const MMSD_DoubleVector& probs, // N
                               MMSD_DoubleFullMatrix& W) const;

    virtual void updateFreedomDegrees(const MMSD_DoubleFullMatrix& W,
                                      const MMSD_DoubleVector& Probs,
                                      const double& maxFD,
                                      MMSD_DoubleVector& nu);
    
    // ----------------------
    // String representation
    // ----------------------
public:
    virtual tString toString() const {
        tString ret=MMSD_Cluster::toString();
        return ret;
    }
    
    
};

#endif