Fix index handling in 1D (#39)

PTsolvers · Aug 28, 2024 · 6e3ae01 · 6e3ae01
1 parent 8d2f553
commit 6e3ae01
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 2 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "Chmy"
 uuid = "33a72cf0-4690-46d7-b987-06506c2248b9"
 authors = ["Ivan Utkin <[email protected]>, Ludovic Raess <[email protected]>, and contributors"]
-version = "0.1.18"
+version = "0.1.19"
 
 [deps]
 Adapt = "79e6a3ab-5dfb-504d-930d-738a2a938a0e"

diff --git a/src/BoundaryConditions/batch.jl b/src/BoundaryConditions/batch.jl
@@ -156,13 +156,16 @@ end
 
 bc!(arch::Architecture, grid::SG, f_bc::Vararg{FieldAndBC}; kwargs...) = bc!(arch, grid, batch(grid, f_bc...; kwargs...))
 
+@inline to_bc_index(::SG, I)    = I .- 1
+@inline to_bc_index(::SG{1}, I) = ()
+
 # batched kernels
 @kernel function bc_kernel!(side::Side, dim::Dim,
                             grid::SG{N},
                             fields::NTuple{K,Field},
                             conditions::NTuple{K,FieldBoundaryCondition}) where {N,K}
     J = @index(Global, NTuple)
-    I = J .- 1
+    I = to_bc_index(grid, J)
     ntuple(Val(K)) do ifield
         Base.@_inline_meta
         @inbounds begin