---

BFDiagonalMatrix.hpp

#ifndef __BFDiagonalMatrix_H__
#define __BFDiagonalMatrix_H__



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

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

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

AFSymMatrix**           d;                                      // N size diagonal

public:
        
BFDiagonalMatrix(int asize, int nsub);                          // empty
BFDiagonalMatrix(int asize, int nsub, float value);                             
BFDiagonalMatrix(const BFDiagonalMatrix& v);
BFDiagonalMatrix operator=(const BFDiagonalMatrix& v);

~BFDiagonalMatrix(){ delete d; }

// 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 N; }

AFSymMatrix** getDblocks(){ return d; }

// base 0
void setBlock0(int i, AFSymMatrix* value){ d[i]= value; }
AFSymMatrix* getBlock0(int i )const { return d[i]; }
// base 1
void setBlock(int i, AFSymMatrix* value){ d[i-1]= value; }
AFSymMatrix* getBlock(int i )const { return d[i-1]; }

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

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