def test_boundary_interpolation_2d():
"""test boundary interpolation for 2d fields"""
grid = CartesianGrid([[0.1, 0.3], [-2, 3]], [3, 3])
field = ScalarField.random_normal(grid)
# test boundary interpolation
bndry_val = np.random.randn(3)
for bndry in grid._iter_boundaries():
val = field.get_boundary_values(*bndry, bc={"value": bndry_val})
np.testing.assert_allclose(val, bndry_val)
def test_boundary_interpolation_2d():
"""test boundary interpolation for 2d fields"""
grid = CartesianGrid([[0.1, 0.3], [-2, 3]], [3, 3])
field = ScalarField.random_normal(grid)
# test boundary interpolation
bndry_val = np.random.randn(3)
for bndry in grid._iter_boundaries():
val = field.get_boundary_values(*bndry, bc={"value": bndry_val})
np.testing.assert_allclose(val, bndry_val)
# boundary conditions have already been enforced
ev = field.make_get_boundary_values(*bndry)
out = ev()
np.testing.assert_allclose(out, bndry_val)
ev(data_full=field._data_full, out=out)
np.testing.assert_allclose(out, bndry_val)
