From 2c80880bd5240d6ed57feade2787e8de8cb7e341 Mon Sep 17 00:00:00 2001
From: Wenqing Wang <wenqing.wang@ufz.de>
Date: Thu, 19 Jul 2018 10:54:19 +0200
Subject: [PATCH] [BGRa] Removed one std::pow
---
MaterialLib/SolidModels/CreepBGRa.cpp | 7 ++++---
1 file changed, 4 insertions(+), 3 deletions(-)
diff --git a/MaterialLib/SolidModels/CreepBGRa.cpp b/MaterialLib/SolidModels/CreepBGRa.cpp
index 598d3998dae..166ca3c3d13 100644
--- a/MaterialLib/SolidModels/CreepBGRa.cpp
+++ b/MaterialLib/SolidModels/CreepBGRa.cpp
@@ -30,7 +30,8 @@ CreepBGRa<DisplacementDim>::integrateStress(
KelvinVector const& eps_prev, KelvinVector const& eps,
KelvinVector const& sigma_prev,
typename MechanicsBase<DisplacementDim>::MaterialStateVariables const&
- /*material_state_variables*/, double const T)
+ /*material_state_variables*/,
+ double const T)
{
using Invariants = MathLib::KelvinVector::Invariants<KelvinVectorSize>;
@@ -58,8 +59,8 @@ CreepBGRa<DisplacementDim>::integrateStress(
jacobian =
KelvinMatrix::Identity() +
pow_norm_s_n1_n_minus_one_2b_G *
- (deviatoric_matrix +
- (_n - 1) * std::pow(norm_s_n1, -2) * s_n1 * s_n1.transpose());
+ (deviatoric_matrix + ((_n - 1) / (norm_s_n1 * norm_s_n1)) *
+ s_n1 * s_n1.transpose());
};
auto const update_residual = [&](ResidualVectorType& r) {
--
GitLab