solver_structure.hpp

Go to the documentation of this file.

#pragma once

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

#include <cmath>
#include <string>
#include <fstream>
#include <sstream>
#include <algorithm>
#include <iostream>
#include <stdlib.h>
#include <stdio.h>

#include "fluid_model.hpp"
#include "numerics_structure.hpp"
#include "variable_structure.hpp"
#include "../../Common/include/gauss_structure.hpp"
#include "../../Common/include/element_structure.hpp"
#include "../../Common/include/geometry_structure.hpp"
#include "../../Common/include/config_structure.hpp"
#include "../../Common/include/matrix_structure.hpp"
#include "../../Common/include/vector_structure.hpp"
#include "../../Common/include/linear_solvers_structure.hpp"
#include "../../Common/include/grid_movement_structure.hpp"

using namespace std;

class CSolver {
protected:
  unsigned short IterLinSolver;  
  unsigned short nVar,          
  nPrimVar,                     
  nPrimVarGrad,                 
  nSecondaryVar,                     
  nSecondaryVarGrad,                 
  nVarGrad,                 
  nDim;                          
  unsigned long nPoint;          
  unsigned long nPointDomain;   
  su2double Max_Delta_Time,  
  Min_Delta_Time;          
  su2double *Residual_RMS,  
  *Residual_Max,        
  *Residual,            
  *Residual_i,          
  *Residual_j;          
  unsigned long *Point_Max; 
  su2double **Point_Max_Coord; 
  su2double *Solution,    
  *Solution_i,        
  *Solution_j;        
  su2double *Vector,  
  *Vector_i,      
  *Vector_j;      
  su2double *Res_Conv,  
  *Res_Visc,        
  *Res_Sour,        
  *Res_Conv_i,      
  *Res_Visc_i,      
  *Res_Conv_j,      
  *Res_Visc_j;      
  su2double **Jacobian_i,  
  **Jacobian_j;          
  su2double **Jacobian_ii,  
  **Jacobian_ij,        
  **Jacobian_ji,        
  **Jacobian_jj;        
  su2double **Smatrix,  
  **Cvector;       
  unsigned short nOutputVariables;  
public:
  
  CSysVector LinSysSol;    
  CSysVector LinSysRes;    
  CSysVector LinSysAux;    
  CSysMatrix Jacobian; 
  CSysMatrix StiffMatrix; 
  CSysVector OutputVariables;    
  string* OutputHeadingNames; /*< \brief vector of strings to store the headings for the exra variables */
  
  CVariable** node;  
  CVariable* node_infty; 
  CSolver(void);
  
  virtual ~CSolver(void);
  
  void SetIterLinSolver(unsigned short val_iterlinsolver);
  
  virtual void Set_MPI_Solution(CGeometry *geometry, CConfig *config);
  
  virtual void Set_MPI_Primitive(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief Set number of linear solver iterations.
  //   * \param[in] val_iterlinsolver - Number of linear iterations.
  //   */
  //  virtual void Set_MPI_Secondary(CGeometry *geometry, CConfig *config);
  
  void SetResidual_RMS(CGeometry *geometry, CConfig *config);
  
  virtual void Set_MPI_Solution_Old(CGeometry *geometry, CConfig *config);
  
  virtual void Set_MPI_Solution_DispOnly(CGeometry *geometry, CConfig *config);
  
  virtual void Set_MPI_Solution_Pred(CGeometry *geometry, CConfig *config);
  
  virtual void Set_MPI_Solution_Pred_Old(CGeometry *geometry, CConfig *config);
  
  virtual void Set_MPI_Solution_Limiter(CGeometry *geometry, CConfig *config);
  
  virtual void Set_MPI_Primitive_Limiter(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief Impose the send-receive boundary condition.
  //   * \param[in] geometry - Geometrical definition of the problem.
  //   * \param[in] config - Definition of the particular problem.
  //   */
  //  virtual void Set_MPI_Secondary_Limiter(CGeometry *geometry, CConfig *config);
  
  virtual void SetNondimensionalization(CGeometry *geometry, CConfig *config, unsigned short iMesh);
  
  virtual CFluidModel* GetFluidModel(void);
  
  unsigned short GetIterLinSolver(void);
  
  su2double GetMax_Delta_Time(void);
  
  su2double GetMin_Delta_Time(void);
  
  virtual su2double GetMax_Delta_Time(unsigned short val_Species);
  
  virtual su2double GetMin_Delta_Time(unsigned short val_Species);
  
  unsigned short GetnVar(void);
  
  unsigned short GetnPrimVar(void);
  
  unsigned short GetnPrimVarGrad(void);
  
  unsigned short GetnSecondaryVar(void);
  
  unsigned short GetnSecondaryVarGrad(void);
  
  unsigned short GetnOutputVariables(void);
  
  virtual void SetResidual_DualTime(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                                    unsigned short iRKStep, unsigned short iMesh, unsigned short RunTime_EqSystem);
  
  void SetRes_RMS(unsigned short val_var, su2double val_residual);
  
  void AddRes_RMS(unsigned short val_var, su2double val_residual);
  
  su2double GetRes_RMS(unsigned short val_var);
  
  void SetRes_Max(unsigned short val_var, su2double val_residual, unsigned long val_point);
  
  void AddRes_Max(unsigned short val_var, su2double val_residual, unsigned long val_point, su2double* val_coord);
  
  su2double GetRes_Max(unsigned short val_var);
  
  virtual su2double GetRes_FEM(unsigned short val_var);
  
  unsigned long GetPoint_Max(unsigned short val_var);
  
  su2double* GetPoint_Max_Coord(unsigned short val_var);
  
  void SetGrid_Movement_Residual(CGeometry *geometry, CConfig *config);
  
  void SetAuxVar_Gradient_GG(CGeometry *geometry, CConfig *config);
  
  void SetAuxVar_Gradient_LS(CGeometry *geometry, CConfig *config);
  
  void SetAuxVar_Surface_Gradient(CGeometry *geometry, CConfig *config);
  
  void SetSolution_Gradient_GG(CGeometry *geometry, CConfig *config);
  
  void SetSolution_Gradient_LS(CGeometry *geometry, CConfig *config);
  
  virtual void Set_MPI_Solution_Gradient(CGeometry *geometry, CConfig *config);
  
  void SetGridVel_Gradient(CGeometry *geometry, CConfig *config);
  
  void SetSolution_Limiter(CGeometry *geometry, CConfig *config);
  
  virtual void SetPrimitive_Limiter(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief A virtual member.
  //   * \param[in] geometry - Geometrical definition of the problem.
  //   * \param[in] config - Definition of the particular problem.
  //   */
  //  virtual void SetSecondary_Limiter(CGeometry *geometry, CConfig *config);
  
  void SetPressureLaplacian(CGeometry *geometry, CConfig *config, su2double *PressureLaplacian);
  
  void Set_OldSolution(CGeometry *geometry);
  
  void Restart_OldGeometry(CGeometry *geometry, CConfig *config);
  
  virtual void SetTime_Step(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                            unsigned short iMesh, unsigned long Iteration);
  
  virtual void Postprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                              unsigned short iMesh);
  
  virtual void Postprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, CNumerics **numerics,
                              unsigned short iMesh);
  virtual void Centered_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                                 CConfig *config, unsigned short iMesh, unsigned short iRKStep);
  
  virtual void Upwind_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                               CConfig *config, unsigned short iMesh);
  
  virtual void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  virtual void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, CNumerics **numerics, unsigned short iMesh, unsigned long Iteration, unsigned short RunTime_EqSystem, bool Output);
  
  virtual void SetUndivided_Laplacian(CGeometry *geometry, CConfig *config);
  
  virtual void Set_MPI_ActDisk(CSolver **solver_container, CGeometry *geometry, CConfig *config);
  
  virtual void Set_MPI_Nearfield(CGeometry *geometry, CConfig *config);
  
  virtual void Set_MPI_Interface(CGeometry *geometry, CConfig *config);

  virtual void Set_MPI_Undivided_Laplacian(CGeometry *geometry, CConfig *config);
  
  virtual void SetMax_Eigenvalue(CGeometry *geometry, CConfig *config);
  
  virtual void Set_MPI_MaxEigenvalue(CGeometry *geometry, CConfig *config);
  
  virtual void SetDissipation_Switch(CGeometry *geometry, CConfig *config);
  
  virtual void Set_MPI_Dissipation_Switch(CGeometry *geometry, CConfig *config);
  
  
  virtual void Compute_ComboObj(CConfig *config);
  
  virtual void BC_Euler_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                             unsigned short val_marker);
  
  virtual void BC_Clamped(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                          unsigned short val_marker);
  
  virtual void BC_Clamped_Post(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                               unsigned short val_marker);
  
  virtual void BC_Normal_Displacement(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                                      unsigned short val_marker);
  
  
  virtual void BC_Normal_Load(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                              unsigned short val_marker);
  
  virtual void BC_Dir_Load(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                           unsigned short val_marker);
  
  virtual void BC_Sine_Load(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                            unsigned short val_marker);
  
  virtual void BC_Pressure(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                           unsigned short val_marker);
  
  virtual void BC_Interface_Boundary(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config, unsigned short val_marker);
  
  virtual void BC_NearField_Boundary(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config, unsigned short val_marker);
  
  virtual void BC_Fluid_Interface(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config);

  virtual void BC_ActDisk_Inlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                                CConfig *config, unsigned short val_marker);
  
  virtual void BC_ActDisk_Outlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                                 CConfig *config, unsigned short val_marker);
  
  virtual void BC_ActDisk(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                          CConfig *config, unsigned short val_marker, bool inlet_surface);
  
  virtual void BC_Isothermal_Wall(CGeometry *geometry,
                                  CSolver **solver_container,
                                  CNumerics *conv_numerics,
                                  CNumerics *visc_numerics,
                                  CConfig *config,
                                  unsigned short val_marker);
  
  virtual void BC_HeatFlux_Wall(CGeometry *geometry, CSolver **solver_container,
                                CNumerics *conv_numerics,
                                CNumerics *visc_numerics, CConfig *config,
                                unsigned short val_marker);
  
  virtual void BC_Dirichlet(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                            unsigned short val_marker);
  
  virtual void BC_Neumann(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                          unsigned short val_marker);
  
  virtual void BC_Far_Field(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                            unsigned short val_marker);
  
  virtual void BC_Sym_Plane(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  
  virtual void BC_Riemann(CGeometry *geometry, CSolver **solver_container,
                          CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  virtual void BC_NonReflecting(CGeometry *geometry, CSolver **solver_container,
                                CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  virtual void BC_Inlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                        CConfig *config, unsigned short val_marker);
  
  virtual void BC_Supersonic_Inlet(CGeometry *geometry, CSolver **solver_container,
                                   CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  virtual void BC_Supersonic_Outlet(CGeometry *geometry, CSolver **solver_container,
                                    CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  virtual void BC_Custom(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                         CConfig *config, unsigned short val_marker);
  
  virtual void BC_Outlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                         CConfig *config, unsigned short val_marker);
  
  virtual void BC_Engine_Inflow(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  virtual void BC_Engine_Exhaust(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  virtual void BC_Dielec(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                         CConfig *config, unsigned short val_marker);
  
  virtual void BC_Electrode(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                            CConfig *config, unsigned short val_marker);
  virtual void Mixing_Process(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short val_Marker);
  
  virtual void MixedOut_Average (su2double val_init_pressure, su2double *val_Averaged_Flux, su2double *val_normal, su2double *pressure_mix, su2double *density_mix);
  
  virtual void MixedOut_Root_Function(su2double *pressure, su2double *val_Averaged_Flux, su2double *val_normal, su2double *valfunc, su2double *density);
  
  virtual void Boundary_Fourier(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short val_Marker, vector<std::complex<su2double> > &c4k,signed long &nboundaryvertex);
  
  virtual void Boundary_Fourier(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short val_Marker, vector<std::complex<su2double> >& c2k,vector<std::complex<su2double> >& c3k,signed long& nboundaryvertex);
  
  virtual void SetExtAveragedValue(CSolver *solver_container, unsigned short intMarker,  unsigned short extMarker);
  
  virtual su2double GetAveragedDensity(unsigned short valMarker);
  
  virtual su2double GetAveragedPressure(unsigned short valMarker);
  
  virtual su2double GetAveragedEnthalpy(unsigned short valMarker);
  
  virtual su2double* GetAveragedVelocity(unsigned short valMarker);
  
  virtual su2double* GetAveragedGridVelocity(unsigned short valMarker);
  
  virtual su2double GetAveragedEntropy(unsigned short valMarker);
  
  virtual su2double GetAveragedTotTemperature(unsigned short valMarker);
  
  virtual su2double GetAveragedTotPressure(unsigned short valMarker);
  
  virtual su2double GetFlowAngle(unsigned short valMarker);
  
  virtual su2double GetAveragedMach(unsigned short valMarker);
  
  virtual su2double GetAveragedNormalMach(unsigned short valMarker);
  
  virtual su2double GetMassFlow(unsigned short valMarker);
  
  virtual su2double GetTotalPressureLoss(unsigned short inMarkerTP);
  
  virtual su2double GetKineticEnergyLoss(unsigned short inMarkerTP);
  
  virtual su2double GetTotalTotalEfficiency(unsigned short inMarkerTP);
  
  virtual su2double GetTotalStaticEfficiency(unsigned short inMarkerTP);
  
  virtual su2double GetEulerianWork(unsigned short inMarkerTP);
  
  virtual su2double GetTotalEnthalpyIn(unsigned short inMarkerTP);
  
  virtual su2double GetFlowAngleIn(unsigned short inMarkerTP);
  
  virtual su2double GetFlowAngleOut(unsigned short inMarkerTP);
  
  virtual su2double GetMassFlowIn(unsigned short inMarkerTP);
  
  virtual su2double GetMassFlowOut(unsigned short inMarkerTP);
  
  virtual su2double GetMachIn(unsigned short inMarkerTP);
  
  virtual su2double GetMachOut(unsigned short inMarkerTP);
  
  virtual su2double GetNormalMachIn(unsigned short inMarkerTP);
  
  virtual su2double GetNormalMachOut(unsigned short inMarkerTP);
  
  virtual su2double GetEnthalpyOut(unsigned short inMarkerTP);
  
  virtual su2double GetVelocityOutIs(unsigned short inMarkerTP);
  
  virtual su2double GetPressureOut(unsigned short inMarkerTP);
  
  virtual su2double GetPressureRatio(unsigned short inMarkerTP);
  
  virtual su2double GetAveragedNormalVelocity(unsigned short valMarker);
  
  virtual su2double GetAveragedTangVelocity(unsigned short valMarker);
  
  virtual su2double GetSlidingState(unsigned short val_marker, unsigned long val_vertex, unsigned short val_state);
  
  virtual void SetSlidingState(unsigned short val_marker, unsigned long val_vertex, unsigned short val_state, su2double component);

  
  virtual void ExplicitRK_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                                    unsigned short iRKStep);
  
  virtual void ExplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  virtual void ImplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  virtual void ImplicitNewmark_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  virtual void ImplicitNewmark_Update(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  
  virtual void ImplicitNewmark_Relaxation(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  virtual void GeneralizedAlpha_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  virtual void GeneralizedAlpha_UpdateDisp(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  virtual void GeneralizedAlpha_UpdateSolution(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  virtual void GeneralizedAlpha_UpdateLoads(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  virtual void Compute_Residual(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                                unsigned short iMesh);
  
  virtual void Pressure_Forces(CGeometry *geometry, CConfig *config);
  
  virtual void TurboPerformance(CSolver *solver,  CConfig *config, unsigned short inMarker,  unsigned short outMarker, unsigned short Kind_TurboPerf, unsigned short inMarkerTP);
  
  virtual void StoreTurboPerformance(CSolver *solver, unsigned short inMarkerTP);
  
  virtual void Momentum_Forces(CGeometry *geometry, CConfig *config);
  
  virtual void Inviscid_DeltaForces(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  virtual void Friction_Forces(CGeometry *geometry, CConfig *config);
  
  virtual void Viscous_DeltaForces(CGeometry *geometry, CConfig *config);
  
  virtual void Wave_Strength(CGeometry *geometry, CConfig *config);
  
  virtual void SetPrimitive_Gradient_GG(CGeometry *geometry, CConfig *config);
  
  virtual void SetPrimitive_Gradient_LS(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief A virtual member.
  //   * \param[in] geometry - Geometrical definition of the problem.
  //   * \param[in] config - Definition of the particular problem.
  //   */
  //  virtual void SetSecondary_Gradient_GG(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief A virtual member.
  //   * \param[in] geometry - Geometrical definition of the problem.
  //   * \param[in] config - Definition of the particular problem.
  //   */
  //  virtual void SetSecondary_Gradient_LS(CGeometry *geometry, CConfig *config);
  
  virtual void Set_MPI_Primitive_Gradient(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief A virtual member.
  //   * \param[in] geometry - Geometrical definition of the problem.
  //   * \param[in] config - Definition of the particular problem.
  //   */
  //  virtual void Set_MPI_Secondary_Gradient(CGeometry *geometry, CConfig *config);
  
  virtual void SetPrimitive_Limiter_MPI(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief A virtual member.
  //   * \param[in] geometry - Geometrical definition of the problem.
  //   * \param[in] config - Definition of the particular problem.
  //   */
  //  virtual void SetSecondary_Limiter_MPI(CGeometry *geometry, CConfig *config);
  
  virtual void SetPreconditioner(CConfig *config, unsigned long iPoint);
  
  virtual void Viscous_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                                CConfig *config, unsigned short iMesh, unsigned short iRKStep);
  
  virtual void AddStiffMatrix(su2double **StiffMatrix_Elem, unsigned long Point_0, unsigned long Point_1, unsigned long Point_2, unsigned long Point_3 );
  
  virtual void Source_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CNumerics *second_numerics,
                               CConfig *config, unsigned short iMesh);
  
  virtual void Source_Template(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                               CConfig *config, unsigned short iMesh);
  
  virtual void SetCSensitivity(unsigned short val_marker, unsigned long val_vertex, su2double val_sensitivity);
  
  virtual void SetForceProj_Vector(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  virtual void SetIntBoundary_Jump(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  virtual void SetTotal_CD(su2double val_Total_CD);
  
  virtual void SetTotal_CL(su2double val_Total_CL);
  
  virtual void SetTotal_NetCThrust(su2double val_Total_NetCThrust);
  
  virtual void SetTotal_Power(su2double val_Total_Power);
  
  virtual void SetTotal_CD_SolidSurf(su2double val_Total_CD_SolidSurf);
  
  virtual void SetTotal_ReverseFlow(su2double val_ReverseFlow);
  
  virtual void SetTotal_MFR(su2double val_Total_MFR);
  
  virtual void SetTotal_Prop_Eff(su2double val_Total_Prop_Eff);
  
  virtual void SetTotal_ByPassProp_Eff(su2double val_Total_ByPassProp_Eff);
  
  virtual void SetTotal_Adiab_Eff(su2double val_Total_Adiab_Eff);
  
  virtual void SetTotal_Poly_Eff(su2double val_Total_Poly_Eff);
  
  virtual void SetTotal_IDC(su2double val_Total_IDC);
  
  virtual void SetTotal_IDC_Mach(su2double val_Total_IDC_Mach);
  
  virtual void SetTotal_IDR(su2double val_Total_IDR);
  
  virtual void SetTotal_DC60(su2double val_Total_DC60);
  
  virtual void SetTotal_Custom(su2double val_Total_Custom, su2double val_coeff);
  
  virtual void AddTotal_Custom(su2double val_Total_Custom, su2double val_coeff);
  
  virtual void SetTotal_CT(su2double val_Total_CT);
  
  virtual void SetTotal_CQ(su2double val_Total_CQ);
  
  virtual void SetTotal_HeatFlux(su2double val_Total_Heat);
  
  virtual void SetTotal_MaxHeatFlux(su2double val_Total_MaxHeat);
  
  virtual void SetDistance(CGeometry *geometry, CConfig *config);
  
  virtual void Inviscid_Sensitivity(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config);
  
  virtual void Smooth_Sensitivity(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config);
  
  virtual void Viscous_Sensitivity(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config);
  
  virtual su2double GetCL_Inv(unsigned short val_marker);
  
  virtual su2double GetCL_Visc(unsigned short val_marker);
  
  virtual su2double GetCMz_Inv(unsigned short val_marker);
  
  virtual su2double GetCMz_Visc(unsigned short val_marker);
  
  virtual su2double GetSurface_CL(unsigned short val_marker);
  
  virtual su2double GetSurface_CD(unsigned short val_marker);
  
  virtual su2double GetSurface_CSF(unsigned short val_marker);
  
  virtual su2double GetSurface_CEff(unsigned short val_marker);
  
  virtual su2double GetSurface_CFx(unsigned short val_marker);
  
  virtual su2double GetSurface_CFy(unsigned short val_marker);
  
  virtual su2double GetSurface_CFz(unsigned short val_marker);
  
  virtual su2double GetSurface_CMx(unsigned short val_marker);
  
  virtual su2double GetSurface_CMy(unsigned short val_marker);
  
  virtual su2double GetSurface_CMz(unsigned short val_marker);
  
  virtual su2double GetSurface_CL_Inv(unsigned short val_marker);
  
  virtual su2double GetSurface_CD_Inv(unsigned short val_marker);
  
  virtual su2double GetSurface_CSF_Inv(unsigned short val_marker);
  
  virtual su2double GetSurface_CEff_Inv(unsigned short val_marker);
  
  virtual su2double GetSurface_CFx_Inv(unsigned short val_marker);
  
  virtual su2double GetSurface_CFy_Inv(unsigned short val_marker);
  
  virtual su2double GetSurface_CFz_Inv(unsigned short val_marker);
  
  virtual su2double GetSurface_CMx_Inv(unsigned short val_marker);
  
  virtual su2double GetSurface_CMy_Inv(unsigned short val_marker);
  
  virtual su2double GetSurface_CMz_Inv(unsigned short val_marker);
  
  virtual su2double GetSurface_CL_Visc(unsigned short val_marker);
  
  virtual su2double GetSurface_CD_Visc(unsigned short val_marker);
  
  virtual su2double GetSurface_CSF_Visc(unsigned short val_marker);
  
  virtual su2double GetSurface_CEff_Visc(unsigned short val_marker);
  
  virtual su2double GetSurface_CFx_Visc(unsigned short val_marker);
  
  virtual su2double GetSurface_CFy_Visc(unsigned short val_marker);
  
  virtual su2double GetSurface_CFz_Visc(unsigned short val_marker);
  
  virtual su2double GetSurface_CMx_Visc(unsigned short val_marker);
  
  virtual su2double GetSurface_CMy_Visc(unsigned short val_marker);
  
  virtual su2double GetSurface_CMz_Visc(unsigned short val_marker);
  
  virtual su2double GetSurface_CL_Mnt(unsigned short val_marker);
  
  virtual su2double GetSurface_CD_Mnt(unsigned short val_marker);
  
  virtual su2double GetSurface_CSF_Mnt(unsigned short val_marker);
  
  virtual su2double GetSurface_CEff_Mnt(unsigned short val_marker);
  
  virtual su2double GetSurface_CFx_Mnt(unsigned short val_marker);
  
  virtual su2double GetSurface_CFy_Mnt(unsigned short val_marker);
  
  virtual su2double GetSurface_CFz_Mnt(unsigned short val_marker);
  
  virtual su2double GetSurface_CMx_Mnt(unsigned short val_marker);
  
  virtual su2double GetSurface_CMy_Mnt(unsigned short val_marker);
  
  virtual su2double GetSurface_CMz_Mnt(unsigned short val_marker);
  
  virtual su2double GetCSF_Visc(unsigned short val_marker);
  
  virtual su2double GetCD_Inv(unsigned short val_marker);
  
  virtual su2double GetInflow_MassFlow(unsigned short val_marker);
  
  virtual void GetSurface_Properties(CGeometry *geometry, CNumerics *conv_numerics,
                                           CNumerics *visc_numerics, CConfig *config, unsigned short iMesh, bool Output);

  virtual void GetPower_Properties(CGeometry *geometry, CConfig *config, unsigned short iMesh, bool Output);

  virtual void GetSurface_Distortion(CGeometry *geometry, CConfig *config, unsigned short iMesh, bool Output);
  
  virtual void SetFarfield_AoA(CGeometry *geometry, CSolver **solver_container,
                               CConfig *config, unsigned short iMesh, bool Output);
  
  virtual void SetActDisk_BCThrust(CGeometry *geometry, CSolver **solver_container,
                                   CConfig *config, unsigned short iMesh, bool Output);

  virtual su2double GetExhaust_MassFlow(unsigned short val_marker);
  
  virtual su2double GetInflow_Pressure(unsigned short val_marker);
  
  virtual su2double GetInflow_Mach(unsigned short val_marker);
  
  virtual su2double GetCSF_Inv(unsigned short val_marker);
  
  virtual su2double GetCEff_Inv(unsigned short val_marker);
  
  virtual su2double GetSurface_HF_Visc(unsigned short val_marker);
  
  virtual su2double GetSurface_MaxHF_Visc(unsigned short val_marker);
  
  virtual su2double GetCD_Visc(unsigned short val_marker);
  
  virtual void SetTotal_ComboObj(su2double ComboObj);
  
  virtual su2double GetTotal_ComboObj(void);
  
  virtual su2double GetTotal_CSF(void);
  
  virtual su2double GetTotal_CEff(void);
  
  virtual su2double GetTotal_CT(void);
  
  virtual su2double GetTotal_CQ(void);
  
  virtual su2double GetTotal_HeatFlux(void);
  
  virtual su2double GetTotal_MaxHeatFlux(void);
  
  virtual su2double Get_PressureDrag(void);
  
  virtual su2double Get_ViscDrag(void);
  
  virtual su2double GetTotal_CMerit(void);
  
  virtual su2double GetTotal_CEquivArea(void);
  
  virtual su2double GetTotal_AeroCD(void);
  
  virtual su2double GetTotal_RadialDistortion(void);
  
  virtual su2double GetTotal_CircumferentialDistortion(void);

  virtual su2double GetTotal_CpDiff(void);
  
  virtual su2double GetTotal_HeatFluxDiff(void);

  virtual su2double GetTotal_CFEA(void);
  
  virtual su2double GetTotal_CNearFieldOF(void);
  
  virtual void AddTotal_ComboObj(su2double val_obj);
  
  virtual void SetTotal_CEquivArea(su2double val_cequivarea);
  
  virtual void SetTotal_AeroCD(su2double val_aerocd);
  
  virtual void SetTotal_RadialDistortion(su2double val_distortion);
  
  virtual void SetTotal_CircumferentialDistortion(su2double val_distortion);

  virtual void SetTotal_CpDiff(su2double val_pressure);
  
  virtual void SetTotal_HeatFluxDiff(su2double val_heat);
  
  virtual void SetTotal_CFEA(su2double val_cfea);
  
  virtual void SetTotal_CNearFieldOF(su2double val_cnearfieldpress);
  
  virtual su2double GetTotal_CL(void);

  virtual su2double GetTotal_CD(void);
  
  virtual su2double GetTotal_NetCThrust(void);
  
  virtual su2double GetTotal_Power(void);
  
  virtual su2double GetTotal_CD_SolidSurf(void);
  
  virtual su2double GetTotal_ReverseFlow(void);
  
  virtual su2double GetTotal_MFR(void);
  
  virtual su2double GetTotal_Prop_Eff(void);
  
  virtual su2double GetTotal_ByPassProp_Eff(void);
  
  virtual su2double GetTotal_Adiab_Eff(void);
  
  virtual su2double GetTotal_Poly_Eff(void);
  
  virtual su2double GetTotal_IDC(void);
  
  virtual su2double GetTotal_IDC_Mach(void);
  
  virtual su2double GetTotal_IDR(void);
  
  virtual su2double GetTotal_DC60(void);
  
  virtual su2double GetTotal_Custom(void);

  virtual su2double GetTotal_CMx(void);
  
  virtual su2double GetTotal_CMy(void);
  
  virtual su2double GetTotal_CMz(void);
  
  virtual su2double GetTotal_CFx(void);
  
  virtual su2double GetTotal_CFy(void);
  
  virtual su2double GetTotal_CFz(void);
  
  virtual su2double GetTotal_CWave(void);
  
  virtual su2double GetTotal_CHeat(void);
  
  virtual su2double GetAllBound_CL_Inv(void);
  
  virtual su2double GetAllBound_CD_Inv(void);
  
  virtual su2double GetAllBound_CSF_Inv(void);
  
  virtual su2double GetAllBound_CEff_Inv(void);
  
  virtual su2double GetAllBound_CMx_Inv(void);
  virtual su2double GetAllBound_CMy_Inv(void);
  
  virtual su2double GetAllBound_CMz_Inv(void);
  
  virtual su2double GetAllBound_CFx_Inv(void);
  
  virtual su2double GetAllBound_CFy_Inv(void);
  
  virtual su2double GetAllBound_CFz_Inv(void);
  
  virtual su2double GetAllBound_CL_Visc(void);
  
  virtual su2double GetAllBound_CD_Visc(void);
  
  virtual su2double GetAllBound_CSF_Visc(void);
  
  virtual su2double GetAllBound_CEff_Visc(void);
  
  virtual su2double GetAllBound_CMx_Visc(void);
  
  virtual su2double GetAllBound_CMy_Visc(void);
  
  virtual su2double GetAllBound_CMz_Visc(void);
  
  virtual su2double GetAllBound_CFx_Visc(void);
  
  virtual su2double GetAllBound_CFy_Visc(void);
  
  virtual su2double GetAllBound_CFz_Visc(void);
  
  virtual su2double GetAllBound_CL_Mnt(void);
  
  virtual su2double GetAllBound_CD_Mnt(void);
  
  virtual su2double GetAllBound_CSF_Mnt(void);
  
  virtual su2double GetAllBound_CEff_Mnt(void);
  
  virtual su2double GetAllBound_CMx_Mnt(void);
  virtual su2double GetAllBound_CMy_Mnt(void);
  
  virtual su2double GetAllBound_CMz_Mnt(void);
  
  virtual su2double GetAllBound_CFx_Mnt(void);
  
  virtual su2double GetAllBound_CFy_Mnt(void);
  
  virtual su2double GetAllBound_CFz_Mnt(void);
  
  virtual su2double GetCPressure(unsigned short val_marker, unsigned long val_vertex);
  
  virtual su2double GetCPressureTarget(unsigned short val_marker, unsigned long val_vertex);
  
  virtual void SetCPressureTarget(unsigned short val_marker, unsigned long val_vertex, su2double val_pressure);
  
  virtual void SetCharacPrimVar(unsigned short val_marker, unsigned long val_vertex, unsigned short val_var, su2double val_value);

  virtual su2double *GetDonorPrimVar(unsigned short val_marker, unsigned long val_vertex);

  virtual void SetDonorPrimVar(unsigned short val_marker, unsigned long val_vertex, unsigned short val_var, su2double val_value);

  virtual void SetDonorAdjVar(unsigned short val_marker, unsigned long val_vertex, unsigned short val_var, su2double val_value);
  
  virtual su2double GetDonorPrimVar(unsigned short val_marker, unsigned long val_vertex, unsigned short val_var);
  
  virtual su2double *GetDonorAdjVar(unsigned short val_marker, unsigned long val_vertex);
  
  virtual su2double GetDonorAdjVar(unsigned short val_marker, unsigned long val_vertex, unsigned short val_var);
  
  virtual unsigned long GetDonorGlobalIndex(unsigned short val_marker, unsigned long val_vertex);
  
  virtual void SetDonorGlobalIndex(unsigned short val_marker, unsigned long val_vertex, unsigned long val_index);
  
  virtual su2double *GetCharacPrimVar(unsigned short val_marker, unsigned long val_vertex);
  
  virtual su2double GetActDisk_DeltaP(unsigned short val_marker, unsigned long val_vertex);
  
  virtual void SetActDisk_DeltaP(unsigned short val_marker, unsigned long val_vertex, su2double val_deltap);
  
  virtual su2double GetActDisk_DeltaT(unsigned short val_marker, unsigned long val_vertex);
  
  virtual void SetActDisk_DeltaT(unsigned short val_marker, unsigned long val_vertex, su2double val_deltat);
  
  virtual su2double GetInlet_Ttotal(unsigned short val_marker, unsigned long val_vertex);
  
  virtual su2double GetInlet_Ptotal(unsigned short val_marker, unsigned long val_vertex);
  
  virtual su2double GetInlet_FlowDir(unsigned short val_marker, unsigned long val_vertex, unsigned short val_dim);
  
  virtual void SetInlet_Ttotal(unsigned short val_marker, unsigned long val_vertex, su2double val_ttotal);
  
  virtual void SetInlet_Ptotal(unsigned short val_marker, unsigned long val_vertex, su2double val_ptotal);
  
  virtual void SetInlet_FlowDir(unsigned short val_marker, unsigned long val_vertex, unsigned short val_dim, su2double val_flowdir);

  virtual su2double GetCSkinFriction(unsigned short val_marker, unsigned long val_vertex, unsigned short val_dim);
  
  virtual su2double GetHeatFlux(unsigned short val_marker, unsigned long val_vertex);
  
  virtual su2double GetHeatFluxTarget(unsigned short val_marker, unsigned long val_vertex);
  
  virtual void SetHeatFluxTarget(unsigned short val_marker, unsigned long val_vertex, su2double val_heat);
  
  virtual su2double GetYPlus(unsigned short val_marker, unsigned long val_vertex);
  
  virtual su2double GetStrainMag_Max(void);
  
  virtual su2double GetOmega_Max(void);
  
  virtual void SetStrainMag_Max(su2double val_strainmag_max);
  
  virtual void SetOmega_Max(su2double val_omega_max);
  
  virtual su2double GetPsiRho_Inf(void);
  
  virtual su2double* GetPsiRhos_Inf(void);
  
  virtual su2double GetPsiE_Inf(void);
  
  virtual su2double GetPhi_Inf(unsigned short val_dim);
  
  virtual su2double GetTotal_Sens_Geo(void);
  
  virtual su2double GetTotal_Sens_Mach(void);
  
  virtual su2double GetTotal_Sens_AoA(void);
  
  virtual su2double GetTotal_Sens_Press(void);
  
  virtual su2double GetTotal_Sens_Temp(void);
  
  virtual su2double GetTotal_Sens_BPress(void);
  
  virtual su2double GetDensity_Inf(void);
  
  virtual su2double GetDensity_Inf(unsigned short val_var);
  
  virtual su2double GetModVelocity_Inf(void);
  
  virtual su2double GetDensity_Energy_Inf(void);
  
  virtual su2double GetPressure_Inf(void);
  
  virtual su2double GetDensity_Velocity_Inf(unsigned short val_dim);
  
  virtual su2double GetDensity_Velocity_Inf(unsigned short val_dim, unsigned short val_var);
  
  virtual su2double GetVelocity_Inf(unsigned short val_dim);
  
  virtual su2double *GetVelocity_Inf(void);
  
  virtual su2double GetViscosity_Inf(void);
  
  virtual su2double GetTke_Inf(void);
  
  virtual su2double GetCSensitivity(unsigned short val_marker, unsigned long val_vertex);
 
  virtual su2double* GetConstants();
  
  virtual su2double GetOneD_TotalPress(void);
  
  virtual void SetOneD_TotalPress(su2double AveragePressure);
  
  virtual su2double GetOneD_Mach(void);
  
  virtual void SetOneD_Mach(su2double AverageMach);
  
  virtual su2double GetOneD_Temp(void);
  
  virtual void SetOneD_Temp(su2double AverageTemperature);
  
  virtual su2double GetOneD_MassFlowRate(void);
  
  virtual void SetOneD_MassFlowRate(su2double MassFlowRate);
  
  virtual su2double GetOneD_FluxAvgPress(void);
  
  virtual void SetOneD_FluxAvgPress(su2double PressureRef);
  virtual su2double GetOneD_FluxAvgDensity(void);
  
  virtual void SetOneD_FluxAvgDensity(su2double DensityRef);
  
  virtual su2double GetOneD_FluxAvgVelocity(void);
  
  virtual void SetOneD_FluxAvgVelocity(su2double VelocityRef);
  
  virtual su2double GetOneD_FluxAvgEntalpy(void);
  virtual void SetOneD_FluxAvgEntalpy(su2double EnthalpyRef);

  virtual void SetFEA_Load(CSolver ***flow_solution, CGeometry **fea_geometry,
                           CGeometry **flow_geometry, CConfig *fea_config,
                           CConfig *flow_config, CNumerics *fea_numerics);
  
  virtual void SetFEA_Load_Int(CSolver ***flow_solution, CGeometry **fea_geometry,
                               CGeometry **flow_geometry, CConfig *fea_config,
                               CConfig *flow_config, CNumerics *fea_numerics);
  
  virtual void Copy_Zone_Solution(CSolver ***solver1_solution,
                                  CGeometry **solver1_geometry,
                                  CConfig *solver1_config,
                                  CSolver ***solver2_solution,
                                  CGeometry **solver2_geometry,
                                  CConfig *solver2_config);
  
  virtual void SetInitialCondition(CGeometry **geometry,
                                   CSolver ***solver_container,
                                   CConfig *config, unsigned long ExtIter);
  
  virtual void ResetInitialCondition(CGeometry **geometry,
                                     CSolver ***solver_container,
                                     CConfig *config, unsigned long ExtIter);
  
  virtual void SetFlow_Displacement(CGeometry **flow_geometry,
                                    CVolumetricMovement *flow_grid_movement,
                                    CConfig *flow_config, CConfig *fea_config,
                                    CGeometry **fea_geometry,
                                    CSolver ***fea_solution);
  
  virtual void SetFlow_Displacement_Int(CGeometry **flow_geometry,
                                        CVolumetricMovement *flow_grid_movement,
                                        CConfig *flow_config, CConfig *fea_config,
                                        CGeometry **fea_geometry,
                                        CSolver ***fea_solution);
  
  virtual void PredictStruct_Displacement(CGeometry **fea_geometry,
                                          CConfig *fea_config,
                                          CSolver ***fea_solution);
  
  virtual void ComputeAitken_Coefficient(CGeometry **fea_geometry,
                                         CConfig *fea_config,
                                         CSolver ***fea_solution,
                                         unsigned long iFSIIter);
  
  
  virtual void SetAitken_Relaxation(CGeometry **fea_geometry,
                                    CConfig *fea_config,
                                    CSolver ***fea_solution);
  
  virtual void Update_StructSolution(CGeometry **fea_geometry,
                                     CConfig *fea_config,
                                     CSolver ***fea_solution);
  
  virtual void LoadRestart(CGeometry **geometry, CSolver ***solver,
                           CConfig *config, int val_iter);
  
  virtual void LoadRestart_FSI(CGeometry *geometry, CSolver ***solver,
                               CConfig *config, int val_iter);
  
  void Gauss_Elimination(su2double** A, su2double* rhs, unsigned short nVar);
  
  void Aeroelastic(CSurfaceMovement *surface_movement, CGeometry *geometry, CConfig *config, unsigned long ExtIter);
  
  void SetUpTypicalSectionWingModel(vector<vector<su2double> >& PHI, vector<su2double>& w, CConfig *config);
  
  void SolveTypicalSectionWingModel(CGeometry *geometry, su2double Cl, su2double Cm, CConfig *config, unsigned short val_Marker, vector<su2double>& displacements);
  
  virtual void RegisterSolution(CGeometry *geometry, CConfig *config);
  
  virtual void RegisterOutput(CGeometry *geometry, CConfig *config);
  
  virtual void SetAdjoint_Output(CGeometry *geometry, CConfig *config);
  
  virtual void ExtractAdjoint_Solution(CGeometry *geometry,  CConfig *config);
  
  virtual void RegisterObj_Func(CConfig *config);
  
  virtual void SetSurface_Sensitivity(CGeometry *geometry, CConfig* config);
  
  virtual void SetSensitivity(CGeometry *geometry, CConfig *config);
  
  virtual void SetAdj_ObjFunc(CGeometry *geometry, CConfig* config);
  
  virtual void SetFSI_ConvValue(unsigned short val_index, su2double val_criteria);
  
  virtual su2double GetFSI_ConvValue(unsigned short val_index);
  
  virtual void Set_Prestretch(CGeometry *geometry, CConfig *config);
  
  virtual void Compute_StiffMatrix(CGeometry *geometry, CSolver **solver_container, CNumerics **numerics, CConfig *config);
  
  virtual void Compute_StiffMatrix_NodalStressRes(CGeometry *geometry, CSolver **solver_container, CNumerics **numerics, CConfig *config);
  
  
  virtual void Compute_MassMatrix(CGeometry *geometry, CSolver **solver_container, CNumerics **numerics, CConfig *config);
  
  virtual void Compute_NodalStressRes(CGeometry *geometry, CSolver **solver_container, CNumerics **numerics, CConfig *config);
  
  virtual void Compute_NodalStress(CGeometry *geometry, CSolver **solver_container, CNumerics **numerics, CConfig *config);
  
  virtual void Compute_DeadLoad(CGeometry *geometry, CSolver **solver_container, CNumerics **numerics, CConfig *config);
  
  virtual void Solve_System(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  
  virtual void Initialize_SystemMatrix(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  virtual void Compute_IntegrationConstants(CConfig *config);
  
  virtual su2double GetWAitken_Dyn(void);
  
  virtual su2double GetWAitken_Dyn_tn1(void);
  
  virtual void SetWAitken_Dyn(su2double waitk);
  
  virtual void SetWAitken_Dyn_tn1(su2double waitk_tn1);
  
  virtual void SetLoad_Increment(su2double val_loadIncrement);
  
  virtual unsigned long SetPrimitive_Variables(CSolver **solver_container, CConfig *config, bool Output);
  
  virtual void SetPressure_Inf(su2double p_inf);
  
  virtual void SetTemperature_Inf(su2double t_inf);
  
  virtual void SetRecording(CGeometry *geometry, CConfig *config, unsigned short kind_recording);
  
  virtual void RegisterVariables(CGeometry *geometry, CConfig *config, bool reset = false);
  
  virtual void ExtractAdjoint_Variables(CGeometry *geometry, CConfig *config);
  
  virtual void SetFreeStream_Solution(CConfig *config);
};

class CBaselineSolver : public CSolver {
public:
  
  CBaselineSolver(void);
  
  CBaselineSolver(CGeometry *geometry, CConfig *config, unsigned short nVar, vector<string> field_names);
  
  CBaselineSolver(CGeometry *geometry, CConfig *config, unsigned short iMesh);
  
  void Set_MPI_Solution(CGeometry *geometry, CConfig *config);
  
  void LoadRestart(CGeometry **geometry, CSolver ***solver, CConfig *config, int val_iter);
  
  void LoadRestart_FSI(CGeometry *geometry, CSolver ***solver, CConfig *config, int val_iter);
  
  virtual ~CBaselineSolver(void);
  
};

class CEulerSolver : public CSolver {
protected:
  
  su2double
  Mach_Inf,  
  Density_Inf,  
  Energy_Inf,      
  Temperature_Inf,      
  Pressure_Inf,    
  *Velocity_Inf;    
  su2double
  *CD_Inv,  
  *CL_Inv,      
  *CSF_Inv,    
  *CMx_Inv,      
  *CMy_Inv,      
  *CMz_Inv,      
  *CFx_Inv,      
  *CFy_Inv,      
  *CFz_Inv,      
  *Surface_CL_Inv, 
  *Surface_CD_Inv, 
  *Surface_CSF_Inv, 
  *Surface_CEff_Inv, 
  *Surface_CFx_Inv,   
  *Surface_CFy_Inv,   
  *Surface_CFz_Inv,   
  *Surface_CMx_Inv,   
  *Surface_CMy_Inv,   
  *Surface_CMz_Inv,   
  *CEff_Inv,        
  *CMerit_Inv,        
  *CT_Inv,      
  *CQ_Inv,      
  *CEquivArea_Inv,        
  *CNearFieldOF_Inv,        
  *CD_Mnt,  
  *CL_Mnt,      
  *CSF_Mnt,    
  *CMx_Mnt,      
  *CMy_Mnt,      
  *CMz_Mnt,      
  *CFx_Mnt,      
  *CFy_Mnt,      
  *CFz_Mnt,      
  *Surface_CL_Mnt, 
  *Surface_CD_Mnt, 
  *Surface_CSF_Mnt, 
  *Surface_CEff_Mnt, 
  *Surface_CFx_Mnt,   
  *Surface_CFy_Mnt,   
  *Surface_CFz_Mnt,   
  *Surface_CMx_Mnt,   
  *Surface_CMy_Mnt,   
  *Surface_CMz_Mnt,   
  *CEff_Mnt,        
  *CMerit_Mnt,        
  *CT_Mnt,      
  *CQ_Mnt,      
  *CEquivArea_Mnt,        
  **CPressure,    
  **CPressureTarget,    
  **HeatFlux,    
  **HeatFluxTarget,    
  **YPlus,    
  ***CharacPrimVar,    
  ***DonorPrimVar,    
  *ForceInviscid,    
  *MomentInviscid,  
  *ForceMomentum,    
  *MomentMomentum;  
  su2double *Inflow_MassFlow,  
  *Exhaust_MassFlow,  
  *Inflow_Pressure,  
  *Inflow_Mach,  
  *Inflow_Area,  
  *Exhaust_Area,  
  *Exhaust_Pressure,  
  *Exhaust_Temperature,  
  Inflow_MassFlow_Total,  
  Exhaust_MassFlow_Total,  
  Inflow_Pressure_Total,  
  Inflow_Mach_Total,  
  InverseDesign;  
  unsigned long **DonorGlobalIndex;    
  su2double **ActDisk_DeltaP,    
  **ActDisk_DeltaT;    
  su2double **Inlet_Ptotal,    
  **Inlet_Ttotal,    
  ***Inlet_FlowDir;    
  su2double
  AllBound_CD_Inv,  
  AllBound_CL_Inv,      
  AllBound_CSF_Inv,      
  AllBound_CMx_Inv,      
  AllBound_CMy_Inv,      
  AllBound_CMz_Inv,      
  AllBound_CFx_Inv,      
  AllBound_CFy_Inv,      
  AllBound_CFz_Inv,      
  AllBound_CEff_Inv,      
  AllBound_CMerit_Inv,      
  AllBound_CT_Inv,      
  AllBound_CQ_Inv,      
  AllBound_CEquivArea_Inv,      
  AllBound_CNearFieldOF_Inv;      
  su2double
  AllBound_CD_Mnt,  
   AllBound_CL_Mnt,      
   AllBound_CSF_Mnt,      
   AllBound_CMx_Mnt,      
   AllBound_CMy_Mnt,      
   AllBound_CMz_Mnt,      
   AllBound_CFx_Mnt,      
   AllBound_CFy_Mnt,      
   AllBound_CFz_Mnt,      
   AllBound_CEff_Mnt,      
   AllBound_CMerit_Mnt,      
   AllBound_CT_Mnt,      
   AllBound_CQ_Mnt;      
  su2double
  OneD_TotalPress, 
  OneD_Mach, 
  OneD_Temp, 
  OneD_PressureRef, 
  OneD_MassFlowRate, 
  OneD_DensityRef, 
  OneD_EnthalpyRef, 
  OneD_VelocityRef, 
  Total_ComboObj, 
  AoA_Prev, 
  Total_CD, 
  Total_CL,    
  Total_CL_Prev,    
  Total_CD_SolidSurf, 
  Total_CD_Prev, 
  Total_NetCThrust, 
  Total_Power, 
  Total_ReverseFlow, 
  Total_IDC,        
  Total_IDC_Mach,        
  Total_IDR,        
  Total_DC60,        
  Total_MFR,     
  Total_Prop_Eff,     
  Total_ByPassProp_Eff,     
  Total_Adiab_Eff,     
  Total_Poly_Eff,     
  Total_NetCThrust_Prev,    
  Total_BCThrust_Prev,    
  Total_Custom,        
  Total_CSF,    
  Total_CMx,      
  Total_CMy,      
  Total_CMz,      
  Total_CFx,      
  Total_CFy,      
  Total_CFz,      
  Total_CEff,      
  Total_CMerit,      
  Total_CT,    
  Total_CQ,    
  Total_Heat,    
  Total_MaxHeat, 
  Total_AeroCD,      
  Total_RadialDistortion,      
  Total_CircumferentialDistortion,      
  Total_CEquivArea,      
  Total_CNearFieldOF,      
  Total_CpDiff,      
  Total_HeatFluxDiff,      
  Total_MassFlowRate;     
  su2double *Surface_CL,   
  *Surface_CD,          
  *Surface_CSF,     
  *Surface_CEff,     
  *Surface_CFx,            
  *Surface_CFy,            
  *Surface_CFz,            
  *Surface_CMx,            
  *Surface_CMy,            
  *Surface_CMz,            
  *Surface_HF_Visc,            
  *Surface_MaxHF_Visc;         
  su2double *iPoint_UndLapl,  
  *jPoint_UndLapl;      
  su2double *SecondaryVar_i,  
  *SecondaryVar_j;      
  su2double *PrimVar_i,  
  *PrimVar_j;      
  su2double **LowMach_Precontioner; 
  unsigned long nMarker,        
  *nVertex;       
  bool space_centered,  
  euler_implicit,      
  least_squares;        
  su2double Gamma;                  
  su2double Gamma_Minus_One;        
  su2double *Primitive,    
  *Primitive_i,        
  *Primitive_j;        
  su2double *Secondary,    
  *Secondary_i,        
  *Secondary_j;        
  su2double Cauchy_Value,  
  Cauchy_Func;      
  unsigned short Cauchy_Counter;  
  su2double *Cauchy_Serie;      
  su2double Old_Func,  
  New_Func;      
  su2double AoA_old;  
  unsigned long AoA_Counter;
  unsigned long BCThrust_Counter;

  CFluidModel  *FluidModel;  
  su2double **AveragedVelocity,
            **AveragedNormal,
            **AveragedGridVel,
            **AveragedFlux,
            **TotalFlux,
            *TotalArea,
            *AveragedNormalVelocity,
            *ExtAveragedNormalVelocity,
            *AveragedTangVelocity,
            *ExtAveragedTangVelocity,
            *AveragedTangGridVelocity,
            *AveragedMach,
            *AveragedNormalMach,
            *AveragedTangMach,
            *AveragedEnthalpy,
            *AveragedPressure,
            *AveragedTotTemperature,
            *AveragedTotPressure,
            *ExtAveragedPressure,
            *ExtAveragedTotTemperature,
            *ExtAveragedTotPressure,
            *AveragedDensity,
            *ExtAveragedDensity,
            *AveragedSoundSpeed,
            *AveragedEntropy,
            *MassFlow,
            *FlowAngle;
  
  su2double *TotalStaticEfficiency,
            *TotalTotalEfficiency,
            *KineticEnergyLoss,
            *TotalPressureLoss,
            *MassFlowIn,
            *MassFlowOut,
            *FlowAngleIn,
            *FlowAngleOut,
            *EulerianWork,
            *TotalEnthalpyIn,
            *PressureRatio,
            *PressureOut,
            *EnthalpyOut,
            *MachIn,
            *MachOut,
            *NormalMachIn,
            *NormalMachOut,
            *VelocityOutIs;

  /* Sliding meshes variables */

  su2double ***SlidingState;

public:
  
  
  CEulerSolver(void);
  
  CEulerSolver(CGeometry *geometry, CConfig *config, unsigned short iMesh);
  
  virtual ~CEulerSolver(void);
  
  void Set_MPI_Solution(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Gradient(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Old(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Limiter(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Primitive_Limiter(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief Impose the send-receive boundary condition.
  //   * \param[in] geometry - Geometrical definition of the problem.
  //   * \param[in] config - Definition of the particular problem.
  //   */
  //  void Set_MPI_Secondary_Limiter(CGeometry *geometry, CConfig *config);
  
  void SetNondimensionalization(CGeometry *geometry, CConfig *config, unsigned short iMesh);
  
  CFluidModel* GetFluidModel(void);
  
  
  su2double GetDensity_Inf(void);
  
  su2double GetModVelocity_Inf(void);
  
  su2double GetDensity_Energy_Inf(void);
  
  su2double GetPressure_Inf(void);
  
  su2double GetDensity_Velocity_Inf(unsigned short val_dim);
  
  su2double GetVelocity_Inf(unsigned short val_dim);
  
  su2double *GetVelocity_Inf(void);
  
  void SetTime_Step(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                    unsigned short iMesh, unsigned long Iteration);
  
  void Centered_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                         CConfig *config, unsigned short iMesh, unsigned short iRKStep);
  
  void Upwind_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                       CConfig *config, unsigned short iMesh);
  
  void ComputeConsExtrapolation(CConfig *config);
  void Source_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CNumerics *second_numerics,
                       CConfig *config, unsigned short iMesh);
  
  void Source_Template(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                       CConfig *config, unsigned short iMesh);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  void Postprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh);
  
  unsigned long SetPrimitive_Variables(CSolver **solver_container, CConfig *config, bool Output);
  
  void SetDissipation_Switch(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Dissipation_Switch(CGeometry *geometry, CConfig *config);
  
  void SetPrimitive_Gradient_GG(CGeometry *geometry, CConfig *config);
  
  void SetPrimitive_Gradient_LS(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief Compute the gradient of the primitive variables using Green-Gauss method,
  //   *        and stores the result in the <i>Gradient_Primitive</i> variable.
  //   * \param[in] geometry - Geometrical definition of the problem.
  //   * \param[in] config - Definition of the particular problem.
  //   */
  //  void SetSecondary_Gradient_GG(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief Compute the gradient of the primitive variables using a Least-Squares method,
  //   *        and stores the result in the <i>Gradient_Primitive</i> variable.
  //   * \param[in] geometry - Geometrical definition of the problem.
  //   * \param[in] config - Definition of the particular problem.
  //   */
  //  void SetSecondary_Gradient_LS(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Primitive_Gradient(CGeometry *geometry, CConfig *config);
  
  void SetPrimitive_Limiter(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief Compute the gradient of the primitive variables using a Least-Squares method,
  //   *        and stores the result in the <i>Gradient_Primitive</i> variable.
  //   * \param[in] geometry - Geometrical definition of the problem.
  //   * \param[in] config - Definition of the particular problem.
  //   */
  //  void Set_MPI_Secondary_Gradient(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief Compute the limiter of the primitive variables.
  //   * \param[in] geometry - Geometrical definition of the problem.
  //   * \param[in] config - Definition of the particular problem.
  //   */
  //  void SetSecondary_Limiter(CGeometry *geometry, CConfig *config);
  
  void SetPreconditioner(CConfig *config, unsigned long iPoint);
  
  void SetUndivided_Laplacian(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Undivided_Laplacian(CGeometry *geometry, CConfig *config);
  
  void SetMax_Eigenvalue(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_MaxEigenvalue(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_ActDisk(CSolver **solver_container, CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Nearfield(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Interface(CGeometry *geometry, CConfig *config);

  void Compute_ComboObj(CConfig *config);
  
  void BC_Euler_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                     unsigned short val_marker);
  
  void BC_Far_Field(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                    CConfig *config, unsigned short val_marker);
  
  void BC_Sym_Plane(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  void BC_Fluid_Interface(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config);
    
  void BC_ActDisk_Inlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                        CConfig *config, unsigned short val_marker);
  
  void BC_ActDisk_Outlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                         CConfig *config, unsigned short val_marker);
  
  void BC_ActDisk(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                  CConfig *config, unsigned short val_marker, bool inlet_surface);
  
  void BC_Interface_Boundary(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                             CConfig *config, unsigned short val_marker);
  
  void BC_NearField_Boundary(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                             CConfig *config, unsigned short val_marker);
  
  void BC_Dirichlet(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                    unsigned short val_marker);
  
  void BC_Riemann(CGeometry *geometry, CSolver **solver_container,
                  CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  void BC_NonReflecting(CGeometry *geometry, CSolver **solver_container,
                        CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  
  void BC_Inlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                CConfig *config, unsigned short val_marker);
  
  void BC_Supersonic_Inlet(CGeometry *geometry, CSolver **solver_container,
                           CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  void BC_Supersonic_Outlet(CGeometry *geometry, CSolver **solver_container,
                            CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  void BC_Custom(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                 CConfig *config, unsigned short val_marker);
  
  void BC_Outlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                 CConfig *config, unsigned short val_marker);
  
  void BC_Engine_Inflow(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                        CConfig *config, unsigned short val_marker);
  
  void BC_Engine_Exhaust(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                         CConfig *config, unsigned short val_marker);
  
  void Mixing_Process(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short val_Marker);
  
  void MixedOut_Average (su2double val_init_pressure, su2double *val_Averaged_Flux, su2double *val_normal, su2double *pressure_mix, su2double *density_mix);
  
  void MixedOut_Root_Function(su2double *pressure, su2double *val_Averaged_Flux, su2double *val_normal, su2double *valfunc, su2double *density);
  
  void Boundary_Fourier(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short val_Marker, vector<std::complex<su2double> > &c4k,signed long &nboundaryvertex);
  
  void Boundary_Fourier(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short val_Marker, vector<std::complex<su2double> >& c2k,vector<std::complex<su2double> >& c3k,signed long& nboundaryvertex);
  
  void SetExtAveragedValue(CSolver *solver_container, unsigned short intMarker,  unsigned short extMarker);
  
  su2double GetAveragedDensity(unsigned short valMarker);
  
  su2double GetAveragedTotPressure(unsigned short valMarker);
  
  su2double GetTotalPressureLoss(unsigned short inMarkerTP);
  
  su2double GetKineticEnergyLoss(unsigned short inMarkerTP);
  
  su2double GetTotalTotalEfficiency(unsigned short inMarkerTP);
  
  su2double GetTotalStaticEfficiency(unsigned short inMarkerTP);
  
  su2double GetEulerianWork(unsigned short inMarkerTP);
  
  su2double GetTotalEnthalpyIn(unsigned short inMarkerTP);
  
  su2double GetFlowAngleIn(unsigned short inMarkerTP);
  
  su2double GetFlowAngleOut(unsigned short inMarkerTP);
  
  su2double GetMassFlowIn(unsigned short inMarkerTP);
  
  su2double GetMassFlowOut(unsigned short inMarkerTP);
  
  su2double GetMachIn(unsigned short inMarkerTP);
  
  su2double GetMachOut(unsigned short inMarkerTP);
  
  su2double GetNormalMachIn(unsigned short inMarkerTP);
  
  su2double GetNormalMachOut(unsigned short inMarkerTP);
  
  su2double GetEnthalpyOut(unsigned short inMarkerTP);
  
  su2double GetVelocityOutIs(unsigned short inMarkerTP);
  
  su2double GetPressureOut(unsigned short inMarkerTP);
  
  su2double GetPressureRatio(unsigned short inMarkerTP);
  
  su2double GetAveragedTotTemperature(unsigned short valMarker);
  
  su2double GetAveragedPressure(unsigned short valMarker);
  
  su2double GetMassFlow(unsigned short valMarker);
  
  su2double GetFlowAngle(unsigned short valMarker);
  
  su2double GetAveragedMach(unsigned short valMarker);
  
  su2double GetAveragedNormalMach(unsigned short valMarker);
  
  su2double GetAveragedEnthalpy(unsigned short valMarker);
  
  su2double GetAveragedEntropy(unsigned short valMarker);
  
  su2double* GetAveragedVelocity(unsigned short valMarker);
  
  su2double* GetAveragedGridVelocity(unsigned short valMarker);
  
  su2double GetAveragedNormalVelocity(unsigned short valMarker);
  
  su2double GetAveragedTangVelocity(unsigned short valMarker);
  
  static bool Compareval(std::vector<su2double> a,std::vector<su2double> b);
  
  void ExplicitRK_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                            unsigned short iRKStep);
  
  void GetSurface_Properties(CGeometry *geometry, CNumerics *conv_numerics,
                                   CNumerics *visc_numerics, CConfig *config, unsigned short iMesh, bool Output);

  void GetPower_Properties(CGeometry *geometry, CConfig *config, unsigned short iMesh, bool Output);
  
  void GetSurface_Distortion(CGeometry *geometry, CConfig *config, unsigned short iMesh, bool Output);
  
  void SetActDisk_BCThrust(CGeometry *geometry, CSolver **solver_container,
                           CConfig *config, unsigned short iMesh, bool Output);
  
  void SetFarfield_AoA(CGeometry *geometry, CSolver **solver_container,
                       CConfig *config, unsigned short iMesh, bool Output);
  
  void ExplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void ImplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void Pressure_Forces(CGeometry *geometry, CConfig *config);
  
  void Momentum_Forces(CGeometry *geometry, CConfig *config);
  
  void TurboPerformance(CSolver *solver,  CConfig *config, unsigned short inMarker,  unsigned short outMarker, unsigned short Kind_TurboPerf , unsigned short inMarkerTP );
  
  void StoreTurboPerformance(CSolver *solver,  unsigned short inMarkerTP );
  
  su2double GetSlidingState(unsigned short val_marker, unsigned long val_vertex, unsigned short val_state);

  su2double GetCL_Inv(unsigned short val_marker);
  
  su2double GetCMz_Inv(unsigned short val_marker);
  
  su2double GetSurface_CL(unsigned short val_marker);
  
  su2double GetSurface_CD(unsigned short val_marker);
  
  su2double GetSurface_CSF(unsigned short val_marker);
  
  su2double GetSurface_CEff(unsigned short val_marker);
  
  su2double GetSurface_CFx(unsigned short val_marker);
  
  su2double GetSurface_CFy(unsigned short val_marker);
  
  su2double GetSurface_CFz(unsigned short val_marker);
  
  su2double GetSurface_CMx(unsigned short val_marker);
  
  su2double GetSurface_CMy(unsigned short val_marker);
  
  su2double GetSurface_CMz(unsigned short val_marker);
  
  su2double GetSurface_CL_Inv(unsigned short val_marker);
  
  su2double GetSurface_CD_Inv(unsigned short val_marker);
  
  su2double GetSurface_CSF_Inv(unsigned short val_marker);
  
  su2double GetSurface_CEff_Inv(unsigned short val_marker);
  
  su2double GetSurface_CFx_Inv(unsigned short val_marker);
  
  su2double GetSurface_CFy_Inv(unsigned short val_marker);
  
  su2double GetSurface_CFz_Inv(unsigned short val_marker);
  
  su2double GetSurface_CMx_Inv(unsigned short val_marker);
  
  su2double GetSurface_CMy_Inv(unsigned short val_marker);
  
  su2double GetSurface_CMz_Inv(unsigned short val_marker);
  
  su2double GetSurface_CL_Mnt(unsigned short val_marker);
  
  su2double GetSurface_CD_Mnt(unsigned short val_marker);
  
  su2double GetSurface_CSF_Mnt(unsigned short val_marker);
  
  su2double GetSurface_CEff_Mnt(unsigned short val_marker);
  
  su2double GetSurface_CFx_Mnt(unsigned short val_marker);
  
  su2double GetSurface_CFy_Mnt(unsigned short val_marker);
  
  su2double GetSurface_CFz_Mnt(unsigned short val_marker);
  
  su2double GetSurface_CMx_Mnt(unsigned short val_marker);
  
  su2double GetSurface_CMy_Mnt(unsigned short val_marker);
  
  su2double GetSurface_CMz_Mnt(unsigned short val_marker);
  
  su2double GetCD_Inv(unsigned short val_marker);
  
  su2double GetInflow_MassFlow(unsigned short val_marker);
  
  su2double GetExhaust_MassFlow(unsigned short val_marker);
  
  su2double GetInflow_Pressure(unsigned short val_marker);
  
  su2double GetInflow_Mach(unsigned short val_marker);
  
  su2double GetCSF_Inv(unsigned short val_marker);
  
  su2double GetCEff_Inv(unsigned short val_marker);
  
  su2double GetTotal_CSF(void);
  
  su2double GetTotal_CEff(void);
  
  su2double GetTotal_CEquivArea(void);
  
  su2double GetTotal_AeroCD(void);

  su2double GetTotal_RadialDistortion(void);

  su2double GetTotal_CircumferentialDistortion(void);

  su2double GetTotal_CpDiff(void);
  
  su2double GetTotal_HeatFluxDiff(void);
  
  su2double GetTotal_CNearFieldOF(void);
  
  void AddTotal_ComboObj(su2double val_obj);
  
  void SetTotal_CEquivArea(su2double val_cequivarea);
  
  void SetTotal_AeroCD(su2double val_aerocd);
  
  void SetTotal_RadialDistortion(su2double val_distortion);

  void SetTotal_CircumferentialDistortion(su2double val_distortion);
  
  void SetTotal_CpDiff(su2double val_pressure);
  
  void SetTotal_HeatFluxDiff(su2double val_heat);
  
  void SetTotal_CNearFieldOF(su2double val_cnearfieldpress);
  
  void SetTotal_ComboObj(su2double ComboObj);
  
  su2double GetTotal_ComboObj(void);
  
  su2double GetTotal_CL(void);

  su2double GetTotal_CD(void);
  
  su2double GetTotal_NetCThrust(void);
  
  su2double GetTotal_Power(void);
  
  su2double GetTotal_CD_SolidSurf(void);
  
  su2double GetTotal_ReverseFlow(void);
  
  su2double GetTotal_MFR(void);
  
  su2double GetTotal_Prop_Eff(void);
  
  su2double GetTotal_ByPassProp_Eff(void);
  
  su2double GetTotal_Adiab_Eff(void);
  
  su2double GetTotal_Poly_Eff(void);
  
  su2double GetTotal_IDC(void);
  
  su2double GetTotal_IDC_Mach(void);
  
  su2double GetTotal_IDR(void);
  
  su2double GetTotal_DC60(void);
  
  su2double GetTotal_Custom(void);

  su2double GetTotal_CMx(void);
  
  su2double GetTotal_CMy(void);
  
  su2double GetTotal_CMz(void);
  
  su2double GetTotal_CFx(void);
  
  su2double GetTotal_CFy(void);
  
  su2double GetTotal_CFz(void);
  
  su2double GetTotal_CT(void);
  
  void SetTotal_CT(su2double val_Total_CT);
  
  su2double GetTotal_CQ(void);
  
  su2double GetTotal_HeatFlux(void);
  
  su2double GetTotal_MaxHeatFlux(void);
  
  void SetTotal_CQ(su2double val_Total_CQ);
  
  void SetTotal_HeatFlux(su2double val_Total_Heat);
  
  void SetTotal_MaxHeatFlux(su2double val_Total_MaxHeat);
  
  su2double GetTotal_CMerit(void);
  
  void SetTotal_CD(su2double val_Total_CD);
  
  void SetTotal_CL(su2double val_Total_CL);

  void SetTotal_NetCThrust(su2double val_Total_NetCThrust);
  
  void SetTotal_Power(su2double val_Total_Power);
  
  void SetTotal_CD_SolidSurf(su2double val_Total_CD_SolidSurf);
  
  void SetTotal_ReverseFlow(su2double val_ReverseFlow);
  
  void SetTotal_MFR(su2double val_Total_MFR);
  
  void SetTotal_Prop_Eff(su2double val_Total_Prop_Eff);
  
  void SetTotal_ByPassProp_Eff(su2double val_Total_ByPassProp_Eff);
  
  void SetTotal_Adiab_Eff(su2double val_Total_Adiab_Eff);
  
  void SetTotal_Poly_Eff(su2double val_Total_Poly_Eff);
  
  void SetTotal_IDC(su2double val_Total_IDC);
  
  void SetTotal_IDC_Mach(su2double val_Total_IDC_Mach);
  
  void SetTotal_IDR(su2double val_Total_IDR);
  
  void SetTotal_DC60(su2double val_Total_DC60);
  
  void SetTotal_Custom(su2double val_Total_Custom, su2double val_coeff);
  
  void AddTotal_Custom(su2double val_Total_Custom, su2double val_coeff);

  su2double GetAllBound_CL_Inv(void);
  
  su2double GetAllBound_CD_Inv(void);
  
  su2double GetAllBound_CSF_Inv(void);
  
  su2double GetAllBound_CEff_Inv(void);
  
  su2double GetAllBound_CMx_Inv(void);
  
  su2double GetAllBound_CMy_Inv(void);
  
  su2double GetAllBound_CMz_Inv(void);
  
  su2double GetAllBound_CFx_Inv(void);
  
  su2double GetAllBound_CFy_Inv(void);
  
  su2double GetAllBound_CFz_Inv(void);
  
  su2double GetAllBound_CL_Mnt(void);
  
  su2double GetAllBound_CD_Mnt(void);
  
  su2double GetAllBound_CSF_Mnt(void);
  
  su2double GetAllBound_CEff_Mnt(void);
  
  su2double GetAllBound_CMx_Mnt(void);
  
  su2double GetAllBound_CMy_Mnt(void);
  
  su2double GetAllBound_CMz_Mnt(void);
  
  su2double GetAllBound_CFx_Mnt(void);
  
  su2double GetAllBound_CFy_Mnt(void);
  
  su2double GetAllBound_CFz_Mnt(void);
  
  su2double GetCPressure(unsigned short val_marker, unsigned long val_vertex);
  
  su2double GetCPressureTarget(unsigned short val_marker, unsigned long val_vertex);
  
  void SetCPressureTarget(unsigned short val_marker, unsigned long val_vertex, su2double val_pressure);
  
  su2double *GetCharacPrimVar(unsigned short val_marker, unsigned long val_vertex);
  
  void SetCharacPrimVar(unsigned short val_marker, unsigned long val_vertex, unsigned short val_var, su2double val_value);
  
  su2double *GetDonorPrimVar(unsigned short val_marker, unsigned long val_vertex);
  
  void SetDonorPrimVar(unsigned short val_marker, unsigned long val_vertex, unsigned short val_var, su2double val_value);
  
  su2double GetDonorPrimVar(unsigned short val_marker, unsigned long val_vertex, unsigned short val_var);
  
  unsigned long GetDonorGlobalIndex(unsigned short val_marker, unsigned long val_vertex);
  
  void SetDonorGlobalIndex(unsigned short val_marker, unsigned long val_vertex, unsigned long val_index);
  
  su2double GetActDisk_DeltaP(unsigned short val_marker, unsigned long val_vertex);
  
  void SetActDisk_DeltaP(unsigned short val_marker, unsigned long val_vertex, su2double val_deltap);
  
  su2double GetActDisk_DeltaT(unsigned short val_marker, unsigned long val_vertex);
  
  void SetActDisk_DeltaT(unsigned short val_marker, unsigned long val_vertex, su2double val_deltat);
  
  su2double GetInlet_Ttotal(unsigned short val_marker, unsigned long val_vertex);
  
  su2double GetInlet_Ptotal(unsigned short val_marker, unsigned long val_vertex);
  
  su2double GetInlet_FlowDir(unsigned short val_marker, unsigned long val_vertex, unsigned short val_dim);
  
  void SetInlet_Ttotal(unsigned short val_marker, unsigned long val_vertex, su2double val_ttotal);
  
  void SetInlet_Ptotal(unsigned short val_marker, unsigned long val_vertex, su2double val_ptotal);
  
  void SetInlet_FlowDir(unsigned short val_marker, unsigned long val_vertex, unsigned short val_dim, su2double val_flowdir);

  su2double GetOneD_TotalPress(void);
  
  void SetOneD_TotalPress(su2double AveragePressure);
  
  su2double GetOneD_Mach(void);
  
  void SetOneD_Mach(su2double AverageMach);
  
  su2double GetOneD_Temp(void);
  
  void SetOneD_Temp(su2double AverageTemperature);
  
  su2double GetOneD_MassFlowRate(void);
  
  void SetOneD_MassFlowRate(su2double MassFlowRate);
  
  su2double GetOneD_FluxAvgPress(void);
  
  void SetOneD_FluxAvgPress(su2double PressureRef);
  
  su2double GetOneD_FluxAvgDensity(void);
  
  void SetOneD_FluxAvgDensity(su2double DensityRef);
  
  su2double GetOneD_FluxAvgVelocity(void);
  
  void SetOneD_FluxAvgVelocity(su2double VelocityRef);
  
  su2double GetOneD_FluxAvgEntalpy(void);
  
  void SetOneD_FluxAvgEntalpy(su2double EnthalpyRef);
  
  void SetResidual_DualTime(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                            unsigned short iRKStep, unsigned short iMesh, unsigned short RunTime_EqSystem);
  
  void SetFlow_Displacement(CGeometry **flow_geometry, CVolumetricMovement *flow_grid_movement, CConfig *flow_config, CConfig *fea_config,
                            CGeometry **fea_geometry, CSolver ***fea_solution);
  
  void SetFlow_Displacement_Int(CGeometry **flow_geometry, CVolumetricMovement *flow_grid_movement, CConfig *flow_config, CConfig *fea_config,
                                CGeometry **fea_geometry, CSolver ***fea_solution);
  
  void LoadRestart(CGeometry **geometry, CSolver ***solver, CConfig *config, int val_iter);
  
  void SetSlidingState(unsigned short val_marker, unsigned long val_vertex, unsigned short val_state, su2double component);
    
  void SetInitialCondition(CGeometry **geometry, CSolver ***solver_container, CConfig *config, unsigned long ExtIter);
  
  void SetPressure_Inf(su2double p_inf);
  
  void SetTemperature_Inf(su2double t_inf);
  
  void SetFreeStream_Solution(CConfig *config);
};

  

class CIncEulerSolver : public CSolver {
protected:
  
  su2double
  Density_Inf,  
  Temperature_Inf,      
  Pressure_Inf,    
  *Velocity_Inf;    
  su2double
  *CD_Inv,  
  *CL_Inv,      
  *CSF_Inv,    
  *CMx_Inv,      
  *CMy_Inv,      
  *CMz_Inv,      
  *CFx_Inv,      
  *CFy_Inv,      
  *CFz_Inv,      
  *Surface_CL_Inv, 
  *Surface_CD_Inv, 
  *Surface_CSF_Inv, 
  *Surface_CEff_Inv, 
  *Surface_CFx_Inv,   
  *Surface_CFy_Inv,   
  *Surface_CFz_Inv,   
  *Surface_CMx_Inv,   
  *Surface_CMy_Inv,   
  *Surface_CMz_Inv,   
  *CEff_Inv,        
  *CMerit_Inv,        
  *CT_Inv,      
  *CQ_Inv,      
  *CD_Mnt,  
  *CL_Mnt,      
  *CSF_Mnt,    
  *CMx_Mnt,      
  *CMy_Mnt,      
  *CMz_Mnt,      
  *CFx_Mnt,      
  *CFy_Mnt,      
  *CFz_Mnt,      
  *Surface_CL_Mnt, 
  *Surface_CD_Mnt, 
  *Surface_CSF_Mnt, 
  *Surface_CEff_Mnt, 
  *Surface_CFx_Mnt,   
  *Surface_CFy_Mnt,   
  *Surface_CFz_Mnt,   
  *Surface_CMx_Mnt,   
  *Surface_CMy_Mnt,   
  *Surface_CMz_Mnt,   
  *CEff_Mnt,        
  *CMerit_Mnt,        
  *CT_Mnt,      
  *CQ_Mnt,      
  **CPressure,    
  **CPressureTarget,    
  **HeatFlux,    
  **HeatFluxTarget,    
  **YPlus,    
  ***CharacPrimVar,    
  *ForceInviscid,    
  *MomentInviscid,  
  *ForceMomentum,    
  *MomentMomentum,  
  InverseDesign;  
  su2double
  AllBound_CD_Inv,  
  AllBound_CL_Inv,      
  AllBound_CSF_Inv,      
  AllBound_CMx_Inv,      
  AllBound_CMy_Inv,      
  AllBound_CMz_Inv,      
  AllBound_CFx_Inv,      
  AllBound_CFy_Inv,      
  AllBound_CFz_Inv,      
  AllBound_CEff_Inv,      
  AllBound_CMerit_Inv,      
  AllBound_CT_Inv,      
  AllBound_CQ_Inv;      
  su2double
  AllBound_CD_Mnt,  
  AllBound_CL_Mnt,      
  AllBound_CSF_Mnt,      
  AllBound_CMx_Mnt,      
  AllBound_CMy_Mnt,      
  AllBound_CMz_Mnt,      
  AllBound_CFx_Mnt,      
  AllBound_CFy_Mnt,      
  AllBound_CFz_Mnt,      
  AllBound_CEff_Mnt,      
  AllBound_CMerit_Mnt,      
  AllBound_CT_Mnt,      
  AllBound_CQ_Mnt;      
  su2double
  OneD_TotalPress, 
  OneD_Mach, 
  OneD_Temp, 
  OneD_PressureRef, 
  OneD_MassFlowRate, 
  OneD_DensityRef, 
  OneD_EnthalpyRef, 
  OneD_VelocityRef, 
  AoA_Prev, 
  Total_CL_Prev,    
  Total_CD_Prev, 
  Total_ComboObj, 
  Total_CD, 
  Total_CL,    
  Total_CSF,    
  Total_CMx,      
  Total_CMy,      
  Total_CMz,      
  Total_CFx,      
  Total_CFy,      
  Total_CFz,      
  Total_CEff,      
  Total_CMerit,      
  Total_CT,    
  Total_CQ,    
  Total_Heat,    
  Total_MaxHeat, 
  Total_CpDiff,      
  Total_HeatFluxDiff,      
  Total_MassFlowRate;     
  su2double *Surface_CL,   
  *Surface_CD,          
  *Surface_CSF,     
  *Surface_CEff,     
  *Surface_CFx,            
  *Surface_CFy,            
  *Surface_CFz,            
  *Surface_CMx,            
  *Surface_CMy,            
  *Surface_CMz;            
  su2double *iPoint_UndLapl,  
  *jPoint_UndLapl;      
  su2double *SecondaryVar_i,  
  *SecondaryVar_j;      
  su2double *PrimVar_i,  
  *PrimVar_j;      
  unsigned long nMarker,        
  *nVertex;       
  bool space_centered,  
  euler_implicit,      
  least_squares;        
  su2double Gamma;                  
  su2double Gamma_Minus_One;        
  su2double *Primitive,    
  *Primitive_i,        
  *Primitive_j;        
  su2double Cauchy_Value,  
  Cauchy_Func;      
  unsigned short Cauchy_Counter;  
  su2double *Cauchy_Serie;      
  su2double Old_Func,  
  New_Func;      
  su2double AoA_old;  
  unsigned long AoA_Counter;
  
  CFluidModel  *FluidModel;  
public:
  
  CIncEulerSolver(void);
  
  CIncEulerSolver(CGeometry *geometry, CConfig *config, unsigned short iMesh);
  
  virtual ~CIncEulerSolver(void);
  
  void Set_MPI_Solution(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Gradient(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Old(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Limiter(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Primitive_Limiter(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief Impose the send-receive boundary condition.
  //   * \param[in] geometry - Geometrical definition of the problem.
  //   * \param[in] config - Definition of the particular problem.
  //   */
  //  void Set_MPI_Secondary_Limiter(CGeometry *geometry, CConfig *config);
  
  void SetNondimensionalization(CGeometry *geometry, CConfig *config, unsigned short iMesh);
  
  CFluidModel* GetFluidModel(void);
  
  
  su2double GetDensity_Inf(void);
  
  su2double GetModVelocity_Inf(void);
  
  su2double GetPressure_Inf(void);
  
  su2double GetDensity_Velocity_Inf(unsigned short val_dim);
  
  su2double GetVelocity_Inf(unsigned short val_dim);
  
  su2double *GetVelocity_Inf(void);
  
  void SetTime_Step(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                    unsigned short iMesh, unsigned long Iteration);
  
  void Centered_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                         CConfig *config, unsigned short iMesh, unsigned short iRKStep);
  
  void Upwind_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                       CConfig *config, unsigned short iMesh);

  void Source_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CNumerics *second_numerics,
                       CConfig *config, unsigned short iMesh);
  
  void Source_Template(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                       CConfig *config, unsigned short iMesh);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  void Postprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh);
  
  unsigned long SetPrimitive_Variables(CSolver **solver_container, CConfig *config, bool Output);
  
  void SetDissipation_Switch(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Dissipation_Switch(CGeometry *geometry, CConfig *config);
  
  void SetPrimitive_Gradient_GG(CGeometry *geometry, CConfig *config);
  
  void SetPrimitive_Gradient_LS(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief Compute the gradient of the primitive variables using Green-Gauss method,
  //   *        and stores the result in the <i>Gradient_Primitive</i> variable.
  //   * \param[in] geometry - Geometrical definition of the problem.
  //   * \param[in] config - Definition of the particular problem.
  //   */
  //  void SetSecondary_Gradient_GG(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief Compute the gradient of the primitive variables using a Least-Squares method,
  //   *        and stores the result in the <i>Gradient_Primitive</i> variable.
  //   * \param[in] geometry - Geometrical definition of the problem.
  //   * \param[in] config - Definition of the particular problem.
  //   */
  //  void SetSecondary_Gradient_LS(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Primitive_Gradient(CGeometry *geometry, CConfig *config);
  
  void SetPrimitive_Limiter(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief Compute the gradient of the primitive variables using a Least-Squares method,
  //   *        and stores the result in the <i>Gradient_Primitive</i> variable.
  //   * \param[in] geometry - Geometrical definition of the problem.
  //   * \param[in] config - Definition of the particular problem.
  //   */
  //  void Set_MPI_Secondary_Gradient(CGeometry *geometry, CConfig *config);
  
  //  /*!
  //   * \brief Compute the limiter of the primitive variables.
  //   * \param[in] geometry - Geometrical definition of the problem.
  //   * \param[in] config - Definition of the particular problem.
  //   */
  //  void SetSecondary_Limiter(CGeometry *geometry, CConfig *config);

  void SetUndivided_Laplacian(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Undivided_Laplacian(CGeometry *geometry, CConfig *config);
  
  void SetMax_Eigenvalue(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_MaxEigenvalue(CGeometry *geometry, CConfig *config);
  
  void BC_Euler_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                     unsigned short val_marker);
  
  void BC_Far_Field(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                    CConfig *config, unsigned short val_marker);
  
  void BC_Sym_Plane(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  void BC_Inlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                CConfig *config, unsigned short val_marker);
  
  void BC_Custom(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                 CConfig *config, unsigned short val_marker);
  
  void BC_Outlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                 CConfig *config, unsigned short val_marker);
  
  static bool Compareval(std::vector<su2double> a,std::vector<su2double> b);
  
  void ExplicitRK_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                            unsigned short iRKStep);

  void SetFarfield_AoA(CGeometry *geometry, CSolver **solver_container,
                       CConfig *config, unsigned short iMesh, bool Output);
  
  void ExplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void Pressure_Forces(CGeometry *geometry, CConfig *config);

  void Momentum_Forces(CGeometry *geometry, CConfig *config);

  void ImplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);

  su2double GetCLift_Inv(unsigned short val_marker);
  
  su2double GetCMz_Inv(unsigned short val_marker);
  
  su2double GetSurface_CL(unsigned short val_marker);
  
  su2double GetSurface_CD(unsigned short val_marker);
  
  su2double GetSurface_CSF(unsigned short val_marker);
  
  su2double GetSurface_CEff(unsigned short val_marker);
  
  su2double GetSurface_CFx(unsigned short val_marker);
  
  su2double GetSurface_CFy(unsigned short val_marker);
  
  su2double GetSurface_CFz(unsigned short val_marker);
  
  su2double GetSurface_CMx(unsigned short val_marker);
  
  su2double GetSurface_CMy(unsigned short val_marker);
  
  su2double GetSurface_CMz(unsigned short val_marker);
  
  su2double GetSurface_CL_Inv(unsigned short val_marker);
  
  su2double GetSurface_CD_Inv(unsigned short val_marker);
  
  su2double GetSurface_CSF_Inv(unsigned short val_marker);
  
  su2double GetSurface_CEff_Inv(unsigned short val_marker);
  
  su2double GetSurface_CFx_Inv(unsigned short val_marker);
  
  su2double GetSurface_CFy_Inv(unsigned short val_marker);
  
  su2double GetSurface_CFz_Inv(unsigned short val_marker);
  
  su2double GetSurface_CMx_Inv(unsigned short val_marker);
  
  su2double GetSurface_CMy_Inv(unsigned short val_marker);
  
  su2double GetSurface_CMz_Inv(unsigned short val_marker);

  su2double GetSurface_CL_Mnt(unsigned short val_marker);

  su2double GetSurface_CD_Mnt(unsigned short val_marker);

  su2double GetSurface_CSF_Mnt(unsigned short val_marker);

  su2double GetSurface_CEff_Mnt(unsigned short val_marker);

  su2double GetSurface_CFx_Mnt(unsigned short val_marker);

  su2double GetSurface_CFy_Mnt(unsigned short val_marker);

  su2double GetSurface_CFz_Mnt(unsigned short val_marker);

  su2double GetSurface_CMx_Mnt(unsigned short val_marker);

  su2double GetSurface_CMy_Mnt(unsigned short val_marker);

  su2double GetSurface_CMz_Mnt(unsigned short val_marker);
  su2double GetCD_Inv(unsigned short val_marker);
  
  su2double GetCSF_Inv(unsigned short val_marker);
  
  su2double GetCEff_Inv(unsigned short val_marker);
  
  su2double GetTotal_CSF(void);
  
  su2double GetTotal_CEff(void);
  
  su2double GetTotal_CpDiff(void);
  
  su2double GetTotal_HeatFluxDiff(void);
  
  void SetTotal_CpDiff(su2double val_pressure);
  
  void SetTotal_HeatFluxDiff(su2double val_heat);
  
  void SetTotal_CLift(su2double val_Total_CLift);
  
  su2double GetTotal_CL(void);

  void SetTotal_ComboObj(su2double ComboObj);

  su2double GetTotal_ComboObj(void);

  su2double GetTotal_CD(void);
  
  su2double GetTotal_CMx(void);
  
  su2double GetTotal_CMy(void);
  
  su2double GetTotal_CMz(void);
  
  su2double GetTotal_CFx(void);
  
  su2double GetTotal_CFy(void);
  
  su2double GetTotal_CFz(void);
  
  su2double GetTotal_CT(void);
  
  void SetTotal_CT(su2double val_Total_CT);
  
  su2double GetTotal_CQ(void);
  
  su2double GetTotal_HeatFlux(void);
  
  su2double GetTotal_MaxHeatFlux(void);
  
  void SetTotal_CQ(su2double val_Total_CQ);
  
  void SetTotal_HeatFlux(su2double val_Total_Heat);
  
  void SetTotal_MaxHeatFlux(su2double val_Total_MaxHeat);
  
  su2double GetTotal_CMerit(void);
  
  void SetTotal_CD(su2double val_Total_CDrag);
  
  su2double GetAllBound_CL_Inv(void);
  
  su2double GetAllBound_CD_Inv(void);
  
  su2double GetAllBound_CSF_Inv(void);
  
  su2double GetAllBound_CEff_Inv(void);
  
  su2double GetAllBound_CMx_Inv(void);
  
  su2double GetAllBound_CMy_Inv(void);
  
  su2double GetAllBound_CMz_Inv(void);
  
  su2double GetAllBound_CFx_Inv(void);
  
  su2double GetAllBound_CFy_Inv(void);
  
  su2double GetAllBound_CFz_Inv(void);
  
  su2double GetAllBound_CL_Mnt(void);

  su2double GetAllBound_CD_Mnt(void);

  su2double GetAllBound_CSF_Mnt(void);

  su2double GetAllBound_CEff_Mnt(void);

  su2double GetAllBound_CMx_Mnt(void);

  su2double GetAllBound_CMy_Mnt(void);

  su2double GetAllBound_CMz_Mnt(void);

  su2double GetAllBound_CFx_Mnt(void);

  su2double GetAllBound_CFy_Mnt(void);

  su2double GetAllBound_CFz_Mnt(void);

  su2double GetCPressure(unsigned short val_marker, unsigned long val_vertex);
  
  su2double GetCPressureTarget(unsigned short val_marker, unsigned long val_vertex);
  
  void SetCPressureTarget(unsigned short val_marker, unsigned long val_vertex, su2double val_pressure);
  
  su2double *GetCharacPrimVar(unsigned short val_marker, unsigned long val_vertex);
  
  su2double GetOneD_TotalPress(void);
  
  void SetOneD_TotalPress(su2double AveragePressure);
  
  su2double GetOneD_Mach(void);
  
  void SetOneD_Mach(su2double AverageMach);
  
  su2double GetOneD_Temp(void);
  
  void SetOneD_Temp(su2double AverageTemperature);
  
  su2double GetOneD_MassFlowRate(void);
  
  void SetOneD_MassFlowRate(su2double MassFlowRate);
  
  su2double GetOneD_FluxAvgPress(void);
  
  void SetOneD_FluxAvgPress(su2double PressureRef);
  
  su2double GetOneD_FluxAvgDensity(void);
  
  void SetOneD_FluxAvgDensity(su2double DensityRef);
  
  su2double GetOneD_FluxAvgVelocity(void);
  
  void SetOneD_FluxAvgVelocity(su2double VelocityRef);
  
  su2double GetOneD_FluxAvgEntalpy(void);
  
  void SetOneD_FluxAvgEntalpy(su2double EnthalpyRef);
  
  void SetResidual_DualTime(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                            unsigned short iRKStep, unsigned short iMesh, unsigned short RunTime_EqSystem);
  
  void SetFlow_Displacement(CGeometry **flow_geometry, CVolumetricMovement *flow_grid_movement, CConfig *flow_config, CConfig *fea_config,
                            CGeometry **fea_geometry, CSolver ***fea_solution);
  
  void SetFlow_Displacement_Int(CGeometry **flow_geometry, CVolumetricMovement *flow_grid_movement, CConfig *flow_config, CConfig *fea_config,
                                CGeometry **fea_geometry, CSolver ***fea_solution);
  
  void LoadRestart(CGeometry **geometry, CSolver ***solver, CConfig *config, int val_iter);
  
  void SetInitialCondition(CGeometry **geometry, CSolver ***solver_container, CConfig *config, unsigned long ExtIter);
  
  void SetPressure_Inf(su2double p_inf);
  
  void SetTemperature_Inf(su2double t_inf);
  
  void SetFreeStream_Solution(CConfig *config);
};

class CNSSolver : public CEulerSolver {
private:
  su2double Viscosity_Inf;  
  su2double Tke_Inf;  
  su2double Prandtl_Lam,   
  Prandtl_Turb;         
  su2double *CD_Visc,  
  *CL_Visc,    
  *CSF_Visc,    
  *CMx_Visc,      
  *CMy_Visc,      
  *CMz_Visc,      
  *CFx_Visc,      
  *CFy_Visc,      
  *CFz_Visc,      
  *Surface_CL_Visc,
  *Surface_CD_Visc,
  *Surface_CSF_Visc,
  *Surface_CEff_Visc,
  *Surface_CFx_Visc,  
  *Surface_CFy_Visc,  
  *Surface_CFz_Visc,  
  *Surface_CMx_Visc,  
  *Surface_CMy_Visc,  
  *Surface_CMz_Visc,  
  *CEff_Visc,      
  *CMerit_Visc,      
  *CT_Visc,    
  *CQ_Visc,    
  *HF_Visc,    
  *MaxHF_Visc, 
  ***CSkinFriction;  
  su2double *ForceViscous,  
  *MomentViscous;      
  su2double AllBound_CD_Visc, 
  AllBound_CL_Visc,    
  AllBound_CSF_Visc,    
  AllBound_CMx_Visc,      
  AllBound_CMy_Visc,      
  AllBound_CMz_Visc,      
  AllBound_CEff_Visc,      
  AllBound_CFx_Visc,      
  AllBound_CFy_Visc,      
  AllBound_CFz_Visc,      
  AllBound_CMerit_Visc,      
  AllBound_CT_Visc,    
  AllBound_CQ_Visc,    
  AllBound_HF_Visc,    
  AllBound_MaxHF_Visc; 
  su2double StrainMag_Max, Omega_Max; 
public:
  
  CNSSolver(void);
  
  CNSSolver(CGeometry *geometry, CConfig *config, unsigned short iMesh);
  
  ~CNSSolver(void);
  
  su2double GetSurface_CL_Visc(unsigned short val_marker);
  
  su2double GetSurface_CD_Visc(unsigned short val_marker);
  
  su2double GetSurface_CSF_Visc(unsigned short val_marker);
  
  su2double GetSurface_CEff_Visc(unsigned short val_marker);
  
  su2double GetSurface_CFx_Visc(unsigned short val_marker);
  
  su2double GetSurface_CFy_Visc(unsigned short val_marker);
  
  su2double GetSurface_CFz_Visc(unsigned short val_marker);
  
  su2double GetSurface_CMx_Visc(unsigned short val_marker);
  
  su2double GetSurface_CMy_Visc(unsigned short val_marker);
  
  su2double GetSurface_CMz_Visc(unsigned short val_marker);
  
  su2double GetAllBound_CL_Visc(void);
  
  su2double GetAllBound_CD_Visc(void);
  
  su2double GetAllBound_CSF_Visc(void);
  
  su2double GetAllBound_CEff_Visc(void);
  
  su2double GetAllBound_CMx_Visc(void);
  
  su2double GetAllBound_CMy_Visc(void);
  
  su2double GetAllBound_CMz_Visc(void);
  
  su2double GetAllBound_CFx_Visc(void);
  
  su2double GetAllBound_CFy_Visc(void);
  
  su2double GetAllBound_CFz_Visc(void);
  
  su2double GetViscosity_Inf(void);
  
  su2double GetTke_Inf(void);
  
  void SetTime_Step(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                    unsigned short iMesh, unsigned long Iteration);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  unsigned long SetPrimitive_Variables(CSolver **solver_container, CConfig *config, bool Output);
  
  void BC_HeatFlux_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  void BC_Isothermal_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                          unsigned short val_marker);
  
  void Friction_Forces(CGeometry *geometry, CConfig *config);
  
  su2double GetSurface_HF_Visc(unsigned short val_marker);
  
  su2double GetSurface_MaxHF_Visc(unsigned short val_marker);
  
  su2double GetCL_Visc(unsigned short val_marker);
  
  su2double GetCMz_Visc(unsigned short val_marker);
  
  su2double GetCSF_Visc(unsigned short val_marker);
  
  su2double GetCD_Visc(unsigned short val_marker);
  
  void Viscous_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                        CConfig *config, unsigned short iMesh, unsigned short iRKStep);
  
  su2double GetCSkinFriction(unsigned short val_marker, unsigned long val_vertex, unsigned short val_dim);
  
  su2double GetHeatFlux(unsigned short val_marker, unsigned long val_vertex);
  
  su2double GetHeatFluxTarget(unsigned short val_marker, unsigned long val_vertex);
  
  void SetHeatFluxTarget(unsigned short val_marker, unsigned long val_vertex, su2double val_heat);
  
  su2double GetYPlus(unsigned short val_marker, unsigned long val_vertex);
  
  su2double GetOmega_Max(void);
  
  su2double GetStrainMag_Max(void);
  
  void SetStrainMag_Max(su2double val_strainmag_max);
  
  void SetOmega_Max(su2double val_omega_max);
  
};

class CIncNSSolver : public CIncEulerSolver {
private:
  su2double Viscosity_Inf;  
  su2double Tke_Inf;  
  su2double Prandtl_Lam,   
  Prandtl_Turb;         
  su2double *CD_Visc,  
  *CL_Visc,    
  *CSF_Visc,    
  *CMx_Visc,      
  *CMy_Visc,      
  *CMz_Visc,      
  *CFx_Visc,      
  *CFy_Visc,      
  *CFz_Visc,      
  *Surface_CL_Visc,
  *Surface_CD_Visc,
  *Surface_CSF_Visc,
  *Surface_CEff_Visc,
  *Surface_CFx_Visc,  
  *Surface_CFy_Visc,  
  *Surface_CFz_Visc,  
  *Surface_CMx_Visc,  
  *Surface_CMy_Visc,  
  *Surface_CMz_Visc,  
  *CEff_Visc,      
  *CMerit_Visc,      
  *CT_Visc,    
  *CQ_Visc,    
  *HF,    
  *MaxHF_Visc, 
  ***CSkinFriction;  
  su2double *ForceViscous,  
  *MomentViscous;      
  su2double AllBound_CD_Visc, 
  AllBound_CL_Visc,    
  AllBound_CSF_Visc,    
  AllBound_CMx_Visc,      
  AllBound_CMy_Visc,      
  AllBound_CMz_Visc,      
  AllBound_CEff_Visc,      
  AllBound_CFx_Visc,      
  AllBound_CFy_Visc,      
  AllBound_CFz_Visc,      
  AllBound_CMerit_Visc,      
  AllBound_CT_Visc,    
  AllBound_CQ_Visc,    
  AllBound_HF_Visc,    
  AllBound_MaxHF_Visc; 
  su2double StrainMag_Max, Omega_Max; 
public:
  
  CIncNSSolver(void);
    
  CIncNSSolver(CGeometry *geometry, CConfig *config, unsigned short iMesh);
    
  ~CIncNSSolver(void);
    
  su2double GetSurface_CL_Visc(unsigned short val_marker);
  
  su2double GetSurface_CD_Visc(unsigned short val_marker);
    
  su2double GetSurface_CSF_Visc(unsigned short val_marker);
    
  su2double GetSurface_CEff_Visc(unsigned short val_marker);
    
  su2double GetSurface_CFx_Visc(unsigned short val_marker);
  
  su2double GetSurface_CFy_Visc(unsigned short val_marker);
    
  su2double GetSurface_CFz_Visc(unsigned short val_marker);
    
  su2double GetSurface_CMx_Visc(unsigned short val_marker);
  
  su2double GetSurface_CMy_Visc(unsigned short val_marker);
  
  su2double GetSurface_CMz_Visc(unsigned short val_marker);
  
  su2double GetAllBound_CL_Visc(void);
  
  su2double GetAllBound_CD_Visc(void);
  
  su2double GetAllBound_CSF_Visc(void);
  
  su2double GetAllBound_CEff_Visc(void);
  
  su2double GetAllBound_CMx_Visc(void);
  
  su2double GetAllBound_CMy_Visc(void);

  su2double GetAllBound_CMz_Visc(void);
  
  su2double GetAllBound_CFx_Visc(void);
  
  su2double GetAllBound_CFy_Visc(void);
  
  su2double GetAllBound_CFz_Visc(void);
  
  su2double GetViscosity_Inf(void);
  
  su2double GetTke_Inf(void);
  
  void SetTime_Step(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                    unsigned short iMesh, unsigned long Iteration);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  unsigned long SetPrimitive_Variables(CSolver **solver_container, CConfig *config, bool Output);
  
  void BC_HeatFlux_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  void Friction_Forces(CGeometry *geometry, CConfig *config);
  
  su2double GetCL_Visc(unsigned short val_marker);
  
  su2double GetCMz_Visc(unsigned short val_marker);
  
  su2double GetCSF_Visc(unsigned short val_marker);
  
  su2double GetCD_Visc(unsigned short val_marker);
  
  void Viscous_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                        CConfig *config, unsigned short iMesh, unsigned short iRKStep);
  
  su2double GetCSkinFriction(unsigned short val_marker, unsigned long val_vertex, unsigned short val_dim);
  
  su2double GetHeatFlux(unsigned short val_marker, unsigned long val_vertex);
  
  su2double GetHeatFluxTarget(unsigned short val_marker, unsigned long val_vertex);
  
  void SetHeatFluxTarget(unsigned short val_marker, unsigned long val_vertex, su2double val_heat);
  
  su2double GetYPlus(unsigned short val_marker, unsigned long val_vertex);
  
  su2double GetOmega_Max(void);
  
  su2double GetStrainMag_Max(void);
  
  void SetStrainMag_Max(su2double val_strainmag_max);
  
  void SetOmega_Max(su2double val_omega_max);
  
};

class CTurbSolver : public CSolver {
protected:
  su2double *FlowPrimVar_i,  
  *FlowPrimVar_j,         
  *lowerlimit,            
  *upperlimit;            
  su2double Gamma;           
  su2double Gamma_Minus_One; 
public:
  
  CTurbSolver(void);
  
  virtual ~CTurbSolver(void);
  
  CTurbSolver(CConfig *config);
  
  void Set_MPI_Solution(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Old(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Gradient(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Limiter(CGeometry *geometry, CConfig *config);
  
  void Upwind_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                       unsigned short iMesh);
  
  void Viscous_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                        CConfig *config, unsigned short iMesh, unsigned short iRKStep);
  
  void BC_Sym_Plane(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  void BC_Euler_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                     unsigned short val_marker);
  
  void ImplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void SetResidual_DualTime(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                            unsigned short iRKStep, unsigned short iMesh, unsigned short RunTime_EqSystem);
  
  void LoadRestart(CGeometry **geometry, CSolver ***solver, CConfig *config, int val_iter);
  
};

class CTurbSASolver: public CTurbSolver {
private:
  su2double nu_tilde_Inf, nu_tilde_Engine, nu_tilde_ActDisk;
  
public:
  CTurbSASolver(void);
  
  CTurbSASolver(CGeometry *geometry, CConfig *config, unsigned short iMesh, CFluidModel* FluidModel);
  
  ~CTurbSASolver(void);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  void Postprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                      unsigned short iMesh);
  
  void Source_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CNumerics *second_numerics,
                       CConfig *config, unsigned short iMesh);
  
  void Source_Template(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                       CConfig *config, unsigned short iMesh);
  
  void BC_HeatFlux_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                        unsigned short val_marker);
  
  void BC_Isothermal_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                          unsigned short val_marker);
  
  void BC_Far_Field(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                    unsigned short val_marker);
  
  void BC_Inlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                unsigned short val_marker);
  
  void BC_Outlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                 unsigned short val_marker);
  
  void BC_Engine_Inflow(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                        CConfig *config, unsigned short val_marker);
  
  void BC_Engine_Exhaust(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                         CConfig *config, unsigned short val_marker);
  
  void BC_Interface_Boundary(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                             CConfig *config, unsigned short val_marker);
  
  void BC_NearField_Boundary(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                             CConfig *config, unsigned short val_marker);
  
  void BC_ActDisk_Inlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                        CConfig *config, unsigned short val_marker);
  
  void BC_ActDisk_Outlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                         CConfig *config, unsigned short val_marker);
  
  void BC_ActDisk(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                  CConfig *config, unsigned short val_marker, bool inlet_surface);

  void SetFreeStream_Solution(CConfig *config);
  
};

class CTurbSSTSolver: public CTurbSolver {
private:
  su2double *constants,  
  kine_Inf,           
  omega_Inf;          
public:
  CTurbSSTSolver(void);
  
  CTurbSSTSolver(CGeometry *geometry, CConfig *config, unsigned short iMesh);
  
  ~CTurbSSTSolver(void);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  void Postprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                      unsigned short iMesh);
  
  void Source_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CNumerics *second_numerics,
                       CConfig *config, unsigned short iMesh);
  
  void Source_Template(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                       CConfig *config, unsigned short iMesh);
  
  void BC_HeatFlux_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                        unsigned short val_marker);
  
  void BC_Isothermal_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                          unsigned short val_marker);
  
  void BC_Far_Field(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                    unsigned short val_marker);
  
  void BC_Inlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                unsigned short val_marker);
  
  void BC_Outlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                 unsigned short val_marker);
  
  su2double* GetConstants();
  
  void SetFreeStream_Solution(CConfig *config);
  
};

class CTransLMSolver: public CTurbSolver {
private:
  su2double Intermittency_Inf, REth_Inf;
public:
  CTransLMSolver(void);
  
  CTransLMSolver(CGeometry *geometry, CConfig *config, unsigned short iMesh);
  
  ~CTransLMSolver(void);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  void Postprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                      unsigned short iMesh);
  
  void Upwind_Residual(CGeometry *geometry, CSolver **solver_container,
                       CNumerics *numerics, CConfig *config, unsigned short iMesh);
  
  void Viscous_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                        CConfig *config, unsigned short iMesh, unsigned short iRKStep);
  
  void Source_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CNumerics *second_numerics,
                       CConfig *config, unsigned short iMesh);
  
  void Source_Template(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                       CConfig *config, unsigned short iMesh);
  
  void BC_HeatFlux_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                        unsigned short val_marker);
  
  void BC_Far_Field(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                    unsigned short val_marker);
  
  void BC_Inlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                unsigned short val_marker);
  
  void BC_Outlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                 unsigned short val_marker);
  
  void BC_Sym_Plane(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                    unsigned short val_marker);
  void ImplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  // Another set of matrix structures for the Lm equations
  CSysMatrix JacobianItmc; 
  su2double *LinSysSolItmc;    
  su2double *LinSysResItmc;    
  su2double *rhsItmc;    
  CSysMatrix JacobianReth; 
  su2double *LinSysSolReth;    
  su2double *LinSysResReth;    
  su2double *rhsReth;    
};

class CAdjEulerSolver : public CSolver {
protected:
  su2double PsiRho_Inf,  
  PsiE_Inf,      
  *Phi_Inf;      
  su2double *Sens_Mach, 
  *Sens_AoA,      
  *Sens_Geo,      
  *Sens_Press,      
  *Sens_Temp,      
  *Sens_BPress,     
  **CSensitivity,    
  ***DonorAdjVar;    
  su2double Total_Sens_Mach;  
  su2double Total_Sens_AoA;    
  su2double Total_Sens_Geo;    
  su2double Total_Sens_Press;    
  su2double Total_Sens_Temp;    
  su2double Total_Sens_BPress;    
  su2double *iPoint_UndLapl,  
  *jPoint_UndLapl;      
  bool space_centered;  
  su2double **Jacobian_Axisymmetric; 
  unsigned long nMarker;        
  su2double Gamma;                  
  su2double Gamma_Minus_One;        
  su2double *FlowPrimVar_i,  
  *FlowPrimVar_j;        
  unsigned long **DonorGlobalIndex;    
  su2double pnorm,
  Area_Monitored; 
  unsigned long AoA_Counter;
  su2double ACoeff, ACoeff_inc, ACoeff_old;
  bool Update_ACoeff;
  
public:
  
  CAdjEulerSolver(void);
  
  CAdjEulerSolver(CGeometry *geometry, CConfig *config, unsigned short iMesh);
  
  virtual ~CAdjEulerSolver(void);
  
  void SetTime_Step(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                    unsigned short iMesh, unsigned long Iteration);
  
  void Set_MPI_Solution(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Gradient(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Limiter(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Old(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Nearfield(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Interface(CGeometry *geometry, CConfig *config);

  void Set_MPI_ActDisk(CSolver **solver_container, CGeometry *geometry, CConfig *config);

  void SetForceProj_Vector(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void SetIntBoundary_Jump(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  su2double GetPsiRho_Inf(void);
  
  su2double GetPsiE_Inf(void);
  
  su2double GetPhi_Inf(unsigned short val_dim);
  
  void Centered_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                         unsigned short iMesh, unsigned short iRKStep);
  
  void Upwind_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                       unsigned short iMesh);
  
  void Source_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CNumerics *second_numerics,
                       CConfig *config, unsigned short iMesh);
  
  void Source_Template(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                       CConfig *config, unsigned short iMesh);
  
  void SetUndivided_Laplacian(CGeometry *geometry, CConfig *config);
  
  su2double *GetDonorAdjVar(unsigned short val_marker, unsigned long val_vertex);
  
  void SetDonorAdjVar(unsigned short val_marker, unsigned long val_vertex, unsigned short val_var, su2double val_value);
  
  su2double GetDonorAdjVar(unsigned short val_marker, unsigned long val_vertex, unsigned short val_var);
  
  unsigned long GetDonorGlobalIndex(unsigned short val_marker, unsigned long val_vertex);
  
  void SetDonorGlobalIndex(unsigned short val_marker, unsigned long val_vertex, unsigned long val_index);

  void Set_MPI_Undivided_Laplacian(CGeometry *geometry, CConfig *config);
  
  void SetDissipation_Switch(CGeometry *geometry, CConfig *config);
  
  void SetFarfield_AoA(CGeometry *geometry, CSolver **solver_container,
                       CConfig *config, unsigned short iMesh, bool Output);
  
  void Set_MPI_Dissipation_Switch(CGeometry *geometry, CConfig *config);
  
  void BC_Euler_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                     unsigned short val_marker);
  
  void BC_Interface_Boundary(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                             CConfig *config, unsigned short val_marker);
  
  void BC_NearField_Boundary(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                             CConfig *config, unsigned short val_marker);
  
  void BC_ActDisk_Inlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                        CConfig *config, unsigned short val_marker);
  
  void BC_ActDisk_Outlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                         CConfig *config, unsigned short val_marker);
  
  void BC_ActDisk(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                  CConfig *config, unsigned short val_marker, bool inlet_surface);

  void BC_Sym_Plane(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                    unsigned short val_marker);
  
  void BC_Far_Field(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                    unsigned short val_marker);
  
  void BC_Inlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                unsigned short val_marker);
  
  
  void BC_Supersonic_Inlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                           unsigned short val_marker);
  
  void BC_Supersonic_Outlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                            unsigned short val_marker);
  
  void BC_Outlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                 unsigned short val_marker);
  
  void BC_Engine_Inflow(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                        CConfig *config, unsigned short val_marker);
  
  void BC_Engine_Exhaust(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics,
                         CConfig *config, unsigned short val_marker);
  
  void ExplicitRK_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                            unsigned short iRKStep);
  
  void ExplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void ImplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  void Inviscid_Sensitivity(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config);
  
  void Smooth_Sensitivity(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config);
  
  su2double GetCSensitivity(unsigned short val_marker, unsigned long val_vertex);
  
  void SetCSensitivity(unsigned short val_marker, unsigned long val_vertex, su2double val_sensitivity);
  
  su2double GetTotal_Sens_Geo(void);
  
  su2double GetTotal_Sens_Mach(void);
  
  su2double GetTotal_Sens_AoA(void);
  
  su2double GetTotal_Sens_Press(void);
  
  su2double GetTotal_Sens_Temp(void);
  
  su2double GetTotal_Sens_BPress(void);
  
  void SetResidual_DualTime(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                            unsigned short iRKStep, unsigned short iMesh, unsigned short RunTime_EqSystem);
  
  void SetInitialCondition(CGeometry **geometry, CSolver ***solver_container, CConfig *config, unsigned long ExtIter);
  
  
};

class CAdjIncEulerSolver : public CSolver {
protected:
  su2double PsiRho_Inf,  
  PsiE_Inf,      
  *Phi_Inf;      
  su2double *Sens_Mach, 
  *Sens_AoA,      
  *Sens_Geo,      
  *Sens_Press,      
  *Sens_Temp,      
  *Sens_BPress,     
  **CSensitivity;    
  su2double Total_Sens_Mach;  
  su2double Total_Sens_AoA;    
  su2double Total_Sens_Geo;    
  su2double Total_Sens_Press;    
  su2double Total_Sens_Temp;    
  su2double Total_Sens_BPress;    
  su2double *iPoint_UndLapl,  
  *jPoint_UndLapl;      
  bool space_centered;  
  su2double **Jacobian_Axisymmetric; 
  unsigned long nMarker;        
  su2double Gamma;                  
  su2double Gamma_Minus_One;        
  su2double *FlowPrimVar_i,  
  *FlowPrimVar_j;        
  su2double pnorm,
  Area_Monitored; 
  unsigned long AoA_Counter;
  su2double ACoeff, ACoeff_inc, ACoeff_old;
  bool Update_ACoeff;

public:
  
  CAdjIncEulerSolver(void);
  
  CAdjIncEulerSolver(CGeometry *geometry, CConfig *config, unsigned short iMesh);
  
  virtual ~CAdjIncEulerSolver(void);
  
  void SetTime_Step(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                    unsigned short iMesh, unsigned long Iteration);
  
  void Set_MPI_Solution(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Gradient(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Limiter(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Old(CGeometry *geometry, CConfig *config);

  void SetForceProj_Vector(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void SetIntBoundary_Jump(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  su2double GetPsiRho_Inf(void);
  
  su2double GetPsiE_Inf(void);
  
  su2double GetPhi_Inf(unsigned short val_dim);
  
  void Centered_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                         unsigned short iMesh, unsigned short iRKStep);
  
  void Upwind_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                       unsigned short iMesh);
  
  void Source_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CNumerics *second_numerics,
                       CConfig *config, unsigned short iMesh);
  
  void Source_Template(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                       CConfig *config, unsigned short iMesh);
  
  void SetUndivided_Laplacian(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Undivided_Laplacian(CGeometry *geometry, CConfig *config);
  
  void SetDissipation_Switch(CGeometry *geometry, CConfig *config);

  void Set_MPI_Dissipation_Switch(CGeometry *geometry, CConfig *config);
  
  void BC_Euler_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                     unsigned short val_marker);
  
  void BC_Interface_Boundary(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config);
  using CSolver::BC_Interface_Boundary;
  
  void BC_NearField_Boundary(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config);
  using CSolver::BC_NearField_Boundary;
  
  void BC_Sym_Plane(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                    unsigned short val_marker);
  
  void BC_Far_Field(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                    unsigned short val_marker);
  
  void BC_Inlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                unsigned short val_marker);
  
  void BC_Outlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                 unsigned short val_marker);

  void ExplicitRK_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                            unsigned short iRKStep);
  
  void ExplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void ImplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  void Inviscid_Sensitivity(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config);
  
  void Smooth_Sensitivity(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config);
  
  su2double GetCSensitivity(unsigned short val_marker, unsigned long val_vertex);
  
  void SetCSensitivity(unsigned short val_marker, unsigned long val_vertex, su2double val_sensitivity);
  
  su2double GetTotal_Sens_Geo(void);
  
  su2double GetTotal_Sens_Mach(void);
  
  su2double GetTotal_Sens_AoA(void);
  
  su2double GetTotal_Sens_Press(void);
  
  su2double GetTotal_Sens_Temp(void);
  
  su2double GetTotal_Sens_BPress(void);
  
  void SetResidual_DualTime(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                            unsigned short iRKStep, unsigned short iMesh, unsigned short RunTime_EqSystem);
  
  void SetInitialCondition(CGeometry **geometry, CSolver ***solver_container, CConfig *config, unsigned long ExtIter);
  
  
};

class CAdjNSSolver : public CAdjEulerSolver {
public:
  
  CAdjNSSolver(void);
  
  CAdjNSSolver(CGeometry *geometry, CConfig *config, unsigned short iMesh);
  
  ~CAdjNSSolver(void);
  
  void SetTime_Step(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                    unsigned short iMesh, unsigned long Iteration);
  
  
  void BC_HeatFlux_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                        unsigned short val_marker);
  
  void BC_Isothermal_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  void Viscous_Sensitivity(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config);
  
  void Viscous_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                        CConfig *config, unsigned short iMesh, unsigned short iRKStep);
  
  void Source_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CNumerics *second_numerics,
                       CConfig *config, unsigned short iMesh);
  
};

class CAdjIncNSSolver : public CAdjIncEulerSolver {
public:
  
  CAdjIncNSSolver(void);
  
  CAdjIncNSSolver(CGeometry *geometry, CConfig *config, unsigned short iMesh);
  
  ~CAdjIncNSSolver(void);
  
  void SetTime_Step(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                    unsigned short iMesh, unsigned long Iteration);
  
  
  void BC_HeatFlux_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                        unsigned short val_marker);
  
  void BC_Isothermal_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  void Viscous_Sensitivity(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config);
  
  void Viscous_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                        CConfig *config, unsigned short iMesh, unsigned short iRKStep);
  
  void Source_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CNumerics *second_numerics,
                       CConfig *config, unsigned short iMesh);
  
};

class CAdjTurbSolver : public CSolver {
private:
  su2double PsiNu_Inf,  
  *FlowSolution_i,  
  *FlowSolution_j;  
  su2double Gamma;                  
  su2double Gamma_Minus_One;        
public:
  
  CAdjTurbSolver(void);
  
  CAdjTurbSolver(CGeometry *geometry, CConfig *config, unsigned short iMesh);
  
  void Set_MPI_Solution(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Old(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Gradient(CGeometry *geometry, CConfig *config);
  
  virtual ~CAdjTurbSolver(void);
  
  void BC_HeatFlux_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                        unsigned short val_marker);
  
  void BC_Isothermal_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                          unsigned short val_marker);
  
  void BC_Far_Field(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                    unsigned short val_marker);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  void Upwind_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                       unsigned short iMesh);
  
  void Viscous_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                        unsigned short iMesh, unsigned short iRKStep);
  
  void Source_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CNumerics *second_numerics,
                       CConfig *config, unsigned short iMesh);
  
  void ImplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
};

class CPoissonSolver : public CSolver {
private:
  su2double *Source_Vector;      
  su2double **StiffMatrix_Elem; 
  su2double **StiffMatrix_Node;  
public:
  
  CPoissonSolver(void);
  
  CPoissonSolver(CGeometry *geometry, CConfig *config);
  
  void Copy_Zone_Solution(CSolver ***solver1_solution, CGeometry **solver1_geometry, CConfig *solver1_config, CSolver ***solver2_solution, CGeometry **solver2_geometry, CConfig *solver2_config);
  
  ~CPoissonSolver(void);
  
  void Viscous_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                        unsigned short iMesh, unsigned short iRKStep);
  
  void AddStiffMatrix(su2double **StiffMatrix_Elem, unsigned long Point_0, unsigned long Point_1, unsigned long Point_2, unsigned long Point_3);
  
  void Compute_Residual(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                        unsigned short iMesh);
  
  void BC_Dirichlet(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short val_marker);
  
  void BC_Neumann(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config, unsigned short val_marker);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  void Source_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CNumerics *second_numerics,
                       CConfig *config, unsigned short iMesh);
  
  void Source_Template(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                       CConfig *config, unsigned short iMesh);
  
  void ImplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
};

class CWaveSolver : public CSolver {
private:
  su2double *CWave;  
  su2double AllBound_CWave;  
  su2double Total_CWave; 
  CSysMatrix StiffMatrixSpace; 
  CSysMatrix StiffMatrixTime;  
  su2double **StiffMatrix_Elem,      
  **StiffMatrix_Node;              
public:
  
  CWaveSolver(void);
  
  CWaveSolver(CGeometry *geometry, CConfig *config);
  
  ~CWaveSolver(void);
  
  void Viscous_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                        unsigned short iMesh, unsigned short iRKStep);
  
  void BC_Euler_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                     unsigned short val_marker);
  
  void BC_Far_Field(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                    unsigned short val_marker);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  void Source_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CNumerics *second_numerics,
                       CConfig *config, unsigned short iMesh);
  
  void Source_Template(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                       CConfig *config, unsigned short iMesh);
  
  void ImplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void SetResidual_DualTime(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                            unsigned short iRKStep, unsigned short iMesh, unsigned short RunTime_EqSystem);
  
  void LoadRestart(CGeometry **geometry, CSolver ***solver, CConfig *config, int val_iter);
  
  void Wave_Strength(CGeometry *geometry, CConfig *config);
  
  void SetSpace_Matrix(CGeometry *geometry,
                       CConfig   *config);
  
  su2double GetTotal_CWave(void);
  
};

class CHeatSolver : public CSolver {
private:
  su2double *CHeat;       
  su2double Total_CHeat; 
  CSysMatrix StiffMatrixSpace; 
  CSysMatrix StiffMatrixTime;   
  su2double **StiffMatrix_Elem; 
  su2double **StiffMatrix_Node;   
public:
  
  CHeatSolver(void);
  
  CHeatSolver(CGeometry *geometry, CConfig *config);
  
  ~CHeatSolver(void);
  
  void Viscous_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                        unsigned short iMesh, unsigned short iRKStep);
  
  void BC_HeatFlux_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                        unsigned short val_marker);
  
  void BC_Isothermal_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  void Source_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CNumerics *second_numerics,
                       CConfig *config, unsigned short iMesh);
  
  void ImplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void SetResidual_DualTime(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                            unsigned short iRKStep, unsigned short iMesh, unsigned short RunTime_EqSystem);
  
  su2double GetTotal_CHeat(void);
  
};

class CFEM_ElasticitySolver : public CSolver {
private:
  
  su2double  Total_CFEA;        
  unsigned long nElement;
  unsigned short nMarker;
  
  su2double *GradN_X,
  *GradN_x;
  
  su2double **Jacobian_c_ij;      
  su2double **Jacobian_s_ij;      
  su2double **Jacobian_k_ij;      
  su2double **MassMatrix_ij;      
  su2double *Res_Stress_i;      
  su2double *Res_Ext_Surf;      
  su2double *Res_Time_Cont;      
  su2double *Res_FSI_Cont;      
  su2double *Res_Dead_Load;      
  su2double *solutionPredictor;    
  su2double *Solution_Interm;    
  su2double *SolRest;      
  su2double *nodeReactions;      
  su2double *normalVertex;      
  su2double **stressTensor;      
  su2double **mZeros_Aux;      
  su2double **mId_Aux;        
  su2double a_dt[9];          
  su2double Conv_Ref[3];        
  su2double Conv_Check[3];      
  su2double FSI_Conv[2];        
  su2double loadIncrement;      
  su2double WAitken_Dyn;        
  su2double WAitken_Dyn_tn1;      
  CSysMatrix MassMatrix;       
  CSysVector TimeRes_Aux;      
  CSysVector TimeRes;        
  CSysVector LinSysReact;      
public:
  
  CElement*** element_container;   
  CFEM_ElasticitySolver(void);
  
  CFEM_ElasticitySolver(CGeometry *geometry, CConfig *config);
  
  virtual ~CFEM_ElasticitySolver(void);
  
  void Set_MPI_Solution(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Old(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_DispOnly(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Pred(CGeometry *geometry, CConfig *config);
  
  void Set_MPI_Solution_Pred_Old(CGeometry *geometry, CConfig *config);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, CNumerics **numerics, unsigned short iMesh, unsigned long Iteration, unsigned short RunTime_EqSystem, bool Output);
  
  
  void SetInitialCondition(CGeometry **geometry, CSolver ***solver_container, CConfig *config, unsigned long ExtIter);
  
  void ResetInitialCondition(CGeometry **geometry, CSolver ***solver_container, CConfig *config, unsigned long ExtIter);
  
  void SetTime_Step(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                    unsigned short iMesh, unsigned long Iteration);
  
  void Compute_StiffMatrix(CGeometry *geometry, CSolver **solver_container, CNumerics **numerics, CConfig *config);
  
  void Compute_StiffMatrix_NodalStressRes(CGeometry *geometry, CSolver **solver_container, CNumerics **numerics, CConfig *config);
  
  void Compute_MassMatrix(CGeometry *geometry, CSolver **solver_container, CNumerics **numerics, CConfig *config);
  
  void Compute_NodalStressRes(CGeometry *geometry, CSolver **solver_container, CNumerics **numerics, CConfig *config);
  
  void Compute_NodalStress(CGeometry *geometry, CSolver **solver_container, CNumerics **numerics, CConfig *config);
  
  void Compute_DeadLoad(CGeometry *geometry, CSolver **solver_container, CNumerics **numerics, CConfig *config);
  
  void Initialize_SystemMatrix(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void Compute_IntegrationConstants(CConfig *config);
  
  void BC_Clamped(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config, unsigned short val_marker);
  
  void BC_Clamped_Post(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                       unsigned short val_marker);
  
  void BC_Normal_Displacement(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                              unsigned short val_marker);
  
  
  void BC_Normal_Load(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                      unsigned short val_marker);
  
  void BC_Dir_Load(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                   unsigned short val_marker);
  
  void BC_Sine_Load(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                    unsigned short val_marker);
  
  
  void BC_Pressure(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                   unsigned short val_marker);
  
  void ImplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void ImplicitNewmark_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void ImplicitNewmark_Update(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void ImplicitNewmark_Relaxation(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void GeneralizedAlpha_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void GeneralizedAlpha_UpdateDisp(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void GeneralizedAlpha_UpdateSolution(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void GeneralizedAlpha_UpdateLoads(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void Postprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config,  CNumerics **numerics,
                      unsigned short iMesh);
  
  void Solve_System(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  su2double GetRes_FEM(unsigned short val_var);
  
  su2double GetTotal_CFEA(void);
  
  void SetTotal_CFEA(su2double val_cfea);
  
  void SetFEA_Load(CSolver ***flow_solution, CGeometry **fea_geometry, CGeometry **flow_geometry, CConfig *fea_config, CConfig *flow_config, CNumerics *fea_numerics);
  
  void SetFEA_Load_Int(CSolver ***flow_solution, CGeometry **fea_geometry, CGeometry **flow_geometry, CConfig *fea_config, CConfig *flow_config, CNumerics *fea_numerics);
  
  void PredictStruct_Displacement(CGeometry **fea_geometry,
                                  CConfig *fea_config,
                                  CSolver ***fea_solution);
  
  void ComputeAitken_Coefficient(CGeometry **fea_geometry,
                                 CConfig *fea_config,
                                 CSolver ***fea_solution,
                                 unsigned long iFSIIter);
  
  void SetAitken_Relaxation(CGeometry **fea_geometry,
                            CConfig *fea_config,
                            CSolver ***fea_solution);
  
  void Update_StructSolution(CGeometry **fea_geometry,
                             CConfig *fea_config,
                             CSolver ***fea_solution);
  
  void SetFSI_ConvValue(unsigned short val_index, su2double val_criteria);
  
  su2double GetFSI_ConvValue(unsigned short val_index);
  
  su2double GetWAitken_Dyn(void);
  
  su2double GetWAitken_Dyn_tn1(void);
  
  void SetWAitken_Dyn(su2double waitk);
  
  void SetWAitken_Dyn_tn1(su2double waitk_tn1);
  
  void SetLoad_Increment(su2double val_loadIncrement);
  
  void Set_Prestretch(CGeometry *geometry, CConfig *config);
};

class CTemplateSolver : public CSolver {
private:
  
public:
  
  CTemplateSolver(void);
  
  CTemplateSolver(CGeometry *geometry, CConfig *config);
  
  ~CTemplateSolver(void);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
  
  void SetTime_Step(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                    unsigned short iMesh, unsigned long Iteration);
  
  void Centered_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                         CConfig *config, unsigned short iMesh, unsigned short iRKStep);
  
  void Upwind_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                       CConfig *config, unsigned short iMesh);
  
  void Source_Residual(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CNumerics *second_numerics,
                       CConfig *config, unsigned short iMesh);
  
  void Source_Template(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                       CConfig *config, unsigned short iMesh);
  
  void BC_Euler_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config,
                     unsigned short val_marker);
  
  void BC_HeatFlux_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                        unsigned short val_marker);
  
  void BC_Far_Field(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                    unsigned short val_marker);
  
  void BC_Inlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                unsigned short val_marker);
  
  void BC_Outlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                 unsigned short val_marker);
  
  void BC_Sym_Plane(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                    unsigned short val_marker);
  
  void BC_Custom(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                 CConfig *config, unsigned short val_marker);
  
  void ExplicitRK_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                            unsigned short iRKStep);
  
  void ExplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
  void ImplicitEuler_Iteration(CGeometry *geometry, CSolver **solver_container, CConfig *config);
  
};

class CDiscAdjSolver : public CSolver {
private:
  unsigned short KindDirect_Solver;
  CSolver *direct_solver;
  su2double *Sens_Mach, 
  *Sens_AoA,      
  *Sens_Geo,      
  *Sens_Press,      
  *Sens_Temp,      
  **CSensitivity;  
  su2double Total_Sens_Mach;  
  su2double Total_Sens_AoA;    
  su2double Total_Sens_Geo;    
  su2double Total_Sens_Press;    
  su2double Total_Sens_Temp;    
  su2double Total_Sens_BPress;    
  su2double ObjFunc_Value;        
  su2double Mach, Alpha, Beta, Pressure, Temperature;
  unsigned long nMarker;        
public:
  
  CDiscAdjSolver(void);
  
  CDiscAdjSolver(CGeometry *geometry, CConfig *config);
  
  CDiscAdjSolver(CGeometry *geometry, CConfig *config, CSolver* solver, unsigned short Kind_Solver, unsigned short iMesh);
  
  ~CDiscAdjSolver(void);
  
  void RegisterSolution(CGeometry *geometry, CConfig *config);
  
  void RegisterOutput(CGeometry *geometry, CConfig *config);
  
  void SetAdjoint_Output(CGeometry *geometry, CConfig *config);
  
  void ExtractAdjoint_Solution(CGeometry *geometry, CConfig *config);
  
  void RegisterObj_Func(CConfig *config);
  
  void SetSurface_Sensitivity(CGeometry *geometry, CConfig* config);
  
  void SetSensitivity(CGeometry *geometry, CConfig *config);
  
  void SetAdj_ObjFunc(CGeometry *geometry, CConfig* config);
  
  
  su2double GetTotal_Sens_Geo(void);
  
  su2double GetTotal_Sens_Mach(void);
  
  su2double GetTotal_Sens_AoA(void);
  
  su2double GetTotal_Sens_Press(void);
  
  su2double GetTotal_Sens_Temp(void);
  
  su2double GetTotal_Sens_BPress(void);
  
  su2double GetCSensitivity(unsigned short val_marker, unsigned long val_vertex);
  
  void SetRecording(CGeometry *geometry, CConfig *config, unsigned short kind_recording);
  
  void RegisterVariables(CGeometry *geometry, CConfig *config, bool reset = false);
  
  void ExtractAdjoint_Variables(CGeometry *geometry, CConfig *config);
  
  void Preprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output);
};
#include "solver_structure.inl"