1c4762a1bSJed Brown 2c4762a1bSJed Brown static char help[] = "Tests PetscSortIntWithPermutation().\n\n"; 3c4762a1bSJed Brown 4c4762a1bSJed Brown #include <petscsys.h> 5c4762a1bSJed Brown 6c4762a1bSJed Brown int main(int argc,char **argv) 7c4762a1bSJed Brown { 8c4762a1bSJed Brown PetscInt i; 9c4762a1bSJed Brown PetscInt x[] = {39, 9, 39, 39, 29},index[5]; 10c4762a1bSJed Brown PetscInt x2[] = {39, 9, 19, 39, 29, 39, 29, 39},index2[8]; 11c4762a1bSJed Brown 12*327415f7SBarry Smith PetscFunctionBeginUser; 139566063dSJacob Faibussowitsch PetscCall(PetscInitialize(&argc,&argv,(char*)0,help)); 149566063dSJacob Faibussowitsch PetscCall(PetscPrintf(PETSC_COMM_SELF,"1st test\n")); 15c4762a1bSJed Brown for (i=0; i<5; i++) index[i] = i; 169566063dSJacob Faibussowitsch PetscCall(PetscSortIntWithPermutation(5, x, index)); 179566063dSJacob 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]])); 18c4762a1bSJed Brown 199566063dSJacob Faibussowitsch PetscCall(PetscPrintf(PETSC_COMM_SELF,"\n2nd test\n")); 20c4762a1bSJed Brown for (i=0; i<8; i++) index2[i] = i; 219566063dSJacob Faibussowitsch PetscCall(PetscSortIntWithPermutation(8, x2, index2)); 229566063dSJacob 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]])); 239566063dSJacob Faibussowitsch PetscCall(PetscFinalize()); 24b122ec5aSJacob Faibussowitsch return 0; 25c4762a1bSJed Brown } 26c4762a1bSJed Brown 27c4762a1bSJed Brown /*TEST 28c4762a1bSJed Brown 29c4762a1bSJed Brown test: 30c4762a1bSJed Brown 31c4762a1bSJed Brown TEST*/ 32