[petsc-dev] Weird interaction between VecDuplicate and VecSetBlockSize

[Lisandro Dalcin]

On 9 April 2012 21:18, Barry Smith <bsmith at mcs.anl.gov> wrote:
>
> On Mar 30, 2012, at 12:44 PM, Lisandro Dalcin wrote:
>
>> On 30 March 2012 19:34, Barry Smith <bsmith at mcs.anl.gov> wrote:
>>>
>>> On Mar 30, 2012, at 11:07 AM, Lisandro Dalcin wrote:
>>>
>>>> On 30 March 2012 19:00, Lisandro Dalcin <dalcinl at gmail.com> wrote:
>>
>> Are you sure? How are we going to support the usage below?:
>>
>> VecCreateSeq(comm,n,&vec)
>> VecSetBlockSize(vec,bs)
>
>    This is not longer supported. If you wish to set the blocksize of a vector you cannot use the convenience routines VecCreateSeq() or VecCreateMPI(). Nor can you can a Vec block size later.
>

Not that I'm asking for it, but have you considered
VecCreateSeq(comm,bs,n,&vec) and VecCreateMPI(comm,&bs,n,N,&vec) ?
VecCreateXXXWithArray() was changed to take the block size...



-- 
Lisandro Dalcin
---------------
CIMEC (INTEC/CONICET-UNL)
Predio CONICET-Santa Fe
Colectora RN 168 Km 472, Paraje El Pozo
3000 Santa Fe, Argentina
Tel: +54-342-4511594 (ext 1011)
Tel/Fax: +54-342-4511169