[petsc-users] More on BlockSize in release 3.3

[Matthew Knepley]

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.

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