[petsc-users] More on BlockSize in release 3.3

[Aldo Bonfiglioli]

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.

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>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ex14f.F
Type: text/x-fortran
Size: 6102 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130103/d5c043d4/attachment.bin>