10   PetscInt *btupprev, *btup;  in main()  local
15   PetscCall(PetscMalloc3(maxdim + 1, &btup, maxdim + 1, &btupprev, maxdim, &gtup));  in main()
24         PetscCall(PetscDTIndexToBary(d + 1, n, k, btup));  in main()
26 …PetscCheck(btup[j] >= 0, PETSC_COMM_SELF, PETSC_ERR_PLIB, "PetscDTIndexToBary, d = %" PetscInt_FMT…  in main()
27           sum += btup[j];  in main()
30         PetscCall(PetscDTBaryToIndex(d + 1, n, btup, &kchk));  in main()
34           while (j >= 0 && btup[j] == btupprev[j]) j--;  in main()
36 …PetscCheck(btup[j] >= btupprev[j], PETSC_COMM_SELF, PETSC_ERR_PLIB, "PetscDTIndexToBary, d = %" Pe…  in main()
37         } else PetscCall(PetscArraycpy(btupprev, btup, d + 1));  in main()
41 …for (j = 0; j < d; j++) PetscCheck(gtup[j] == btup[d - 1 - j], PETSC_COMM_SELF, PETSC_ERR_PLIB, "P…  in main()
45   PetscCall(PetscFree3(btup, btupprev, gtup));  in main()