From 99eb77b3dc96e4ad38483fe77afddafc38480589 Mon Sep 17 00:00:00 2001 From: Philip Eisenlohr Date: Wed, 28 Mar 2007 13:48:33 +0000 Subject: [PATCH] mandel fun at end of hypela2 --- trunk/mpie_cpfem_marc.f90 | 10 ++++++---- 1 file changed, 6 insertions(+), 4 deletions(-) diff --git a/trunk/mpie_cpfem_marc.f90 b/trunk/mpie_cpfem_marc.f90 index c3becf113..3de2d8c8c 100644 --- a/trunk/mpie_cpfem_marc.f90 +++ b/trunk/mpie_cpfem_marc.f90 @@ -177,10 +177,12 @@ ! Marc: 11, 22, 33, 12, 23, 13 s(1:ngens)=invnrmMandel(1:ngens)*CPFEM_stress_all(1:ngens, nn, cp_en) d(1:ngens,1:ngens)=CPFEM_jaco_old(1:ngens,1:ngens, nn, cp_en) - forall(i=1:ngens) d(1:ngens,i)=d(1:ngens,i)*invnrmMandel(1:ngens) - d(1:ngens,1:ngens)=transpose(d(1:ngens,1:ngens)) - forall(i=1:ngens) d(1:ngens,i)=d(1:ngens,i)*nrmMandel(1:ngens) - d(1:ngens,1:ngens)=transpose(d(1:ngens,1:ngens)) + forall(i=1:ngens) d(1:ngens,i)=d(1:ngens,i)*nrmMandel(1:ngens) + forall(i=1:ngens) d(i,1:ngens)=d(i,1:ngens)*invnrmMandel(1:ngens) + + !d(1:ngens,1:ngens)=transpose(d(1:ngens,1:ngens)) + !forall(i=1:ngens) d(1:ngens,i)=d(1:ngens,i)*invnrmMandel(1:ngens) + !d(1:ngens,1:ngens)=transpose(d(1:ngens,1:ngens)) return end