mandel fun at end of hypela2

This commit is contained in:
Philip Eisenlohr 2007-03-28 13:48:33 +00:00
parent 9b96eb9984
commit 99eb77b3dc
1 changed files with 6 additions and 4 deletions

View File

@ -177,10 +177,12 @@
! Marc: 11, 22, 33, 12, 23, 13 ! Marc: 11, 22, 33, 12, 23, 13
s(1:ngens)=invnrmMandel(1:ngens)*CPFEM_stress_all(1:ngens, nn, cp_en) 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) 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) 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(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 return
end end