Load the package:

Create a numerical region of a disk:

Extract the bounding box:

Create a numerical region of a full region in two dimensions with a bounding box:

Create a numerical region of an empty region in two dimensions:

Create a numerical region of a disk and extract the predicates and the predicate variables:

Create a numerical region of a full region in two dimensions:

Extract the ElementMesh:

Create an ElementMesh with ToElementMesh:

The ElementMesh field is now populated:

The boundary ElementMesh is now also populated:

Clear the meshes and constraints:

Create a numerical region of a full region in two dimensions:

Initially, boundary and element meshes are not populated:

Create a boundary mesh:

The boundary mesh is now populated:

Create a mesh:

Define a NumericalRegion based on the mesh:

Retrieve the predicate and the symbolic region of the NumericalRegion:

Retrieve the bounding box and variables:

The meshes are the same:

There is no boundary mesh associated:

Create a numerical region of a full region in two dimensions:

Set up variable and solution data:

Initialize the coefficients and boundary conditions that have a discontinuity:

Inspect the constraints introduced by the PDE coefficients and boundary conditions:

Visualize the boundary mesh with the autogenerated internal boundary along the discontinuity:

Visualize the full mesh: