From 170f895676a38e6650b879f10b56f7e379b8602c Mon Sep 17 00:00:00 2001
From: Andreas Kloeckner <inform@tiker.net>
Date: Tue, 8 Jul 2014 07:51:03 -0500
Subject: [PATCH] Tweak node/vertex consistency tolerance
---
meshmode/mesh/__init__.py | 12 +++++++++---
1 file changed, 9 insertions(+), 3 deletions(-)
diff --git a/meshmode/mesh/__init__.py b/meshmode/mesh/__init__.py
index f17d112..7d19b1a 100644
--- a/meshmode/mesh/__init__.py
+++ b/meshmode/mesh/__init__.py
@@ -24,7 +24,7 @@ THE SOFTWARE.
import numpy as np
import modepy as mp
-#import numpy.linalg as la
+import numpy.linalg as la
from pytools import Record
__doc__ = """
@@ -297,9 +297,15 @@ def _test_node_vertex_consistency_simplex(mesh, mgrp):
np.sum((map_vertices - grp_vertices)**2, axis=0),
axis=-1))
- tol = 1e2 * np.finfo(per_element_vertex_errors.dtype).eps
+ tol = 1e3 * np.finfo(per_element_vertex_errors.dtype).eps
- assert np.max(per_element_vertex_errors) < tol, np.max(per_element_vertex_errors)
+ from meshmode.mesh.processing import find_bounding_box
+
+ bbox_min, bbox_max = find_bounding_box(mesh)
+ size = la.norm(bbox_max-bbox_min)
+
+ assert np.max(per_element_vertex_errors) < tol*size, \
+ np.max(per_element_vertex_errors)
return True
--
GitLab