grid_adaptation_structure.hpp

Go to the documentation of this file.

#pragma once

#include "./mpi_structure.hpp"

#include <cmath>
#include <iostream>
#include <cstdlib>
#include <fstream>

#include "geometry_structure.hpp"
#include "config_structure.hpp"

using namespace std;

class CGridAdaptation {
protected:
    unsigned long nPoint_new,   
    nElem_new;                  
    unsigned short nDim,    
    nVar;                   
    su2double **ConsVar_Sol,    
    **ConsVar_Res,          
    **ConsVar_Adapt;        
    su2double **AdjVar_Sol, 
    **AdjVar_Res,           
    **AdjVar_Adapt;         
    su2double **LinVar_Sol, 
    **LinVar_Res,           
    **LinVar_Adapt;         
    su2double **Gradient,       
    **Gradient_Flow,        
    **Gradient_Adj;         
    su2double *Index;           
public:

    CGridAdaptation(CGeometry *geometry, CConfig *config);

    ~CGridAdaptation(void);
    
    void GetFlowSolution(CGeometry *geometry, CConfig *config);
    
    void GetFlowResidual(CGeometry *geometry, CConfig *config);
    
    void GetAdjSolution(CGeometry *geometry, CConfig *config);
    
    void GetAdjResidual(CGeometry *geometry, CConfig *config);
    
    void SetComplete_Refinement(CGeometry *geometry, unsigned short strength);
    
    void SetNo_Refinement(CGeometry *geometry, unsigned short strength);
    
    void SetWake_Refinement(CGeometry *geometry, unsigned short strength);
    
    void SetSupShock_Refinement(CGeometry *geometry, CConfig *config);
    
    void SetNearField_Refinement(CGeometry *geometry, CConfig *config);
    
    void SetHomothetic_Adaptation2D(CGeometry *geometry, CPhysicalGeometry *geo_adapt, CConfig *config);
    
    void SetHomothetic_Adaptation3D(CGeometry *geometry, CPhysicalGeometry *geo_adapt, CConfig *config);
        
    long CheckTriangleCode(bool *AdaptCode);
    
    long CheckRectCode(bool *AdaptCode);
    
    long CheckRectExtCode(bool *AdaptCode);
    
    long CheckTetraCode(bool *AdaptCode);
    
    long CheckHexaCode(bool *AdaptCode);
    
    long CheckPyramCode(bool *AdaptCode);
    
    void TriangleDivision(long code, long *nodes, long *edges, long **Division, long *nPart);
    
    void RectDivision(long code, long *nodes, long **Division, long *nPart);
    
    void RectExtDivision(long code, long *nodes, long **Division, long *nPart);
    
    void TetraDivision(long code, long *nodes, long *edges, long **Division, long *nPart);
    
    void HexaDivision(long code, long *nodes, long **Division, long *nPart);
    
    void PyramDivision(long code, long *nodes, long **Division, long *nPart);
    
    void SetIndicator_Flow(CGeometry *geometry, CConfig *config, unsigned short strength);
    
    void SetIndicator_Adj(CGeometry *geometry, CConfig *config, unsigned short strength);
    
    void SetIndicator_FlowAdj(CGeometry *geometry, CConfig *config);
    
    void SetIndicator_Robust(CGeometry *geometry, CConfig *config);
    
    void SetIndicator_Computable(CGeometry *geometry, CConfig *config);
    
    void SetIndicator_Computable_Robust(CGeometry *geometry, CConfig *config);
    
    void SetRestart_FlowSolution(CConfig *config, CPhysicalGeometry *geo_adapt, string mesh_flowfilename);
    
    void SetRestart_AdjSolution(CConfig *config, CPhysicalGeometry *geo_adapt, string mesh_adjfilename);
    
    void SetRestart_LinSolution(CConfig *config, CPhysicalGeometry *geo_adapt, string mesh_linfilename);
    
    void SetSensorElem(CGeometry *geometry, CConfig *config, unsigned long max_elem);

};

#include "grid_adaptation_structure.inl"