little cleanup and add new option to config template

Author: pcarruscag

SU2_CFD/src/solvers/CSolver.cpp

@@ -385,12 +385,8 @@ void CSolver::InitiatePeriodicComms(CGeometry *geometry,
 
   /*--- Allocate buffers for matrices that need rotation. ---*/
 
-  su2double **jacBlock = new su2double*[ICOUNT];
-  su2double **rotBlock = new su2double*[ICOUNT];
-  for (iVar = 0; iVar < ICOUNT; iVar++) {
-    jacBlock[iVar] = new su2double[JCOUNT];
-    rotBlock[iVar] = new su2double[JCOUNT];
-  }
+  su2activematrix jacBlock(ICOUNT,JCOUNT);
+  su2activematrix rotBlock(ICOUNT,JCOUNT);
 
   /*--- Check to make sure we have created a large enough buffer
    for these comms during preprocessing. It will be reallocated whenever
@@ -1207,13 +1203,6 @@ void CSolver::InitiatePeriodicComms(CGeometry *geometry,
   delete [] rotPrim_i;
   delete [] rotPrim_j;
 
-  for (iVar = 0; iVar < ICOUNT; iVar++) {
-    delete [] jacBlock[iVar];
-    delete [] rotBlock[iVar];
-  }
-  delete [] jacBlock;
-  delete [] rotBlock;
-
 }
 
 void CSolver::CompletePeriodicComms(CGeometry *geometry,

config_template.cfg

@@ -1102,6 +1102,9 @@ TIME_DISCRE_ADJFLOW= EULER_IMPLICIT
 % Relaxation coefficient (also for discrete adjoint problems)
 RELAXATION_FACTOR_ADJOINT= 1.0
 %
+% Enable (if != 0) quasi-Newton acceleration/stabilization of discrete adjoints
+QUASI_NEWTON_NUM_SAMPLES= 20
+%
 % Reduction factor of the CFL coefficient in the adjoint problem
 CFL_REDUCTION_ADJFLOW= 0.8
 %