From 65a0a921c78927c4ec369ccba73f25cd05436a23 Mon Sep 17 00:00:00 2001
From: Matt Smith <mjsmith6@illinois.edu>
Date: Wed, 2 Mar 2022 20:20:42 -0600
Subject: [PATCH] Make the result of `dt_geometric_factors` in 1D consistent
with other dimensions (#225)
* return a radius, not a diameter
* add comment
---
grudge/dt_utils.py | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/grudge/dt_utils.py b/grudge/dt_utils.py
index 4c477216..b41c683c 100644
--- a/grudge/dt_utils.py
+++ b/grudge/dt_utils.py
@@ -266,7 +266,8 @@ def dt_geometric_factors(
)
if dcoll.dim == 1:
- return freeze(cell_vols)
+ # Inscribed "circle" radius is half the cell size
+ return freeze(cell_vols/2)
dd_face = DOFDesc("all_faces", dd.discretization_tag)
face_discr = dcoll.discr_from_dd(dd_face)
--
GitLab