python/tests/test-3-basis.py

1 # Copyright (c) 2017-2026, Lawrence Livermore National Security, LLC and other CEED contributors
2 # All Rights Reserved. See the top-level LICENSE and NOTICE files for details.
4 # SPDX-License-Identifier: BSD-2-Clause
16 import check
20 # -------------------------------------------------------------------------------
22 # -------------------------------------------------------------------------------
28         result *= math.tanh(x[d] - center)
40 # -------------------------------------------------------------------------------
42 # -------------------------------------------------------------------------------
58         stdout, stderr, ref_stdout = check.output(capsys)
62 # -------------------------------------------------------------------------------
64 # -------------------------------------------------------------------------------
78 # -------------------------------------------------------------------------------
80 # -------------------------------------------------------------------------------
96                                         (2**(dim - d))) // (2**(dim - d - 1)) else -1
135                 assert math.fabs(u[i] - fx) < 1E-4
137 # -------------------------------------------------------------------------------
139 # -------------------------------------------------------------------------------
157                                         (2**(dim - d))) // (2**(dim - d - 1)) else -1
189         # Check if 1' * G * u = u' * (G' * 1)
196         assert math.fabs(sum1 - sum2) < 10. * TOL
198 # -------------------------------------------------------------------------------
199 # Test creation and destruction of a 2D Simplex non-tensor H1 basis
200 # -------------------------------------------------------------------------------
220 # -------------------------------------------------------------------------------
221 # Test integration with a 2D Simplex non-tensor H1 basis
222 # -------------------------------------------------------------------------------
255     # Check values at quadrature points
260         assert math.fabs(sum - 17. / 24.) < 10. * TOL
262 # -------------------------------------------------------------------------------
263 # Test grad with a 2D Simplex non-tensor H1 basis
264 # -------------------------------------------------------------------------------
296     # Check values at quadrature points
300             assert math.fabs(out_array[0 * Q + i] - value) < 10. * TOL
303             assert math.fabs(out_array[1 * Q + i] - value) < 10. * TOL
305 # -------------------------------------------------------------------------------
306 # Test interpolation with a 2D Quad non-tensor H(div) basis
307 # -------------------------------------------------------------------------------
332     # Check values at quadrature points
336             assert math.fabs(out_array[1 * Q + i] - 1.) < 10. * TOL
338 # -------------------------------------------------------------------------------
339 # Test interpolation with a 2D Simplex non-tensor H(curl) basis
340 # -------------------------------------------------------------------------------
369     # Check values at quadrature points
372             assert math.fabs(out_array[0 * Q + i] - 1.) < 10. * TOL
375     in_array[0] = -1.
383     # Check values at quadrature points
386             assert math.fabs(out_array[1 * Q + i] - 1.) < 10. * TOL
388 # -------------------------------------------------------------------------------