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
15 // The example has no dependencies, and is designed to be self-contained. For
20 // line argument (-ceed).
29 // ----------------------------------------------------------------------------
31 // ----------------------------------------------------------------------------
32 fn main() -> libceed::Result<()> {  in main()
39 fn example_1_vector(options: opt::Opt) -> libceed::Result<()> {  in example_1_vector()
48         test,  in example_1_vector()
58         if test {  in example_1_vector()
71         println!("    Ceed specification [-c] : {}", ceed_spec);  in example_1_vector()
72         println!("    Mesh dimension     [-d] : {}", dim);  in example_1_vector()
73         println!("    Mesh degree        [-m] : {}", mesh_degree);  in example_1_vector()
74         println!("    Solution degree    [-p] : {}", solution_degree);  in example_1_vector()
75         println!("    Num. 1D quadr. pts [-q] : {}", num_qpts);  in example_1_vector()
76         println!("    Approx. # unknowns [-s] : {}", problem_size);  in example_1_vector()
78             "    QFunction source   [-g] : {}",  in example_1_vector()
137     // -- QFunction from user closure  in example_1_vector()
151                             - jacobian[i + q * 1] * jacobian[i + q * 2])  in example_1_vector()
162                                 - jacobian[i + q * 5] * jacobian[i + q * 7])  in example_1_vector()
163                             - jacobian[i + q * 1]  in example_1_vector()
165                                     - jacobian[i + q * 5] * jacobian[i + q * 6])  in example_1_vector()
168                                     - jacobian[i + q * 4] * jacobian[i + q * 6]))  in example_1_vector()
184     // -- QFunction from gallery  in example_1_vector()
189     // -- QFunction for use with Operator  in example_1_vector()
217     // -- QFunction from user closure  in example_1_vector()
235     // -- QFunction from gallery  in example_1_vector()
237     // -- QFunction for use with Operator  in example_1_vector()
279             volume - exact_volume  in example_1_vector()
284         _ => 1E-5,  in example_1_vector()
286     let error = (volume - exact_volume).abs();  in example_1_vector()
291                 "Volume error too large - expected: {:.12e}, actual: {:.12e}",  in example_1_vector()
299 // ----------------------------------------------------------------------------
301 // ----------------------------------------------------------------------------
302 #[cfg(test)]
306     #[test]
314             problem_size_requested: -1,  in example_1_vector_1d()
315             test: true,  in example_1_vector_1d()
322     #[test]
330             problem_size_requested: -1,  in example_1_vector_2d()
331             test: true,  in example_1_vector_2d()
338     #[test]
346             problem_size_requested: -1,  in example_1_vector_3d()
347             test: true,  in example_1_vector_3d()
354     #[test]
362             problem_size_requested: -1,  in example_1_vector_1d_gallery()
363             test: true,  in example_1_vector_1d_gallery()
370     #[test]
378             problem_size_requested: -1,  in example_1_vector_2d_gallery()
379             test: true,  in example_1_vector_2d_gallery()
386     #[test]
394             problem_size_requested: -1,  in example_1_vector_3d_gallery()
395             test: true,  in example_1_vector_3d_gallery()
403 // ----------------------------------------------------------------------------