On Thu, May 7, 2020 at 11:15 AM Lisandro Dalcin <[email protected]> wrote:
> On Thu, 7 May 2020 at 17:53, Hammond, Glenn E via petsc-users < > [email protected]> wrote: > >> Jed, >> >> We call MatMeshToCellGraph() to generate the dual matrix. This function >> relies upon ParMETIS. Can you point me to a similar function that does not >> require ParMETIS ( e.g. for PTScotch)? >> >> > Maybe a starting point: > > DMPlexCreateFromCellList(comm, ..., &dm); > PetscSectionCreate(comm, &cellPartSection); > DMPlexGetPartitioner(dm, &partitioner); > PetscPartitionerSetType(partitioner, PETSCPARTITIONERPTSCOTCH); > PetscPartitionerDMPlexPartition(partitioner, dm, NULL, cellPartSection, > &cellPart); > > One minor annoyance is that the first call will need the vertex > coordinates. Matthew, any better way? > I can just fix it to allow NULL for coordinates. Matt > -- > Lisandro Dalcin > ============ > Research Scientist > Extreme Computing Research Center (ECRC) > King Abdullah University of Science and Technology (KAUST) > http://ecrc.kaust.edu.sa/ > -- 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/>
