---

FloatSMatrix.hpp

#ifndef __FloatSMatrix_H__
#define __FloatSMatrix_H__



#include <stdio.h>
#include "FloatVector.hpp"

#include "AFSymMatrix.hpp"


class FloatSMatrix : public AFSymMatrix
{
private:

float**         data;

public:
        
FloatSMatrix(int n, float value=0.0);
FloatSMatrix(const FloatSMatrix& m);
FloatSMatrix operator=(const FloatSMatrix& m);

virtual ~FloatSMatrix();

virtual AFSymMatrix* copyAnn();
virtual AFloatMatrix* copy();
virtual AFloatMatrix* t(AFloatMatrix* result=0);

// base 0
virtual void set0(int i, int j, float value){ data[i][j]= value; }
virtual float get0(int i, int j){ return data[i][j]; }

// base 1
virtual void set(int i, int j, float value){ data[i-1][j-1]= value; }
virtual float get(int i, int j){ return data[i-1][j-1]; }

// computations
virtual void setAll(float value=0.0);

virtual float det();
virtual float trace();

virtual float norme2();
virtual float normeF();         // Frobenuis norme= trace(At*A)
virtual float norme1();
virtual float normeInf();

virtual float sum();
virtual float sum2();

virtual float minimum();
virtual float maximum();
virtual float sigma();

// math vectorial op
FloatSMatrix sqrt();
FloatSMatrix pow(float pow);
FloatSMatrix sqr();
FloatSMatrix exp();
FloatSMatrix log();
FloatSMatrix log(float base);

FloatSMatrix sin();
FloatSMatrix cos();
FloatSMatrix tan();     

// 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+=(FloatSMatrix& m);
void operator-=(FloatSMatrix& m);

virtual void add(AFloatMatrix& m);
virtual void subst(AFloatMatrix& m);

// Krylov Spaces
FloatVector** arnoldi(FloatSMatrix* h); // H(n-1,n-1)

// u= A*v including implicit transposition
virtual FloatVector* mult_Av(FloatVector* v, FloatVector* result=0);
virtual FloatVector* mult_ATv(FloatVector* v, FloatVector* result=0);
// 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
FloatVector* mult_Av(FloatVector* v, int col0, int row0, bool incremental=false, FloatVector* result=0);
FloatVector* mult_ATv(FloatVector* v, int col0, int row0, bool incremental=false, FloatVector* result=0);
// u= v*A
FloatVector* mult_vA(FloatVector* v, int row0, int col0, bool incremental=false, FloatVector* result=0)
                        { return mult_ATv(v, col0, row0, incremental, result); }
FloatVector* mult_vAT(FloatVector* vt, int row0, int col0, bool incremental=false, FloatVector* result=0)
                        { return mult_Av(vt, col0, row0, incremental, result); }

// condition number Estimation
float cond1();                          // to verify, may need more accuracy

// direct inversion U, L
AFloatMatrix* inverseOfU0(AFloatMatrix* uinv=0);                // Uinv = U^-1
AFloatMatrix* inverseOfU1(AFloatMatrix* uinv=0);        // U1inv = U1^-1

// incomplete factorisation
virtual AFloatMatrix* ilu0();
virtual AFloatMatrix* milu0();
virtual AFloatMatrix* rilu0(float omega=0.95);
virtual AFloatMatrix* iluth(double threshold=0.001);
virtual AFloatMatrix* ilut(int p, double eps=0.001);            // Threshold is related to magitude in row
virtual AFloatMatrix* icf();
virtual AFloatMatrix* ildlt();

// approximate inversion
virtual AFloatMatrix* apinv(int iterOuter, int iterInner);              // column based iterative inversion

// friends
friend FloatVectorT operator*(FloatVectorT& vt, FloatSMatrix& m);       
friend FloatVector operator*(FloatSMatrix& m, FloatVector& v);  

virtual void output(){ printf("Full "); AFSymMatrix::output(); }
};

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