---

BSp1FloatSMatrix.hpp

#ifndef __BSp1FloatSMatrix_H__
#define __BSp1FloatSMatrix_H__


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

#include "AFSymMatrix.hpp"

class BSp1FloatSMatrix 
{
protected:

int                     width;                  // n blocks

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

int                     nItems, maxItems;

int*                    ia;
int*                    ja;
AFSymMatrix**   data;

public:
        
BSp1FloatSMatrix(int w, int nmax);
BSp1FloatSMatrix(const BSp1FloatSMatrix& m);
BSp1FloatSMatrix operator=(const BSp1FloatSMatrix& m);

virtual ~BSp1FloatSMatrix();

virtual BSp1FloatSMatrix* copy();

// 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 width; }
int getBlockHeight(){ return width; }
int size(){ return width*width; }


virtual int numberOfItems(){return nItems;}
virtual int getMaxItems(){return maxItems;}

virtual void addItem(int x,int y, AFSymMatrix* value=0);

virtual bool getItem(int index, int* x, int* y, AFSymMatrix** value);
virtual bool setItem(int index, AFSymMatrix* value);

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

// base 1
void setBlock(int i, int j, AFSymMatrix* value){ setBlock0(i-1,j-1, value); }
AFSymMatrix* getBlock(int i, int j){ return getBlock0(i-1,j-1); }

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

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

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

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

virtual float norme2();
virtual float sum();
virtual float sigma(){return 0.0;}

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

// 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);


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