Hi Matthew, After setting DMSetUseNatural to true and calling DMPlexGetGlobalToNatural, I call PestcSFView right away, it gives segmentation fault. I have also tried DMGetNaturalSF, it also gives segmentation fault when calling PetscSFView. I use PETSC_VIEWER_STDOUT_WORLD as PetscViewer
Thanks, Josh Matthew Knepley <[email protected]> 於 2023年6月9日 週五 下午1:04寫道： > On Fri, Jun 9, 2023 at 1:46 PM YuSh Lo <[email protected]> wrote: > >> Hi Barry, >> >> Is there any way to use the mapping generated by DMPlexDistribute along >> with AO? >> > > For Plex, if you turn on > > https://petsc.org/main/manualpages/DM/DMSetUseNatural/ > > before DMPlexDistribute(), it will compute and store a GlobalToNatural > map. This can be > used to map vectors back and forth, but you can extract the SF > > DMPlexGetGlobalToNaturalSF > <https://petsc.org/main/manualpages/DMPlex/DMPlexGetGlobalToNaturalSF/> > > and use that to remap your IS, by extracting the indices. > > THanks, > > Matt > > >> Thanks, >> Josh >> >> >> Barry Smith <[email protected]> 於 2023年6月9日 週五 上午10:42寫道： >> >>> >>> You might be looking for >>> https://petsc.org/release/manualpages/AO/AO/#ao >>> >>> >>> On Jun 9, 2023, at 11:02 AM, Mark Adams <[email protected]> wrote: >>> >>> An IS is just an array of integers. We need your context. >>> Is this question for sparse matrices? If so look at the documentation on >>> the AIJ matrix construction and the global vertex numbering system. >>> >>> Mark >>> >>> On Thu, Jun 8, 2023 at 1:15 PM YuSh Lo <[email protected]> wrote: >>> >>>> Hi, >>>> >>>> I have an IS that contains some vertex that is in natural numbering. >>>> How do I map them to global numbering without being distributed? >>>> >>>> Thanks, >>>> Josh >>>> >>> >>> > > -- > 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/> >
