[petsc-users] array still valid after DMDAVecRestoreArrayF90()?

[Jed Brown]

No, you should only Get an array once.
On Jul 16, 2012 8:15 PM, "Panruo Wu" <pwu at mymail.mines.edu> wrote:

> How about in the inner block use a different array name
> like array2?
>
> Panruo Wu
>
> On Sun, Jul 15, 2012 at 10:39 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:
>
>> On Sun, Jul 15, 2012 at 11:20 PM, Panruo Wu <pwu at mymail.mines.edu> wrote:
>>
>>> Hello,
>>>
>>> For some reason, it would be convenient for me to do something
>>> like
>>>
>>> call DMDAVecGetArrayF90(da, global, array, ierr)
>>>
>>> call DMDAVecGetArrayF90(da, global, array, ierr)
>>> ! read-only access to array here.
>>> call DMDAVecRestoreArrayF90(da, global, array,ierr)
>>> ! read&write access to array here.
>>> call DMDAVecRestoreArrayF90(da, global, array, ierr)
>>>
>>
>> No, it is not okay to lie about the array status.
>>
>>
>>>
>>> Will this piece of code do what I expect?  I wrote a simply
>>> program which gave positive answer but I want to know if
>>> it's guaranteed to work or only happened to be working.
>>>
>>> Thanks!
>>> Panruo Wu
>>>
>>> On Sat, Jul 14, 2012 at 1:30 PM, Panruo Wu <pwu at mymail.mines.edu> wrote:
>>>
>>>> Thank you Barry!
>>>>
>>>> Panruo
>>>>
>>>>
>>>> On Sat, Jul 7, 2012 at 3:10 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:
>>>>
>>>>>
>>>>>    So long as you have the same boundary types and the same array
>>>>> sizes in the i and j direction they give the same distribution.
>>>>>
>>>>>     Barry
>>>>> On Jul 7, 2012, at 3:58 PM, Panruo Wu wrote:
>>>>>
>>>>> > Hello,
>>>>> >
>>>>> > If I create 2 DAs with (almost) identical parameters except DA name
>>>>> > and dof like:
>>>>> >
>>>>> > call DMDACreate2d(PETSC_COMM_WORLD, DMDA_BOUNDARY_GHOSTED,
>>>>> > DMDA_BOUNDARY_GHOSTED, &
>>>>> >   stype, M, N, m, n, dof1, s &
>>>>> >   PETSC_NULL_INTEGER, PETSC_NULL_INTEGER, &
>>>>> >   da1, ierr)
>>>>> >
>>>>> >
>>>>> >
>>>>> > call DMDACreate2d(PETSC_COMM_WORLD, DMDA_BOUNDARY_GHOSTED,
>>>>> > DMDA_BOUNDARY_GHOSTED, &
>>>>> >   stype,  M, N, m, n, dof2, s &
>>>>> >   PETSC_NULL_INTEGER, PETSC_NULL_INTEGER, &
>>>>> >   da2, ierr)
>>>>> >
>>>>> >
>>>>> > my question is, will the two DAs have the same distribution scheme?
>>>>> > Specifically,
>>>>> > will the DMDAGetCorners() give the same results when querying da1 &
>>>>> da2?
>>>>> >
>>>>> > Thanks,
>>>>> > Panruo Wu
>>>>>
>>>>>
>>>>
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120716/318ce972/attachment.html>