<div dir="ltr"><div dir="ltr"><div>Good evening,</div><div><br></div><div>I'm trying to make a parallel assembling of two networks in PETSc using power2.c as a template.</div><div>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. <br></div><div><br></div><div>Running with this new setup power2.c gives:</div><div><br></div><div>[1]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>[1]PETSC ERROR: Invalid argument<br>[1]PETSC ERROR: Global vertex 9 on rank 1 was unclaimed<br>[1]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.<br>[1]PETSC ERROR: Petsc Release Version 3.10.2, unknown <br>[1]PETSC ERROR: ./prova2 on a arch-linux2-c-debug named SGI-W-02 by sgilab Thu Nov 22 19:11:40 2018<br>[1]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=0 --with-fc=0 --download-f2cblaslapack --download-mpich<br>[1]PETSC ERROR: #1 DMPlexBuildFromCellList_Parallel_Internal() line 2632 in /home/sgilab/petsc/src/dm/impls/plex/plexcreate.c<br>[1]PETSC ERROR: #2 DMPlexCreateFromCellListParallel() line 2776 in /home/sgilab/petsc/src/dm/impls/plex/plexcreate.c<br>[1]PETSC ERROR: #3 DMNetworkLayoutSetUp_Parallel() line 328 in /home/sgilab/petsc/src/snes/examples/tutorials/network/power2/prova2.c<br>[1]PETSC ERROR: #4 main() line 476 in /home/sgilab/petsc/src/snes/examples/tutorials/network/power2/prova2.c<br>[1]PETSC ERROR: PETSc Option Table entries:<br>[1]PETSC ERROR: -ksp_type gmres<br>[1]PETSC ERROR: -pc_type bjacobi<br>[1]PETSC ERROR: -snes_atol 1e-8<br>[1]PETSC ERROR: -snes_converged_reason<br>[1]PETSC ERROR: -snes_linesearch_type basic<br>[1]PETSC ERROR: -snes_rtol 1e-20<br>[1]PETSC ERROR: -snes_type newtonls<br>[1]PETSC ERROR: -sub_pc_factor_mat_ordering_type qmd<br>[1]PETSC ERROR: -sub_pc_type lu<br>[1]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint@mcs.anl.gov----------<br></div><div><br></div><div>I was wondering whether I am using DMPlexCreateFromCellListParallel in the right context and for the right purpose.<br></div><div><br></div><div>Do you have any suggestions? What do you think?</div><div><br></div><div>Thank you in advance,</div><div><br></div><div>Andrea<br></div></div></div>