From 27dafa185b81eeddebf58a0ea3ae93b59a109437 Mon Sep 17 00:00:00 2001 From: Jean-Paul Pelteret Date: Mon, 23 Apr 2018 09:54:45 +0200 Subject: [PATCH] Use function from Physics module instead of hardcoded equivalent Since the addition of fixes to the AD number support, all of the functions in the Physics module now work with all of our supported AD number types. Previously, this line would not compile for nested Sacado Rad-DFad types. --- .../cook_membrane.cc | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/Quasi_static_Finite_strain_Compressible_Elasticity/cook_membrane.cc b/Quasi_static_Finite_strain_Compressible_Elasticity/cook_membrane.cc index 54a5d7a..59de4a4 100644 --- a/Quasi_static_Finite_strain_Compressible_Elasticity/cook_membrane.cc +++ b/Quasi_static_Finite_strain_Compressible_Elasticity/cook_membrane.cc @@ -1896,8 +1896,7 @@ Point grid_y_transform (const Point &pt_in) const Tensor<2,dim,ADNumberType> &grad_u = scratch.solution_grads_u_total[q_point]; const Tensor<2,dim,ADNumberType> F = Physics::Elasticity::Kinematics::F(grad_u); const ADNumberType det_F = determinant(F); -// const Tensor<2,dim,ADNumberType> F_bar = Physics::Elasticity::Kinematics::F_iso(F); - const Tensor<2,dim,ADNumberType> F_bar = ADNumberType(std::pow(determinant(F),-1.0/dim))*F; + const Tensor<2,dim,ADNumberType> F_bar = Physics::Elasticity::Kinematics::F_iso(F); const SymmetricTensor<2,dim,ADNumberType> b_bar = Physics::Elasticity::Kinematics::b(F_bar); const Tensor<2,dim,ADNumberType> F_inv = invert(F); Assert(det_F > ADNumberType(0.0), ExcInternalError()); -- 2.39.5