Re: [petsc-users] PetscSF Fortran interface

2023-01-12 Thread Junchao Zhang
How about this? PetscInt, pointer :: remoteOffsets(:) call PetscSFCreateRemoteOffsetsf90(distributionSF, section_filt_l, leafSection, remoteoffsets, ierr ) call PetscSFCreateSectionSFF90(distributionSF, section_filt_l, remoteoffsets, leafSection, distributionSF_dof, ierr) call

Re: [petsc-users] PetscSF Fortran interface

2023-01-12 Thread Nicholas Arnold-Medabalimi
That is what I tried initially however, I get a segmentation fault. I can confirm it's due to the remote offsets because if I try and output remoteoffsets between the Distribute Section and Create Section it throws the same segmentation fault. Thanks for the help Nicholas On Thu, Jan 12, 2023 at

Re: [petsc-users] PetscSF Fortran interface

2023-01-12 Thread Junchao Zhang
On Thu, Jan 12, 2023 at 1:28 PM Nicholas Arnold-Medabalimi < narno...@umich.edu> wrote: > Hi Junchao > > Going back to this merge request. I'm not sure I follow exactly the usage > described in the commit history. I have a c prototype of what I am trying > to do which is > >

Re: [petsc-users] PetscSF Fortran interface

2023-01-12 Thread Nicholas Arnold-Medabalimi
Hi Junchao Going back to this merge request. I'm not sure I follow exactly the usage described in the commit history. I have a c prototype of what I am trying to do which is PetscSectionCreate(PETSC_COMM_WORLD, ); PetscSFDistributeSection(redistributionSF, filteredSection_local, ,

Re: [petsc-users] PetscSF Fortran interface

2023-01-11 Thread Nicholas Arnold-Medabalimi
Hi Junchao I hope I didn't make making any errors; I created the merge request following those instructions and commented @ you. https://gitlab.com/petsc/petsc/-/merge_requests/5969 Sincerely Nicholas On Wed, Jan 11, 2023 at 9:16 AM Junchao Zhang wrote: > Hi, Nicholas, > See >

Re: [petsc-users] PetscSF Fortran interface

2023-01-11 Thread Junchao Zhang
Hi, Nicholas, See https://petsc.org/release/developers/contributing/#starting-a-new-feature-branch on how to contribute. You will need to create your own fork of petsc, then create a feature branch and your code, and then ask for a merge request to merge your branch to the petsc repo.

Re: [petsc-users] PetscSF Fortran interface

2023-01-10 Thread Nicholas Arnold-Medabalimi
Hi Junchao Apologies for not seeing that. Usually, the fortran90-specific functions have notes on the original C version, and I also can't see PetscSFCreateSectionSFF90 on the function list on the doc site. Thanks so much, and I saw your notes on the merge request. I don't suppose

Re: [petsc-users] PetscSF Fortran interface

2023-01-10 Thread Junchao Zhang
Hi, Nicholas, It seems we have implemented it, but with another name, PetscSFCreateSectionSFF90, see https://gitlab.com/petsc/petsc/-/merge_requests/5386 Try it to see if it works! --Junchao Zhang On Tue, Jan 10, 2023 at 11:45 AM Nicholas Arnold-Medabalimi < narno...@umich.edu> wrote: >

Re: [petsc-users] PetscSF Fortran interface

2023-01-10 Thread Nicholas Arnold-Medabalimi
Er to be honest I still can't get my stub to compile properly, and I don't know how to go about making a merge request. But here is what I am attempting right now. Let me know how best to proceed Its not exactly clear to me how to setup up the remote offset properly. in

Re: [petsc-users] PetscSF Fortran interface

2023-01-10 Thread Junchao Zhang
Hi, Nicholas, Could you make a merge request to PETSc and then our Fortran experts can comment on your MR? Thanks. --Junchao Zhang On Tue, Jan 10, 2023 at 11:10 AM Nicholas Arnold-Medabalimi < narno...@umich.edu> wrote: > Hi Junchao > > I think I'm almost there, but I could use some

Re: [petsc-users] PetscSF Fortran interface

2023-01-10 Thread Nicholas Arnold-Medabalimi
Hi Junchao I think I'm almost there, but I could use some insight into how to use the PETSC_F90_2PTR_PROTO and F90Array1dAccess for the remoteOffset parameter input so if another function comes up, I can add it myself without wasting your time. I am very grateful for your help and time.

Re: [petsc-users] PetscSF Fortran interface

2023-01-09 Thread Nicholas Arnold-Medabalimi
Hi Junchao Thanks again for your help in November. I've been using the your merge request branch quite heavily. Would it be possible to add a petscsfcreatesectionsf interface as well? I'm trying to write it myself using your commits as a guide but I have been struggling with handling the section

Re: [petsc-users] PetscSF Fortran interface

2022-11-19 Thread Junchao Zhang
On Sat, Nov 19, 2022 at 8:05 PM Nicholas Arnold-Medabalimi < narno...@umich.edu> wrote: > Hi > > Thanks, this is awesome. Thanks for the very prompt fix. Just one > question: will the array outputs on the fortran side copies (and need to be > deallocated) or direct access to the dmplex? > Direct

Re: [petsc-users] PetscSF Fortran interface

2022-11-19 Thread Nicholas Arnold-Medabalimi
Hi Thanks, this is awesome. Thanks for the very prompt fix. Just one question: will the array outputs on the fortran side copies (and need to be deallocated) or direct access to the dmplex? Sincerely Nicholas On Sat, Nov 19, 2022 at 8:21 PM Junchao Zhang wrote: > Hi, Nicholas, > See this

Re: [petsc-users] PetscSF Fortran interface

2022-11-19 Thread Junchao Zhang
Hi, Nicholas, See this MR, https://gitlab.com/petsc/petsc/-/merge_requests/5860 It is in testing, but you can try branch jczhang/add-petscsf-fortran to see if it works for you. Thanks. --Junchao Zhang On Sat, Nov 19, 2022 at 4:16 PM Nicholas Arnold-Medabalimi < narno...@umich.edu> wrote:

Re: [petsc-users] PetscSF Fortran interface

2022-11-19 Thread Nicholas Arnold-Medabalimi
Hi Junchao Thanks. I was wondering if there is any update on this. I may write a small interface for those two routines myself in the interim but I'd appreciate any insight you have. Sincerely Nicholas On Wed, Nov 16, 2022 at 10:39 PM Junchao Zhang wrote: > Hi, Nicholas, > I will have a

Re: [petsc-users] PetscSF Fortran interface

2022-11-16 Thread Junchao Zhang
Hi, Nicholas, I will have a look and get back to you. Thanks. --Junchao Zhang On Wed, Nov 16, 2022 at 9:27 PM Nicholas Arnold-Medabalimi < narno...@umich.edu> wrote: > Hi Petsc Users > > I'm in the process of adding some Petsc for mesh management into an > existing Fortran Solver. It has