[petsc-dev] PetscSF regression in maint

[Lisandro Dalcin]

Trivial Python script creating a periodic 1d plex with 8 segments. The
script runs fine with tag v3.7.2, however it fails with current maint.

[dalcinl at kw14821 test-dmplex]$ cat mk-plex-1d.py
import petsc4py, sys
petsc4py.init(sys.argv)

from petsc4py import PETSc
import numpy as np

nel = 8

if PETSc.COMM_WORLD.rank == 0:
    a = np.arange(0, nel, dtype='i')
    b = np.arange(1, nel+1, dtype='i')
    b[-1] = a[0]
    elems = np.vstack((a, b)).T.copy()
    coord = np.linspace(0, 1, nel+1)[:-1].reshape(-1, 1)
else:
    elems = np.zeros([0,2], dtype='i')
    coord = np.zeros([0,1])
dm = PETSc.DMPlex().createFromCellList(1, elems, coord)

dm.view()
dm.setAdjacencyUseCone(True)
dm.setAdjacencyUseClosure(False)
dm.distribute(overlap=1)
dm.view()


[dalcinl at kw14821 test-dmplex]$ mpiexec -n 2 python mk-plex-1d.py
DM Object: 2 MPI processes
  type: plex
DM_0x84000004_0 in 1 dimensions:
  0-cells: 8 0
  1-cells: 8 0
Labels:
  depth: 2 strata of sizes (8, 8)
Traceback (most recent call last):
  File "mk-plex-1d.py", line 23, in <module>
    dm.distribute(overlap=1)
  File "PETSc/DMPlex.pyx", line 465, in
petsc4py.PETSc.DMPlex.distribute (src/petsc4py.PETSc.c:226540)
petsc4py.PETSc.Error: error code 671769091
[0] DMPlexDistribute() line 1603 in
/home/devel/petsc/3.7/src/dm/impls/plex/plexdistribute.c
[0] DMPlexPartitionLabelInvert() line 1742 in
/home/devel/petsc/3.7/src/dm/impls/plex/plexpartition.c
[0] DMPlexDistributeData() line 950 in
/home/devel/petsc/3.7/src/dm/impls/plex/plexdistribute.c
[0] PetscSFDistributeSection() line 1972 in
/home/devel/petsc/3.7/src/vec/is/utils/vsectionis.c
[0] PetscSFCreateEmbeddedSF() line 874 in
/home/devel/petsc/3.7/src/vec/is/sf/interface/sf.c
[0] PetscSFBcastEnd() line 1001 in
/home/devel/petsc/3.7/src/vec/is/sf/interface/sf.c
[0] PetscSFBcastEnd_Basic() line 945 in
/home/devel/petsc/3.7/src/vec/is/sf/impls/basic/sfbasic.c
[0] PetscSFBasicGetPackInUse() line 804 in
/home/devel/petsc/3.7/src/vec/is/sf/impls/basic/sfbasic.c
[0] MPIPetsc_Type_compare() line 106 in
/home/devel/petsc/3.7/src/vec/is/sf/interface/sftype.c
[0] MPIPetsc_Type_free() line 28 in
/home/devel/petsc/3.7/src/vec/is/sf/interface/sftype.c
Traceback (most recent call last):
  File "mk-plex-1d.py", line 23, in <module>
    dm.distribute(overlap=1)
  File "PETSc/DMPlex.pyx", line 465, in
petsc4py.PETSc.DMPlex.distribute (src/petsc4py.PETSc.c:226540)
petsc4py.PETSc.Error: error code 470442499
[1] DMPlexDistribute() line 1603 in
/home/devel/petsc/3.7/src/dm/impls/plex/plexdistribute.c
[1] DMPlexPartitionLabelInvert() line 1742 in
/home/devel/petsc/3.7/src/dm/impls/plex/plexpartition.c
[1] DMPlexDistributeData() line 950 in
/home/devel/petsc/3.7/src/dm/impls/plex/plexdistribute.c
[1] PetscSFDistributeSection() line 1972 in
/home/devel/petsc/3.7/src/vec/is/utils/vsectionis.c
[1] PetscSFCreateEmbeddedSF() line 874 in
/home/devel/petsc/3.7/src/vec/is/sf/interface/sf.c
[1] PetscSFBcastEnd() line 1001 in
/home/devel/petsc/3.7/src/vec/is/sf/interface/sf.c
[1] PetscSFBcastEnd_Basic() line 945 in
/home/devel/petsc/3.7/src/vec/is/sf/impls/basic/sfbasic.c
[1] PetscSFBasicGetPackInUse() line 804 in
/home/devel/petsc/3.7/src/vec/is/sf/impls/basic/sfbasic.c
[1] MPIPetsc_Type_compare() line 106 in
/home/devel/petsc/3.7/src/vec/is/sf/interface/sftype.c
[1] MPIPetsc_Type_free() line 28 in
/home/devel/petsc/3.7/src/vec/is/sf/interface/sftype.c

===================================================================================
=   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
=   PID 31698 RUNNING AT kw14821
=   EXIT CODE: 1
=   CLEANING UP REMAINING PROCESSES
=   YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
===================================================================================


-- 
Lisandro Dalcin
============
Research Scientist
Computer, Electrical and Mathematical Sciences & Engineering (CEMSE)
Extreme Computing Research Center (ECRC)
King Abdullah University of Science and Technology (KAUST)
http://ecrc.kaust.edu.sa/

4700 King Abdullah University of Science and Technology
al-Khawarizmi Bldg (Bldg 1), Office # 0109
Thuwal 23955-6900, Kingdom of Saudi Arabia
http://www.kaust.edu.sa

Office Phone: +966 12 808-0459