[petsc-users] More on BlockSize in release 3.3

Matthew Knepley knepley at gmail.com
Thu Jan 3 09:07:15 CST 2013


On Thu, Jan 3, 2013 at 8:57 AM, Matthew Knepley <knepley at gmail.com> wrote:

> On Thu, Jan 3, 2013 at 6:50 AM, Aldo Bonfiglioli <
> aldo.bonfiglioli at unibas.it> wrote:
>
>> Dear developers,
>> I have modified
>> src/vec/vec/examples/tutorials/ex14f.F
>>
>> to test 4 different combinations of
>> VecSetValues(Blocked) calls on the
>> local/global representation of the ghosted vector.
>>
>> The runtime option -job now controls the behavior.
>>
>> Things work except when using VecSetValuesBlocked
>> on the LOCAL representation of the ghosted Vector.
>>
>> Is this maybe because the local representation (lx)
>> does not seem to inherit the blocksize from the global vector (gx) it
>> comes from?
>> With 3.3 I cannot any longer set it explicitly using VecSetBlockSize.
>>
>> The modified code is enclosed.
>>
>
> The problem is that VecDuplicate() does not preserve the local block size.
> Tracking it down.
>

Pushed a fix to 3.3 which will come out with the next patch release, or
pull it (and its in dev).

    Matt


>     Matt
>
>
>> Regards,
>> Aldo
>>
>>
>>
>> --
>> Dr. Aldo Bonfiglioli
>> Associate professor of Fluid Flow Machinery
>> Scuola di Ingegneria
>> Universita' della Basilicata
>> V.le dell'Ateneo lucano, 10 85100 Potenza ITALY
>> tel:+39.0971.205203 fax:+39.0971.205215
>>
>>
>> Publications list <http://publicationslist.org/aldo.bonfiglioli>
>>
>
>
>
> --
> 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
>



-- 
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130103/f76962bb/attachment.html>


More information about the petsc-users mailing list