fixes for VecPlaceArray/VecResetArray

Lisandro Dalcin dalcinl at gmail.com
Wed Dec 3 17:37:13 CST 2008


On Wed, Dec 3, 2008 at 8:08 PM, Matthew Knepley <knepley at gmail.com> wrote:
> On Wed, Dec 3, 2008 at 5:04 PM, Lisandro Dalcin <dalcinl at gmail.com> wrote:
>>
>> I've pushed some fixes for VecPlaceArray/VecResetArray, trying to the
>> pair a bit more smarter. However, the I'm not sure what's the use case
>> in VecPlaceArray for accepting NULL for the array. Just in case, I did
>> not changed this, but if NULL is actually pased, then the call is now
>> non-op. Is this fine? I'm missing something??
>
> I think you would like to reset the array back to NULL sometimes, to clear
> out
> memory of an old array it no longer owns. Is this right?
>

Not sure if I understood you comment, but the companion call
VecResetArray() restores the original array. Before my fix, the Vec
was left actually pointing to a NULL.



>
>>
>> --
>> Lisandro Dalcín
>> ---------------
>> Centro Internacional de Métodos Computacionales en Ingeniería (CIMEC)
>> Instituto de Desarrollo Tecnológico para la Industria Química (INTEC)
>> Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)
>> PTLC - Güemes 3450, (3000) Santa Fe, Argentina
>> Tel/Fax: +54-(0)342-451.1594
>
> --
> What most experimenters take for granted before they begin their experiments
> is infinitely more interesting than any results to which their experiments
> lead.
> -- Norbert Wiener
>



-- 
Lisandro Dalcín
---------------
Centro Internacional de Métodos Computacionales en Ingeniería (CIMEC)
Instituto de Desarrollo Tecnológico para la Industria Química (INTEC)
Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)
PTLC - Güemes 3450, (3000) Santa Fe, Argentina
Tel/Fax: +54-(0)342-451.1594




More information about the petsc-dev mailing list