[petsc-dev] Weird interaction between VecDuplicate and VecSetBlockSize

Barry Smith bsmith at mcs.anl.gov
Mon Apr 9 13:18:32 CDT 2012


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. 

   Barry

   This prevents the mysterious changing of the block size in a layout shared by several vectors.

> 
> We cannot certainly have a long-living vec with bs=-1, right?
> 
> 
> -- 
> 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




More information about the petsc-dev mailing list