fix multizone memory errors

Author: pcarruscag

SU2_CFD/include/drivers/CDriverBase.hpp

@@ -62,8 +62,8 @@ class CDriverBase {
       *nInst,            /*!< \brief Total number of instances in the problem (per zone). */
       **interface_types; /*!< \brief Type of coupling between the distinct (physical) zones. */
 
-  CConfig* driver_config; /*!< \brief Definition of the driver configuration. */
-  COutput* driver_output; /*!< \brief Definition of the driver output. */
+  CConfig* driver_config = nullptr; /*!< \brief Definition of the driver configuration. */
+  COutput* driver_output = nullptr; /*!< \brief Definition of the driver output. */
 
   CConfig** config_container;           /*!< \brief Definition of the particular problem. */
   COutput** output_container;           /*!< \brief Pointer to the COutput class. */
@@ -74,8 +74,8 @@ class CDriverBase {
   CSurfaceMovement** surface_movement;  /*!< \brief Surface movement classes of the problem. */
   CVolumetricMovement*** grid_movement; /*!< \brief Volume grid movement classes of the problem. */
 
-  CConfig* main_config;     /*!< \brief Reference to base (i.e. ZONE 0) configuration (used in driver API). */
-  CGeometry* main_geometry; /*!< \brief Reference to base (i.e. ZONE, INST, MESH 0) geometry (used in driver API). */
+  CConfig* main_config = nullptr;     /*!< \brief Reference to base (i.e. ZONE 0) configuration (used in driver API). */
+  CGeometry* main_geometry = nullptr; /*!< \brief Reference to base (i.e. ZONE, INST, MESH 0) geometry (used in driver API). */
 
  public:
   /*!

SU2_CFD/src/drivers/CDriverBase.cpp

@@ -55,12 +55,6 @@ void CDriverBase::SetContainers_Null() {
   for (iZone = 0; iZone < nZone; iZone++) {
     nInst[iZone] = 1;
   }
-
-  driver_config = nullptr;
-  driver_output = nullptr;
-
-  main_config = nullptr;
-  main_geometry = nullptr;
 }
 
 void CDriverBase::CommonPostprocessing() {

SU2_DEF/src/drivers/CDeformationDriver.cpp

@@ -49,10 +49,6 @@ CDeformationDriver::CDeformationDriver(char* confFile, SU2_Comm MPICommunicator)
   rank = SU2_MPI::GetRank();
   size = SU2_MPI::GetSize();
 
-  /*--- Initialize containers. --- */
-
-  SetContainers_Null();
-
   /*--- Preprocessing of the config files. ---*/
 
   Input_Preprocessing();
@@ -102,6 +98,10 @@ void CDeformationDriver::Input_Preprocessing() {
 
   nZone = driver_config->GetnZone();
 
+  /*--- Initialize containers. --- */
+
+  SetContainers_Null();
+
   /*--- Loop over all zones to initialize the various classes. In most
    cases, nZone is equal to one. This represents the solution of a partial
    differential equation on a single block, unstructured mesh. ---*/

SU2_DEF/src/drivers/CDiscAdjDeformationDriver.cpp

@@ -55,10 +55,6 @@ CDiscAdjDeformationDriver::CDiscAdjDeformationDriver(char* confFile, SU2_Comm MP
   rank = SU2_MPI::GetRank();
   size = SU2_MPI::GetSize();
 
-  /*--- Initialize containers. --- */
-
-  SetContainers_Null();
-
   /*--- Preprocessing of the config files. ---*/
 
   Input_Preprocessing();
@@ -109,6 +105,10 @@ void CDiscAdjDeformationDriver::Input_Preprocessing() {
 
   nZone = driver_config->GetnZone();
 
+  /*--- Initialize containers. --- */
+
+  SetContainers_Null();
+
   /*--- Loop over all zones to initialize the various classes. In most
    * cases, nZone is equal to one. This represents the solution of a partial
    * differential equation on a single block, unstructured mesh. ---*/

TestCases/incomp_navierstokes/streamwise_periodic/chtPinArray_2d/configMaster.cfg

@@ -69,7 +69,7 @@ DEFORM_LINEAR_SOLVER_PREC= ILU
 DEFORM_LINEAR_SOLVER_ERROR= 1E-14
 %
 DEFORM_NONLINEAR_ITER= 1
-DEFORM_LINEAR_SOLVER_ITER= 1000
+DEFORM_LINEAR_SOLVER_ITER= 10
 %
 DEFORM_CONSOLE_OUTPUT= YES
 DEFORM_STIFFNESS_TYPE= WALL_DISTANCE