gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

MESH: Change summation order in restrict_block().

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2013-12-19 10:18:06 -02:00
parent 2a8eee95c0
commit b21b7f99d1
1 changed files with 12 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
src/mesh.F90
View File

@ -1318,21 +1318,21 @@ module mesh
! !
#if NDIMS == 2 #if NDIMS == 2
pparent%u(1:nv,is:it,js:jt, 1 ) = & pparent%u(1:nv,is:it,js:jt, 1 ) = &
2.5d-01 * ((pchild%u(1:nv,il:iu:2,jl:ju:2,1) & 2.50d-01 * ((pchild%u(1:nv,il:iu:2,jl:ju:2, 1 ) &
+ pchild%u(1:nv,ip:iu:2,jp:ju:2,1)) & + pchild%u(1:nv,ip:iu:2,jp:ju:2, 1 )) &
+ (pchild%u(1:nv,ip:iu:2,jl:ju:2,1) & + (pchild%u(1:nv,il:iu:2,jp:ju:2, 1 ) &
+ pchild%u(1:nv,il:iu:2,jp:ju:2,1))) + pchild%u(1:nv,ip:iu:2,jl:ju:2, 1 )))
#endif /* NDIMS == 2 */ #endif /* NDIMS == 2 */
#if NDIMS == 3 #if NDIMS == 3
pparent%u(1:nv,is:it,js:jt,ks:kt) = & pparent%u(1:nv,is:it,js:jt,ks:kt) = &
1.25d-01 * ((pchild%u(1:nv,il:iu:2,jl:ju:2,kl:ku:2) & 1.25d-01 * (((pchild%u(1:nv,il:iu:2,jl:ju:2,kl:ku:2) &
+ pchild%u(1:nv,ip:iu:2,jp:ju:2,kp:ku:2)) & + pchild%u(1:nv,ip:iu:2,jp:ju:2,kp:ku:2)) &
+ (pchild%u(1:nv,il:iu:2,jl:ju:2,kp:ku:2) & + (pchild%u(1:nv,il:iu:2,jl:ju:2,kp:ku:2) &
+ pchild%u(1:nv,ip:iu:2,jp:ju:2,kl:ku:2)) & + pchild%u(1:nv,ip:iu:2,jp:ju:2,kl:ku:2))) &
+ (pchild%u(1:nv,ip:iu:2,jl:ju:2,kl:ku:2) & + ((pchild%u(1:nv,il:iu:2,jp:ju:2,kp:ku:2) &
+ pchild%u(1:nv,il:iu:2,jp:ju:2,kp:ku:2)) & + pchild%u(1:nv,ip:iu:2,jl:ju:2,kl:ku:2)) &
+ (pchild%u(1:nv,ip:iu:2,jl:ju:2,kp:ku:2) & + (pchild%u(1:nv,il:iu:2,jp:ju:2,kl:ku:2) &
+ pchild%u(1:nv,il:iu:2,jp:ju:2,kl:ku:2))) + pchild%u(1:nv,ip:iu:2,jl:ju:2,kp:ku:2))))
#endif /* NDIMS == 3 */ #endif /* NDIMS == 3 */
end do ! p = 1, nchild end do ! p = 1, nchild