From a97be463aab83838623ec696bed4bba684fcf386 Mon Sep 17 00:00:00 2001
From: Sean Engelstad <sean.engelstad1@gmail.com>
Date: Tue, 8 Oct 2024 19:03:05 -0400
Subject: [PATCH] fix introduced shear loads bug

---
 src/constitutive/TACSGPBladeStiffenedShellConstitutive.cpp | 7 ++-----
 1 file changed, 2 insertions(+), 5 deletions(-)

diff --git a/src/constitutive/TACSGPBladeStiffenedShellConstitutive.cpp b/src/constitutive/TACSGPBladeStiffenedShellConstitutive.cpp
index 5ea989e3..ae605434 100644
--- a/src/constitutive/TACSGPBladeStiffenedShellConstitutive.cpp
+++ b/src/constitutive/TACSGPBladeStiffenedShellConstitutive.cpp
@@ -2150,9 +2150,7 @@ TACSGPBladeStiffenedShellConstitutive::computeCriticalGlobalShearLoadSens(
     TacsScalar shear_geom_sens[2];
     ksAggregationSens(
         shear_geom, 2, this->ksWeight, shear_geom_sens);
-    TacsScalar sg2_sens = shear_geom_sens[1];
-    TacsScalar c_rho0_sens = -2.0 * sg2_sens / rho_0;
-
+    TacsScalar c_rho0_sens = shear_geom_sens[1] * shear_geom[1] * -2.0 / rho_0;
 
     // compute the overall sensitivities
     *D11sens += N12sens * N12crit * 0.25 / D11;
@@ -2352,8 +2350,7 @@ TACSGPBladeStiffenedShellConstitutive::computeCriticalLocalShearLoadSens(
     TacsScalar shear_geom_sens[2];
     ksAggregationSens(
         shear_geom, 2, this->ksWeight, shear_geom_sens);
-    TacsScalar sg2_sens = shear_geom_sens[1];
-    TacsScalar c_rho0_sens = -2.0 * sg2_sens / rho_0;    
+    TacsScalar c_rho0_sens = shear_geom_sens[1] * shear_geom[1] * -2.0 / rho_0; 
 
     // compute the overall sensitivities
     *D11sens += N12sens * N12crit * 0.25 / D11;