19   for (CeedInt dim = 1; dim <= 3; dim++) {  in main()  local
20     CeedInt num_qpts = CeedIntPow(q, dim);  in main()
22     CeedVectorCreate(ceed, num_qpts * dim * dim, &dx);  in main()
24     CeedVectorCreate(ceed, num_qpts * dim * num_comp, &du);  in main()
26 … CeedScalar dx_array[num_qpts * dim * dim], w_array[num_qpts], du_array[num_qpts * dim * num_comp];  in main()
31       for (CeedInt d = 0; d < dim; d++) {  in main()
32         for (CeedInt g = 0; g < dim; g++) {  in main()
34             dx_array[i + (g * dim + d) * num_qpts] = d == g;  in main()
39         for (CeedInt g = 0; g < dim; g++) {  in main()
50     CeedVectorCreate(ceed, num_qpts * dim * (dim + 1) / 2, &q_data);  in main()
52     CeedVectorCreate(ceed, num_qpts * dim * num_comp, &dv);  in main()
56     snprintf(name_setup, sizeof name_setup, "Poisson%" CeedInt_FMT "DBuild", dim);  in main()
65     snprintf(name_apply, sizeof name_apply, "Vector3Poisson%" CeedInt_FMT "DApply", dim);  in main()
82           for (CeedInt g = 0; g < dim; g++) sum += v_array[i + (g * num_comp + c) * num_qpts];  in main()
84         if (fabs(sum - dim * (c + 1)) > 10 * CEED_EPSILON) {  in main()
86 …dInt_FMT "D volume error in component %" CeedInt_FMT ": %f != %f\n", dim, c, sum, dim * (c + 1.0));  in main()