xref: /petsc/src/sys/tests/ex7.c (revision 732aec7a18f2199fb53bb9a2f3aef439a834ce31)
1c4762a1bSJed Brown static char help[] = "Tests PetscSortIntWithPermutation().\n\n";
2c4762a1bSJed Brown 
3c4762a1bSJed Brown #include <petscsys.h>
4c4762a1bSJed Brown 
main(int argc,char ** argv)5d71ae5a4SJacob Faibussowitsch int main(int argc, char **argv)
6d71ae5a4SJacob Faibussowitsch {
7c4762a1bSJed Brown   PetscInt i;
8c4762a1bSJed Brown   PetscInt x[]  = {39, 9, 39, 39, 29}, index[5];
9c4762a1bSJed Brown   PetscInt x2[] = {39, 9, 19, 39, 29, 39, 29, 39}, index2[8];
10c4762a1bSJed Brown 
11327415f7SBarry Smith   PetscFunctionBeginUser;
12*c8025a54SPierre Jolivet   PetscCall(PetscInitialize(&argc, &argv, NULL, help));
139566063dSJacob Faibussowitsch   PetscCall(PetscPrintf(PETSC_COMM_SELF, "1st test\n"));
14c4762a1bSJed Brown   for (i = 0; i < 5; i++) index[i] = i;
159566063dSJacob Faibussowitsch   PetscCall(PetscSortIntWithPermutation(5, x, index));
169566063dSJacob Faibussowitsch   for (i = 0; i < 5; i++) PetscCall(PetscPrintf(PETSC_COMM_SELF, " %" PetscInt_FMT "     %" PetscInt_FMT "     %" PetscInt_FMT "\n", x[i], index[i], x[index[i]]));
17c4762a1bSJed Brown 
189566063dSJacob Faibussowitsch   PetscCall(PetscPrintf(PETSC_COMM_SELF, "\n2nd test\n"));
19c4762a1bSJed Brown   for (i = 0; i < 8; i++) index2[i] = i;
209566063dSJacob Faibussowitsch   PetscCall(PetscSortIntWithPermutation(8, x2, index2));
219566063dSJacob Faibussowitsch   for (i = 0; i < 8; i++) PetscCall(PetscPrintf(PETSC_COMM_SELF, " %" PetscInt_FMT "     %" PetscInt_FMT "     %" PetscInt_FMT "\n", x2[i], index2[i], x2[index2[i]]));
229566063dSJacob Faibussowitsch   PetscCall(PetscFinalize());
23b122ec5aSJacob Faibussowitsch   return 0;
24c4762a1bSJed Brown }
25c4762a1bSJed Brown 
26c4762a1bSJed Brown /*TEST
27c4762a1bSJed Brown 
28c4762a1bSJed Brown    test:
29c4762a1bSJed Brown 
30c4762a1bSJed Brown TEST*/
31