---

FIdentity.hpp

#ifndef __FIdentity_H__
#define __FIdentity_H__


#include <stdio.h>
#include <math.h>

#include "FloatVector.hpp"

#include "AFSymMatrix.hpp"


class FIdentity : public AFSymMatrix
{
public:
        
FIdentity(int n):AFSymMatrix(n){}
FIdentity(const FIdentity& m):AFSymMatrix(m){}
FIdentity operator=(const FIdentity& m){ if(this!=NULL)width= m.width; return *this; }

virtual ~FIdentity(){}

virtual AFloatMatrix* copy(){ return new FIdentity(*this); }
virtual AFloatMatrix* t(AFloatMatrix* result=0);

// base 0
virtual void set0(int i, int j, float value){}
virtual float get0(int i, int j){ return (i==j ? 1.0 : 0.0); }

// base 1
virtual void set(int i, int j, float value){}
virtual float get(int i, int j){ return (i==j ? 1.0 : 0.0); }

// computations
virtual void setAll(float value=0.0){}

virtual float det(){return 1.0;}
virtual float trace(){return getWidth();}

virtual float norme2(){return getWidth();}
virtual float sum(){return getWidth();}
virtual float sum2(){return getWidth();}

virtual float minimum(){return 0.0;}
virtual float maximum(){return 1.0;}
virtual float sigma();


// scalar single operations
void operator+=(float value){}
void operator-=(float value){}
void operator*=(float value){}
void operator/=(float value){}

virtual void add(float value){}
virtual void subst(float value){}
virtual void mult(float value){}
virtual void div(float value){}

// vector single operations
void operator+=(FIdentity& m){}
void operator-=(FIdentity& m){}

virtual void add(AFloatMatrix& m){}
virtual void subst(AFloatMatrix& m){}


// u= A*v including implicit transposition
virtual FloatVector* mult_Av(FloatVector* v, FloatVector* result=0);
virtual FloatVector* mult_ATv(FloatVector* v, FloatVector* result=0){ return mult_Av(v, result); }
// u= v*A
virtual FloatVector* mult_vA(FloatVector* v, FloatVector* result=0){ return mult_ATv(v, result); }
virtual FloatVector* mult_vAT(FloatVector* vt, FloatVector* result=0){ return mult_Av(vt, result); }

// sub vector usage
// u= A*v
virtual FloatVector* mult_Av(FloatVector* v, int col0, int row0, bool incremental=false, FloatVector* result=0);
virtual FloatVector* mult_ATv(FloatVector* v, int col0, int row0, bool incremental=false, FloatVector* result=0)
                { return mult_Av(v, col0, row0, incremental, result); }
        
// u= v*A
virtual FloatVector* mult_vA(FloatVector* v, int row0, int col0, bool incremental=false, FloatVector* result=0)
                { return mult_ATv(v, col0, row0, incremental, result); }
virtual FloatVector* mult_vAT(FloatVector* vt, int row0, int col0, bool incremental=false, FloatVector* result=0)
                { return mult_Av(vt, col0, row0, incremental, result); }

virtual void output(){ printf("I(%d,%d)\n", getWidth(), getHeight()); } 
virtual void output(FILE* file){ fprintf(file, "I(%d,%d)\n", getWidth(), getHeight()); } 
};

#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.