Re: [petsc-users] Fortran alternative for DMDAGetElements?

2023-06-30 Thread Matthew Knepley
On Fri, Jun 30, 2023 at 10:48 AM Ngoc Mai Monica Huynh < ngocmaimonica.hu...@unipv.it> wrote: > Hi, > yes, it runs and now also my code. > I moved it from the extension .F to .F90. > (With my older codes the extension .F still works fine, but not with this > one) > Yes, you need .F90 to properly

Re: [petsc-users] Fortran alternative for DMDAGetElements?

2023-06-30 Thread Ngoc Mai Monica Huynh
Hi, yes, it runs and now also my code. I moved it from the extension .F to .F90. (With my older codes the extension .F still works fine, but not with this one) Thank you for the patience and support! Monica > On 30 Jun 2023, at 15:54, Barry Smith wrote: > > > I glued your code fragme

Re: [petsc-users] Fortran alternative for DMDAGetElements?

2023-06-30 Thread Barry Smith
I glued your code fragment into a stand-alone program and it runs fine for me on 16 ranks. Does this simple program run for you? program main #include use petsc implicit none integer ierr MPI_Comm comm DM da3d ISLocalToGlobalMapping map PetscInt

Re: [petsc-users] Fortran alternative for DMDAGetElements?

2023-06-30 Thread Ngoc Mai Monica Huynh
Yes, I have the #include at the top of the code. Thank you very much for your help. I’ll let you know if I have any improvements from my side. Looking forward to hearing from you. Thanks, Monica > On 30 Jun 2023, at 15:08, Matthew Knepley wrote: > > On Fri, Jun 30, 2023 at 8:38 AM Ngoc

Re: [petsc-users] Fortran alternative for DMDAGetElements?

2023-06-30 Thread Matthew Knepley
On Fri, Jun 30, 2023 at 8:38 AM Ngoc Mai Monica Huynh < ngocmaimonica.hu...@unipv.it> wrote: > Yes, it compiles and run correctly > Okay, then we try to alter that example until it looks like your test. One thing is the #include at the top. Do you have that in your code? If Fortran does not find

Re: [petsc-users] Fortran alternative for DMDAGetElements?

2023-06-30 Thread Ngoc Mai Monica Huynh
Yes, it compiles and run correctly Monica > On 30 Jun 2023, at 12:50, Matthew Knepley wrote: > > On Fri, Jun 30, 2023 at 6:47 AM Ngoc Mai Monica Huynh > mailto:ngocmaimonica.hu...@unipv.it>> wrote: > Hi, > > I have no problem now in compiling, thank you for providing the Fortran > int

Re: [petsc-users] Fortran alternative for DMDAGetElements?

2023-06-30 Thread Matthew Knepley
On Fri, Jun 30, 2023 at 6:47 AM Ngoc Mai Monica Huynh < ngocmaimonica.hu...@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 > relate

Re: [petsc-users] Fortran alternative for DMDAGetElements?

2023-06-30 Thread Ngoc Mai Monica Huynh
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: --

Re: [petsc-users] Fortran alternative for DMDAGetElements?

2023-06-29 Thread Matthew Knepley
On Thu, Jun 29, 2023 at 3:05 PM Ngoc Mai Monica Huynh < ngocmaimonica.hu...@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 ha > scritto: > >  > >

Re: [petsc-users] Fortran alternative for DMDAGetElements?

2023-06-29 Thread Ngoc Mai Monica Huynh
Thank you. Does this mean that DMDARestoreElements() is supported as well now? Monica > Il giorno 29 giu 2023, alle ore 20:17, Barry Smith ha > scritto: > >  > > The code is ready in the branch > barry/2023-06-29/add-dmdagetelements-fortran > https://gitlab.com/petsc/petsc/-/merge_reque

Re: [petsc-users] Fortran alternative for DMDAGetElements?

2023-06-29 Thread Barry Smith
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 > wrote: > > That would be amazing, thank you very much! > Monica > >> On 29 Jun 2023, at

Re: [petsc-users] Fortran alternative for DMDAGetElements?

2023-06-29 Thread Ngoc Mai Monica Huynh
That would be amazing, thank you very much! Monica > On 29 Jun 2023, at 18:38, Barry Smith wrote: > > > I can provide the Fortran interface this afternoon. > > Barry > > >> On Jun 29, 2023, at 10:48 AM, Ngoc Mai Monica Huynh >> wrote: >> >> Hi everyone, >> >> I would need to

Re: [petsc-users] Fortran alternative for DMDAGetElements?

2023-06-29 Thread Barry Smith
I can provide the Fortran interface this afternoon. Barry > On Jun 29, 2023, at 10:48 AM, Ngoc Mai Monica Huynh > 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 t

[petsc-users] Fortran alternative for DMDAGetElements?

2023-06-29 Thread Ngoc Mai Monica Huynh
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