corrected error in deMandeling of jacoby matrix
This commit is contained in:
parent
09c47f15e3
commit
6f3f5ee23c
|
@ -178,7 +178,7 @@
|
||||||
cp_en = mesh_FEasCP('elem', n(1))
|
cp_en = mesh_FEasCP('elem', n(1))
|
||||||
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,1)*invnrmMandel(1:ngens)
|
forall(i=1:ngens) d(1:ngens,i)=d(1:ngens,i)*invnrmMandel(1:ngens)
|
||||||
return
|
return
|
||||||
|
|
||||||
END SUBROUTINE
|
END SUBROUTINE
|
||||||
|
|
Loading…
Reference in New Issue