From 39e4145f5cca40b024715f15e35aa89d4fe4d217 Mon Sep 17 00:00:00 2001 From: Shivam Gupta Date: Fri, 10 Apr 2015 17:49:27 -0500 Subject: [PATCH] Fixed bug where elements was not being used to generate connectivity, and also, cleared velements when a midpoint had already been created when refining to make purpose of velements be only to store hanging node connectivity --- meshmode/mesh/refinement.py | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/meshmode/mesh/refinement.py b/meshmode/mesh/refinement.py index cd2ce945..9cb9a21d 100644 --- a/meshmode/mesh/refinement.py +++ b/meshmode/mesh/refinement.py @@ -467,7 +467,9 @@ class Refiner(object): verts.append(vertices_idx) vertices_idx += 1 else: - verts.append(self.pair_map[vertex_pair].midpoint) + cur_midpoint = self.pair_map[vertex_pair].midpoint + verts.append(cur_midpoint) + self.vertex_to_ray[cur_midpoint][cur_pmap.ray].value.velements = [None, None] # }}} @@ -674,13 +676,11 @@ class Refiner(object): for ray in self.rays: curnode = ray.first while curnode is not None: - ''' if len(curnode.value.elements) >= 2: if curnode.value.elements[0] is not None: element_to_element[curnode.value.elements[0]].update(curnode.value.elements) if curnode.value.elements[1] is not None: element_to_element[curnode.value.elements[1]].update(curnode.value.elements) - ''' if len(curnode.value.velements) >= 2: if curnode.value.velements[0] is not None: element_to_element[curnode.value.velements[0]].update(curnode.value.velements) -- GitLab