18   CeedInt             p = 3, q = 4, dim = 2, num_comp = 2;  in main()  local
23   CeedScalar          assembled_values[num_comp * num_comp * num_dofs * num_dofs];  in main()
24   CeedScalar          assembled_true[num_comp * num_comp * num_dofs * num_dofs];  in main()
41   CeedVectorCreate(ceed, num_comp * num_dofs, &u);  in main()
42   CeedVectorCreate(ceed, num_comp * num_dofs, &v);  in main()
56 …CeedElemRestrictionCreate(ceed, num_elem, p * p, num_comp, num_dofs, num_comp * num_dofs, CEED_MEM…  in main()
64   CeedBasisCreateTensorH1Lagrange(ceed, dim, num_comp, p, q, CEED_GAUSS, &basis_u);  in main()
74   CeedQFunctionAddInput(qf_mass, "u", num_comp, CEED_EVAL_INTERP);  in main()
75   CeedQFunctionAddOutput(qf_mass, "v", num_comp, CEED_EVAL_INTERP);  in main()
97   for (CeedInt k = 0; k < num_comp * num_comp * num_dofs * num_dofs; k++) {  in main()
107 …for (CeedInt k = 0; k < num_entries; ++k) assembled_values[rows[k] * num_comp * num_dofs + cols[k]…  in main()
115   for (CeedInt j = 0; j < num_dofs * num_comp; j++) {  in main()
132 …for (CeedInt i = 0; i < num_dofs * num_comp; i++) assembled_true[i * num_dofs * num_comp + j] = v_…  in main()
137   for (CeedInt i = 0; i < num_comp * num_dofs; i++) {  in main()
138     for (CeedInt j = 0; j < num_comp * num_dofs; j++) {  in main()
139 …if (fabs(assembled_values[i * num_dofs * num_comp + j] - assembled_true[i * num_dofs * num_comp + …  in main()
141 …CeedInt_FMT "] Error in assembly: %f != %f\n", i, j, assembled_values[i * num_dofs * num_comp + j],  in main()
142                assembled_true[i * num_dofs * num_comp + j]);  in main()