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()
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()
186         let name = format!("Mass{}DBuild", dim);  in example_1_vector()
189     // -- QFunction for use with Operator  in example_1_vector()
217     // -- QFunction from user closure  in example_1_vector()
221         for c in 0..ncomp_u {  in example_1_vector()
223                 .skip(c * q)  in example_1_vector()
224                 .zip(u.iter().skip(c * q).zip(qdata.iter()))  in example_1_vector()
235     // -- QFunction from gallery  in example_1_vector()
237     // -- QFunction for use with Operator  in example_1_vector()
259     for c in 0..ncomp_u {  in example_1_vector()
261         u.view_mut()?.iter_mut().skip(c * q).take(q).for_each(|u| {  in example_1_vector()
262             *u = (c + 1) as libceed::Scalar;  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()
290             message: format!(  in example_1_vector()
291                 "Volume error too large - expected: {:.12e}, actual: {:.12e}",  in example_1_vector()
299 // ----------------------------------------------------------------------------
301 // ----------------------------------------------------------------------------
314             problem_size_requested: -1,  in example_1_vector_1d()
330             problem_size_requested: -1,  in example_1_vector_2d()
346             problem_size_requested: -1,  in example_1_vector_3d()
362             problem_size_requested: -1,  in example_1_vector_1d_gallery()
378             problem_size_requested: -1,  in example_1_vector_2d_gallery()
394             problem_size_requested: -1,  in example_1_vector_3d_gallery()
403 // ----------------------------------------------------------------------------