---

G.hpp

#ifndef __G_H__
#define __G_H__



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

class RealSMatrix;

class G
{
public:

double          c;
double          s;

int                     i;
int                     j;

public:
        
G(double x1, double x2);
G(RealSMatrix& m, int i, int j);
G(const G& m);
G operator=(const G& m);
~G(){}

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

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

double computeGx(RealVector& x, int i, int j, bool xUpdate=false);              // computed u,sigma via x
double computeGx(double* x, int i, int j);
double computeGx(double x1, double x2);

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

double getCos(){ return c; }
double getSin(){ return s; }
double getTheta(){ return atan(s/c); }                          

// Givens transformations
// ---------------------------
RealVector* Gx(RealVector& v, int i, int j, RealVector* dest=0);                // v= G*x
double* Gx(double* v, int i, int j, double* dest=0);

RealSMatrix* GA(RealSMatrix& A, int i, int j, RealSMatrix* dest=0);             // M= Gt*A
RealSMatrix* AG(RealSMatrix& A, int i, int j, RealSMatrix* dest=0);             // M= A*G


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