fixes and update ref file

Author: pcarruscag

SU2_CFD/include/solvers/CDiscAdjFEASolver.hpp

@@ -92,7 +92,7 @@ class CDiscAdjFEASolver final : public CSolver {
       }
     }
 
-    void UpdateTotal() {
+    void Store() {
       /*--- Clears the old values such that on the next time step the total is
        * incremented instead of updated. ---*/
       for (auto i = 0u; i < size; ++i) OldSens[i] = 0.0;

SU2_CFD/src/solvers/CDiscAdjFEASolver.cpp

@@ -346,21 +346,24 @@ void CDiscAdjFEASolver::Preprocessing(CGeometry *geometry, CSolver **solver_cont
 
 void CDiscAdjFEASolver::SetSensitivity(CGeometry *geometry, CConfig *config, CSolver*){
 
-  E.UpdateTotal();
-  Nu.UpdateTotal();
-  Rho.UpdateTotal();
-  Rho_DL.UpdateTotal();
-  if (de_effects) EField.UpdateTotal();
-  if (fea_dv) DV.UpdateTotal();
+  const bool time_domain = config->GetTime_Domain();
+
+  /*--- Store the final material sensitivities for the time step to increment them on the next time step. ---*/
+  if (time_domain) {
+    E.Store();
+    Nu.Store();
+    Rho.Store();
+    Rho_DL.Store();
+    if (de_effects) EField.Store();
+    if (fea_dv) DV.Store();
+  }
 
   /*--- Extract the topology optimization density sensitivities. ---*/
 
   direct_solver->ExtractAdjoint_Variables(geometry, config);
 
   /*--- Extract the geometric sensitivities ---*/
 
-  const bool time_domain = config->GetTime_Domain();
-
   for (unsigned long iPoint = 0; iPoint < nPoint; iPoint++) {
 
     auto Coord = geometry->nodes->GetCoord(iPoint);

TestCases/disc_adj_fsi/dyn_fsi/grad_dv.opt.ref

@@ -1,9 +1,9 @@
 INDEX	GRAD
-0	-3.461460667601001e-03
-1	-1.841786311588662e-03
-2	-7.915536257748965e-04
-3	-2.739622082729719e-04
-4	-2.734869133461108e-04
-5	-7.881162428890194e-04
-6	-1.828978290516676e-03
-7	-3.427219398258311e-03
+0	-3.461460667601000e-03
+1	-1.841786311588663e-03
+2	-7.915536257748967e-04
+3	-2.739622082729717e-04
+4	-2.734869133461104e-04
+5	-7.881162428890206e-04
+6	-1.828978290516677e-03
+7	-3.427219398258316e-03