---

BFTriDiag3_1.hpp

#ifndef __BFTriDiag3_1_H__
#define __BFTriDiag3_1_H__


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

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

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

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

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

~BFTriDiag3_1(){ 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); }

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

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

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

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

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

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

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

inline void setBlockL(int i, FDiagonalMatrix* value){ low[i-1]= value; }
inline FDiagonalMatrix* 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); }

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

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.