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

Commit

Permalink
a bracket
Browse files Browse the repository at this point in the history
  • Loading branch information
xtalax committed Jan 5, 2022
1 parent 85b8be8 commit 5a5792c
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion src/derivative_operators/derivative_operator.jl
Original file line number Diff line number Diff line change
Expand Up @@ -184,7 +184,7 @@ function CompleteCenteredDifference{N}(derivative_order::Int,

stencil_coefs = convert(SVector{stencil_length, T}, (1/dx^derivative_order) * calculate_weights(derivative_order, zero(T), dummy_x))
_low_boundary_coefs = SVector{boundary_stencil_length, T}[convert(SVector{boundary_stencil_length, T}, (1/dx^derivative_order) * calculate_weights(derivative_order, oneunit(T)*x0, left_boundary_x)) for x0 in L_boundary_deriv_spots]
low_boundary_coefs = convert(SVector{boundary_point_count},vcat(_low_boundary_coefs)
low_boundary_coefs = convert(SVector{boundary_point_count},vcat(_low_boundary_coefs))

# _high_boundary_coefs = SVector{boundary_stencil_length, T}[convert(SVector{boundary_stencil_length, T}, (1/dx^derivative_order) * calculate_weights(derivative_order, oneunit(T)*x0, reverse(right_boundary_x))) for x0 in R_boundary_deriv_spots]
high_boundary_coefs = convert(SVector{boundary_point_count},reverse(map(reverse, _low_boundary_coefs*(-1)^derivative_order)))
Expand Down

0 comments on commit 5a5792c

Please sign in to comment.