294 …tscDTUniformTensorQuadrature(PetscInt dim, PetscInt num_comp, PetscInt num_points, PetscReal a, Pe…  in PetscDTUniformTensorQuadrature()  argument
297 …total_points = dim > 1 ? (dim > 2 ? (num_points * num_points * num_points) : (num_points * num_poi…  in PetscDTUniformTensorQuadrature()
317       PetscReal step = (b - a) / num_points;  in PetscDTUniformTensorQuadrature()
319       for (PetscInt i = 0; i < num_points; i++) {  in PetscDTUniformTensorQuadrature()
326       PetscCall(PetscMalloc1(num_points, &coords_1d));  in PetscDTUniformTensorQuadrature()
327       PetscReal step = (b - a) / num_points;  in PetscDTUniformTensorQuadrature()
329       for (PetscInt i = 0; i < num_points; i++) coords_1d[i] = step * (i + 0.5) + a;  in PetscDTUniformTensorQuadrature()
330       for (PetscInt i = 0; i < num_points; i++) {  in PetscDTUniformTensorQuadrature()
331         for (PetscInt j = 0; j < num_points; j++) {  in PetscDTUniformTensorQuadrature()
332           coords[(i * num_points + j) * dim + 0] = coords_1d[i];  in PetscDTUniformTensorQuadrature()
333           coords[(i * num_points + j) * dim + 1] = coords_1d[j];  in PetscDTUniformTensorQuadrature()
334 …        for (PetscInt c = 0; c < num_comp; c++) weights[(i * num_points + j) * num_comp + c] = 1.0;  in PetscDTUniformTensorQuadrature()
341       PetscCall(PetscMalloc1(num_points, &coords_1d));  in PetscDTUniformTensorQuadrature()
342       PetscReal step = (b - a) / num_points;  in PetscDTUniformTensorQuadrature()
344       for (PetscInt i = 0; i < num_points; i++) coords_1d[i] = step * (i + 0.5) + a;  in PetscDTUniformTensorQuadrature()
345       for (PetscInt i = 0; i < num_points; i++) {  in PetscDTUniformTensorQuadrature()
346         for (PetscInt j = 0; j < num_points; j++) {  in PetscDTUniformTensorQuadrature()
347           for (PetscInt k = 0; k < num_points; k++) {  in PetscDTUniformTensorQuadrature()
348             coords[((i * num_points + j) * num_points + k) * dim + 0] = coords_1d[i];  in PetscDTUniformTensorQuadrature()
349             coords[((i * num_points + j) * num_points + k) * dim + 1] = coords_1d[j];  in PetscDTUniformTensorQuadrature()
350             coords[((i * num_points + j) * num_points + k) * dim + 2] = coords_1d[k];  in PetscDTUniformTensorQuadrature()
351 …for (PetscInt c = 0; c < num_comp; c++) weights[((i * num_points + j) * num_points + k) * num_comp…  in PetscDTUniformTensorQuadrature()
364   PetscCall(PetscQuadratureSetOrder(*q, num_points));  in PetscDTUniformTensorQuadrature()