[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: memory bleeding

Ok! thanks! It runs so good!!

I supposed it when I read yours fisrts mails but I tried it last day and I didn't write the code correctly. When I have readed this mail I've tried again and I've found the error!

jordi


En/na Matthew Knepley ha escrit: 
The code you have there will create a matrix for every loop iteration,
and will only destroy one at the end. The new MatMatMult() interface
has an argument which allows the user to pass in a matrix.

   Matt

On 5/10/06, *Jordi Marcé Nogué* <jordi.marce@xxxxxxx <mailto:jordi.marce@xxxxxxx>> wrote:

    Yes, of course I destroy the matrix.... The scheme of my code when I
    create a diagonal Mass-lumping matrix is the code below.

    In MatSeqAIJSetPreallocation(M_aux1,6,PETSC_NULL); the number is 6
    because in a general coordinates it's possible this matrix changes in a
    6x6 full matrix.

    "element3D *p_element = lfiber[i].getElement(e)"  is a internal
    procedure to obtain information about the element

    ------------------------------------------------------------------

    Mat M_aux1, M_aux2;

    MatCreateSeqAIJ(PETSC_COMM_SELF, 6, 6, 6, PETSC_NULL, &M_aux1);
    MatCreateSeqAIJ(PETSC_COMM_SELF, 6, 6, 6, PETSC_NULL, &M_aux2);

    MatSeqAIJSetPreallocation(M_aux1,6,PETSC_NULL);
    MatSeqAIJSetPreallocation(M_aux2,6,PETSC_NULL);

    MatSetFromOptions(M_aux1);
    MatSetOption(M_aux1, MAT_SYMMETRIC);
    MatSetOption(M_aux1, MAT_IGNORE_ZERO_ENTRIES);

    MatSetFromOptions(M_aux2);
    MatSetOption(M_aux2, MAT_SYMMETRIC);
    MatSetOption(M_aux2, MAT_IGNORE_ZERO_ENTRIES);


    for(uint32_t i=0; i<nfiber;i++)
    {
         for(uint32_t e=0; e<lfiber[i].nelements;e++)
         {

          MatZeroentries(M_aux1);
          MatZeroentries(M_aux2);

          element3D *p_element = lfiber[i].getElement(e);

          p_element->updateTs();

          MatSetValue(M_aux1, 0, 0, p_element->L0 * value.pho / 6.,
    INSERT_VALUES);
          MatSetValue(M_aux1, 1, 1, p_element->L0 * value.pho / 6.,
    INSERT_VALUES);
          MatSetValue(M_aux1, 2, 2, p_element->L0 * value.pho / 6.,
    INSERT_VALUES);
          MatSetValue(M_aux1, 3, 3, p_element->L0 * value.pho / 6.,
    INSERT_VALUES);
          MatSetValue(M_aux1, 4, 4, p_element->L0 * value.pho / 6.,
    INSERT_VALUES);
          MatSetValue(M_aux1, 5, 5, p_element->L0 * value.pho / 6.,
    INSERT_VALUES);


          MatAssemblyBegin(M_aux1,MAT_FINAL_ASSEMBLY);
          MatAssemblyEnd(M_aux1,MAT_FINAL_ASSEMBLY);


          MatMatMult(M_aux1,p_element->T,&M_aux2);
              MatMatMult(p_element->TT,M_aux2,&M_aux1);

          }

          // here I work with the matrix M_aux2, but in this point the
    memory
          // is constant. The code doesn't waste memory
    }

    MatDestroy(M_aux1);
    MatDestroy(M_aux2);

    ------------------------


    Thanks,

    best regards


    --
    Jordi Marcé-Nogué
    Dept. Resistència de Materials i Estructures a l'Enginyeria
    Universitat Politècnica de Catalunya (UPC)

    Edifici T45 - despatx 137
    ETSEIAT (Terrassa)

    phone: +34 937 398 728
    mail: jordi.marce@xxxxxxx <mailto:jordi.marce@xxxxxxx>




--
"Failure has a thousand explanations. Success doesn't need one" -- Sir Alec Guiness 


--
Jordi Marcé-Nogué
Dept. Resistència de Materials i Estructures a l'Enginyeria
Universitat Politècnica de Catalunya (UPC)

Edifici T45 - despatx 137
ETSEIAT (Terrassa)

phone: +34 937 398 728
mail: jordi.marce@xxxxxxx