[petsc-users] issue with DMPlexCreateFromCellListParallel

[Matthew Knepley]

On Thu, Nov 22, 2018 at 1:28 PM Andrea Ceresoli via petsc-users <
petsc-users at mcs.anl.gov> wrote:

> Good evening,
>
> I'm trying to make a parallel assembling of two networks in PETSc using
> power2.c as a template.
> Instead of reading data only from rank0 process, I want 2 processes (rank0
> and rank1) to read their own data (each process reads case9.m). As far as
> my understanding to make this work I modified the function
> DMNetworkLayoutSetUp when the function DMPlexCreateFromCellList is called
> and replaced it with DMPlexCreateFromCellListParallel when the size is
> greater than 1.
>
> Running with this new setup power2.c gives:
>
> [1]PETSC ERROR: --------------------- Error Message
> --------------------------------------------------------------
> [1]PETSC ERROR: Invalid argument
> [1]PETSC ERROR: Global vertex 9 on rank 1 was unclaimed
> [1]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html
> for trouble shooting.
> [1]PETSC ERROR: Petsc Release Version 3.10.2, unknown
> [1]PETSC ERROR: ./prova2 on a arch-linux2-c-debug named SGI-W-02 by sgilab
> Thu Nov 22 19:11:40 2018
> [1]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=0 --with-fc=0
> --download-f2cblaslapack --download-mpich
> [1]PETSC ERROR: #1 DMPlexBuildFromCellList_Parallel_Internal() line 2632
> in /home/sgilab/petsc/src/dm/impls/plex/plexcreate.c
> [1]PETSC ERROR: #2 DMPlexCreateFromCellListParallel() line 2776 in
> /home/sgilab/petsc/src/dm/impls/plex/plexcreate.c
> [1]PETSC ERROR: #3 DMNetworkLayoutSetUp_Parallel() line 328 in
> /home/sgilab/petsc/src/snes/examples/tutorials/network/power2/prova2.c
> [1]PETSC ERROR: #4 main() line 476 in
> /home/sgilab/petsc/src/snes/examples/tutorials/network/power2/prova2.c
> [1]PETSC ERROR: PETSc Option Table entries:
> [1]PETSC ERROR: -ksp_type gmres
> [1]PETSC ERROR: -pc_type bjacobi
> [1]PETSC ERROR: -snes_atol 1e-8
> [1]PETSC ERROR: -snes_converged_reason
> [1]PETSC ERROR: -snes_linesearch_type basic
> [1]PETSC ERROR: -snes_rtol 1e-20
> [1]PETSC ERROR: -snes_type newtonls
> [1]PETSC ERROR: -sub_pc_factor_mat_ordering_type qmd
> [1]PETSC ERROR: -sub_pc_type lu
> [1]PETSC ERROR: ----------------End of Error Message -------send entire
> error message to petsc-maint at mcs.anl.gov----------
>
> I was wondering whether I am using DMPlexCreateFromCellListParallel in the
> right context and for the right purpose.
>
> Do you have any suggestions? What do you think?
>

It looks like your input is wrong. Please have a look at

  src/dm/impls/plex/examples/tests/ex18.c

It has all our tests for parallel loading.

  Thanks,

     Matt


> Thank you in advance,
>
> Andrea
>


-- 
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/20181122/14a3c1f4/attachment-0001.html>