pzanalysiserror.h

Go to the documentation of this file.

#ifndef TPZANALYSISERRORH
#define TPZANALYSISERRORH

#include <fstream>        // for ostream, operator<<
#include "pzanalysis.h"   // for TPZAnalysis
#include "pzcompel.h"     // for TPZCompElSide, TPZCompEl (ptr only)
#include "pzmanvector.h"  // for TPZManVector
#include "pzreal.h"       // for REAL
#include "pzstack.h"      // for TPZStack
class TPZCompMesh;

class TPZAnalysisError : public TPZAnalysis {
                
                TPZManVector<int64_t> fElIndexes;
                TPZManVector<REAL> fElErrors;
                TPZStack<TPZCompElSide> fSingular;
                REAL fTotalError;
                REAL fAdmissibleError;
                
                REAL fEtaAdmissible;
                int64_t fNIterations;
                
                public :
                TPZAnalysisError(TPZCompMesh *mesh,std::ostream &out);
                ~TPZAnalysisError() {};
                
                void SetAdaptivityParameters(REAL EtaAdmissible, int64_t NIterations);
                
                void hp_Adaptive_Mesh_Design(std::ostream &out,REAL &EtaAdmissible);
                
                void GetSingularElements(TPZStack<TPZCompElSide> &listel);
                
                void ZoomInSingularity(REAL csi, TPZCompElSide elside, REAL singularity_order = 0.9);
                
                void HPAdapt(REAL CurrentEtaAdmissible, std::ostream &out);
                
                REAL MaximLocalError();
                
                REAL h_Parameter(TPZCompEl *cel);
                
                void MathematicaPlot();
                
                void EvaluateError(REAL CurrentEtaAdmissible, bool store_error, std::ostream &out);
                
private:
                void PlotLocal(int64_t iter, REAL CurrentEtaAdmissible, std::ostream &out);
                
                void ExpandConnected(TPZStack<TPZCompElSide> &singel);          
};

#endif