1 2 static char help[] = "Tests MatMult() on MatLoad() matrix \n\n"; 3 4 #include <petscmat.h> 5 6 int main(int argc,char **args) 7 { 8 Mat A; 9 Vec x,b; 10 PetscErrorCode ierr; 11 PetscViewer fd; /* viewer */ 12 char file[PETSC_MAX_PATH_LEN]; /* input file name */ 13 PetscBool flg; 14 15 ierr = PetscInitialize(&argc,&args,(char*)0,help);if (ierr) return ierr; 16 /* Determine file from which we read the matrix A */ 17 CHKERRQ(PetscOptionsGetString(NULL,NULL,"-f",file,sizeof(file),&flg)); 18 PetscCheck(flg,PETSC_COMM_WORLD,PETSC_ERR_USER,"Must indicate binary file with the -f option"); 19 20 /* Load matrix A */ 21 CHKERRQ(PetscViewerBinaryOpen(PETSC_COMM_WORLD,file,FILE_MODE_READ,&fd)); 22 CHKERRQ(MatCreate(PETSC_COMM_WORLD,&A)); 23 CHKERRQ(MatLoad(A,fd)); 24 flg = PETSC_FALSE; 25 CHKERRQ(VecCreate(PETSC_COMM_WORLD,&x)); 26 CHKERRQ(PetscOptionsGetString(NULL,NULL,"-vec",file,sizeof(file),&flg)); 27 if (flg) { 28 if (file[0] == '0') { 29 PetscInt m; 30 PetscScalar one = 1.0; 31 CHKERRQ(PetscInfo(0,"Using vector of ones for RHS\n")); 32 CHKERRQ(MatGetLocalSize(A,&m,NULL)); 33 CHKERRQ(VecSetSizes(x,m,PETSC_DECIDE)); 34 CHKERRQ(VecSetFromOptions(x)); 35 CHKERRQ(VecSet(x,one)); 36 } 37 } else { 38 CHKERRQ(VecLoad(x,fd)); 39 CHKERRQ(PetscViewerDestroy(&fd)); 40 } 41 CHKERRQ(VecDuplicate(x,&b)); 42 CHKERRQ(MatMult(A,x,b)); 43 44 /* Print (for testing only) */ 45 CHKERRQ(MatView(A,0)); 46 CHKERRQ(VecView(b,0)); 47 /* Free data structures */ 48 CHKERRQ(MatDestroy(&A)); 49 CHKERRQ(VecDestroy(&x)); 50 CHKERRQ(VecDestroy(&b)); 51 ierr = PetscFinalize(); 52 return ierr; 53 } 54