1 /// @file
2 /// Test viewing of mass matrix operator
3 /// \test Test viewing of mass matrix operator
4 #include <ceed.h>
5 #include <math.h>
6 #include <stdio.h>
7 #include <stdlib.h>
8
9 #include "t500-operator.h"
10
main(int argc,char ** argv)11 int main(int argc, char **argv) {
12 Ceed ceed;
13 CeedElemRestriction elem_restriction_x, elem_restriction_u, elem_restriction_q_data;
14 CeedBasis basis_x, basis_u;
15 CeedQFunction qf_setup, qf_mass;
16 CeedOperator op_setup, op_mass;
17 CeedVector q_data;
18 CeedInt num_elem = 15, p = 5, q = 8;
19 CeedInt num_nodes_x = num_elem + 1, num_nodes_u = num_elem * (p - 1) + 1;
20 CeedInt ind_x[num_elem * 2], ind_u[num_elem * p];
21
22 CeedInit(argv[1], &ceed);
23
24 CeedVectorCreate(ceed, num_elem * q, &q_data);
25
26 // Restrictions
27 for (CeedInt i = 0; i < num_elem; i++) {
28 ind_x[2 * i + 0] = i;
29 ind_x[2 * i + 1] = i + 1;
30 }
31 CeedElemRestrictionCreate(ceed, num_elem, 2, 1, 1, num_nodes_x, CEED_MEM_HOST, CEED_USE_POINTER, ind_x, &elem_restriction_x);
32
33 for (CeedInt i = 0; i < num_elem; i++) {
34 for (CeedInt j = 0; j < p; j++) {
35 ind_u[p * i + j] = 2 * (i * (p - 1) + j);
36 }
37 }
38 CeedElemRestrictionCreate(ceed, num_elem, p, 2, 1, 2 * num_nodes_u, CEED_MEM_HOST, CEED_USE_POINTER, ind_u, &elem_restriction_u);
39
40 CeedInt strides_q_data[3] = {1, q, q};
41 CeedElemRestrictionCreateStrided(ceed, num_elem, q, 1, q * num_elem, strides_q_data, &elem_restriction_q_data);
42
43 // Bases
44 CeedBasisCreateTensorH1Lagrange(ceed, 1, 1, 2, q, CEED_GAUSS, &basis_x);
45 CeedBasisCreateTensorH1Lagrange(ceed, 1, 2, p, q, CEED_GAUSS, &basis_u);
46
47 // QFunctions
48 CeedQFunctionCreateInterior(ceed, 1, setup, setup_loc, &qf_setup);
49 CeedQFunctionAddInput(qf_setup, "weight", 1, CEED_EVAL_WEIGHT);
50 CeedQFunctionAddInput(qf_setup, "dx", 1 * 1, CEED_EVAL_GRAD);
51 CeedQFunctionAddOutput(qf_setup, "rho", 1, CEED_EVAL_NONE);
52
53 CeedQFunctionCreateInterior(ceed, 1, mass, mass_loc, &qf_mass);
54 CeedQFunctionAddInput(qf_mass, "rho", 1, CEED_EVAL_NONE);
55 CeedQFunctionAddInput(qf_mass, "u", 2, CEED_EVAL_INTERP);
56 CeedQFunctionAddOutput(qf_mass, "v", 2, CEED_EVAL_INTERP);
57
58 // Operators
59 CeedOperatorCreate(ceed, qf_setup, CEED_QFUNCTION_NONE, CEED_QFUNCTION_NONE, &op_setup);
60 CeedOperatorSetField(op_setup, "weight", CEED_ELEMRESTRICTION_NONE, basis_x, CEED_VECTOR_NONE);
61 CeedOperatorSetField(op_setup, "dx", elem_restriction_x, basis_x, CEED_VECTOR_ACTIVE);
62 CeedOperatorSetField(op_setup, "rho", elem_restriction_q_data, CEED_BASIS_NONE, CEED_VECTOR_ACTIVE);
63
64 CeedOperatorCreate(ceed, qf_mass, CEED_QFUNCTION_NONE, CEED_QFUNCTION_NONE, &op_mass);
65 CeedOperatorSetField(op_mass, "rho", elem_restriction_q_data, CEED_BASIS_NONE, q_data);
66 CeedOperatorSetField(op_mass, "u", elem_restriction_u, basis_u, CEED_VECTOR_ACTIVE);
67 CeedOperatorSetField(op_mass, "v", elem_restriction_u, basis_u, CEED_VECTOR_ACTIVE);
68
69 CeedOperatorSetName(op_setup, "setup");
70 CeedOperatorViewTerse(op_setup, stdout);
71 CeedOperatorView(op_setup, stdout);
72
73 // Check tabs and CeedObject functionality
74 {
75 CeedOperator op_mass_copy = NULL;
76
77 CeedOperatorReferenceCopy(op_mass, &op_mass_copy);
78 CeedOperatorSetName(op_mass_copy, "mass");
79 CeedOperatorSetNumViewTabs(op_mass_copy, 1);
80 CeedOperatorViewTerse(op_mass_copy, stdout);
81 CeedObjectView((CeedObject)op_mass_copy, stdout);
82 CeedObjectDestroy((CeedObject *)&op_mass_copy);
83 }
84
85 CeedVectorDestroy(&q_data);
86 CeedElemRestrictionDestroy(&elem_restriction_u);
87 CeedElemRestrictionDestroy(&elem_restriction_x);
88 CeedElemRestrictionDestroy(&elem_restriction_q_data);
89 CeedBasisDestroy(&basis_u);
90 CeedBasisDestroy(&basis_x);
91 CeedQFunctionDestroy(&qf_setup);
92 CeedQFunctionDestroy(&qf_mass);
93 CeedOperatorDestroy(&op_setup);
94 CeedOperatorDestroy(&op_mass);
95 CeedDestroy(&ceed);
96 return 0;
97 }
98