small fix

Author: pcarruscag

SU2_CFD/include/solvers/CFEASolver.hpp

@@ -527,10 +527,7 @@ class CFEASolver : public CSolver {
    * \param[in] numerics - Implementation of numerical method.
    * \param[in] of_comp_mode - Mode to compute just the objective function.
    */
-  inline void Postprocessing(CGeometry *geometry,
-                             CConfig *config,
-                             CNumerics **numerics,
-                             bool of_comp_mode) final;
+  void Postprocessing(CGeometry *geometry, CConfig *config, CNumerics **numerics, bool of_comp_mode) final;
 
   /*!
    * \brief Routine to solve the Jacobian-Residual linearized system.

SU2_CFD/src/solvers/CFEASolver.cpp

@@ -1916,14 +1916,12 @@ void CFEASolver::Postprocessing(CGeometry *geometry, CConfig *config, CNumerics
   /*--- Compute the objective function. ---*/
 
   const auto kindObjFunc = config->GetKind_ObjFunc();
-
-  if (((kindObjFunc == REFERENCE_GEOMETRY) || (kindObjFunc == REFERENCE_NODE)) &&
-      ((config->GetDV_FEA() == YOUNG_MODULUS) || (config->GetDV_FEA() == DENSITY_VAL))) {
-
-    Stiffness_Penalty(geometry, numerics, config);
-  }
+  const bool penalty = ((kindObjFunc == REFERENCE_GEOMETRY) || (kindObjFunc == REFERENCE_NODE)) &&
+                       ((config->GetDV_FEA() == YOUNG_MODULUS) || (config->GetDV_FEA() == DENSITY_VAL));
 
   if (of_comp_mode) {
+    if (penalty) Stiffness_Penalty(geometry, numerics, config);
+
     switch (kindObjFunc) {
       case REFERENCE_GEOMETRY: Compute_OFRefGeom(geometry, config); break;
       case REFERENCE_NODE:     Compute_OFRefNode(geometry, config); break;
@@ -1939,6 +1937,7 @@ void CFEASolver::Postprocessing(CGeometry *geometry, CConfig *config, CNumerics
     Compute_NodalStress(geometry, numerics, config);
 
     /*--- Compute functions for monitoring and output. ---*/
+    if (penalty) Stiffness_Penalty(geometry, numerics, config);
     Compute_OFRefNode(geometry, config);
     Compute_OFCompliance(geometry, config);
     if (config->GetRefGeom()) Compute_OFRefGeom(geometry, config);