---

I.hpp

#ifndef __I_I__
#define __I_I__



#include <stdio.h>
#include <string.h>
#include "RealVector.hpp"
#include "RealSMatrix.hpp"

class I
{
private:

int                     width;

public:
        
I(int n){width=n;}
I(const I& m){width= m.width;}
I operator=(const I& m){ if(this!=NULL)width= m.width; return *this; }
~I(){}

I* copy(){ return new I(*this); }

I t(){ return *this; }                                  // The Matrix is symmetric
I inverse(){ return *this; }                            // The Matrix is orthogonal

int size(){ return width*width; }

int getWidth()const { return width; }
int getHeight()const { return width; }

// base 0
void set0(int i, int j, double value){}
double get0(int i, int j)const { return (i==j ? 1.0 : 0.0); }

// base 1
void set(int i, int j, double value){}
double get(int i, int j)const { return (i==j ? 1.0 : 0.0); }


// Identity transformations
// ---------------------------
RealVector* Ix(RealVector& v, RealVector* dest=0)               // v= I*x
        {
        if( dest == 0 ) dest= v.copy();
        else *dest= v;
        return dest;
        }

double* Ix(double* v, double* dest=0)
        {
        if( dest == 0 ) dest= new double[width];
        memcpy(dest, v, width*sizeof(double));
        return dest;
        }

RealSMatrix* IA(RealSMatrix& A, RealSMatrix* dest=0)            // M= I*A
        {
        if( dest == 0 ) dest= A.copy();
        else *dest= A;
        return dest;
        }

        
RealSMatrix* AI(RealSMatrix& A, RealSMatrix* dest=0)            // M= A*I
        {
        if( dest == 0 ) dest= A.copy();
        else *dest= A;
        return dest;
        }

void output(){printf("I(%d)\n", width); RealSMatrix ident(*this); ident.output();}
void output(FILE* file){fprintf(file, "I(%d)\n", width); RealSMatrix ident(*this); ident.output(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.