Link to the source code for this example:

Open vector_fe_ex10 in new tab.


The console output of the program:

***************************************************************
* Running Example vector_fe_ex10:
*   ./example-opt dim=2 element_type=TRI6 -pc_type lu
***************************************************************

 Mesh Information:
  elem_dimensions()={2}
  elem_default_orders()={2}
  supported_nodal_order()=2
  spatial_dimension()=2
  n_nodes()=961
    n_local_nodes()=961
  n_elem()=450
    n_local_elem()=450
    n_active_elem()=450
  n_subdomains()=1
  n_elemsets()=0
  n_partitions()=1
  n_processors()=1
  n_threads()=1
  processor_id()=0
  is_prepared()=true
  is_replicated()=true

 EquationSystems
  n_systems()=1
   System #0, "GradDiv"
    Type "LinearImplicit"
    Variables="u"
    Finite Element Types="RAVIART_THOMAS"
    Approximation Orders="FIRST"
    n_dofs()=705
    n_local_dofs()=705
    max(n_local_dofs())=705
    n_constrained_dofs()=0
    n_local_constrained_dofs()=0
    max(local unconstrained dofs)=705
    n_vectors()=1
    n_matrices()=1
    DofMap Sparsity
      Average  On-Processor Bandwidth <= 4.82979
      Average Off-Processor Bandwidth <= 0
      Maximum  On-Processor Bandwidth <= 5
      Maximum Off-Processor Bandwidth <= 0
    DofMap Constraints
      Number of DoF Constraints = 0

~~ Vector field (u) ~~
L2 error is: 0.171261
HDiv semi-norm error is: 0.872515
HDiv error is: 0.889164

***************************************************************
* Done Running Example vector_fe_ex10:
*   ./example-opt dim=2 element_type=TRI6 -pc_type lu
***************************************************************
***************************************************************
* Running Example vector_fe_ex10:
*   ./example-opt dim=2 element_type=TRI7 -pc_type lu
***************************************************************

 Mesh Information:
  elem_dimensions()={2}
  elem_default_orders()={3}
  supported_nodal_order()=3
  spatial_dimension()=2
  n_nodes()=1411
    n_local_nodes()=1411
  n_elem()=450
    n_local_elem()=450
    n_active_elem()=450
  n_subdomains()=1
  n_elemsets()=0
  n_partitions()=1
  n_processors()=1
  n_threads()=1
  processor_id()=0
  is_prepared()=true
  is_replicated()=true

 EquationSystems
  n_systems()=1
   System #0, "GradDiv"
    Type "LinearImplicit"
    Variables="u"
    Finite Element Types="RAVIART_THOMAS"
    Approximation Orders="FIRST"
    n_dofs()=705
    n_local_dofs()=705
    max(n_local_dofs())=705
    n_constrained_dofs()=0
    n_local_constrained_dofs()=0
    max(local unconstrained dofs)=705
    n_vectors()=1
    n_matrices()=1
    DofMap Sparsity
      Average  On-Processor Bandwidth <= 4.82979
      Average Off-Processor Bandwidth <= 0
      Maximum  On-Processor Bandwidth <= 5
      Maximum Off-Processor Bandwidth <= 0
    DofMap Constraints
      Number of DoF Constraints = 0

~~ Vector field (u) ~~
L2 error is: 0.171261
HDiv semi-norm error is: 0.872515
HDiv error is: 0.889164

***************************************************************
* Done Running Example vector_fe_ex10:
*   ./example-opt dim=2 element_type=TRI7 -pc_type lu
***************************************************************
***************************************************************
* Running Example vector_fe_ex10:
*   ./example-opt dim=2 element_type=QUAD8 -pc_type lu
***************************************************************

 Mesh Information:
  elem_dimensions()={2}
  elem_default_orders()={2}
  supported_nodal_order()=2
  spatial_dimension()=2
  n_nodes()=736
    n_local_nodes()=736
  n_elem()=225
    n_local_elem()=225
    n_active_elem()=225
  n_subdomains()=1
  n_elemsets()=0
  n_partitions()=1
  n_processors()=1
  n_threads()=1
  processor_id()=0
  is_prepared()=true
  is_replicated()=true

 EquationSystems
  n_systems()=1
   System #0, "GradDiv"
    Type "LinearImplicit"
    Variables="u"
    Finite Element Types="RAVIART_THOMAS"
    Approximation Orders="FIRST"
    n_dofs()=480
    n_local_dofs()=480
    max(n_local_dofs())=480
    n_constrained_dofs()=0
    n_local_constrained_dofs()=0
    max(local unconstrained dofs)=480
    n_vectors()=1
    n_matrices()=1
    DofMap Sparsity
      Average  On-Processor Bandwidth <= 6.625
      Average Off-Processor Bandwidth <= 0
      Maximum  On-Processor Bandwidth <= 7
      Maximum Off-Processor Bandwidth <= 0
    DofMap Constraints
      Number of DoF Constraints = 0

~~ Vector field (u) ~~
L2 error is: 0.171869
HDiv semi-norm error is: 1.06744
HDiv error is: 1.08118

***************************************************************
* Done Running Example vector_fe_ex10:
*   ./example-opt dim=2 element_type=QUAD8 -pc_type lu
***************************************************************
***************************************************************
* Running Example vector_fe_ex10:
*   ./example-opt dim=2 element_type=QUAD9 -pc_type lu
***************************************************************

 Mesh Information:
  elem_dimensions()={2}
  elem_default_orders()={2}
  supported_nodal_order()=2
  spatial_dimension()=2
  n_nodes()=961
    n_local_nodes()=961
  n_elem()=225
    n_local_elem()=225
    n_active_elem()=225
  n_subdomains()=1
  n_elemsets()=0
  n_partitions()=1
  n_processors()=1
  n_threads()=1
  processor_id()=0
  is_prepared()=true
  is_replicated()=true

 EquationSystems
  n_systems()=1
   System #0, "GradDiv"
    Type "LinearImplicit"
    Variables="u"
    Finite Element Types="RAVIART_THOMAS"
    Approximation Orders="FIRST"
    n_dofs()=480
    n_local_dofs()=480
    max(n_local_dofs())=480
    n_constrained_dofs()=0
    n_local_constrained_dofs()=0
    max(local unconstrained dofs)=480
    n_vectors()=1
    n_matrices()=1
    DofMap Sparsity
      Average  On-Processor Bandwidth <= 6.625
      Average Off-Processor Bandwidth <= 0
      Maximum  On-Processor Bandwidth <= 7
      Maximum Off-Processor Bandwidth <= 0
    DofMap Constraints
      Number of DoF Constraints = 0

~~ Vector field (u) ~~
L2 error is: 0.171869
HDiv semi-norm error is: 1.06744
HDiv error is: 1.08118

***************************************************************
* Done Running Example vector_fe_ex10:
*   ./example-opt dim=2 element_type=QUAD9 -pc_type lu
***************************************************************
***************************************************************
* Running Example vector_fe_ex10:
*   ./example-opt dim=2 order=3 element_type=TRI6 -pc_type lu
***************************************************************

 Mesh Information:
  elem_dimensions()={2}
  elem_default_orders()={2}
  supported_nodal_order()=2
  spatial_dimension()=2
  n_nodes()=961
    n_local_nodes()=961
  n_elem()=450
    n_local_elem()=450
    n_active_elem()=450
  n_subdomains()=1
  n_elemsets()=0
  n_partitions()=1
  n_processors()=1
  n_threads()=1
  processor_id()=0
  is_prepared()=true
  is_replicated()=true

 EquationSystems
  n_systems()=1
   System #0, "GradDiv"
    Type "LinearImplicit"
    Variables="u"
    Finite Element Types="RAVIART_THOMAS"
    Approximation Orders="THIRD"
    n_dofs()=4815
    n_local_dofs()=4815
    max(n_local_dofs())=4815
    n_constrained_dofs()=0
    n_local_constrained_dofs()=0
    max(local unconstrained dofs)=4815
    n_vectors()=1
    n_matrices()=1
    DofMap Sparsity
      Average  On-Processor Bandwidth <= 19.8224
      Average Off-Processor Bandwidth <= 0
      Maximum  On-Processor Bandwidth <= 27
      Maximum Off-Processor Bandwidth <= 0
    DofMap Constraints
      Number of DoF Constraints = 0

~~ Vector field (u) ~~
L2 error is: 0.000476859
HDiv semi-norm error is: 0.00418619
HDiv error is: 0.00421326

***************************************************************
* Done Running Example vector_fe_ex10:
*   ./example-opt dim=2 order=3 element_type=TRI6 -pc_type lu
***************************************************************
***************************************************************
* Running Example vector_fe_ex10:
*   ./example-opt dim=2 order=3 element_type=TRI7 -pc_type lu
***************************************************************

 Mesh Information:
  elem_dimensions()={2}
  elem_default_orders()={3}
  supported_nodal_order()=3
  spatial_dimension()=2
  n_nodes()=1411
    n_local_nodes()=1411
  n_elem()=450
    n_local_elem()=450
    n_active_elem()=450
  n_subdomains()=1
  n_elemsets()=0
  n_partitions()=1
  n_processors()=1
  n_threads()=1
  processor_id()=0
  is_prepared()=true
  is_replicated()=true

 EquationSystems
  n_systems()=1
   System #0, "GradDiv"
    Type "LinearImplicit"
    Variables="u"
    Finite Element Types="RAVIART_THOMAS"
    Approximation Orders="THIRD"
    n_dofs()=4815
    n_local_dofs()=4815
    max(n_local_dofs())=4815
    n_constrained_dofs()=0
    n_local_constrained_dofs()=0
    max(local unconstrained dofs)=4815
    n_vectors()=1
    n_matrices()=1
    DofMap Sparsity
      Average  On-Processor Bandwidth <= 19.8224
      Average Off-Processor Bandwidth <= 0
      Maximum  On-Processor Bandwidth <= 27
      Maximum Off-Processor Bandwidth <= 0
    DofMap Constraints
      Number of DoF Constraints = 0

~~ Vector field (u) ~~
L2 error is: 0.000476859
HDiv semi-norm error is: 0.00418619
HDiv error is: 0.00421326

***************************************************************
* Done Running Example vector_fe_ex10:
*   ./example-opt dim=2 order=3 element_type=TRI7 -pc_type lu
***************************************************************
***************************************************************
* Running Example vector_fe_ex10:
*   ./example-opt dim=2 order=3 element_type=QUAD8 -pc_type lu
***************************************************************

 Mesh Information:
  elem_dimensions()={2}
  elem_default_orders()={2}
  supported_nodal_order()=2
  spatial_dimension()=2
  n_nodes()=736
    n_local_nodes()=736
  n_elem()=225
    n_local_elem()=225
    n_active_elem()=225
  n_subdomains()=1
  n_elemsets()=0
  n_partitions()=1
  n_processors()=1
  n_threads()=1
  processor_id()=0
  is_prepared()=true
  is_replicated()=true

 EquationSystems
  n_systems()=1
   System #0, "GradDiv"
    Type "LinearImplicit"
    Variables="u"
    Finite Element Types="RAVIART_THOMAS"
    Approximation Orders="THIRD"
    n_dofs()=4140
    n_local_dofs()=4140
    max(n_local_dofs())=4140
    n_constrained_dofs()=0
    n_local_constrained_dofs()=0
    max(local unconstrained dofs)=4140
    n_vectors()=1
    n_matrices()=1
    DofMap Sparsity
      Average  On-Processor Bandwidth <= 30.3913
      Average Off-Processor Bandwidth <= 0
      Maximum  On-Processor Bandwidth <= 45
      Maximum Off-Processor Bandwidth <= 0
    DofMap Constraints
      Number of DoF Constraints = 0

~~ Vector field (u) ~~
L2 error is: 0.000327013
HDiv semi-norm error is: 0.00205269
HDiv error is: 0.00207857

***************************************************************
* Done Running Example vector_fe_ex10:
*   ./example-opt dim=2 order=3 element_type=QUAD8 -pc_type lu
***************************************************************
***************************************************************
* Running Example vector_fe_ex10:
*   ./example-opt dim=2 order=3 element_type=QUAD9 -pc_type lu
***************************************************************

 Mesh Information:
  elem_dimensions()={2}
  elem_default_orders()={2}
  supported_nodal_order()=2
  spatial_dimension()=2
  n_nodes()=961
    n_local_nodes()=961
  n_elem()=225
    n_local_elem()=225
    n_active_elem()=225
  n_subdomains()=1
  n_elemsets()=0
  n_partitions()=1
  n_processors()=1
  n_threads()=1
  processor_id()=0
  is_prepared()=true
  is_replicated()=true

 EquationSystems
  n_systems()=1
   System #0, "GradDiv"
    Type "LinearImplicit"
    Variables="u"
    Finite Element Types="RAVIART_THOMAS"
    Approximation Orders="THIRD"
    n_dofs()=4140
    n_local_dofs()=4140
    max(n_local_dofs())=4140
    n_constrained_dofs()=0
    n_local_constrained_dofs()=0
    max(local unconstrained dofs)=4140
    n_vectors()=1
    n_matrices()=1
    DofMap Sparsity
      Average  On-Processor Bandwidth <= 30.3913
      Average Off-Processor Bandwidth <= 0
      Maximum  On-Processor Bandwidth <= 45
      Maximum Off-Processor Bandwidth <= 0
    DofMap Constraints
      Number of DoF Constraints = 0

~~ Vector field (u) ~~
L2 error is: 0.000327013
HDiv semi-norm error is: 0.00205269
HDiv error is: 0.00207857

***************************************************************
* Done Running Example vector_fe_ex10:
*   ./example-opt dim=2 order=3 element_type=QUAD9 -pc_type lu
***************************************************************
***************************************************************
* Running Example vector_fe_ex10:
*   ./example-opt dim=3 element_type=TET14 grid_size=6 -pc_type lu
***************************************************************

 Mesh Information:
  elem_dimensions()={3}
  elem_default_orders()={3}
  supported_nodal_order()=3
  spatial_dimension()=3
  n_nodes()=19045
    n_local_nodes()=19045
  n_elem()=5184
    n_local_elem()=5184
    n_active_elem()=5184
  n_subdomains()=1
  n_elemsets()=0
  n_partitions()=1
  n_processors()=1
  n_threads()=1
  processor_id()=0
  is_prepared()=true
  is_replicated()=true

 EquationSystems
  n_systems()=1
   System #0, "GradDiv"
    Type "LinearImplicit"
    Variables="u"
    Finite Element Types="RAVIART_THOMAS"
    Approximation Orders="FIRST"
    n_dofs()=10800
    n_local_dofs()=10800
    max(n_local_dofs())=10800
    n_constrained_dofs()=0
    n_local_constrained_dofs()=0
    max(local unconstrained dofs)=10800
    n_vectors()=1
    n_matrices()=1
    DofMap Sparsity
      Average  On-Processor Bandwidth <= 6.76
      Average Off-Processor Bandwidth <= 0
      Maximum  On-Processor Bandwidth <= 7
      Maximum Off-Processor Bandwidth <= 0
    DofMap Constraints
      Number of DoF Constraints = 0

~~ Vector field (u) ~~
L2 error is: 0.399985
HDiv semi-norm error is: 1.63
HDiv error is: 1.67836

***************************************************************
* Done Running Example vector_fe_ex10:
*   ./example-opt dim=3 element_type=TET14 grid_size=6 -pc_type lu
***************************************************************
***************************************************************
* Running Example vector_fe_ex10:
*   ./example-opt dim=3 element_type=HEX27 -pc_type lu
***************************************************************

 Mesh Information:
  elem_dimensions()={3}
  elem_default_orders()={2}
  supported_nodal_order()=2
  spatial_dimension()=3
  n_nodes()=29791
    n_local_nodes()=29791
  n_elem()=3375
    n_local_elem()=3375
    n_active_elem()=3375
  n_subdomains()=1
  n_elemsets()=0
  n_partitions()=1
  n_processors()=1
  n_threads()=1
  processor_id()=0
  is_prepared()=true
  is_replicated()=true

 EquationSystems
  n_systems()=1
   System #0, "GradDiv"
    Type "LinearImplicit"
    Variables="u"
    Finite Element Types="RAVIART_THOMAS"
    Approximation Orders="FIRST"
    n_dofs()=10800
    n_local_dofs()=10800
    max(n_local_dofs())=10800
    n_constrained_dofs()=0
    n_local_constrained_dofs()=0
    max(local unconstrained dofs)=10800
    n_vectors()=1
    n_matrices()=1
    DofMap Sparsity
      Average  On-Processor Bandwidth <= 10.375
      Average Off-Processor Bandwidth <= 0
      Maximum  On-Processor Bandwidth <= 11
      Maximum Off-Processor Bandwidth <= 0
    DofMap Constraints
      Number of DoF Constraints = 0

~~ Vector field (u) ~~
L2 error is: 0.243059
HDiv semi-norm error is: 1.50958
HDiv error is: 1.52903

***************************************************************
* Done Running Example vector_fe_ex10:
*   ./example-opt dim=3 element_type=HEX27 -pc_type lu
***************************************************************