search

geodynamics / pylith

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

Only one Neumann BC allowed in PyLith multiphysics #47

Closed baagaard-usgs closed 6 years ago

baagaard-usgs commented 6 years ago

In the current PyLith multiphysics branch (knepley/feature-petsc-fe), we will not allow multiple Neumann BC per field. PETSc currently assumes that there is only one Neumann BC per field and we have not yet identified a good way to allow multiple Neumann BC per field. The workaround is to use the spatial database to create the desired spatial variation in the parameters.

baagaard-usgs commented 6 years ago

This issue is no longer applicable. We have refactored the boundary integration DMPlex routines to allow multiple boundary integrations.