tpzeqnarray.h

Go to the documentation of this file.

#ifndef TPZEQNARRAY_H
#define TPZEQNARRAY_H

#include "pzmatrix.h"
#include "pzstack.h"
#include "pzfmatrix.h"
#include "pzreal.h"
#include "pzmanvector.h"
#include "pzvec.h"

#include <stdlib.h>

#ifdef BORLAND
#include <stdio.h>
#endif

template<class TVar>
class TPZEqnArray {     
                
                enum MSymmetric {  EIsUndefined, EIsSymmetric, EIsNonSymmetric };
                
public:
                static void main();
                
    void SetNonSymmetric();
                
    void SetSymmetric();
                
    int IsSymmetric();
                
    TPZEqnArray();
    ~TPZEqnArray();
    void EqnForward(TPZFMatrix<TVar> & F, DecomposeType dec);
    void EqnBackward(TPZFMatrix<TVar> & U, DecomposeType dec);
                
    void Reset();
    void BeginEquation(int eq);
    void AddTerm(int col, TVar val)
                {
                                fIndex.Push(col);
                                fEqValues.Push(val);
                                fLastTerm++;
                } 
    void EndEquation();
                
    void Read(char * inputfile);
                
    void Write(char * outputfile);
                
    void Print(const char * name, std::ostream & out);
                
    void Write(FILE * outputfile);
                
    void Read(FILE * inputfile);
                
private:
                
                MSymmetric fSymmetric;
                
                
    int fNumEq;
                
    TPZStack < int , 100 > fEqStart;
    TPZStack < int , 100 > fEqNumber;
                
    TPZStack < TVar, 1000 > fEqValues;
                
    TPZStack < int , 1000 > fIndex;
                
    int fLastTerm;
};

#endif //TPZEQNARRAY_H