From d2c52ac1f33b597706dae04aeab59369d3f58f33 Mon Sep 17 00:00:00 2001
From: Shivam Gupta <sgupta72@illinois.edu>
Date: Tue, 21 Apr 2015 01:18:36 -0500
Subject: [PATCH] Fixed issue with connectivity to self
---
meshmode/mesh/refinement.py | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/meshmode/mesh/refinement.py b/meshmode/mesh/refinement.py
index a57d5233..7e9f7c1e 100644
--- a/meshmode/mesh/refinement.py
+++ b/meshmode/mesh/refinement.py
@@ -692,10 +692,12 @@ class Refiner(object):
for ivertex in grp[iel_grp]:
element_to_element[element_index].update(
vertex_to_element[ivertex])
- if self.hanging_vertex_element[ivertex] != -1:
+ if self.hanging_vertex_element[ivertex] != -1 and element_index != self.hanging_vertex_element[ivertex]:
element_to_element[element_index].update([self.hanging_vertex_element[ivertex]])
element_to_element[self.hanging_vertex_element[ivertex]].update([element_index])
element_index += 1
+ for iel, neighbors in enumerate(element_to_element):
+ neighbors.remove(iel)
#print self.ray_elements
'''
for ray in self.rays:
--
GitLab