<div dir="ltr"><div class="gmail_quote"><div dir="ltr">On Thu, Nov 22, 2018 at 1:28 PM Andrea Ceresoli via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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" target="_blank">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></div></blockquote><div><br></div><div>It looks like your input is wrong. Please have a look at</div><div><br></div><div>  src/dm/impls/plex/examples/tests/ex18.c</div><div><br></div><div>It has all our tests for parallel loading.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div dir="ltr"><div>Thank you in advance,</div><div><br></div><div>Andrea<br></div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>