Skip to content

Commit ecb1743

Browse files
committed
fix off-by-one
1 parent d2a2a89 commit ecb1743

File tree

1 file changed

+2
-2
lines changed

1 file changed

+2
-2
lines changed

Src/Base/AMReX_MultiFabUtil.cpp

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -845,12 +845,12 @@ namespace amrex
845845
{
846846
// For node (i,j,k), compute curl using face values that touch this node
847847
// curl_x = dw/dy - dv/dz
848-
curlarr(i,j,k,0) = dxinv[1] * (warr(i-1,j,k-1) - warr(i-1,j-1,k-1))
848+
curlarr(i,j,k,0) = dxinv[1] * (warr(i-1,j,k) - warr(i-1,j-1,k))
849849
- dxinv[2] * (varr(i-1,j-1,k) - varr(i-1,j-1,k-1));
850850

851851
// curl_y = du/dz - dw/dx
852852
curlarr(i,j,k,1) = dxinv[2] * (uarr(i-1,j-1,k) - uarr(i-1,j-1,k-1))
853-
- dxinv[0] * (warr(i,j-1,k-1) - warr(i-1,j-1,k-1));
853+
- dxinv[0] * (warr(i,j-1,k) - warr(i-1,j-1,k));
854854

855855
// curl_z = dv/dx - du/dy
856856
curlarr(i,j,k,2) = dxinv[0] * (varr(i,j-1,k-1) - varr(i-1,j-1,k-1))

0 commit comments

Comments
 (0)