[petsc-users] block size

[Blaise Bourdin]

search for an older thread on this list

possible fixes are to either delete / not create the .info files or to assign different prefix to each vec

B
On Jul 18, 2012, at 12:39 PM, Anton Popov wrote:

> "Comment out one of vectors, your code works fine."
> I wouldn't say "fine" because I need to store both vectors.
> 
>  So what is the workaround?
> 
> - I'm not allowed to mix vectors with different block size in the same file?
> - I'm not allowed to write more then ONE vector in a file?
> 
> 
> On 7/18/12 4:10 PM, Hong Zhang wrote:
>> You use same datafile 'vec.dat' for writing two different vectors,
>> V3 and V1:
>> 
>>  ierr = VecView(V1, view_out);                               
>>  ierr = VecView(V3, view_out);                               
>> 
>> // here, vec.dat holds V3
>> 
>> Then read it in the order
>>  ierr = VecLoad(V1, view_in); 
>> //crash here because reading V3 into V1
>> 
>>  ierr = VecLoad(V3, view_in); 
>> 
>> Comment out one of vectors, your code works fine.
>> 
>> Hong
>> 
>> 
>> 
>> 
>> On Wed, Jul 18, 2012 at 7:54 AM, Anton Popov <popov at uni-mainz.de> wrote:
>> Dear petsc team,
>> 
>> could you please tell me what's wrong with the attached example file?
>> I run it on 4 processors with petsc-3.3-p1.
>> 
>> What could error message "Local size 1000 not compatible with block size 3!" mean?
>> 
>> I've another question related to this issue. What is the real purpose of PetscViewerBinarySkipInfo function?
>> I see no reason to skip creating "info" file, because the file produced by the attached example seems to be correct.
>> 
>> Moreover, similar block size error occurs in our code while reading file with multiple vectors, irrespective whether "info" file exists or not.
>> 
>> Thank you,
>> 
>> Anton
>> 
>> 
> 
> 

-- 
Department of Mathematics and Center for Computation & Technology
Louisiana State University, Baton Rouge, LA 70803, USA
Tel. +1 (225) 578 1612, Fax  +1 (225) 578 4276 http://www.math.lsu.edu/~bourdin







-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120718/070d044c/attachment-0001.html>