[petsc-users] Fortran alternative for DMDAGetElements?

[Matthew Knepley]

On Fri, Jun 30, 2023 at 6:47 AM Ngoc Mai Monica Huynh <
ngocmaimonica.huynh at unipv.it> wrote:

> Hi,
>
> I have no problem now in compiling, thank you for providing the Fortran
> interface.
> I have a follow up question.
> When running the code, I get this error, which I’m pretty sure it is
> related to DMDAGetElements(), since up to that line everything works fine.
>
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
> [0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation,
> probably memory access out of range
> [0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger
> [0]PETSC ERROR: or see https://petsc.org/release/faq/#valgrind and
> https://petsc.org/release/faq/
> [0]PETSC ERROR: ---------------------  Stack Frames
> ------------------------------------
> [0]PETSC ERROR: No error traceback is available, the problem could be in
> the main program.
> --------------------------------------------------------------------------
> MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
> with errorcode 59.
>
> NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
> You may or may not see output from other processes, depending on
> exactly when Open MPI kills them.
> —————————————————————————————————————
>
> The lines of code I’m working on are the following:
>
>       integer ierr
>
>       MPI_Comm comm
>       DM da3d
>       ISLocalToGlobalMapping map
>       PetscInt nel,nen
>       PetscInt, pointer :: e_loc(:)
>
>       call DMDACreate3d(comm,DM_BOUNDARY_NONE,DM_BOUNDARY_NONE,
>      &                  DM_BOUNDARY_NONE,DMDA_STENCIL_BOX,433,41,29,
>      &                  8,2,1,3,1,PETSC_NULL_INTEGER,
>      &                  PETSC_NULL_INTEGER,PETSC_NULL_INTEGER,
>      &                  da3d,ierr)
>       call DMSetMatType(da3d,MATIS,ierr)
>       call DMSetFromOptions(da3d,ierr)
>       call DMDASetElementType(da3d,DMDA_ELEMENT_Q1,ierr)
>       call DMSetUp(da3d,ierr)
>       call DMGetLocalToGlobalMapping(da3d,map,ierr)
>
>       call DMDAGetElements(da3d,nel,nen,e_loc,ierr)
>
> By printing in a dummy way any kind of message before and after
> DMDAGetElements(), I cannot pass over it.
> Unfortunately, I cannot run with the debug option on this machine.
> Am I calling the routine in a wrong way?
>

Does

  src/dm/tutorials/ex11f90.F90

run for you?

  Thanks,

     Matt


> Thanks,
> Monica
>
>
> On 29 Jun 2023, at 21:09, Matthew Knepley <knepley at gmail.com> wrote:
>
> On Thu, Jun 29, 2023 at 3:05 PM Ngoc Mai Monica Huynh <
> ngocmaimonica.huynh at unipv.it> wrote:
>
>> Thank you.
>> Does this mean that DMDARestoreElements() is supported as well now?
>>
>
> Yes.
>
>  Thanks,
>
>     Matt
>
>
>> Monica
>>
>>
>> Il giorno 29 giu 2023, alle ore 20:17, Barry Smith <bsmith at petsc.dev> ha
>> scritto:
>>
>> 
>>
>>   The code is ready in the branch
>> *barry/2023-06-29/add-dmdagetelements-fortran *
>> https://gitlab.com/petsc/petsc/-/merge_requests/6647
>>
>>   Barry
>>
>>
>> On Jun 29, 2023, at 12:41 PM, Ngoc Mai Monica Huynh <
>> ngocmaimonica.huynh at unipv.it> wrote:
>>
>> That would be amazing, thank you very much!
>> Monica
>>
>> On 29 Jun 2023, at 18:38, Barry Smith <bsmith at petsc.dev> wrote:
>>
>>
>> I can provide the Fortran interface this afternoon.
>>
>>  Barry
>>
>>
>> On Jun 29, 2023, at 10:48 AM, Ngoc Mai Monica Huynh <
>> ngocmaimonica.huynh at unipv.it> wrote:
>>
>> Hi everyone,
>>
>> I would need to use the routine DMDAGetElements() in our Fortran code.
>> However, as I read from the manual, there is no Fortran support for this
>> routine.
>> Is there any similar alternative there?
>>
>> Many thanks!
>> Best regards,
>> Monica Huynh
>>
>>
>>
>>
>>
>
> --
> 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
>
> https://www.cse.buffalo.edu/~knepley/
> <http://www.cse.buffalo.edu/~knepley/>
>
>
>

-- 
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

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230630/ec53ae62/attachment.html>