[petsc-users] error in calling VecGetArrayf90()

Barry Smith bsmith at mcs.anl.gov
Wed Jan 5 13:35:24 CST 2011


  I totally lost it. Yes VecGetArrayF90 starts with 1 on each process. 


   Barry

I lost my mind and was thinking of DAVecGetArrayF90()



On Jan 5, 2011, at 1:28 PM, Jed Brown wrote:

> 1:(Iend-Istart) is "local" indexing in the global vector.  Since there us no scatter, there is no local vector and thus no indexing that includes ghosts (which would be consistent use of the term "local").  I think Ethan is correct that the only way to get proper "natural" indexing is to have an inner function.  Am I missing something?
> 
> 
>> On Jan 5, 2011 11:21 AM, "Barry Smith" <bsmith at mcs.anl.gov> wrote:
>> 
>> 
>>  I don't think this is correct. You are suppose to use the local indexing for each process. With the strange index starting at 1 instead of 0.
>> 
>> 
>>   Barry
>> 
>> 
>> On Jan 5, 2011, at 1:03 PM, Ethan Coon wrote:
>> 
>> > On all processors, the array you get is indexed:
>> ...
>> 
> 



More information about the petsc-users mailing list