Collection of improvements

Author: Francisco Palacios

Common/include/geometry_structure.hpp

@@ -591,7 +591,7 @@ class CGeometry {
 	 * \brief A virtual member.
 	 * \param[in] config - Definition of the particular problem.
 	 */
-	void ComputeAirfoil_Section(double *Plane_P0, double *Plane_Normal, unsigned short iSection,
+	void ComputeAirfoil_Section(double *Plane_P0, double *Plane_Normal,
                                       double MinXCoord, double MaxXCoord, double *FlowVariable,
                                       vector<double> &Xcoord_Airfoil, vector<double> &Ycoord_Airfoil,
                                       vector<double> &Zcoord_Airfoil, vector<double> &Variable_Airfoil,

Common/include/grid_adaptation_structure.hpp

@@ -141,13 +141,6 @@ class CGridAdaptation {
 	 */		
 	void SetWake_Refinement(CGeometry *geometry, unsigned short strength);
 
-	/*! 
-	 * \brief Do an adaptation of the computational grid on the two phase problem interphase.
-	 * \param[in] geometry - Geometrical definition of the problem.
-	 * \param[in] strength - Adaptation Strength.	 
-	 */		
-	void SetTwoPhase_Refinement(CGeometry *geometry, unsigned short strength);
-	
 	/*! 
 	 * \brief Do an adaptation of the computational grid on the supersonic shock region.
 	 * \param[in] geometry - Geometrical definition of the problem.

Common/include/grid_movement_structure.hpp

@@ -529,14 +529,6 @@ class CFreeFormDefBox : public CGridMovement {
 	 * \return Value of the Derivative of the Bernstein polynomial.
 	 */		
 	double GetBernsteinDerivative(short val_n, short val_i, double val_t, short val_order);
-	
-  /*!
-	 * \brief The routine computes F(u,v,w)=||X(u,v,w)-(x,y,z)||^2  evaluated at (u,v,w).
-	 * \param[in] val_coord - Parametric coordiates of the target point.
-	 * \param[in] xyz - Cartesians coordinates of the point.
-	 * \return Value of the analytical objective function.
-	 */
-	double GetFFDObjFunc(double *val_coord, double *xyz);
 
 	/*! 
 	 * \brief The routine computes the gradient of F(u,v,w)=||X(u,v,w)-(x,y,z)||^2  evaluated at (u,v,w).

Common/include/option_structure.hpp

@@ -863,7 +863,6 @@ enum ENUM_ADAPT {
   WAKE = 12,			/*!< \brief Do a grid refinement on the wake. */
   SMOOTHING = 14,		/*!< \brief Do a grid smoothing of the geometry. */
   SUPERSONIC_SHOCK = 15,	/*!< \brief Do a grid smoothing. */
-  TWOPHASE = 16,			/*!< \brief Do a grid refinement on the free surface interphase. */
   PERIODIC = 17			/*!< \brief Add the periodic halo cells. */
 };
 static const map<string, ENUM_ADAPT> Adapt_Map = CCreateMap<string, ENUM_ADAPT>
@@ -882,8 +881,7 @@ static const map<string, ENUM_ADAPT> Adapt_Map = CCreateMap<string, ENUM_ADAPT>
 ("WAKE", WAKE)
 ("SMOOTHING", SMOOTHING)
 ("SUPERSONIC_SHOCK", SUPERSONIC_SHOCK)
-("PERIODIC", PERIODIC)
-("TWOPHASE", TWOPHASE);
+("PERIODIC", PERIODIC);
 
 /*!
  * \brief types of input file formats

Common/src/config_structure.cpp

@@ -1111,7 +1111,6 @@ void CConfig::SetConfig_Options(unsigned short val_iZone, unsigned short val_nZo
   /* DESCRIPTION: Evaluate inverse design on the surface  */
   addBoolOption("INV_DESIGN_HEATFLUX", InvDesign_HeatFlux, false);
 
-
   /* CONFIG_CATEGORY: Unsupported options */
   /*--- Options that are experimental and not intended for general use ---*/
 
@@ -1148,6 +1147,15 @@ void CConfig::SetConfig_Options(unsigned short val_iZone, unsigned short val_nZo
   /* DESCRIPTION: Setup for design variables */
   addPythonOption("DEFINITION_DV");
 
+  /* DESCRIPTION: Maximum number of iterations */
+  addPythonOption("OPT_ITERATIONS");
+  
+  /* DESCRIPTION: Requested accuracy */
+  addPythonOption("OPT_ACCURACY");
+  
+  /* DESCRIPTION: Setup for design variables */
+  addPythonOption("BOUND_DV");
+  
   /* DESCRIPTION: Current value of the design variables */
   addPythonOption("DV_VALUE_NEW");
 
@@ -3387,7 +3395,7 @@ void CConfig::SetOutput(unsigned short val_software, unsigned short val_izone) {
 
 	if (val_software == SU2_MSH) {
 		switch (Kind_Adaptation) {
-		case FULL: case WAKE: case TWOPHASE: case FULL_FLOW: case FULL_ADJOINT: case FULL_LINEAR: case SMOOTHING: case SUPERSONIC_SHOCK:
+		case FULL: case WAKE: case FULL_FLOW: case FULL_ADJOINT: case FULL_LINEAR: case SMOOTHING: case SUPERSONIC_SHOCK:
 			break;
 		case GRAD_FLOW:
 			cout << "Read flow solution from: " << Solution_FlowFileName << "." << endl;
@@ -4083,7 +4091,6 @@ void CConfig::SetOutput(unsigned short val_software, unsigned short val_izone) {
       case PERIODIC: cout << "Grid modification to run periodic bc problems." << endl; break;
       case FULL: cout << "Grid adaptation using a complete refinement." << endl; break;
       case WAKE: cout << "Grid adaptation of the wake." << endl; break;
-      case TWOPHASE: cout << "Grid adaptation of the interphase of a free surface flow." << endl; break;
       case FULL_FLOW: cout << "Flow grid adaptation using a complete refinement." << endl; break;
       case FULL_ADJOINT: cout << "Adjoint grid adaptation using a complete refinement." << endl; break;
       case FULL_LINEAR: cout << "Linear grid adaptation using a complete refinement." << endl; break;

Common/src/dual_grid_structure.cpp

@@ -453,7 +453,7 @@ CVertex::~CVertex() {
 }
 
 void CVertex::SetNodes_Coord(double *val_coord_Edge_CG, double *val_coord_FaceElem_CG, double *val_coord_Elem_CG) {
-	double vec_a[3], vec_b[3], Dim_Normal[3];
+  double vec_a[3] = {0.0,0.0,0.0}, vec_b[3] = {0.0,0.0,0.0}, Dim_Normal[3] = {0.0,0.0,0.0};
 	unsigned short iDim;
 
 	for (iDim = 0; iDim < nDim; iDim++) {