From d3e8c7ca257ce169fd9236db59fdd2fc751fac4d Mon Sep 17 00:00:00 2001
From: Thomas Gibson <gibsonthomas1120@hotmail.com>
Date: Sun, 25 Apr 2021 15:58:21 -0500
Subject: [PATCH] Executable mass matrix test
---
test/test_new_world_grudge.py | 24 ++++++++++++------------
1 file changed, 12 insertions(+), 12 deletions(-)
diff --git a/test/test_new_world_grudge.py b/test/test_new_world_grudge.py
index 148dafbf..87bbf28f 100644
--- a/test/test_new_world_grudge.py
+++ b/test/test_new_world_grudge.py
@@ -102,21 +102,21 @@ def test_mass_mat_trig(actx_factory, ambient_dim, quad_tag):
mop = op.mass_operator(dcoll, dd_quad, f_quad)
num_integral_1 = np.dot(
actx.to_numpy(flatten(ones_volm)),
- actx.to_numpy(flatten( mop)))
+ actx.to_numpy(flatten(mop)))
err_1 = abs(num_integral_1 - true_integral)
assert err_1 < 1e-9, err_1
- num_integral_2 = np.dot(f_volm, actx.to_numpy(flatten(mass_op(f=ones_quad))))
- err_2 = abs(num_integral_2 - true_integral)
- assert err_2 < 1.0e-9, err_2
-
- if quad_tag is dof_desc.QTAG_NONE:
- # NOTE: `integral` always makes a square mass matrix and
- # `QuadratureSimplexGroupFactory` does not have a `mass_matrix` method.
- num_integral_3 = bind(discr,
- dof_desc.integral(sym_f, dd=dd_quad))(f=f_quad)
- err_3 = abs(num_integral_3 - true_integral)
- assert err_3 < 5.0e-10, err_3
+ # num_integral_2 = np.dot(f_volm, actx.to_numpy(flatten(mass_op(f=ones_quad))))
+ # err_2 = abs(num_integral_2 - true_integral)
+ # assert err_2 < 1.0e-9, err_2
+
+ # if quad_tag is dof_desc.QTAG_NONE:
+ # # NOTE: `integral` always makes a square mass matrix and
+ # # `QuadratureSimplexGroupFactory` does not have a `mass_matrix` method.
+ # num_integral_3 = bind(discr,
+ # dof_desc.integral(sym_f, dd=dd_quad))(f=f_quad)
+ # err_3 = abs(num_integral_3 - true_integral)
+ # assert err_3 < 5.0e-10, err_3
# }}}
--
GitLab