diff --git a/meshmode/array_context.py b/meshmode/array_context.py index 2e3ae4af97c3f49374c2e661987e60bfb792f401..6118e9d5530c1ee46aff3ff75ed4d0413698420b 100644 --- a/meshmode/array_context.py +++ b/meshmode/array_context.py @@ -162,7 +162,7 @@ class ArrayContext: v = make_zero_and_vars(var_names, params=size_names) domain = v[0].domain() for vname, sname in zip(var_names, size_names): - domain = domain & v[0].le_set(v[vname]) & v[vname].le_set(v[sname]) + domain = domain & v[0].le_set(v[vname]) & v[vname].lt_set(v[sname]) domain_bset, = domain.get_basic_sets() diff --git a/meshmode/discretization/connection/chained.py b/meshmode/discretization/connection/chained.py index 1be269bb7e9a7326a50a134cb4c4fbda18f37968..53e57e2ea02fbc9395244f10cd07f23752bbbeca 100644 --- a/meshmode/discretization/connection/chained.py +++ b/meshmode/discretization/connection/chained.py @@ -202,7 +202,7 @@ def flatten_chained_connection(actx, connection): return connection if not connection.connections: - return make_same_mesh_connection(connection.to_discr, + return make_same_mesh_connection(actx, connection.to_discr, connection.from_discr) # recursively build direct connections