DM+DS: Reorganized DS handling completely do that multiple regions can be supported
- A region is labeled portion of the mesh which has the same approximation space for the solution. Note that this i

DM+DS: Reorganized DS handling completely do that multiple regions can be supported
- A region is labeled portion of the mesh which has the same approximation space for the solution. Note that this is NOT the same as the label indicating the support of a field. Regions are the refinement of the field supports into unique approximation space pieces. Right now, we only support a single strategy, namely one space on the regular mesh and one on the hybrid cells, but we only have to change DMCreateDS() to extend the capability.
- Now instead of using PetscDSSetDiscretization() and DMSetDS(), you use DMSetField() and DMCreateDS(). Also, you now must use DMCopyDisc(), instead of DMSetDS(), which copies both the field and DS structure to another DM. These changes have been made in all examples.
- DMGetCellDS() returns the approximation space for that cell
- DMGetDS() return the default space since most people will not use regions. There is still a lot of work left t generalize all methods to multiple DSes, but all old code still functions.
- DMGetRegion/RegionNumDS() gets spaces by label and region number
- DMCreateDS() is the setup method which makes the regions and corresponding DSes. With some work, this could be converted to a lazy initialization.
- Cleaned up access to DS since now we can have more than one
- Projection was improved to extract the DS for the particular cells it is iterating over
- The DS are set from options during Section creation. The idea is that people using DS would also use automatic Section creation, or they would manage things themselves and call SetFromOptions manually.
- I updated the subDM and superDM constructors, but they only work for a single region

show more ...

DM: Reorganized DM fields
- A field is defined by a discretization (PetscFE or PetscFV) and a region (DMLabel or NULL)
- This will feed into the creation of PetscDS for every different cell system
-

DM: Reorganized DM fields
- A field is defined by a discretization (PetscFE or PetscFV) and a region (DMLabel or NULL)
- This will feed into the creation of PetscDS for every different cell system
- Changed interface for Get/SetField to add the region

show more ...

Merge remote-tracking branch 'origin/knepley/feature-section-label' into knepley/feature-multiple-ds

* origin/knepley/feature-section-label:
Plex: Changed interface for DMPlexCreateSection() - Now

Merge remote-tracking branch 'origin/knepley/feature-section-label' into knepley/feature-multiple-ds

* origin/knepley/feature-section-label:
Plex: Changed interface for DMPlexCreateSection() - Now get dim and numFields from the DM - Can give a label indicating where a field is supported - Updated use in examples
Plex: Move section calculation into another file

show more ...

Merge branch 'knepley/feature-label-object' [PR #1284]

* knepley/feature-label-object:
DMLabel: Make it a PetscObject

Merge remote-tracking branch 'origin/knepley/feature-label-object' into knepley/feature-section-label

* origin/knepley/feature-label-object:
DMLabel: When distributing, do not need a valid Label f

Merge remote-tracking branch 'origin/knepley/feature-label-object' into knepley/feature-section-label

* origin/knepley/feature-label-object:
DMLabel: When distributing, do not need a valid Label from processes which we are sending to
Plex: Cannot assume label "markers" exists
Forest ex2: Fixed DMLabelCreate()
DMLabel: Remove Fortran binding for DMLabelGetName()
TS ex11: Fixed interface for DMLabelCreate()
DMLabel: Make it a PetscObject - Add a CLASSID - Constructor takes an MPI_Comm - Remove interface that is in Object - Update interface in library - Fixed location of destruct in DMPlexCreateSubmeshGeneric_Interpolated() - Added DMLabelReset()

show more ...

Plex: Cannot assume label "markers" exists

Merge branch 'master' into knepley/fix-plex-partition-parallel

* master: (414 commits)
very small updates to users manual
syntax is define(XXX) in examples tests not defined(XXX)
Increase patc

Merge branch 'master' into knepley/fix-plex-partition-parallel

* master: (414 commits)
very small updates to users manual
syntax is define(XXX) in examples tests not defined(XXX)
Increase patchlevel to 3.10.3
petsc4py: updat to latest maint snapshot with SNESTEST updates
Fix link for CUDA example in PETSc from GPUs.html page Minor fixes to installation instructions
Mat: Fixed types in Pardiso interface
Converted example comment to latex
This example crashes (for unknown reasons, probably related to the handling of extended derived types passed as arguments) with the PGI Fortran compiler
Remove outdated requires: from top of several examples
remove outdated test requires: command from top of source file
Fix the documentation for -ksp_gmres_cgs_refinement_type
Commit f20c2d659b168217cb455989f28afe052d71a64c broke lex.py so the users manual was no longer being built. There reverted the changes to lex.py
Plex: Add PetscPartitioner to fort list Suggested-by: Jiaoyan Li <Jiaoyan.Li@inl.gov>
Doc: Fix spelling Suggested-by: Manuel Colera Rico <m.colera@upm.es>
Mat ex88: check in double or quad precision only
Fix error with PetscComplex when PETSc is compiled with C and used in C++
configure: fix typo endwith -> endswith
Sys ex23: Migrated to ts_max_time in yaml file as well.
AO ex5: Fixed missing CHKERRQ();
Mat test ex88: add tests for MatMultAdd and MatMultTranposeAdd
...

Conflicts:
src/dm/impls/plex/examples/tests/ex12.c

show more ...

Plex: Move section calculation into another file

DMLabel: Make it a PetscObject
- Add a CLASSID
- Constructor takes an MPI_Comm
- Remove interface that is in Object
- Update interface in library
- Fixed location of destruct in DMPlexCreateSubmeshGe

DMLabel: Make it a PetscObject
- Add a CLASSID
- Constructor takes an MPI_Comm
- Remove interface that is in Object
- Update interface in library
- Fixed location of destruct in DMPlexCreateSubmeshGeneric_Interpolated()
- Added DMLabelReset()

show more ...

Merge branch 'master' into knepley/feature-snes-patch

Not sure if this has done the right thing because there are no tests
yet. Also, the fortran stubs break (it doesn't know what to do with
SNESCo

Merge branch 'master' into knepley/feature-snes-patch

Not sure if this has done the right thing because there are no tests
yet. Also, the fortran stubs break (it doesn't know what to do with
SNESCompositeType), and I don't know how to fix it.

show more ...

Merged petsc/petsc into master

Merged petsc/petsc into master

Merge branch 'stefano_zampini/fix-p4est-squashed' [PR #1229]

* stefano_zampini/fix-p4est-squashed:
Few fixes to the p4est/p8est implementations of DMForest.

DMPlex tests ex1: add tests for meshes with p4est

move function to check cell shape from ex1 to DMPlex public interface

disable mesh distribution when p4est is used from the initial DM

flag broken

DMPlex tests ex1: add tests for meshes with p4est

move function to check cell shape from ex1 to DMPlex public interface

disable mesh distribution when p4est is used from the initial DM

flag broken label tests
flag broken DMPlexDistributeOverlap with non-conforming PLEX

disable long runtime tests

show more ...

DMPlex: customize tikzscale parameter from command line

DMBoundary: temporary fix to plex section boundary issue

DMPlexCreateSection isn't handling the case where
there are more dofs than components on a point properly.
There seems to be an ad hoc API th

DMBoundary: temporary fix to plex section boundary issue

DMPlexCreateSection isn't handling the case where
there are more dofs than components on a point properly.
There seems to be an ad hoc API that if you pass numcomps == 0
to DMAddBoundary, it constraints all components. That's not good, but
I documented it and it avoided the plex problem until we can fix it.

show more ...

DMPlex: extend CheckSymmetry to non-conforming meshes

DMPlex: fix bug in DMPlexCreatePointNumbering

DMPlex should not make any assumption on how cells,vertices,faces and edges are layout in the DAG

DMPlex: fix deadlock

PetscSectionVecView is collective!

DMCreateInterpolation_Plex: according to the docs, the scaling argument is optional

DMPlexCreateSpectralClosurePermutation: fix for 3d

minor fixes

Merge branch 'master' into haplav/feature-mataij-load-from-hdf5

because I need the test loop fix

Merge branch 'haplav/fix-plex-interpolate-interface-cones' [PR #1227]

* haplav/fix-plex-interpolate-interface-cones:
DMPlexIntepolate: fix orientation of faces in interface
Fixes orientation of face

Merge branch 'haplav/fix-plex-interpolate-interface-cones' [PR #1227]

* haplav/fix-plex-interpolate-interface-cones:
DMPlexIntepolate: fix orientation of faces in interface
Fixes orientation of faces in interface (i.e. order and orientations of the corresponding cones) so that they conform across ranks.

show more ...

Merge branch 'knepley/feature-plex-label-viz' [PR #1222]

* knepley/feature-plex-label-viz:
Plex: Added DMPlexCreateLabelField() to view labels with VTK
Plex: Automatically view label using -dm_label

Merge branch 'knepley/feature-plex-label-viz' [PR #1222]

* knepley/feature-plex-label-viz:
Plex: Added DMPlexCreateLabelField() to view labels with VTK
Plex: Automatically view label using -dm_label_view

show more ...