From 3d16443e10490ea43193b40941a5beb478ad34c3 Mon Sep 17 00:00:00 2001 From: Cory Mikida Date: Fri, 8 May 2020 16:02:02 -0500 Subject: [PATCH] Small changes to get SBP-DG cases running with Grudge --- meshmode/discretization/__init__.py | 3 +-- meshmode/mesh/generation.py | 9 +++++---- 2 files changed, 6 insertions(+), 6 deletions(-) diff --git a/meshmode/discretization/__init__.py b/meshmode/discretization/__init__.py index 6a4fcda7..70a93798 100644 --- a/meshmode/discretization/__init__.py +++ b/meshmode/discretization/__init__.py @@ -356,9 +356,8 @@ class Discretization(object): "stride:auto,stride:auto,stride:auto") return knl - result = self.empty(dtype=self.real_dtype, extra_dims=(self.ambient_dim,)) - with cl.CommandQueue(self.cl_context) as queue: + result = self.empty(queue=queue, dtype=self.real_dtype, extra_dims=(self.ambient_dim,)) for grp in self.groups: if grp.nelements == 0: continue diff --git a/meshmode/mesh/generation.py b/meshmode/mesh/generation.py index 3540dadd..de5d9581 100644 --- a/meshmode/mesh/generation.py +++ b/meshmode/mesh/generation.py @@ -658,7 +658,7 @@ def generate_urchin(order, m, n, est_rel_interp_tolerance, min_rad=0.2): # {{{ generate_box_mesh def generate_box_mesh(axis_coords, order=1, coord_dtype=np.float64, - group_factory=None): + group_factory=None, boundary_tags=[]): """Create a semi-structured mesh. :param axis_coords: a tuple with a number of entries corresponding @@ -776,14 +776,15 @@ def generate_box_mesh(axis_coords, order=1, coord_dtype=np.float64, from meshmode.mesh import Mesh return Mesh(vertices, [grp], - is_conforming=True) + is_conforming=True, boundary_tags=boundary_tags) # }}} # {{{ generate_regular_rect_mesh -def generate_regular_rect_mesh(a=(0, 0), b=(1, 1), n=(5, 5), order=1): +def generate_regular_rect_mesh(a=(0, 0), b=(1, 1), n=(5, 5), order=1, + boundary_tags=[]): """Create a semi-structured rectangular mesh. :param a: the lower left hand point of the rectangle @@ -797,7 +798,7 @@ def generate_regular_rect_mesh(a=(0, 0), b=(1, 1), n=(5, 5), order=1): axis_coords = [np.linspace(a_i, b_i, n_i) for a_i, b_i, n_i in zip(a, b, n)] - return generate_box_mesh(axis_coords, order=order) + return generate_box_mesh(axis_coords, order=order, boundary_tags=boundary_tags) # }}} -- GitLab