[petsc-users] [petsc-dev] Getting 2d array with updated ghost values from DM global vector

[Barry Smith]

On Jul 3, 2012, at 3:08 AM, Blaise Bourdin wrote:

> 
> On Jul 3, 2012, at 4:10 AM, Barry Smith wrote:
> 
>> 
>> Blaise,
>> 
>>  I don't understand why the patch does anything:
>> 
>> -  *ierr = VecRestoreArray(*v,0);if (*ierr) return;
>> +  PetscScalar *fa;
>> +  *ierr = F90Array1dAccess(a,PETSC_SCALAR,(void**)&fa PETSC_F90_2PTR_PARAM(ptrd));
>> +  *ierr = VecRestoreArray(*v,&fa);if (*ierr) return;
>> *ierr = F90Array1dDestroy(&a,PETSC_SCALAR PETSC_F90_2PTR_PARAM(ptrd));
>> 
>> All that passing &fa into VecRestoreArray() does is cause fa to be zeroed. Why would that have any affect on anything?
> 
> 
> Not sure either, I quite don't understand this code, but I noticed that the logic of VecRestoreArrayF90 was different from that of DMDAVecRestoreArrayF90
> 
> src/vec/vec/interface/f90-custom/zvectorf90.c:33
> PetscScalar *fa;
> *__ierr = F90Array1dAccess(ptr,PETSC_SCALAR,(void**)&fa PETSC_F90_2PTR_PARAM(ptrd));if (*__ierr) return;
> *__ierr = F90Array1dDestroy(ptr,PETSC_SCALAR PETSC_F90_2PTR_PARAM(ptrd));if (*__ierr) return;

     It could be the above line is important; but the Accesser and restore array are not.

     I'll have Satish apply the patch.

    Thanks

     Barry


> *__ierr = VecRestoreArray(*x,&fa);
> 
> Why aren't the calls to F90Array1dAccess and F90Array1dDestroy necessary in the context of DMDAVecGetArrayF90?
> 
> Blaise
> 
> 
> 
>> 
>> Thanks
>> 
>>   Barry
>> 
>> On Jul 2, 2012, at 10:10 AM, Blaise Bourdin wrote:
>> 
>>> Hi,
>>> 
>>> There appears to be a bug in DMDAVecRestoreArrayF90. It is probably only triggered when the intel compilers. gfortran and intel seem to have very different internal implementations of fortran90 allocatable arrays. 
>>> 
>>> Developers, can you check if the attached patch makes sense? It will not fix the case of a 3d da with dof>1 since F90Array4dAccess is not implemented. Other than that, it seems to fix ex11f90 under linux and mac OS
>>> 
>>> <DMDAVecRestoreArrayF90.patch>
>>> 
>>> Blaise
>>> 
>>> 
>>> 
>>> On Jul 2, 2012, at 5:41 PM, TAY wee-beng wrote:
>>> 
>>>> On 2/7/2012 2:49 PM, Matthew Knepley wrote:
>>>>> On Mon, Jul 2, 2012 at 3:58 AM, TAY wee-beng <zonexo at gmail.com> wrote:
>>>>> Hi,
>>>>> 
>>>>> I have used DMDACreate2d for my code and then use:
>>>>> 
>>>>> call DMLocalToGlobalBegin(da,b_rhs_semi_local,INSERT_VALUES,b_rhs_semi_global,ierr)
>>>>> 
>>>>> call DMLocalToGlobalEnd(da,b_rhs_semi_local,INSERT_VALUES,b_rhs_semi_global,ierr)
>>>>> 
>>>>> to construct the global DM vector b_rhs_semi_global
>>>>> 
>>>>> Now I want to get the values with ghost values in a 2d array locally which is declared as:
>>>>> 
>>>>> real(8), allocatable :: array2d(:,:)
>>>>> 
>>>>> I guess I should use DMDAGetGhostCorners to get the corressponding indices and allocate it.  But what should I do next? How can I use something like VecGetArrayF90 to get to the pointer to access the local vector?
>>>>> 
>>>>> I can't use DMDAVecGetArrayF90/DMDAVecRestoreArrayF90 since I'm using intel fortran and they can't work. I can't use gfortran at the moment since I've problems with HYPRE with           gfortran in 3D.
>>>>> 
>>>>> Are you certain of this? That used to be true, but the current version should work for any F90.
>>>>> 
>>>>> Matt
>>>> 
>>>> I just tested 3.3-p1 and it still doesn't work (example ex11f90 in dm). Is there a chance petsc-dev can work?
>>>>> 
>>>>> Thanks
>>>>> 
>>>>> -- 
>>>>> Yours sincerely,
>>>>> 
>>>>> TAY wee-beng
>>>>> 
>>>>> 
>>>>> 
>>>>> 
>>>>> -- 
>>>>> What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
>>>>> -- Norbert Wiener
>>>> 
>>>> 
>>> 
>>> -- 
>>> 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
>>> 
>>> 
>>> 
>>> 
>>> 
>>> 
>>> 
>> 
> 
> -- 
> 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
> 
> 
> 
> 
> 
> 
>