From d5d71d2321c4ae766d5f1c6cc1e50dbd0881a3b7 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Sun, 25 Feb 2018 13:42:02 +0100 Subject: [PATCH] math_transpose33 has no benefit here --- src/spectral_utilities.f90 | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/src/spectral_utilities.f90 b/src/spectral_utilities.f90 index c0863cebb..1a277ef63 100644 --- a/src/spectral_utilities.f90 +++ b/src/spectral_utilities.f90 @@ -935,7 +935,6 @@ subroutine utilities_constitutiveResponse(P,P_av,C_volAvg,C_minmaxAvg,& debug_reset, & debug_info use math, only: & - math_transpose33, & math_rotate_forward33, & math_det33 use mesh, only: & @@ -994,10 +993,10 @@ subroutine utilities_constitutiveResponse(P,P_av,C_volAvg,C_minmaxAvg,& call MPI_Allreduce(MPI_IN_PLACE,P_av,9,MPI_DOUBLE,MPI_SUM,PETSC_COMM_WORLD,ierr) if (debugRotation) & write(6,'(/,a,/,3(3(2x,f12.4,1x)/))',advance='no') ' Piola--Kirchhoff stress (lab) / MPa =',& - math_transpose33(P_av)*1.e-6_pReal + transpose(P_av)*1.e-6_pReal P_av = math_rotate_forward33(P_av,rotation_BC) write(6,'(/,a,/,3(3(2x,f12.4,1x)/))',advance='no') ' Piola--Kirchhoff stress / MPa =',& - math_transpose33(P_av)*1.e-6_pReal + transpose(P_av)*1.e-6_pReal flush(6) max_dPdF = 0.0_pReal