pzskylmatpar.h

Go to the documentation of this file.

#ifndef TSKYLPARMATH
#define TSKYLPARMATH

#include "pzmatrix.h"
#include "pzvec.h"
#include "pzmanvector.h"
#include "pzskylmat.h"


#include <stdlib.h>
#include <signal.h>
#include <time.h>

//parallel libraries
#include <pthread.h>

#include <set>

template<class TVar>
class TPZSkylParMatrix : public TPZSkylMatrix<TVar>
{
public:
                static int main();
                static int main_nada();
                
                TPZSkylParMatrix();
                TPZSkylParMatrix(const int64_t dim);
                TPZSkylParMatrix(const int64_t dim, const TPZVec<int64_t> &skyline,int NumThreads);
                TPZSkylParMatrix(const TPZSkylParMatrix<TVar> &A);
                
                CLONEDEF(TPZSkylParMatrix)
    //    : TPZMatrix(A.Dim(), A.Dim()), fElem(0), fStorage(0) {Copy(A); }
                
                
                virtual ~TPZSkylParMatrix();
                
                int Decompose_Cholesky(std::list<int64_t> &singular) override;
                int Decompose_Cholesky() override;
                
                int Decompose_LDLt(std::list<int64_t> &singular) override;
                int Decompose_LDLt() override;
                
                void SetSkyline(const TPZVec<int64_t> &skyline);

    public:
int ClassId() const override;

private:
                
                static void *ParallelLDLt(void *t);
                static void *ParallelLDLt2(void *t);
                static void *ParallelCholesky(void *t);
                
                void ColumnToWork(int64_t &lcol, int64_t &lprevcol);
                void ColumnToWork(int64_t &lcol);
                void DecomposeColumnCholesky(int64_t lcol, int64_t lprevcol);
                void DecomposeColumnLDLt(int64_t lcol, int64_t lprevcol);
                void DecomposeColumnLDLt2(int64_t lcol);
                void PrintState();
                
public:
                TPZVec<int> fDec;  
                TPZVec<int64_t> fSkyline;
                
private:
                int64_t fEqDec, fNthreads;
                int * fThreadUsed;
                std::set<int64_t> fColUsed;
                int fNumDecomposed;
                bool fCorrectSingular;
                std::list<int64_t> fSingular;
                
#ifdef OOPARLIB
                
                virtual int Unpack( TReceiveStorage *buf );
                static TSaveable *CreateInstance(TReceiveStorage *buf);
                virtual int Pack( TSendStorage *buf ) const;
                virtual std::string ClassName() const   { return( "TPZSkylParMatrix"); }
                virtual int DerivedFrom(const int64_t Classid) const;
                virtual int DerivedFrom(const char *classname) const; // a class with name classname
                
#endif
                
};

#endif