Link to the source code for this example:

The console output of the program:

***************************************************************
* Running Example systems_of_equations_ex9:
*   ./example-opt
***************************************************************

 Mesh Information:
  elem_dimensions()={3}
  spatial_dimension()=3
  n_nodes()=1608
    n_local_nodes()=1608
  n_elem()=1110
    n_local_elem()=1110
    n_active_elem()=1110
  n_subdomains()=2
  n_partitions()=1
  n_processors()=1
  n_threads()=1
  processor_id()=0

 EquationSystems
  n_systems()=1
   System #0, "Elasticity"
    Type "LinearImplicit"
    Variables={ "u" "v" "w" }
    Finite Element Types="LAGRANGE"
    Approximation Orders="FIRST"
    n_dofs()=4824
    n_local_dofs()=4824
    n_constrained_dofs()=504
    n_local_constrained_dofs()=504
    n_vectors()=1
    n_matrices()=1
    DofMap Sparsity
      Average  On-Processor Bandwidth <= 70.9254
      Average Off-Processor Bandwidth <= 0
      Maximum  On-Processor Bandwidth <= 108
      Maximum Off-Processor Bandwidth <= 0
    DofMap Constraints
      Number of DoF Constraints = 504
      Average DoF Constraint Length= 2.35714
      Number of Node Constraints = 132
      Maximum Node Constraint Length= 4
      Average Node Constraint Length= 4

***************************************************************
* Done Running Example systems_of_equations_ex9:
*   ./example-opt
***************************************************************