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
11 // 3D body using matrix-free application of a mass operator.  Arbitrary mesh and
14 // The example has no dependencies, and is designed to be self-contained. For
19 // line argument (-ceed).
28 // ----------------------------------------------------------------------------
30 // ----------------------------------------------------------------------------
31 fn main() -> libceed::Result<()> {  in main()
38 fn example_1(options: opt::Opt) -> libceed::Result<()> {  in example_1()
47         test,  in example_1()
57         if test {  in example_1()
69         println!("    Ceed specification [-c] : {}", ceed_spec);  in example_1()
70         println!("    Mesh dimension     [-d] : {}", dim);  in example_1()
71         println!("    Mesh degree        [-m] : {}", mesh_degree);  in example_1()
72         println!("    Solution degree    [-p] : {}", solution_degree);  in example_1()
73         println!("    Num. 1D quadr. pts [-q] : {}", num_qpts);  in example_1()
74         println!("    Approx. # unknowns [-s] : {}", problem_size);  in example_1()
76             "    QFunction source   [-g] : {}",  in example_1()
133     // -- QFunction from user closure  in example_1()
147                             - jacobian[i + q * 1] * jacobian[i + q * 2])  in example_1()
158                                 - jacobian[i + q * 5] * jacobian[i + q * 7])  in example_1()
159                             - jacobian[i + q * 1]  in example_1()
161                                     - jacobian[i + q * 5] * jacobian[i + q * 6])  in example_1()
164                                     - jacobian[i + q * 4] * jacobian[i + q * 6]))  in example_1()
180     // -- QFunction from gallery  in example_1()
185     // -- QFunction for use with Operator  in example_1()
213     // -- QFunction from user closure  in example_1()
227     // -- QFunction from gallery  in example_1()
229     // -- QFunction for use with Operator  in example_1()
261             volume - exact_volume  in example_1()
266         _ => 1E-5,  in example_1()
268     let error = (volume - exact_volume).abs();  in example_1()
273                 "Volume error too large - expected: {:.12e}, actual: {:.12e}",  in example_1()
281 // ----------------------------------------------------------------------------
283 // ----------------------------------------------------------------------------
284 #[cfg(test)]
288     #[test]
296             problem_size_requested: -1,  in example_1_1d()
297             test: true,  in example_1_1d()
304     #[test]
312             problem_size_requested: -1,  in example_1_2d()
313             test: true,  in example_1_2d()
320     #[test]
328             problem_size_requested: -1,  in example_1_3d()
329             test: true,  in example_1_3d()
336     #[test]
344             problem_size_requested: -1,  in example_1_1d_gallery()
345             test: true,  in example_1_1d_gallery()
352     #[test]
360             problem_size_requested: -1,  in example_1_2d_gallery()
361             test: true,  in example_1_2d_gallery()
368     #[test]
376             problem_size_requested: -1,  in example_1_3d_gallery()
377             test: true,  in example_1_3d_gallery()
385 // ----------------------------------------------------------------------------