From 65a0a921c78927c4ec369ccba73f25cd05436a23 Mon Sep 17 00:00:00 2001 From: Matt Smith 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 4c4772169..b41c683cd 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)