---

BFloatSMatrix.hpp

#ifndef __BFloatSMatrix_H__
#define __BFloatSMatrix_H__



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

class BFloatSMatrix
{
private:

int                             width;                  // n blocks

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

AFSymMatrix**           p1;
AFSymMatrix***          data;                           // sub matrices

public:
        
BFloatSMatrix(int n, int nsub);                 // Empty block Matrix    ///
BFloatSMatrix(int n, int nsub, float value);            
BFloatSMatrix(int n, AFSymMatrix& m);
BFloatSMatrix(const BFloatSMatrix& bm);                 
BFloatSMatrix operator=(const BFloatSMatrix& bm);       

~BFloatSMatrix();                                                       

BFloatSMatrix* 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; }

// On items
// base 0
void setBlock0(int i, int j, AFSymMatrix* value){ data[i][j]= value; }
AFSymMatrix* getBlock0(int i, int j){ return data[i][j]; }

// base 1
void setBlock(int i, int j, AFSymMatrix* value){ data[i-1][j-1]= value; }
AFSymMatrix* getBlock(int i, int j){ return data[i-1][j-1]; }

// On values
// base 0
void set0(int i, int j, float value)
        { 
        int ix= i/widthSubMatrix; 
        int iy= j/widthSubMatrix;
        data[ix][iy]->set0(i%widthSubMatrix, j%widthSubMatrix, value); 
        }
float get0(int i, int j)
        { 
        int ix= i/widthSubMatrix; 
        int iy= j/widthSubMatrix;
        return data[ix][iy]->get0(i%widthSubMatrix, j%widthSubMatrix); }

// base 1
void set(int i, int j, float value)     
        { 
        int ix= (i-1)/widthSubMatrix; 
        int iy= (j-1)/widthSubMatrix;
        if( data[ix][iy] != 0 )
                data[ix][iy]->set0((i-1)%widthSubMatrix, (j-1)%widthSubMatrix, value); 
        }
float get(int i, int j) 
        { 
        int ix= (i-1)/widthSubMatrix; 
        int iy= (j-1)/widthSubMatrix;
        if( data[ix][iy] == 0 ) return 0.0;
        return data[ix][iy]->get0((i-1)%widthSubMatrix, (j-1)%widthSubMatrix); 
        }

// computations
void setAll(float value=0.0);
void setAll(AFSymMatrix& m);

float det();
float trace();

float norme2();
float sum();                                    

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

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