223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

fixed bug in divergenceFDM

This commit is contained in:
Martin Diehl 2012-08-27 15:43:17 +00:00
parent 9214da4465
commit efd1b513cf
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
code/math.f90
View File

@ -3212,7 +3212,7 @@ function math_divergenceFDM(geomdim,order,field)
+ (/1_pInt,1_pInt,1_pInt/) + (/1_pInt,1_pInt,1_pInt/)
do l = 1_pInt, vec_tens do l = 1_pInt, vec_tens
math_divergenceFDM(i+1_pInt,j+1_pInt,k+1_pInt,l) = math_divergenceFDM(i+1_pInt,j+1_pInt,k+1_pInt,l) & math_divergenceFDM(i+1_pInt,j+1_pInt,k+1_pInt,l) = math_divergenceFDM(i+1_pInt,j+1_pInt,k+1_pInt,l) &
+ FDcoefficient(m,order) * & + FDcoefficient(m,order+1_pInt) * &
((field(coordinates(1,1),coordinates(1,2),coordinates(1,3),l,1)- & ((field(coordinates(1,1),coordinates(1,2),coordinates(1,3),l,1)- &
field(coordinates(2,1),coordinates(2,2),coordinates(2,3),l,1))*real(res(1),pReal)/geomdim(1) +& field(coordinates(2,1),coordinates(2,2),coordinates(2,3),l,1))*real(res(1),pReal)/geomdim(1) +&
(field(coordinates(3,1),coordinates(3,2),coordinates(3,3),l,2)- & (field(coordinates(3,1),coordinates(3,2),coordinates(3,3),l,2)- &