---

H.hpp

#ifndef __H_H__
#define __H_H__



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

class RealSMatrix;

class H
{
private:

int                     width;
RealVector              u;

double          sigma;

public:
        
H(int n, double value=0.0):u(0.0){width=n; sigma=0.0;}
H(int n, double* x, bool xUpdate=false);
H(RealVector& x, bool xUpdate=false);
H(const H& m);
H operator=(const H& m);
~H(){}

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

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

double computeHx(RealVector& x, bool xUpdate=false);            // computed u,sigma via x
double computeHx(double* x);

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

double getBeta(){ return 2.0/u.vT_v(); }                // beta= 2/(uT*u)

// base 0
void setU0(int i, double value){ u.set0(i, value); }
double getU0(int i)const { return u.get0(i); }

// base 1
void setU(int i, double value){ u.set0(i, value); }
double getU(int i)const { return u.get0(i); }

double getSigma(){return sigma;}
RealVector& getU(){return u;}

// Householder transformations
// ---------------------------
RealVector* Hx(RealVector& v, RealVector* dest=0);              // v= H*x
double* Hx(double* v, double* dest=0);

RealSMatrix* HA(RealSMatrix& A, RealSMatrix* dest=0);           // M= H*A
RealSMatrix* AH(RealSMatrix& A, RealSMatrix* dest=0);           // M= A*H


// friends
friend RealVector operator*(H& h, RealVector& v);       // v= H*x 
friend RealSMatrix operator*(RealSMatrix& m, H& h);     // M= H*A
friend RealSMatrix operator*(H& h, RealSMatrix& m);     // M= A*H
        
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.