pzreferredcompel.h

Go to the documentation of this file.

#ifndef PZSPECIAL
#define PZSPECIAL

struct TPZElementMatrix;
class TPZCompEl;
#include "pzintel.h"
class TPZGeoEl;
class TPZCompMesh;
#include "pzfmatrix.h"
#include "pzvec.h"
#include "pzmanvector.h"
#include "pzmaterialdata.h"
#include <map>
#include <set>

template<class TCOMPEL>
class TPZReferredCompEl : public TCOMPEL {
public:
                
                TPZReferredCompEl(TPZCompMesh &mesh, TPZGeoEl *gel, int64_t &index);
                
                TPZReferredCompEl();
                
                TPZReferredCompEl(TPZCompMesh &mesh, const TPZReferredCompEl<TCOMPEL> &copy);
                
                TPZReferredCompEl(TPZCompMesh &mesh,
                      const TPZReferredCompEl<TCOMPEL> &copy,
                      std::map<int64_t,int64_t> & gl2lcConMap,
                      std::map<int64_t,int64_t> & gl2lcElMap);
                
                ~TPZReferredCompEl();
                
                virtual void SetCreateFunctions(TPZCompMesh *mesh) override;
                
                TPZCompEl * ReferredElement();
                
//    virtual void ComputeSolution(TPZVec<REAL> &qsi, TPZFMatrix<REAL> &phi, TPZFMatrix<REAL> &dphix,
//                                 const TPZFMatrix<REAL> &axes, TPZSolVec &sol, TPZGradSolVec &dsol);
                
    virtual void ComputeSolution(TPZVec<REAL> &qsi, TPZMaterialData &data) override;
    
                virtual void ComputeSolution(TPZVec<REAL> &qsi,
                                                                                                                                 TPZVec<REAL> &normal,
                                                                                                                                 TPZSolVec &leftsol, TPZGradSolVec &dleftsol,TPZFMatrix<REAL> &leftaxes,
                                                                                                                                 TPZSolVec &rightsol, TPZGradSolVec &drightsol,TPZFMatrix<REAL> &rightaxes) override;
                
                virtual void Print(std::ostream & out = std::cout) const override;
                
    public:
int ClassId() const override;
    public:
protected:
                
                void AppendOtherSolution(TPZVec<REAL> &qsi, TPZSolVec &sol,
                                                                                                                 TPZGradSolVec &dsol,  TPZFMatrix<REAL> &axes);
                
    void AppendOtherSolution(TPZVec<REAL> &qsi, TPZSolVec &sol);
    
                void AppendOtherSolution(TPZVec<REAL> &qsi, TPZSolVec &sol,
                                                                                                                 TPZGradSolVec &dsol,  const TPZFMatrix<REAL> &axes);
                
                void AppendOtherSolution(TPZVec<REAL> &qsi, TPZVec<REAL> &normal,
                                                                                                                 TPZSolVec &leftsol, TPZGradSolVec &dleftsol, TPZFMatrix<REAL> &leftaxes,
                                                                                                                 TPZSolVec &rightsol, TPZGradSolVec &drightsol,TPZFMatrix<REAL> &rightaxes);
};

void AdjustSolutionDerivatives(TPZFMatrix<STATE> &dsolfrom, TPZFMatrix<REAL> &axesfrom,
                               TPZFMatrix<STATE> &dsolto, const TPZFMatrix<REAL> &axesto);

TPZCompEl *CreateReferredDisc(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateReferredPointEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateReferredLinearEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateReferredQuadEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateReferredTriangleEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateReferredCubeEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateReferredPrismEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateReferredPyramEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateReferredTetraEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);

template<class TVar>
void Append(TPZVec<TVar> &u1, TPZVec<TVar> &u2, TPZVec<TVar> &u12);
template<class TVar>
void Append(TPZFMatrix<TVar> &u1, TPZFMatrix<TVar> &u2, TPZFMatrix<TVar> &u12);
bool AreEqual(const TPZVec<REAL> &A, const TPZVec<REAL> &B, REAL tol = 1e-10);

#endif