Lines Matching +full:format +full:- +full:c
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 // area of a 3D body using matrix-free application of a diffusion operator.
16 // The example has no dependencies, and is designed to be self-contained. For
21 // line argument (-ceed).
30 // ----------------------------------------------------------------------------
32 // ----------------------------------------------------------------------------
33 fn main() -> libceed::Result<()> { in main()
40 fn example_2_vector(options: opt::Opt) -> libceed::Result<()> { in example_2_vector()
72 println!(" Ceed specification [-c] : {}", ceed_spec); in example_2_vector()
73 println!(" Mesh dimension [-d] : {}", dim); in example_2_vector()
74 println!(" Mesh degree [-m] : {}", mesh_degree); in example_2_vector()
75 println!(" Solution degree [-p] : {}", solution_degree); in example_2_vector()
76 println!(" Num. 1D quadr. pts [-q] : {}", num_qpts); in example_2_vector()
77 println!(" Approx. # unknowns [-s] : {}", problem_size); in example_2_vector()
79 " QFunction source [-g] : {}", in example_2_vector()
145 // -- QFunction from user closure in example_2_vector()
161 let qw = weights[i] / (j11 * j22 - j21 * j12); in example_2_vector()
164 qdata[i + q * 2] = -qw * (j11 * j12 + j21 * j22); in example_2_vector()
175 - jacobian[i + q * ((j + 1) % 3 + 3 * ((k + 2) % 3))] in example_2_vector()
220 // -- QFunction from gallery in example_2_vector()
222 let name = format!("Poisson{}DBuild", dim); in example_2_vector()
225 // -- QFunction for use with Operator in example_2_vector()
253 // -- QFunction from user closure in example_2_vector()
259 for c in 0..ncomp_u { in example_2_vector()
261 .skip(c * q) in example_2_vector()
262 .zip(ug.iter().skip(c * q).zip(qdata.iter())) in example_2_vector()
273 for c in 0..ncomp_u { in example_2_vector()
274 let du = [ug[i + (c + 0 * ncomp_u) * q], ug[i + (c + 1 * ncomp_u) * q]]; in example_2_vector()
276 vg[i + (c + j * ncomp_u) * q] = in example_2_vector()
290 for c in 0..ncomp_u { in example_2_vector()
292 ug[i + (c + 0 * ncomp_u) * q], in example_2_vector()
293 ug[i + (c + 1 * ncomp_u) * q], in example_2_vector()
294 ug[i + (c + 2 * ncomp_u) * q], in example_2_vector()
297 vg[i + (c + j * ncomp_u) * q] = du[0] * dxdxdxdx_t[0][j] in example_2_vector()
315 // -- QFunction from gallery in example_2_vector()
317 let name = format!("Vector3Poisson{}DApply", dim); in example_2_vector()
320 // -- QFunction for use with Operator in example_2_vector()
343 for c in 0..ncomp_u { in example_2_vector()
347 .skip(c * q) in example_2_vector()
352 * (c + 1) as libceed::Scalar; in example_2_vector()
371 println!("Surface area error : {:.12e}", area - exact_area); in example_2_vector()
374 1 => 1E-5, in example_2_vector()
375 _ => 1E-1, in example_2_vector()
377 let error = (area - exact_area).abs(); in example_2_vector()
381 message: format!( in example_2_vector()
382 "Volume error too large - expected: {:.12e}, actual: {:.12e}", in example_2_vector()
390 // ----------------------------------------------------------------------------
392 // ----------------------------------------------------------------------------
405 problem_size_requested: -1, in example_2_vector_1d()
421 problem_size_requested: -1, in example_2_vector_2d()
437 problem_size_requested: -1, in example_2_vector_3d()
453 problem_size_requested: -1, in example_2_vector_1d_gallery()
469 problem_size_requested: -1, in example_2_vector_2d_gallery()
485 problem_size_requested: -1, in example_2_vector_3d_gallery()
494 // ----------------------------------------------------------------------------