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

[WIP] Fixing Convolutions #104

Merged
merged 1 commit into from
Jun 10, 2019
Merged

[WIP] Fixing Convolutions #104

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
11 changes: 5 additions & 6 deletions src/derivative_operators/convolutions.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -22,7 +22,7 @@ function convolve_interior!(x_temp::AbstractVector{T}, x::AbstractVector{T}, A::
end
end

function convolve_BC_left!(x_temp::AbstractVector{T}, x::AbstractVector, A::DerivativeOperator{T,S}) where {T<:Real,S<:SVector}
function convolve_BC_left!(x_temp::AbstractVector{T}, x::AbstractVector{T}, A::DerivativeOperator{T,S}) where {T<:Real,S<:SVector}
coeffs = A.low_boundary_coefs
Threads.@threads for i in 1 : A.boundary_length
xtempi = coeffs[i][1]*x[1]
Expand All @@ -33,15 +33,14 @@ function convolve_BC_left!(x_temp::AbstractVector{T}, x::AbstractVector, A::Deri
end
end

function convolve_BC_right!(x_temp::AbstractVector{T}, x::AbstractVector, A::DerivativeOperator{T,S}) where {T<:Real,S<:SVector}
coeffs = A.low_boundary_coefs
bc_start = length(x) - A.stencil_length
function convolve_BC_right!(x_temp::AbstractVector{T}, x::AbstractVector{T}, A::DerivativeOperator{T,S}) where {T<:Real,S<:SVector}
coeffs = A.high_boundary_coefs
Threads.@threads for i in 1 : A.boundary_length
xtempi = coeffs[i][end]*x[end]
@inbounds for idx in A.stencil_length:-1:2
@inbounds for idx in A.stencil_length:-1:1
xtempi += coeffs[i][end-idx] * x[end-idx]
end
x_temp[i] = xtempi
x_temp[end-A.boundary_length+i] = xtempi
end
end

Expand Down