---

SingleUnitFastIca.hpp

#ifndef __SingleUnitFastIca__
#define __SingleUnitFastIca__

#include "AContrastFunction.hpp"
#include "RealVector.hpp"
#include "RealSMatrix.hpp"

class SingleUnitFastIca 
{
protected:

AContrastFunction*      contrast;

RealVector*                     wold;
RealVector*                     w;

int                             count;
RealVector*                     evg;                    // cummulate estimation needed in FastICA
RealVector*                     evg_Cinv;               // used in case unsphered data
double                  evgdiff;

double                  epsilon;

public:

SingleUnitFastIca(double eps=1.0e-4);
SingleUnitFastIca(AContrastFunction* func, double eps)
                {contrast=func; w=0; wold=0; evg=0; evgdiff=0.0; count=0; epsilon=eps; }
virtual ~SingleUnitFastIca()
 {
 if( w != 0 )
        { 
        delete w;
        delete wold;

        delete evg;
        }
 }

virtual void set(AContrastFunction*     func){contrast= func;}

virtual RealVector* getW(){return w;}

// learning
virtual void init(int size);

virtual void startLearning()
        { count=0; evg->setAll(0.0); evg_Cinv->setAll(0.0); evgdiff=0.0; } // to use to start a new iteration
virtual void learn(RealVector* sample);
virtual bool endLearning();
virtual bool endLearning(RealSMatrix* C, RealSMatrix* Cinv);

virtual void normalizeW();                      // w/||w||
virtual void normalizeW(RealSMatrix* C);  // w/sqrt(wT*C*w)

// out
virtual double forward(RealVector* src);                // maximise absolute value of kurtosis

virtual void output()
        {
        printf("Single Unit FastICA");
        }                               
};

#endif

Restoreinpaint sourceforge project `C++/Java Image Processing, Restoration, Inpainting Project'.

Bernard De Cuyper: Open Project Leader: Concept, design and development. Bernard De Cuyper & Eddy Fraiha 2002, 2003. Bernard De Cuyper 2004. Open and free, for friendly usage only. Modifications on Belgium ground of this piece of artistic work, by governement institutions or companies, must be notified to Bernard De Cuyper.