---

BFTriDiagonalMatrix.hpp

#ifndef __BFTriDiagonalMatrix_H__
#define __BFTriDiagonalMatrix_H__


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

class BFTriDiagonalMatrix
{
private:
        
int                             N;                              // number of matrices on the diagonal

int                             totalWidth;                     // n values
int                             widthSubMatrix;         // n value in sub matrix

AFSymMatrix**           up;                                     // N-1 size upper diagonal 
AFSymMatrix**           d;                                      // N size diagonal
AFSymMatrix**           low;                                    // N-1 size lower diagonal 

public:
        
BFTriDiagonalMatrix(int asize, int nsub);                               // empty
BFTriDiagonalMatrix(int asize, int nsub, float value);  
BFTriDiagonalMatrix(int asize, int nsub, float lval, float dval, float uval);                   
BFTriDiagonalMatrix(const BFTriDiagonalMatrix& v);
BFTriDiagonalMatrix operator=(const BFTriDiagonalMatrix& v);

~BFTriDiagonalMatrix(){ delete d; delete up; delete low;}

// total values
int getWidth(){ return totalWidth; }
int getHeight(){ return totalWidth; }

// item size
int getItemWidth(){ return widthSubMatrix; }
int getItemHeight(){ return widthSubMatrix; }

// number of blocks
int getBlockWidth(){ return N; }
int getBlockHeight(){ return N; }
inline int size()const { return (3*N-2); }

AFSymMatrix** getUblocks(){ return up; }
AFSymMatrix** getDblocks(){ return d; }
AFSymMatrix** getLblocks(){ return low; }

// base 0
void setBlock0(int i, int j, AFSymMatrix* value);
AFSymMatrix* getBlock0(int i, int j)const ;

inline void setBlockD0(int i, AFSymMatrix* value){ d[i]= value; }
inline AFSymMatrix* getBlockD0(int i)const { return d[i]; }

inline void setBlockU0(int i, AFSymMatrix* value){ up[i]= value; }
inline AFSymMatrix* getBlockU0(int i)const { return up[i]; }

inline void setBlockL0(int i, AFSymMatrix* value){ low[i]= value; }
inline AFSymMatrix* getBlockL0(int i)const { return low[i]; }

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

void setBlock(int i, int j, AFSymMatrix* value);
AFSymMatrix* getBlock(int i, int j)const ;

inline void setBlockD(int i, AFSymMatrix* value){ d[i-1]= value; }
inline AFSymMatrix* getBlockD(int i)const { return d[i-1]; }

inline void setBlockU(int i, AFSymMatrix* value){ up[i-1]= value; }
inline AFSymMatrix* getBlockU(int i)const { return up[i-1]; }

inline void setBlockL(int i, AFSymMatrix* value){ low[i-1]= value; }
inline AFSymMatrix* getBlockL(int i)const { return low[i-1]; }

// base 1
void set(int i, int j, float value);
float get(int i, int j)const ;

// computations
float det();
float trace();

float sum();

float minimum();
float maximum();
float mean(){ return sum()/size(); }

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

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

void output();
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.