Applying penalty stiffness to initial tangent

Author: mhscott

SRC/material/uniaxial/PenaltyMaterial.cpp

@@ -163,6 +163,15 @@ PenaltyMaterial::getTangent(void)
     return 0.0;
 }
 
+double 
+PenaltyMaterial::getInitialTangent(void)
+{
+  if (theMaterial)
+    return theMaterial->getInitialTangent() + penalty;
+  else
+    return 0.0;
+}
+
 double 
 PenaltyMaterial::getDampTangent(void)
 {

SRC/material/uniaxial/PenaltyMaterial.h

@@ -50,7 +50,7 @@ class PenaltyMaterial : public UniaxialMaterial
     double getStress(void);
     double getTangent(void);
     double getDampTangent(void);
-    double getInitialTangent(void) {return theMaterial->getInitialTangent();}
+    double getInitialTangent(void);
 
     int commitState(void);
     int revertToLastCommit(void);