[petsc-users] Petsc memory consumption keep increasing in my loop

[Medane TCHAKOROM]

Hello,

I have problem with a code i'am working on.

To illustrate my problem, here is an example:


int main(int argc, char *argv[])
{

     PetscErrorCode ierr;

     ierr = PetscInitialize(&argc, &argv, (char *)0, NULL);
     if (ierr)
         return ierr;

     int i = 0;
     for (i = 0; i < 1; i++)
     {
         Mat A;
         ierr = MatCreate(PETSC_COMM_WORLD, &A);
         CHKERRQ(ierr);
         ierr = MatSetSizes(A, 16, 16, PETSC_DECIDE,PETSC_DECIDE);
         CHKERRQ(ierr);
         ierr = MatSetFromOptions(A);
         CHKERRQ(ierr);
         ierr = MatSetUp(A);
         CHKERRQ(ierr);


         ierr = MatAssemblyBegin(A, MAT_FINAL_ASSEMBLY);
         CHKERRQ(ierr);
         ierr = MatAssemblyEnd(A, MAT_FINAL_ASSEMBLY);
         CHKERRQ(ierr);



         /// SOME CODE HERE....

         MatDestroy(&A);


     }


     FILE *fPtr;
     fPtr = fopen("petsc_dump_file.txt", "a");
     PetscMallocDump(fPtr);
     fclose(fPtr);

     ierr = PetscFinalize();
     CHKERRQ(ierr);

     return 0;
}



The problem is , in the loop, the memory consumption keep increasing 
till the end of the program.

I checked memory leak with PetscMallocDump, and found out that the 
problem may be due to matrix creation.

I'am new to Petsc and i don't know if i'am doing something wrong. Thanks


Médane