---

SelectiveThomasLS.hpp

#ifndef __SelectiveThomasLS_H__
#define __SelectiveThomasLS_H__



#include "AbsTriDiagonalLS.hpp"


class SelectiveThomasLS : public AbsTriDiagonalLS
{
private:
                                                
RealVector*                     uOut;   
RealVector*                     dOut;
RealVector*                     lOut;   

RealVector*                     y;      

int*                            frozenPixel;    
double                  dfrozen;
double                  dfrozenInv;


protected:
        
void computeLRdecomposition();

void forwardSubstitution();
void backwardSubstitution();

public:
        
SelectiveThomasLS(int asize, double cstFrozen=1.0):AbsTriDiagonalLS(asize)
        {
        uOut= new RealVector(N-1);
        dOut= new RealVector(N);
        lOut= new RealVector(N-1);

        y= new RealVector(N);

        dfrozen= cstFrozen;
        dfrozenInv= 1.0/dfrozen;

        frozenPixel= new int[N];
        memset(frozenPixel, 0, N*sizeof(int));
        }

virtual void initFreezing(int val=0)
        {
        memset(frozenPixel, val, N*sizeof(int));
        }

virtual ~SelectiveThomasLS(){delete uOut; delete dOut; delete lOut; delete y; delete[] frozenPixel;}

virtual void freeze(int item, int val=0)  // freezing used for marking identity, diagonal only values
        {
        frozenPixel[item]= val;
        }

virtual void loadA(RTriDiagonalMatrix* m);
virtual void loadMatrix(RealVector* u1, RealVector* d1, RealVector* l1);
virtual void loadX(RealVector* data);
virtual void loadB(RealVector* data);


virtual bool solve();                                                           // SOLVING to define later

        
virtual void output()
        {
        printf("Thomas ");
        AbsTriDiagonalLS::output();
        }

virtual void output(FILE* file)
        {
        fprintf(file, "Thomas ");
        AbsTriDiagonalLS::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.