---

F_LBiDiag.hpp

#ifndef __F_LBiDiag_H__
#define __F_LBiDiag_H__



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

#include "AFSymMatrix.hpp"

class F_LBiDiag : public AFSymMatrix
{
private:
        
FloatVector**           band;           // band0= Diagonal                      

public:
        
F_LBiDiag(int asize, float val=0.0);
F_LBiDiag(const F_LBiDiag& v);
F_LBiDiag operator=(const F_LBiDiag& v);

virtual ~F_LBiDiag(){delete band[0]; delete band[1]; delete[] band; }

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

void load(F_LBiDiag* m);
void load(FloatVector* d, FloatVector* dnear);

virtual int size()const { return (2*width-1); }

FloatVector* getBand(int i){ return band[i]; }
FloatVector* getLvalues(int i){ return band[1]; }
FloatVector* getDvalues(){ return band[0]; }

virtual void setBand(int i, FloatVector* v){band[i]->load(v);}

// base 0
virtual void set0(int i, int j, float value);
virtual float get0(int i, int j);

inline void setD0(int i, float value){ band[0]->set0(i, value); }
inline float getD0(int i) const{ return band[0]->get0(i); }

virtual void setL0(int i, float value){ band[1]->set0(i, value); }
virtual float getL0(int i)const{ return band[1]->get0(i); }

// base 1
virtual void set(int i, int j, float value){ set0(i-1, j-1, value); }
virtual float get(int i, int j){ return get0(i-1, j-1); }

inline void setD(int i, float value){ band[0]->set(i, value); }
inline float getD(int i) const{ return band[0]->get(i); }

virtual void setL(int i, float value){ band[1]->set(i, value); }
virtual float getL(int i)const{ return band[1]->get(i); }

// computations

virtual void setAll(float value=0.0){ band[0]->setAll(value); band[1]->setAll(value); }

// direct solvers
virtual FloatVector* solve(FloatVector* b, FloatVector* x=0);           // Lx=b

virtual FloatVector* linearSolverL0(FloatVector* b, FloatVector* x=0);          // Lx=b
virtual FloatVector* linearSolverL1(FloatVector* b, FloatVector* x=0);          // L1x=b

// direct inversion L
AFSymMatrix* invert();          // Linv = L^-1

virtual float det(){return 0;}
virtual float trace(){return 0;}

virtual float norme2();
virtual float sum();

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

// scalar single operations
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+=(F_LBiDiag& m);
void operator-=(F_LBiDiag& m);

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

// u= A*v including implicit transposition
FloatVector* mult_Av(FloatVector* v, FloatVector* result=0);
FloatVector* mult_ATv(FloatVector* v, FloatVector* result=0);
// u= v*A
FloatVector* mult_vA(FloatVector* v, FloatVector* result=0)     { return mult_ATv(v, result); }
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); }


virtual void output();
virtual void output(FILE* file);
};

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