Skip to content
This repository was archived by the owner on Jul 19, 2023. It is now read-only.

Commit e5a3da2

Browse files
committed
removed unnessecary boundary point
1 parent bfed2d8 commit e5a3da2

1 file changed

Lines changed: 1 addition & 1 deletion

File tree

src/derivative_operators/derivative_operator.jl

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -221,7 +221,7 @@ function CompleteHalfCenteredDifference(derivative_order::Int,
221221
right_boundary_x = reverse(-boundary_stencil_length:-1)
222222

223223
# ? Is fornberg valid when taking an x0 outside of the stencil i.e at the boundary?
224-
xoffset = 1.5:boundary_stencil_length+0.5
224+
xoffset = 1.5:boundary_stencil_length-0.5
225225

226226
boundary_point_count = div(centered_stencil_length,2) # -1 due to the ghost point
227227
# Because it's a N x (N+2) operator, the last stencil on the sides are the [b,0,x,x,x,x] stencils, not the [0,x,x,x,x,x] stencils, since we're never solving for the derivative at the boundary point.

0 commit comments

Comments
 (0)