1 2 static char help[] = "Test memory leak when duplicating a redundant matrix.\n\n"; 3 4 /* 5 Include "petscmat.h" so that we can use matrices. 6 automatically includes: 7 petscsys.h - base PETSc routines petscvec.h - vectors 8 petscmat.h - matrices 9 petscis.h - index sets petscviewer.h - viewers 10 */ 11 #include <petscmat.h> 12 13 int main(int argc,char **args) 14 { 15 Mat A,Ar,C; 16 PetscViewer fd; /* viewer */ 17 char file[PETSC_MAX_PATH_LEN]; /* input file name */ 18 PetscErrorCode ierr; 19 PetscInt ns=2; 20 PetscMPIInt size; 21 PetscSubcomm subc; 22 PetscBool flg; 23 24 ierr = PetscInitialize(&argc,&args,(char*)0,help);if (ierr) return ierr; 25 /* 26 Determine files from which we read the two linear systems 27 (matrix and right-hand-side vector). 28 */ 29 CHKERRQ(PetscOptionsGetString(NULL,NULL,"-f0",file,sizeof(file),&flg)); 30 PetscCheck(flg,PETSC_COMM_WORLD,PETSC_ERR_USER,"Must indicate binary file with the -f0 option"); 31 CHKERRQ(PetscViewerBinaryOpen(PETSC_COMM_WORLD,file,FILE_MODE_READ,&fd)); 32 CHKERRMPI(MPI_Comm_size(PETSC_COMM_WORLD,&size)); 33 CHKERRQ(PetscPrintf(PETSC_COMM_WORLD,"Reading matrix with %d processors\n",size)); 34 CHKERRQ(MatCreate(PETSC_COMM_WORLD,&A)); 35 CHKERRQ(MatLoad(A,fd)); 36 CHKERRQ(PetscViewerDestroy(&fd)); 37 /* 38 Determines amount of subcomunicators 39 */ 40 CHKERRQ(PetscOptionsGetInt(NULL,NULL,"-nsub",&ns,NULL)); 41 CHKERRQ(PetscPrintf(PETSC_COMM_WORLD,"Splitting in %" PetscInt_FMT " subcommunicators\n",ns)); 42 CHKERRQ(PetscSubcommCreate(PetscObjectComm((PetscObject)A),&subc)); 43 CHKERRQ(PetscSubcommSetNumber(subc,ns)); 44 CHKERRQ(PetscSubcommSetType(subc,PETSC_SUBCOMM_CONTIGUOUS)); 45 CHKERRQ(PetscSubcommSetFromOptions(subc)); 46 CHKERRQ(MatCreateRedundantMatrix(A,0,PetscSubcommChild(subc),MAT_INITIAL_MATRIX,&Ar)); 47 CHKERRQ(PetscPrintf(PETSC_COMM_WORLD,"Copying matrix\n")); 48 CHKERRQ(MatDuplicate(Ar,MAT_COPY_VALUES,&C)); 49 CHKERRQ(MatAXPY(Ar,0.1,C,DIFFERENT_NONZERO_PATTERN)); 50 CHKERRQ(PetscSubcommDestroy(&subc)); 51 52 /* 53 Free memory 54 */ 55 CHKERRQ(MatDestroy(&A)); 56 CHKERRQ(MatDestroy(&Ar)); 57 CHKERRQ(MatDestroy(&C)); 58 ierr = PetscFinalize(); 59 return ierr; 60 } 61 62 /*TEST 63 64 test: 65 nsize: 4 66 requires: !complex double !defined(PETSC_USE_64BIT_INDICES) 67 args: -f0 ${wPETSC_DIR}/share/petsc/datafiles/matrices/ns-real-int32-float64 -malloc_dump 68 69 TEST*/ 70