search

geodynamics / pylith

PyLith is a finite element code for the solution of dynamic and quasi-static tectonic deformation problems.
Other
154 stars 98 forks source link

Error when attempting to use multiple solution subfields with different basis order #98

Closed baagaard-usgs closed 5 years ago

baagaard-usgs commented 5 years ago

PyLith branch: knepley/features-petsc-fe

Test case

cd tests_auto/linearelasticity/nofaults-2d
make export-data
pylith gravity_incompressible.cfg gravity_incompressible_quad.cfg

Error

[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[0]PETSC ERROR: Null argument, when expecting valid pointer
[0]PETSC ERROR: Null Object: Parameter # 1
[0]PETSC ERROR: See https://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
[0]PETSC ERROR: Petsc Development GIT revision: v3.11.2-824-g5f7a1134a9  GIT Date: 2019-05-30 08:47:24 -0400
[0]PETSC ERROR: /tools/common/cig/gcc-7.3.0/bin/mpinemesis on a arch-gcc-7.3.0_debug named arling by brad Thu May 30 16:41:59 2019
[0]PETSC ERROR: Configure options --PETSC_ARCH=arch-gcc-7.3.0_debug --with-debugging=1 --with-clanguage=c --with-mpi-compilers=1 --with-shared-libraries=1 --with-64-bit-points=1 --with-large-file-io=1 --CC=mpicc --CXX=mpicxx --FC= --F77= LDFLAGS= --with-blaslapack-lib="-L/usr/lib/atlas-base -llapack_atlas -llapack -latlas -lblas" --with-lgrind=0 --with-c2html=0 --with-chaco=1 --download-chaco=1 --download-parmetis=1 --download-metis=1 --with-ml=1 --download-ml=1 --download-superlu --with-hdf5=1 --with-zlib=1 --with-hdf5-dir=/tools/common/hdf5-1.10.2/gcc-7.3.0 --download-triangle --with-fc=0
[0]PETSC ERROR: #1 PetscDualSpaceGetDimension() line 615 in /tools/common/petsc-dev/src/dm/dt/dualspace/interface/dualspace.c
[0]PETSC ERROR: #2 PetscDualSpaceGetAllPointsUnion() line 395 in /tools/common/petsc-dev/src/dm/impls/plex/plexproject.c
[0]PETSC ERROR: #3 DMProjectLocal_Generic_Plex() line 583 in /tools/common/petsc-dev/src/dm/impls/plex/plexproject.c
[0]PETSC ERROR: #4 DMProjectFieldLabelLocal_Plex() line 813 in /tools/common/petsc-dev/src/dm/impls/plex/plexproject.c
[0]PETSC ERROR: #5 DMProjectFieldLabelLocal() line 7565 in /tools/common/petsc-dev/src/dm/interface/dm.c
[0]PETSC ERROR: #6 DMPlexInsertBoundaryValuesEssentialField() line 799 in /tools/common/petsc-dev/src/dm/impls/plex/plexfem.c
[0]PETSC ERROR: #7 virtual void pylith::feassemble::ConstraintSpatialDB::setSolution(pylith::topology::Field*, double)() line 151 in /home/brad/src/cig/pylith/libsrc/pylith/feassemble/ConstraintSpatialDB.cc
baagaard-usgs commented 5 years ago

This was associated with using a discontinuous basis.