---

AbsTriDiagonalLS.hpp

#ifndef __AbsTriDiagonalLS_H__
#define __AbsTriDiagonalLS_H__


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

class AbsTriDiagonalLS
{
protected:
        
int                                     N;

RTriDiagonalMatrix*     A;
RealVector*                     x;                                      
 
RealVector*                     b;                                       

public:
        
AbsTriDiagonalLS(int asize)
        {
        N=asize;
        
        x= new RealVector(N);
        b= new RealVector(N);

        A= new RTriDiagonalMatrix(N);
        }

virtual ~AbsTriDiagonalLS(){ delete A; delete x; delete b;}

virtual void initFreezing(int val=0){}
virtual void freeze(int item, int val=0){}  // freezing used later for marking identity diagonal only values 

virtual void loadA(RTriDiagonalMatrix* m){ A->load(m); }
virtual void loadMatrix(RealVector* u1, RealVector* d1, RealVector* l1){ A->load(u1,d1,l1); }
virtual void loadX(RealVector* data){ x->load(data); }
virtual void loadB(RealVector* data){ b->load(data); }

virtual int size(){ return (3*N-2); }

virtual int getWidth(){ return N; }
virtual int getHeight(){ return N; }

virtual RealVector* getX(){ return x; }
virtual RealVector* getB(){ return b; }

virtual RTriDiagonalMatrix* getA(){ return A; }

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

virtual bool solve(RealVector* aX, RealVector* aB)      // A not changed
        {
        loadX(aX);
        loadB(aB);

        return solve();
        }

virtual bool solve(RTriDiagonalMatrix*  aA, RealVector* aX, RealVector* aB)     // new system
        {
        loadA(aA);

        loadX(aX);
        loadB(aB);

        return solve();
        }
        

virtual RealVector* evaluate(RealVector* aX, RealVector* aB)
        {
        solve(aX, aB);

        return getX();
        }

virtual RealVector* evaluate(RTriDiagonalMatrix* aA, RealVector* aX, RealVector* aB)
        {
        solve(aA, aX, aB);

        return getX();
        }
        
virtual void output()
        {
        printf("Linear System\n");
        A->output();
        printf("x: "); x->output();
        printf("b: "); b->output();
        }

virtual void output(FILE* file)
        {
        fprintf(file, "Linear System\n");
        A->output(file);
        fprintf(file, "x: "); x->output(file);
        fprintf(file, "b: "); b->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.