1 static char help[] = "Tests MatLoad() with uneven dimensions set in program\n\n"; 2 3 #include <petscmat.h> 4 5 int main(int argc,char **args) 6 { 7 Mat A; 8 PetscViewer fd; 9 char file[PETSC_MAX_PATH_LEN]; 10 PetscErrorCode ierr; 11 PetscBool flg; 12 PetscMPIInt rank; 13 14 ierr = PetscInitialize(&argc,&args,(char*)0,help);if (ierr) return ierr; 15 CHKERRMPI(MPI_Comm_rank(PETSC_COMM_WORLD,&rank)); 16 17 /* Determine files from which we read the matrix */ 18 CHKERRQ(PetscOptionsGetString(NULL,NULL,"-f",file,sizeof(file),&flg)); 19 PetscCheck(flg,PETSC_COMM_WORLD,PETSC_ERR_USER,"Must indicate binary file with the -f"); 20 21 /* Load matrices */ 22 CHKERRQ(PetscViewerBinaryOpen(PETSC_COMM_WORLD,file,FILE_MODE_READ,&fd)); 23 CHKERRQ(MatCreate(PETSC_COMM_WORLD,&A)); 24 CHKERRQ(MatSetFromOptions(A)); 25 CHKERRQ(MatSetBlockSize(A,2)); 26 if (rank == 0) { 27 CHKERRQ(MatSetSizes(A, 4, PETSC_DETERMINE, PETSC_DETERMINE,PETSC_DETERMINE)); 28 } else { 29 CHKERRQ(MatSetSizes(A, 8, PETSC_DETERMINE, PETSC_DETERMINE,PETSC_DETERMINE)); 30 } 31 CHKERRQ(MatLoad(A,fd)); 32 CHKERRQ(PetscViewerDestroy(&fd)); 33 CHKERRQ(MatDestroy(&A)); 34 ierr = PetscFinalize(); 35 return ierr; 36 } 37 38 /*TEST 39 40 test: 41 nsize: 2 42 args: -mat_type aij -mat_view -f ${wPETSC_DIR}/share/petsc/datafiles/matrices/ns-real-int32-float64 43 requires: double !complex !defined(PETSC_USE_64BIT_INDICES) 44 45 test: 46 suffix: 2 47 nsize: 2 48 args: -mat_type baij -mat_view -f ${wPETSC_DIR}/share/petsc/datafiles/matrices/ns-real-int32-float64 49 requires: double !complex !defined(PETSC_USE_64BIT_INDICES) 50 51 test: 52 suffix: 3 53 nsize: 2 54 args: -mat_type sbaij -mat_view -f ${wPETSC_DIR}/share/petsc/datafiles/matrices/ns-real-int32-float64 55 requires: double !complex !defined(PETSC_USE_64BIT_INDICES) 56 57 TEST*/ 58