6 static PetscErrorCode CheckSymmetry(PetscInt dim, PetscInt order, PetscBool tensor)  in CheckSymmetry()  argument
20 …PetscCall(DMPlexCreateReferenceCell(PETSC_COMM_SELF, DMPolytopeTypeSimpleShape(dim, tensor ? PETSC…  in CheckSymmetry()
35 …nc, &ids, nFunc, &idsCopy, nFunc, &idsCopy2, nFunc * dim, &vals, nFunc * dim, &valsCopy, nFunc * d…  in CheckSymmetry()
46     for (j = 0; j < dim; j++) vals[dim * i + j] = valsCopy2[dim * i + j] = (PetscScalar)points[j];  in CheckSymmetry()
75 …for (l = 0; l < dim; l++) valsCopy[kLocal * dim + l] = vals[(offset + k) * dim + l] * (flip ? flip…  in CheckSymmetry()
83 …Order %" PetscInt_FMT ", Point %" PetscInt_FMT " Symmetry %" PetscInt_FMT, dim, tensor ? "Tensor" …  in CheckSymmetry()
88 …      PetscCall(VecCreateSeqWithArray(PETSC_COMM_SELF, dim, numDofs[depth] * dim, valsCopy, &vec));  in CheckSymmetry()
97 …for (l = 0; l < dim; l++) valsCopy2[(offset + k) * dim + l] = valsCopy[kLocal * dim + l] * (flip ?…  in CheckSymmetry()
100 …, functional %" PetscInt_FMT ": (%" PetscInt_FMT " != %" PetscInt_FMT ")", dim, tensor ? "Tensor" …  in CheckSymmetry()
101         for (l = 0; l < dim; l++) {  in CheckSymmetry()
102 …dim * k + l] == vals[dim * k + l], PETSC_COMM_SELF, PETSC_ERR_PLIB, "Symmetry failure: %" PetscInt…  in CheckSymmetry()
117   PetscInt dim, order, tensor;  in main()  local
122     for (dim = 1; dim <= 3; dim++) {  in main()
123       if (dim == 1 && tensor) continue;  in main()
124 …for (order = 0; order <= (tensor ? 5 : 6); order++) PetscCall(CheckSymmetry(dim, order, tensor ? P…  in main()