variable_structure.hpp

Go to the documentation of this file.

#pragma once

#include "../../Common/include/mpi_structure.hpp"

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

#include "../../Common/include/config_structure.hpp"
#include "fluid_model.hpp"


using namespace std;

class CVariable {
protected:
  
  su2double *Solution,    
  *Solution_Old;      
  bool Non_Physical;      
  su2double *Solution_time_n,  
  *Solution_time_n1;      
  su2double **Gradient;    
  su2double *Limiter;        
  su2double *Solution_Max;    
  su2double *Solution_Min;    
  su2double AuxVar;      
  su2double *Grad_AuxVar;  
  su2double Delta_Time;  
  su2double Max_Lambda,  
  Max_Lambda_Inv,    
  Max_Lambda_Visc,  
  Lambda;        
  su2double Sensor;  
  su2double *Undivided_Laplacian;  
  su2double *Res_TruncError,  
  *Residual_Old,    
  *Residual_Sum;    
  static unsigned short nDim;    
  unsigned short nVar;    
  unsigned short nPrimVar, nPrimVarGrad;    
  unsigned short nSecondaryVar, nSecondaryVarGrad;    
public:
  
  CVariable(void);
  
  CVariable(unsigned short val_nvar, CConfig *config);
  
  CVariable(unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  virtual ~CVariable(void);
  
  void SetSolution(su2double *val_solution);
  
  void SetSolution(unsigned short val_var, su2double val_solution);
  
  void Add_DeltaSolution(unsigned short val_var, su2double val_solution);
  
  void SetNon_Physical(bool val_value);
  
  su2double GetNon_Physical(void);
  
  su2double GetSolution(unsigned short val_var);
  
  su2double GetSolution_Old(unsigned short val_var);
  
  void SetSolution_Old(su2double *val_solution_old);
  
  void SetSolution_Old(unsigned short val_var, su2double val_solution_old);
  
  void Set_OldSolution(void);
  
  void Set_Solution(void);
  
  void Set_Solution_time_n(void);
  
  void Set_Solution_time_n1(void);
  
  void Set_Solution_time_n(su2double* val_sol);
  
  void Set_Solution_time_n1(su2double* val_sol);
  
  void SetVelSolutionZero(void);
  
  void SetVelSolutionVector(su2double *val_vector);
  
  void SetVelSolutionOldZero(void);
  
  void SetVelSolutionOldVector(su2double *val_vector);
  
  void SetSolutionZero(void);
  
  void SetSolutionZero(unsigned short val_var);
  
  void AddSolution(unsigned short val_var, su2double val_solution);
  
  void AddClippedSolution(unsigned short val_var, su2double val_solution,
                          su2double lowerlimit, su2double upperlimit);
  
  void AddConservativeSolution(unsigned short val_var, su2double val_solution,
                               su2double val_density, su2double val_density_old, su2double lowerlimit,
                               su2double upperlimit);
  
  su2double *GetSolution(void);
  
  su2double *GetSolution_Old(void);
  
  su2double *GetSolution_time_n(void);
  
  su2double *GetSolution_time_n1(void);
  
  void SetResidual_Old(su2double *val_residual_old);
  
  void AddResidual_Sum(su2double *val_residual);
  
  void SetResidualSumZero(void);
  
  virtual void SetVel_ResTruncError_Zero(unsigned short iSpecies);
  
  su2double *GetResidual_Sum(void);
  
  su2double *GetResidual_Old(void);
  
  void GetResidual_Sum(su2double *val_residual);
  
  void SetAuxVar(su2double val_auxvar);
  
  su2double GetAuxVar(void);
  
  void SetAuxVarGradientZero(void);
  
  void SetAuxVarGradient(unsigned short val_dim, su2double val_gradient);
  
  void AddAuxVarGradient(unsigned short val_dim, su2double val_value);
  
  void SubtractAuxVarGradient(unsigned short val_dim, su2double val_value);
  
  su2double *GetAuxVarGradient(void);
  
  su2double GetAuxVarGradient(unsigned short val_dim);
  
  void AddRes_TruncError(su2double *val_truncation_error);
  
  void SubtractRes_TruncError(su2double *val_truncation_error);
  
  void SetRes_TruncErrorZero(void);
  
  void SetVal_ResTruncError_Zero(unsigned short val_var);
  
  void SetVel_ResTruncError_Zero(void);
  
  void SetEnergy_ResTruncError_Zero(void);
  
  su2double *GetResTruncError(void);
  
  void GetResTruncError(su2double *val_trunc_error);
  
  void SetGradient(su2double **val_gradient);
  
  void SetGradient(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  void SetGradientZero(void);
  
  void AddGradient(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  void SubtractGradient(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  su2double **GetGradient(void);
  
  su2double GetGradient(unsigned short val_var, unsigned short val_dim);
  
  void SetLimiter(unsigned short val_var, su2double val_limiter);
  
  virtual void SetLimiterPrimitive(unsigned short val_species, unsigned short val_var, su2double val_limiter);
  
  virtual su2double GetLimiterPrimitive(unsigned short val_species, unsigned short val_var);
  
  void SetSolution_Max(unsigned short val_var, su2double val_solution);
  
  void SetSolution_Min(unsigned short val_var, su2double val_solution);
  
  su2double *GetLimiter(void);
  
  su2double GetLimiter(unsigned short val_var);
  
  su2double GetSolution_Max(unsigned short val_var);
  
  su2double GetSolution_Min(unsigned short val_var);
  
  virtual su2double GetPreconditioner_Beta();
  
  virtual void SetPreconditioner_Beta(su2double val_Beta);
  
  virtual su2double* GetWindGust();
  
  virtual void SetWindGust(su2double* val_WindGust);
  
  virtual su2double* GetWindGustDer();
  
  virtual void SetWindGustDer(su2double* val_WindGust);
  
  void SetDelta_Time(su2double val_delta_time);
  
  virtual void SetDelta_Time(su2double val_delta_time, unsigned short iSpecies);
  
  su2double GetDelta_Time(void);
  
  virtual su2double GetDelta_Time(unsigned short iSpecies);
  
  void SetMax_Lambda(su2double val_max_lambda);
  
  void SetMax_Lambda_Inv(su2double val_max_lambda);
  
  virtual void SetMax_Lambda_Inv(su2double val_max_lambda, unsigned short val_species);
  
  void SetMax_Lambda_Visc(su2double val_max_lambda);
  
  virtual void SetMax_Lambda_Visc(su2double val_max_lambda, unsigned short val_species);
  
  void AddMax_Lambda(su2double val_max_lambda);
  
  void AddMax_Lambda_Inv(su2double val_max_lambda);
  
  void AddMax_Lambda_Visc(su2double val_max_lambda);
  
  su2double GetMax_Lambda(void);
  
  su2double GetMax_Lambda_Inv(void);
  
  su2double GetMax_Lambda_Visc(void);
  
  void SetLambda(su2double val_lambda);
  
  virtual void SetLambda(su2double val_lambda, unsigned short val_iSpecies);
  
  void AddLambda(su2double val_lambda);
  
  virtual void AddLambda(su2double val_lambda, unsigned short val_iSpecies);
  
  su2double GetLambda(void);
  
  virtual su2double GetLambda(unsigned short val_iSpecies);
  
  void SetSensor(su2double val_sensor);
  
  virtual void SetSensor(su2double val_sensor, unsigned short iSpecies);
  
  su2double GetSensor(void);
  
  virtual su2double GetSensor(unsigned short iSpecies);
  
  void SetUndivided_Laplacian(unsigned short val_var, su2double val_undivided_laplacian);
  
  void AddUnd_Lapl(su2double *val_und_lapl);
  
  void SubtractUnd_Lapl(su2double *val_und_lapl);
  
  void SubtractUnd_Lapl(unsigned short val_var, su2double val_und_lapl);
  
  void SetUnd_LaplZero(void);
  
  void SetUnd_Lapl(unsigned short val_var, su2double val_und_lapl);
  
  su2double *GetUndivided_Laplacian(void);
  
  su2double GetUndivided_Laplacian(unsigned short val_var);
  
  virtual su2double GetDensity(void);
  
  virtual su2double GetDensity(unsigned short val_iSpecies);
  
  virtual su2double GetMassFraction(unsigned short val_Species);
  
  virtual su2double GetEnergy(void);
  
  virtual su2double *GetForceProj_Vector(void);
  
  virtual su2double *GetObjFuncSource(void);
  
  virtual su2double *GetIntBoundary_Jump(void);
  
  virtual su2double GetEddyViscosity(void);
  
  virtual su2double GetEnthalpy(void);
  
  virtual su2double GetPressure(void);
  
  virtual su2double GetProjVel(su2double *val_vector);
  
  virtual su2double GetProjVel(su2double *val_vector, unsigned short val_species);
  
  virtual su2double GetSoundSpeed(void);
  
  virtual su2double GetBetaInc2(void);
  
  virtual su2double GetTemperature(void);
  
  virtual su2double GetTemperature_ve(void);
  
  virtual su2double GetRhoCv_tr(void);
  
  virtual su2double GetRhoCv_ve(void);
  
  virtual su2double GetVelocity(unsigned short val_dim);
  
  virtual su2double GetVelocity2(void);
  
  virtual su2double GetVelocity2(unsigned short val_species);
  
  virtual su2double GetLaminarViscosity(void);
  
  
  virtual su2double GetLaminarViscosity(unsigned short iSpecies);
  
  virtual su2double* GetDiffusionCoeff(void);
  
  virtual su2double GetThermalConductivity(void);
  
  virtual su2double GetSpecificHeatCp(void);
  
  virtual su2double GetThermalConductivity_ve(void);
  
  virtual void SetGammaSep(su2double gamma_sep);
  
  virtual void SetGammaEff(void);
  
  virtual su2double GetIntermittency();
  
  virtual su2double *GetVorticity(void);
  
  virtual su2double GetStrainMag(void);
  
  virtual void SetForceProj_Vector(su2double *val_ForceProj_Vector);
  
  virtual void SetObjFuncSource(su2double *val_SetObjFuncSource);
  
  virtual void SetIntBoundary_Jump(su2double *val_IntBoundary_Jump);
  
  virtual void SetEddyViscosity(su2double eddy_visc);
  
  virtual void SetEnthalpy(void);
  
  virtual bool SetPrimVar(CConfig *config);
  
  virtual bool SetPrimVar(CFluidModel *FluidModel);
  
  virtual void SetSecondaryVar(CFluidModel *FluidModel);
  
  virtual bool Cons2PrimVar(CConfig *config, su2double *U, su2double *V,
                            su2double *dPdU, su2double *dTdU,
                            su2double *dTvedU);
  virtual void Prim2ConsVar(CConfig *config, su2double *V, su2double *U);
  
  virtual bool SetPrimVar(su2double SharpEdge_Distance, bool check, CConfig *config);
  
  virtual bool SetPrimVar(su2double eddy_visc, su2double turb_ke, CConfig *config);
  
  virtual bool SetPrimVar(su2double eddy_visc, su2double turb_ke, CFluidModel *FluidModel);
  
  virtual bool SetPrimVar(su2double Density_Inf, CConfig *config);
  
  virtual bool SetPrimVar(su2double Density_Inf, su2double Viscosity_Inf, su2double eddy_visc, su2double turb_ke, CConfig *config);
  
  virtual su2double GetPrimitive(unsigned short val_var);
  
  virtual void SetPrimitive(unsigned short val_var, su2double val_prim);
  
  virtual void SetPrimitive(su2double *val_prim);
  
  virtual su2double *GetPrimitive(void);
  
  virtual su2double GetSecondary(unsigned short val_var);
  
  virtual void SetSecondary(unsigned short val_var, su2double val_secondary);
  
  virtual void SetSecondary(su2double *val_secondary);
  
  virtual void SetdPdrho_e(su2double dPdrho_e);
  
  virtual void SetdPde_rho(su2double dPde_rho);
  
  virtual void SetdTdrho_e(su2double dTdrho_e);
  
  virtual void SetdTde_rho(su2double dTde_rho);
  
  virtual void Setdmudrho_T(su2double dmudrho_T);
  
  virtual void SetdmudT_rho(su2double dmudT_rho);
  
  virtual void Setdktdrho_T(su2double dktdrho_T);
  
  virtual void SetdktdT_rho(su2double dktdT_rho);
  
  virtual su2double *GetSecondary(void);
  
  virtual void SetDensity(su2double val_density);
  
  virtual void SetPressure(void);
  
  virtual void SetVelocity(void);
  
  virtual void SetBetaInc2(su2double val_betainc2);
  
  virtual void SetPhi_Old(su2double *val_phi);
  
  virtual bool SetPressure(su2double Gamma);
  
  virtual bool SetPressure(CConfig *config);
  
  virtual bool SetPressure(su2double Gamma, su2double turb_ke);
  
  virtual su2double CalcEve(su2double *V, CConfig *config, unsigned short val_Species);
  
  virtual su2double CalcHs(su2double *V, CConfig *config, unsigned short val_Species);
  
  virtual su2double CalcCvve(su2double val_Tve, CConfig *config, unsigned short val_Species);
  
  virtual void CalcdPdU(su2double *V, CConfig *config, su2double *dPdU);
  
  virtual void CalcdTdU(su2double *V, CConfig *config, su2double *dTdU);
  
  virtual void CalcdTvedU(su2double *V, CConfig *config, su2double *dTdU);
  
  virtual su2double *GetdPdU(void);
  
  virtual su2double *GetdTdU(void);
  
  virtual su2double *GetdTvedU(void);
  
  virtual bool SetDensity(void);
  
  virtual void SetDeltaPressure(su2double *val_velocity, su2double Gamma);
  
  virtual bool SetSoundSpeed(su2double Gamma);
  
  virtual bool SetSoundSpeed(CConfig *config);
  
  virtual bool SetSoundSpeed(void);
  
  virtual bool SetTemperature(su2double Gas_Constant);
  
  virtual bool SetTemperature_ve(su2double val_Tve);
  
  virtual bool SetTemperature(CConfig *config);
  
  virtual void SetPrimitive(CConfig *config);
  
  virtual void SetPrimitive(CConfig *config, su2double *Coord);
  
  virtual void SetWallTemperature(su2double Temperature_Wall);
  
  virtual void SetWallTemperature(su2double* Temperature_Wall);
  
  virtual void SetThermalCoeff(CConfig *config);
  
  virtual void SetStress_FEM(unsigned short iVar, su2double val_stress);
  
  virtual void AddStress_FEM(unsigned short iVar, su2double val_stress);
  
  virtual su2double *GetStress_FEM(void);
  
  virtual void SetVonMises_Stress(su2double val_stress);
  
  virtual su2double GetVonMises_Stress(void);
  
  virtual void SetTraction(unsigned short iVar, unsigned short jVar, su2double val_traction);
  
  virtual void AddTraction(unsigned short iVar, unsigned short jVar, su2double val_traction);
  
  virtual su2double **GetTraction(void);
  
  virtual void Add_SurfaceLoad_Res(su2double *val_surfForce);
  
  virtual su2double *Get_SurfaceLoad_Res(void);
  
  virtual su2double Get_SurfaceLoad_Res(unsigned short iVar);
  
  virtual void Clear_SurfaceLoad_Res(void);
  
  virtual void Set_SurfaceLoad_Res_n(void);
  
  virtual su2double Get_SurfaceLoad_Res_n(unsigned short iVar);
  
  virtual void Add_BodyForces_Res(su2double *val_bodyForce);
  
  virtual su2double *Get_BodyForces_Res(void);
  
  virtual su2double Get_BodyForces_Res(unsigned short iVar);
  
  virtual void Clear_BodyForces_Res(void);
  
  virtual void Set_FlowTraction(su2double *val_flowTraction);
  
  virtual void Add_FlowTraction(su2double *val_flowTraction);
  
  virtual su2double *Get_FlowTraction(void);
  
  virtual su2double Get_FlowTraction(unsigned short iVar);
  
  virtual void Set_FlowTraction_n(void);
  
  virtual su2double Get_FlowTraction_n(unsigned short iVar);
  
  virtual void Clear_FlowTraction(void);
  
  virtual void SetVelocity2(void);
  
  virtual void SetVelocity_Old(su2double *val_velocity);
  
  virtual void SetLaminarViscosity(su2double laminarViscosity);
  
  virtual void SetLaminarViscosity(CConfig *config);
  
  virtual void SetThermalConductivity(su2double thermalConductivity);
  
  virtual void SetThermalConductivity(CConfig *config);
  
  virtual void SetSpecificHeatCp(su2double Cp);
  
  virtual bool SetVorticity(bool val_limiter);
  
  virtual bool SetStrainMag(bool val_limiter);
  
  virtual void SetVelSolutionOldDVector(void);
  
  virtual void SetVelSolutionDVector(void);
  
  virtual void SetGradient_PrimitiveZero(unsigned short val_primvar);
  
  virtual void AddGradient_Primitive(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  virtual void SubtractGradient_Primitive(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  virtual su2double GetGradient_Primitive(unsigned short val_var, unsigned short val_dim);
  
  virtual su2double GetLimiter_Primitive(unsigned short val_var);
  
  virtual void SetGradient_Primitive(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  virtual void SetLimiter_Primitive(unsigned short val_var, su2double val_value);
  
  virtual su2double **GetGradient_Primitive(void);
  
  virtual su2double *GetLimiter_Primitive(void);
  
  virtual void SetGradient_SecondaryZero(unsigned short val_secondaryvar);
  
  virtual void AddGradient_Secondary(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  virtual void SubtractGradient_Secondary(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  virtual su2double GetGradient_Secondary(unsigned short val_var, unsigned short val_dim);
  
  virtual su2double GetLimiter_Secondary(unsigned short val_var);
  
  virtual void SetGradient_Secondary(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  virtual void SetLimiter_Secondary(unsigned short val_var, su2double val_value);
  
  virtual su2double **GetGradient_Secondary(void);
  
  virtual su2double *GetLimiter_Secondary(void);
  
  virtual void SetBlendingFunc(su2double val_viscosity, su2double val_dist, su2double val_density);
  
  virtual su2double GetF1blending(void);
  
  virtual su2double GetF2blending(void);
  
  virtual su2double GetCrossDiff(void) { return 0.0; };
  
  virtual su2double GetmuT(void);
  
  virtual void SetmuT(su2double val_muT);
  
  virtual void AddMax_Lambda_Inv(su2double val_max_lambda, unsigned short iSpecies);
  
  virtual void AddMax_Lambda_Visc(su2double val_max_lambda, unsigned short iSpecies);
  
  virtual void SetHarmonicBalance_Source(unsigned short val_var, su2double val_source);
  
  virtual su2double GetHarmonicBalance_Source(unsigned short val_var);
  
  virtual void SetEddyViscSens(su2double *val_EddyViscSens, unsigned short numTotalVar);
  
  virtual su2double *GetEddyViscSens(void);
  
  virtual void SetSolution_Direct(su2double *val_solution_direct);
  
  virtual su2double *GetSolution_Direct(void);
  
  virtual void SetSolution_time_n(void);
  
  virtual void SetSolution_time_n(unsigned short val_var, su2double val_solution);
  
  virtual void SetSolution_time_n(su2double *val_solution_time_n);
  
  
  virtual void SetSolution_Vel(su2double *val_solution);
  
  virtual void SetSolution_Vel(unsigned short val_var, su2double val_solution_vel);
  
  virtual void SetSolution_Vel_time_n(su2double *val_solution_vel_time_n);
  
  virtual void SetSolution_Vel_time_n(void);
  
  virtual void SetSolution_Vel_time_n(unsigned short val_var, su2double val_solution_vel_time_n);
  
  virtual su2double GetSolution_time_n(unsigned short val_var);
  
  virtual su2double GetSolution_Vel(unsigned short val_var);
  
  virtual su2double *GetSolution_Vel(void);
  
  virtual su2double GetSolution_Vel_time_n(unsigned short val_var);
  
  virtual su2double *GetSolution_Vel_time_n(void);
  
  
  virtual void SetSolution_Accel(su2double *val_solution_accel);
  
  virtual void SetSolution_Accel(unsigned short val_var, su2double val_solution_accel);
  
  virtual void SetSolution_Accel_time_n(su2double *val_solution_accel_time_n);
  
  virtual void SetSolution_Accel_time_n(void);
  
  virtual void SetSolution_Accel_time_n(unsigned short val_var, su2double val_solution_accel_time_n);
  
  virtual su2double GetSolution_Accel(unsigned short val_var);
  
  virtual su2double *GetSolution_Accel(void);
  
  virtual su2double GetSolution_Accel_time_n(unsigned short val_var);
  
  virtual su2double *GetSolution_Accel_time_n(void);
  
  
  virtual void SetSolution_Pred(void);
  
  virtual void SetSolution_Pred(su2double *val_solution_pred);
  
  virtual void SetSolution_Pred(unsigned short val_var, su2double val_solution_pred);
  
  virtual su2double GetSolution_Pred(unsigned short val_var);
  
  virtual su2double *GetSolution_Pred(void);
  
  virtual void SetSolution_Pred_Old(void);
  
  virtual void SetSolution_Pred_Old(su2double *val_solution_pred_Old);
  
  virtual void SetSolution_Pred_Old(unsigned short val_var, su2double val_solution_pred_old);
  
  virtual su2double GetSolution_Pred_Old(unsigned short val_var);
  
  virtual su2double *GetSolution_Pred_Old(void);
  
  virtual void SetPrestretch(unsigned short iVar, su2double val_prestretch);
  
  virtual su2double *GetPrestretch(void);
  
  virtual su2double GetPrestretch(unsigned short iVar);
  
  void RegisterSolution(bool input);
  
  void RegisterSolution_time_n();
  
  void RegisterSolution_time_n1();
  
  void SetAdjointSolution(su2double *adj_sol);
  
  void GetAdjointSolution(su2double *adj_sol);
  
  void SetAdjointSolution_time_n(su2double *adj_sol);
  
  void GetAdjointSolution_time_n(su2double *adj_sol);
  
  void SetAdjointSolution_time_n1(su2double *adj_sol);
  
  void GetAdjointSolution_time_n1(su2double *adj_sol);
  
  virtual void SetSensitivity(unsigned short iDim, su2double val);
  
  virtual su2double GetSensitivity(unsigned short iDim);
  
  virtual void SetDual_Time_Derivative(unsigned short iVar, su2double der);
  
  virtual void SetDual_Time_Derivative_n(unsigned short iVar, su2double der);
  
  virtual su2double GetDual_Time_Derivative(unsigned short iVar);
  
  virtual su2double GetDual_Time_Derivative_n(unsigned short iVar);
};

class CBaselineVariable : public CVariable {
public:
  
  CBaselineVariable(void);
  
  CBaselineVariable(su2double *val_solution, unsigned short val_nvar, CConfig *config);
  
  virtual ~CBaselineVariable(void);
  
};

class CPotentialVariable : public CVariable {
  su2double *Charge_Density;
public:
  
  CPotentialVariable(void);
  
  CPotentialVariable(su2double val_potential, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  ~CPotentialVariable(void);
  
  su2double* GetChargeDensity();
  
  void SetChargeDensity(su2double positive_charge, su2double negative_charge);
  
};

class CWaveVariable : public CVariable {
protected:
  su2double *Solution_Direct;  
public:
  
  CWaveVariable(void);
  
  CWaveVariable(su2double *val_wave, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  ~CWaveVariable(void);
  
  void SetSolution_Direct(su2double *val_solution_direct);
  
  su2double *GetSolution_Direct(void);
  
};

class CHeatVariable : public CVariable {
protected:
  su2double *Solution_Direct;  
public:
  
  CHeatVariable(void);
  
  CHeatVariable(su2double *val_Heat, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  ~CHeatVariable(void);
  
  void SetSolution_Direct(su2double *val_solution_direct);
  
  su2double *GetSolution_Direct(void);
  
};

class CFEM_ElasVariable : public CVariable {
protected:
  
  bool dynamic_analysis;          
  bool fsi_analysis;            
  su2double *Stress;              
  su2double *FlowTraction;          
  su2double *FlowTraction_n;          
  //  su2double *Residual_Int;          /*!< \brief Internal stress term for the calculation of the residual */
  su2double *Residual_Ext_Surf;        
  su2double *Residual_Ext_Surf_n;        
  su2double *Residual_Ext_Body;        
  su2double VonMises_Stress;         
  su2double *Solution_time_n;        
  su2double *Solution_Vel,          
  *Solution_Vel_time_n;          
  su2double *Solution_Accel,          
  *Solution_Accel_time_n;          
  su2double *Solution_Pred,          
  *Solution_Pred_Old;            
  su2double *Prestretch;        
public:
  
  CFEM_ElasVariable(void);
  
  CFEM_ElasVariable(su2double *val_fea, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  ~CFEM_ElasVariable(void);
  
  su2double *GetStress_FEM(void);
  
  void SetStress_FEM(unsigned short iVar, su2double val_stress);
  
  void AddStress_FEM(unsigned short iVar, su2double val_stress);
  
  void Add_SurfaceLoad_Res(su2double *val_surfForce);
  
  su2double *Get_SurfaceLoad_Res(void);
  
  su2double Get_SurfaceLoad_Res(unsigned short iVar);
  
  void Clear_SurfaceLoad_Res(void);
  
  void Set_SurfaceLoad_Res_n(void);
  
  su2double Get_SurfaceLoad_Res_n(unsigned short iVar);
  
  void Add_BodyForces_Res(su2double *val_bodyForce);
  
  void Clear_BodyForces_Res(void);
  
  su2double *Get_BodyForces_Res(void);
  
  su2double Get_BodyForces_Res(unsigned short iVar);
  
  void Set_FlowTraction(su2double *val_flowTraction);
  
  void Add_FlowTraction(su2double *val_flowTraction);
  
  su2double *Get_FlowTraction(void);
  
  su2double Get_FlowTraction(unsigned short iVar);
  
  void Set_FlowTraction_n(void);
  
  su2double Get_FlowTraction_n(unsigned short iVar);
  
  void Clear_FlowTraction(void);
  
  void SetSolution_time_n(void);
  
  void SetSolution_time_n(su2double *val_solution_time_n);
  
  void SetSolution_time_n(unsigned short val_var, su2double val_solution);
  
  void SetSolution_Vel(su2double *val_solution_vel);
  
  void SetSolution_Vel(unsigned short val_var, su2double val_solution_vel);
  
  void SetSolution_Vel_time_n(void);
  
  void SetSolution_Vel_time_n(su2double *val_solution_vel_time_n);
  
  void SetSolution_Vel_time_n(unsigned short val_var, su2double val_solution_vel_time_n);
  
  su2double GetSolution_time_n(unsigned short val_var);
  
  su2double GetSolution_Vel(unsigned short val_var);
  
  su2double *GetSolution_Vel(void);
  
  su2double GetSolution_Vel_time_n(unsigned short val_var);
  
  su2double *GetSolution_Vel_time_n(void);
  
  void SetSolution_Accel(su2double *val_solution_accel);
  
  void SetSolution_Accel(unsigned short val_var, su2double val_solution_accel);
  
  void SetSolution_Accel_time_n(su2double *val_solution_accel_time_n);
  
  void SetSolution_Accel_time_n(void);
  
  void SetSolution_Accel_time_n(unsigned short val_var, su2double val_solution_accel_time_n);
  
  su2double GetSolution_Accel(unsigned short val_var);
  
  su2double *GetSolution_Accel(void);
  
  su2double GetSolution_Accel_time_n(unsigned short val_var);
  
  su2double *GetSolution_Accel_time_n(void);
  
  void SetSolution_Pred(void);
  
  void SetSolution_Pred(su2double *val_solution_pred);
  
  void SetSolution_Pred(unsigned short val_var, su2double val_solution_pred);
  
  su2double GetSolution_Pred(unsigned short val_var);
  
  su2double *GetSolution_Pred(void);
  
  void SetSolution_Pred_Old(void);
  
  void SetSolution_Pred_Old(su2double *val_solution_pred_Old);
  
  void SetSolution_Pred_Old(unsigned short val_var, su2double val_solution_pred_old);
  
  su2double GetSolution_Pred_Old(unsigned short val_var);
  
  su2double *GetSolution_Pred_Old(void);
  
  void SetPrestretch(unsigned short iVar, su2double val_prestretch);
  
  su2double *GetPrestretch(void);
  
  su2double GetPrestretch(unsigned short iVar);
  
  void SetVonMises_Stress(su2double val_stress);
  
  su2double GetVonMises_Stress(void);
  
  
};

class CFEABoundVariable : public CVariable {
protected:
  su2double **Traction;  
public:
  
  CFEABoundVariable(void);
  
  CFEABoundVariable(unsigned short val_nDim, unsigned short val_nvar, unsigned short val_nElBound, CConfig *config);
  
  ~CFEABoundVariable(void);
  
  void SetTraction(unsigned short iVar, unsigned short jVar, su2double val_traction);
  
  void AddTraction(unsigned short iVar, unsigned short jVar, su2double val_traction);
  
  su2double **GetTraction(void);
  
};

class CEulerVariable : public CVariable {
protected:
  su2double  Velocity2;      
  su2double *HB_Source;     
  su2double  Precond_Beta;  
  su2double *WindGust;      
  su2double *WindGustDer;   
  /*--- Primitive variable definition ---*/
  
  su2double *Primitive;  
  su2double **Gradient_Primitive;  
  su2double *Limiter_Primitive;    
  /*--- Secondary variable definition ---*/
  
  su2double *Secondary;            
  su2double **Gradient_Secondary;  
  su2double *Limiter_Secondary;   
public:
  
  CEulerVariable(void);
  
  CEulerVariable(su2double val_density, su2double *val_velocity, su2double val_energy, unsigned short val_nDim,
                 unsigned short val_nvar, CConfig *config);
  
  CEulerVariable(su2double *val_solution, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  virtual ~CEulerVariable(void);
  
  void SetGradient_PrimitiveZero(unsigned short val_primvar);
  
  void AddGradient_Primitive(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  void SubtractGradient_Primitive(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  su2double GetGradient_Primitive(unsigned short val_var, unsigned short val_dim);
  
  su2double GetLimiter_Primitive(unsigned short val_var);
  
  void SetGradient_Primitive(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  void SetLimiter_Primitive(unsigned short val_var, su2double val_value);
  
  su2double **GetGradient_Primitive(void);
  
  su2double *GetLimiter_Primitive(void);
  
  void SetGradient_SecondaryZero(unsigned short val_secondaryvar);
  
  void AddGradient_Secondary(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  void SubtractGradient_Secondary(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  su2double GetGradient_Secondary(unsigned short val_var, unsigned short val_dim);
  
  su2double GetLimiter_Secondary(unsigned short val_var);
  
  void SetGradient_Secondary(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  void SetLimiter_Secondary(unsigned short val_var, su2double val_value);
  
  su2double **GetGradient_Secondary(void);
  
  su2double *GetLimiter_Secondary(void);
  
  void SetdPdrho_e(su2double dPdrho_e);
  
  void SetdPde_rho(su2double dPde_rho);
  
  bool SetPressure(su2double Gamma);
  
  bool SetSoundSpeed(su2double Gamma);
  
  void SetEnthalpy(void);
  
  bool SetPrimVar(CFluidModel *FluidModel);
  
  void SetSecondaryVar(CFluidModel *FluidModel);
  
  su2double GetPrimitive(unsigned short val_var);
  
  void SetPrimitive(unsigned short val_var, su2double val_prim);
  
  void SetPrimitive(su2double *val_prim);
  
  su2double *GetPrimitive(void);
  
  su2double GetSecondary(unsigned short val_var);
  
  void SetSecondary(unsigned short val_var, su2double val_secondary);
  
  void SetSecondary(su2double *val_secondary);
  
  su2double *GetSecondary(void);
  
  bool SetDensity(void);
  
  bool SetTemperature(su2double Gas_Constant);
  
  su2double GetVelocity2(void);
  
  su2double GetPressure(void);
  
  su2double GetSoundSpeed(void);
  
  su2double GetEnthalpy(void);
  
  su2double GetDensity(void);
  
  su2double GetEnergy(void);
  
  su2double GetTemperature(void);
  
  su2double GetVelocity(unsigned short val_dim);
  
  su2double GetProjVel(su2double *val_vector);
  
  void SetVelocity(void);
  
  void SetVelocity_Old(su2double *val_velocity);
  
  void SetHarmonicBalance_Source(unsigned short val_var, su2double val_source);
  
  su2double GetHarmonicBalance_Source(unsigned short val_var);
  
  su2double GetPreconditioner_Beta();
  
  void SetPreconditioner_Beta(su2double val_Beta);
  
  su2double* GetWindGust();
  
  void SetWindGust(su2double* val_WindGust);
  
  su2double* GetWindGustDer();
  
  void SetWindGustDer(su2double* val_WindGust);
};

class CIncEulerVariable : public CVariable {
protected:
  su2double Velocity2;      
  su2double Precond_Beta;  
  su2double *WindGust;           
  su2double *WindGustDer;        
  /*--- Primitive variable definition ---*/
  
  su2double *Primitive;  
  su2double **Gradient_Primitive;  
  su2double *Limiter_Primitive;    
public:
  
  CIncEulerVariable(void);
  
  CIncEulerVariable(su2double val_pressure, su2double *val_velocity, unsigned short val_nDim,
                    unsigned short val_nvar, CConfig *config);
  
  CIncEulerVariable(su2double *val_solution, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  virtual ~CIncEulerVariable(void);
  
  void SetGradient_PrimitiveZero(unsigned short val_primvar);
  
  void AddGradient_Primitive(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  void SubtractGradient_Primitive(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  su2double GetGradient_Primitive(unsigned short val_var, unsigned short val_dim);
  
  su2double GetLimiter_Primitive(unsigned short val_var);
  
  void SetGradient_Primitive(unsigned short val_var, unsigned short val_dim, su2double val_value);
  
  void SetLimiter_Primitive(unsigned short val_var, su2double val_value);
  
  su2double **GetGradient_Primitive(void);
  
  su2double *GetLimiter_Primitive(void);
  
  void SetPressure();
  
  su2double GetPrimitive(unsigned short val_var);
  
  void SetPrimitive(unsigned short val_var, su2double val_prim);
  
  void SetPrimitive(su2double *val_prim);
  
  su2double *GetPrimitive(void);
  
  void SetDensity(su2double val_density);
  
  void SetVelocity(void);
  
  void SetBetaInc2(su2double val_betainc2);
  
  su2double GetVelocity2(void);
  
  su2double GetPressure(void);
  
  su2double GetBetaInc2(void);
  
  su2double GetDensity(void);
  
  su2double GetVelocity(unsigned short val_dim);
  
  su2double GetProjVel(su2double *val_vector);
  
  void SetVelocity_Old(su2double *val_velocity);

  su2double* GetWindGust();
  
  void SetWindGust(su2double* val_WindGust);
  
  su2double* GetWindGustDer();
  
  void SetWindGustDer(su2double* val_WindGust);
  
  bool SetPrimVar(su2double Density_Inf, CConfig *config);
  
};

class CNSVariable : public CEulerVariable {
private:
  su2double Prandtl_Lam;     
  su2double Prandtl_Turb;    
  su2double Temperature_Ref; 
  su2double Viscosity_Ref;   
  su2double Viscosity_Inf;   
  su2double Vorticity[3];    
  su2double StrainMag;       
public:
  
  CNSVariable(void);
  
  CNSVariable(su2double val_density, su2double *val_velocity,
              su2double val_energy, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  CNSVariable(su2double *val_solution, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  ~CNSVariable(void);
  
  void SetLaminarViscosity(su2double laminarViscosity);
  
  void SetThermalConductivity(su2double thermalConductivity);
  
  void SetSpecificHeatCp(su2double Cp);
  
  bool SetVorticity(bool val_limiter);
  
  bool SetStrainMag(bool val_limiter);
  
  void SetEddyViscosity(su2double eddy_visc);
  
  su2double GetLaminarViscosity(void);
  
  su2double GetThermalConductivity(void);
  
  su2double GetEddyViscosity(void);
  
  su2double GetSpecificHeatCp(void);
  
  void SetWallTemperature(su2double temperature_wall);
  
  su2double *GetVorticity(void);
  
  su2double GetStrainMag(void);
  
  void SetdTdrho_e(su2double dTdrho_e);
  
  void SetdTde_rho(su2double dTde_rho);
  
  void Setdmudrho_T(su2double dmudrho_T);
  
  void SetdmudT_rho(su2double dmudT_rho);
  
  void Setdktdrho_T(su2double dktdrho_T);
  
  void SetdktdT_rho(su2double dktdT_rho);
  
  bool SetPrimVar(su2double eddy_visc, su2double turb_ke, CFluidModel *FluidModel);
  using CVariable::SetPrimVar;
  
  void SetSecondaryVar(CFluidModel *FluidModel);
};

class CIncNSVariable : public CIncEulerVariable {
private:
  su2double Prandtl_Lam;     
  su2double Prandtl_Turb;    
  su2double Temperature_Ref; 
  su2double Viscosity_Ref;   
  su2double Viscosity_Inf;   
  su2double Vorticity[3];    
  su2double StrainMag;       
public:
  
  CIncNSVariable(void);
  
  CIncNSVariable(su2double val_pressure, su2double *val_velocity, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  CIncNSVariable(su2double *val_solution, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  ~CIncNSVariable(void);
  
  void SetLaminarViscosity(su2double laminarViscosity);
  
  bool SetVorticity(bool val_limiter);
  
  bool SetStrainMag(bool val_limiter);
  
  void SetEddyViscosity(su2double eddy_visc);
  
  su2double GetLaminarViscosity(void);
  
  su2double GetEddyViscosity(void);
  
  su2double *GetVorticity(void);
  
  su2double GetStrainMag(void);
  
  bool SetPrimVar(su2double Density_Inf, su2double Viscosity_Inf, su2double eddy_visc, su2double turb_ke, CConfig *config);
  using CVariable::SetPrimVar;
  
};

class CTurbVariable : public CVariable {
protected:
  su2double muT;                
  su2double *HB_Source;          
public:
  CTurbVariable(void);
  
  CTurbVariable(unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  virtual ~CTurbVariable(void);
  
  su2double GetmuT();
  
  void SetmuT(su2double val_muT);
};

class CTurbSAVariable : public CTurbVariable {
public:
  CTurbSAVariable(void);
  
  CTurbSAVariable(su2double val_nu_tilde, su2double val_muT, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  ~CTurbSAVariable(void);
  
  void SetHarmonicBalance_Source(unsigned short val_var, su2double val_source);
  
  su2double GetHarmonicBalance_Source(unsigned short val_var);
  
};


class CTurbMLVariable : public CTurbVariable {
public:
  CTurbMLVariable(void);
  
  CTurbMLVariable(su2double val_nu_tilde, su2double val_muT, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  ~CTurbMLVariable(void);
  
  void SetHarmonicBalance_Source(unsigned short val_var, su2double val_source);
  
  su2double GetHarmonicBalance_Source(unsigned short val_var);
  
};

class CTransLMVariable : public CTurbVariable {
protected:
  su2double gamma_sep;
  
public:
  
  CTransLMVariable(void);
  
  CTransLMVariable(su2double val_nu_tilde, su2double val_intermittency, su2double val_REth, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  ~CTransLMVariable(void);
  
  su2double GetIntermittency(void);
  
  void SetGammaSep(su2double gamma_sep_in);
  
  void SetGammaEff(void);
  
};

class CTurbSSTVariable : public CTurbVariable {
protected:
  su2double sigma_om2,
  beta_star;
  su2double F1,    
  F2,            
  CDkw;           
public:
  CTurbSSTVariable(void);
  
  CTurbSSTVariable(su2double val_rho_kine, su2double val_rho_omega, su2double val_muT, unsigned short val_nDim, unsigned short val_nvar,
                   su2double *constants, CConfig *config);
  
  ~CTurbSSTVariable(void);
  
  void SetBlendingFunc(su2double val_viscosity, su2double val_dist, su2double val_density);
  
  su2double GetF1blending(void);
  
  su2double GetF2blending(void);
  
  su2double GetCrossDiff(void);
};


class CAdjEulerVariable : public CVariable {
protected:
  su2double *Psi;    
  su2double *ForceProj_Vector;  
  su2double *ObjFuncSource;    
  su2double *IntBoundary_Jump;  
  su2double *HB_Source;    
  bool incompressible;
public:
  
  CAdjEulerVariable(void);
  
  CAdjEulerVariable(su2double val_psirho, su2double *val_phi, su2double val_psie, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  CAdjEulerVariable(su2double *val_solution, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  virtual ~CAdjEulerVariable(void);
  
  bool SetPrimVar(su2double SharpEdge_Distance, bool check, CConfig *config);
  
  void SetPhi_Old(su2double *val_phi);
  
  void SetForceProj_Vector(su2double *val_ForceProj_Vector);
  
  void SetObjFuncSource(su2double *val_SetObjFuncSource);
  
  void SetIntBoundary_Jump(su2double *val_IntBoundary_Jump);
  
  su2double *GetForceProj_Vector(void);
  
  su2double *GetObjFuncSource(void);
  
  su2double *GetIntBoundary_Jump(void);
  
  void SetHarmonicBalance_Source(unsigned short val_var, su2double val_source);
  
  su2double GetHarmonicBalance_Source(unsigned short val_var);
};

class CAdjIncEulerVariable : public CVariable {
protected:
  su2double *Psi;    
  su2double *ForceProj_Vector;  
  su2double *ObjFuncSource;    
  su2double *IntBoundary_Jump;  
  bool incompressible;
public:
  
  CAdjIncEulerVariable(void);
  
  CAdjIncEulerVariable(su2double val_psirho, su2double *val_phi, su2double val_psie, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  CAdjIncEulerVariable(su2double *val_solution, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  virtual ~CAdjIncEulerVariable(void);
  
  bool SetPrimVar(su2double SharpEdge_Distance, bool check, CConfig *config);
  
  void SetPhi_Old(su2double *val_phi);
  
  void SetForceProj_Vector(su2double *val_ForceProj_Vector);
  
  void SetObjFuncSource(su2double *val_SetObjFuncSource);
  
  void SetIntBoundary_Jump(su2double *val_IntBoundary_Jump);
  
  su2double *GetForceProj_Vector(void);
  
  su2double *GetObjFuncSource(void);
  
  su2double *GetIntBoundary_Jump(void);
  
  void SetTimeSpectral_Source(unsigned short val_var, su2double val_source);
  
  su2double GetTimeSpectral_Source(unsigned short val_var);
};

class CAdjNSVariable : public CAdjEulerVariable {  
private:
  
public:
  
  CAdjNSVariable(void);
  
  CAdjNSVariable(su2double val_psirho, su2double *val_phi, su2double val_psie, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  CAdjNSVariable(su2double *val_solution, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  ~CAdjNSVariable(void);
  
  void SetPhi_Old(su2double *val_phi);
  
  void SetForceProj_Vector(su2double *val_ForceProj_Vector);
  
  su2double *GetForceProj_Vector(void);
  
  void SetVelSolutionOldDVector(void);
  
  void SetVelSolutionDVector(void);
  
};

class CAdjIncNSVariable : public CAdjIncEulerVariable {
private:
  
public:
  
  CAdjIncNSVariable(void);
  
  CAdjIncNSVariable(su2double val_psirho, su2double *val_phi, su2double val_psie, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  CAdjIncNSVariable(su2double *val_solution, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  ~CAdjIncNSVariable(void);
  
  void SetPhi_Old(su2double *val_phi);
  
  void SetForceProj_Vector(su2double *val_ForceProj_Vector);
  
  su2double *GetForceProj_Vector(void);
  
  void SetVelSolutionOldDVector(void);
  
  void SetVelSolutionDVector(void);
  
};

class CAdjTurbVariable : public CVariable {
protected:
  su2double *dmuT_dUTvar;       
  su2double **dRTstar_dUTvar;   
  su2double **dFT_dUTvar;   
  su2double *EddyViscSens;    
public:
  
  CAdjTurbVariable(void);
  
  CAdjTurbVariable(su2double val_psinu_inf, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  ~CAdjTurbVariable(void);
  
  void SetEddyViscSens(su2double *val_EddyViscSens, unsigned short numTotalVar);
  
  su2double *GetEddyViscSens(void);
};

class CTemplateVariable : public CVariable {
public:
  
  CTemplateVariable(void);
  
  CTemplateVariable(su2double val_potential, unsigned short val_nDim, unsigned short val_nvar, CConfig *config);
  
  ~CTemplateVariable(void);
};

class CDiscAdjVariable : public CVariable {
private:
  su2double* Sensitivity; /* Vector holding the derivative of target functional with respect to the coordinates at this node*/
  su2double* Solution_Direct;
  su2double* DualTime_Derivative;
  su2double* DualTime_Derivative_n;
  
public:
  CDiscAdjVariable(void);
  
  ~CDiscAdjVariable(void);
  
  CDiscAdjVariable(su2double *val_solution, unsigned short val_ndim, unsigned short val_nvar, CConfig *config);
  
  void SetSensitivity(unsigned short iDim, su2double val);
  
  su2double GetSensitivity(unsigned short iDim);
  
  void SetDual_Time_Derivative(unsigned short iVar, su2double der);
  
  void SetDual_Time_Derivative_n(unsigned short iVar, su2double der);
  
  su2double GetDual_Time_Derivative(unsigned short iVar);
  
  su2double GetDual_Time_Derivative_n(unsigned short iVar);
  
  void SetSolution_Direct(su2double *sol);
  
  su2double* GetSolution_Direct();
};


#include "variable_structure.inl"