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

Description
Hello!
GeneralBC stalls when I try to use a vector for dx instead of a scalar.
I found this example in ../test/robin.jl
G = GeneralBC([-10.0, 1.0, 6.0], [-10.0, 1.0, 6.0], 1.0, 3)
When I modify it to
dx=ones(24)*0.03
G = GeneralBC([-10.0, 1.0, 6.0], [-10.0, 1.0, 6.0], dx, 3)
I get the error
DimensionMismatch("tried to assign 4×4 array to 1×4 destination") at ../src/derivative_operators/BC_operators.jl:115
I am trying to debug. Happy to collaborate with someone to fix this.
Cris