Skip to content
Snippets Groups Projects
Commit 0a025178 authored by Alexandru Fikl's avatar Alexandru Fikl
Browse files

fix some remaining issues with 1d meshes

parent 8b127637
Branches
Tags
1 merge request!71Fix some remaining issues with 1D meshes
...@@ -43,7 +43,7 @@ def _make_cross_face_batches(queue, ...@@ -43,7 +43,7 @@ def _make_cross_face_batches(queue,
return cl.array.to_device(queue, ary, array_queue=None) return cl.array.to_device(queue, ary, array_queue=None)
from meshmode.discretization.connection.direct import InterpolationBatch from meshmode.discretization.connection.direct import InterpolationBatch
if tgt_bdry_discr.ambient_dim == 1: if tgt_bdry_discr.dim == 0:
yield InterpolationBatch( yield InterpolationBatch(
from_group_index=i_src_grp, from_group_index=i_src_grp,
from_element_indices=to_dev(src_bdry_element_indices), from_element_indices=to_dev(src_bdry_element_indices),
......
...@@ -177,7 +177,9 @@ class QuadratureSimplexElementGroup(SimplexElementGroupBase): ...@@ -177,7 +177,9 @@ class QuadratureSimplexElementGroup(SimplexElementGroupBase):
@memoize_method @memoize_method
def _quadrature_rule(self): def _quadrature_rule(self):
dims = self.mesh_el_group.dim dims = self.mesh_el_group.dim
if dims == 1: if dims == 0:
return mp.Quadrature(np.empty((0, 1)), np.empty((0, 1)))
elif dims == 1:
return mp.LegendreGaussQuadrature(self.order) return mp.LegendreGaussQuadrature(self.order)
else: else:
return mp.XiaoGimbutasSimplexQuadrature(self.order, dims) return mp.XiaoGimbutasSimplexQuadrature(self.order, dims)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment