examples/python/ex3_volume.py

2 # Copyright (c) 2017-2026, Lawrence Livermore National Security, LLC and other CEED contributors.
3 # All Rights Reserved. See the top-level LICENSE and NOTICE files for details.
5 # SPDX-License-Identifier: BSD-2-Clause
14 #     python ex3_volume -c /cpu/self
15 #     python ex3_volume -c /gpu/cuda
50         print(f"    Ceed specification [-c] : {args.ceed}")
51         print(f"    Mesh dimension     [-d] : {dim}")
52         print(f"    Mesh degree        [-m] : {mesh_degree}")
53         print(f"    Solution degree    [-p] : {sol_degree}")
54         print(f"    Num. 1D quadr. pts [-q] : {num_qpts}")
55         print(f"    Approx. # unknowns [-s] : {problem_size}")
56         print(f"    QFunction source   [-g] : {'gallery' if args.gallery else 'user'}")
58     # Check - Gallery not supported
67     # Tensor-product Lagrange basis for mesh coordinates
71     # Tensor-product Lagrange basis for solution
114                               os.path.join(file_dir, "ex3-volume.h:build_mass_diff"))
132                               os.path.join(file_dir, "ex3-volume.h:apply_mass_diff"))
151     u.set_value(1.0)  # Set all entries of u to 1.0
156     # Compute volume by summing all entries in v
159         # Simply sum all values to compute the volume
166         print(f"Volume error         : {volume - exact_volume:.14g}")
168         # Test mode - check if error is within tolerance
169         tol = 200 * libceed.EPSILON if dim == 1 else 1e-5
170         if abs(volume - exact_volume) > tol:
171             print(f"Volume error : {volume - exact_volume:.14g}")