Dirichlet BC don't work in 1D

[boriskaus]

MWE:

grid   = UniformGrid(arch; origin=(0, ), extent=(2, ), dims=(11,))
u = Field(arch, grid, Center())
bc!(arch, grid, u => Dirichlet() )

[luraess]

Seems that I_f computed in

Chmy.jl/src/BoundaryConditions/first_order_boundary_condition.jl

Line 60 in 8ba91ea

I_f = insert_dim(dim, I, itp_halo_index(side, dim, f))

returns a Tuple of size N=2 in 1D which then gets passed to

Chmy.jl/src/BoundaryConditions/first_order_boundary_condition.jl

Line 67 in 8ba91ea

a = value(bc, grid, Center(), dim, I_f...)

as ::Vararg{Integer,N} where N is expected to be N=1 in 1D.

In 2D case, one gets:

(I, dim) = ((0,), Dim{2}())
I_f = (0, 0)

while in 1D case:

(I, dim) = ((0,), Dim{1}())
I_f = (0, 0)

where I_f has 2 elements while N=1.