Contact Requirements Traceability Matrix

[frictionless-mortar-3d]
  type = 'Exodiff'
  input = 'frictionless-mortar-3d.i'
  min_ad_size = 100
  exodiff = 'frictionless-mortar-3d_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'The system shall solve a 3D frictionless bouncing block problem with mortar constraint'
[]

[frictionless-mortar-3d-test-derivative-trimming]
  type = 'RunApp'
  input = 'frictionless-mortar-3d-test-derivative-trimming.i'
  requirement = 'The system shall solve a 3D frictionless bouncing block problem with mortar constraints where the primary suface is composed of a single element and the secondary side is composed of first order faces with a required derivative container size of less than 50'
  ad_indexing_type = 'global'
[]

[frictionless-mortar-3d-action]
  type = 'CSVDiff'
  input = 'frictionless-mortar-3d-action.i'
  min_ad_size = 100
  csvdiff = 'frictionless-mortar-3d-action_out_contact-pressure_0002.csv'
  mesh_mode = 'REPLICATED'
  requirement = 'The system shall solve a 3D frictionless bouncing block problem with mortar constraint using the contact action'
  rel_err = 1.0e-5
  abs_zero = 1.0e-5
[]

[frictionless-mortar-3d-action-correct-edge-dropping]
  type = 'CSVDiff'
  input = 'frictionless-mortar-3d-action.i'
  min_ad_size = 100
  csvdiff = 'frictionless-mortar-3d-action-correct-edge-dropping_out_contact-pressure_0002.csv'
  mesh_mode = 'REPLICATED'
  requirement = 'The system shall solve a 3D frictionless bouncing block problem with mortar constraint using the contact action and selecting the temporary flag correct edge dropping'
  cli_args = 'Contact/mortar/correct_edge_dropping=true Outputs/file_base=frictionless-mortar-3d-action-correct-edge-dropping_out'
  rel_err = 1.0e-5
  abs_zero = 1.0e-5
[]

[frictionless-mortar-3d-debug-mesh]
  type = 'Exodiff'
  input = 'frictionless-mortar-3d.i'
  min_ad_size = 100
  exodiff = 'mortar_segment_mesh.e'
  mesh_mode = 'REPLICATED'
  cli_args = 'Executioner/end_time=0.5 Executioner/dtmin=0.5 Constraints/normal_lm/debug_mesh=true Constraints/normal_x/debug_mesh=true Constraints/normal_y/debug_mesh=true Constraints/normal_z/debug_mesh=true'
  requirement = 'The system shall solve a 3D frictionless bouncing block problem with mortar constraint and output the mortar segment mesh for debugging purposes'
  max_parallel = 1
[]

[frictionless-mortar-3d-debug-mesh-geometry]
  type = 'Exodiff'
  input = 'frictionless-mortar-3d.i'
  min_ad_size = 100
  exodiff = 'nodal_geometry_only.e'
  mesh_mode = 'REPLICATED'
  cli_args = 'Executioner/end_time=0.5 Executioner/dtmin=0.5 Constraints/normal_lm/debug_mesh=true Constraints/normal_x/debug_mesh=true Constraints/normal_y/debug_mesh=true Constraints/normal_z/debug_mesh=true Outputs/exodus=false Problem/solve=false'
  requirement = 'The system shall solve a 3D frictionless bouncing block problem with mortar constraint and output the normal and tangent vector fields for debugging purposes'
  max_parallel = 1
[]

[frictional-mortar-3d]
  type = 'Exodiff'
  input = 'frictional-mortar-3d.i'
  min_ad_size = 100
  exodiff = 'frictional-mortar-3d_out.e'
  requirement = 'The system shall solve a 3D frictional bouncing block problem with mortar constraint using nodal-attached geometry'
  rel_err = 1.0e-5
  abs_zero = 1.0e-5
  max_parallel = 8
[]

[frictional-mortar-3d-action]
  type = 'CSVDiff'
  input = 'frictional-mortar-3d-action.i'
  min_ad_size = 100
  csvdiff = 'frictional-mortar-3d-action_out_frictional-pressure_0001.csv frictional-mortar-3d-action_out_frictional-pressure-3d_0001.csv frictional-mortar-3d-action_out_tangent_x_0001.csv frictional-mortar-3d-action_out_tangent_y_0001.csv'
  requirement = 'The system shall solve a 3D frictional bouncing block problem with mortar constraints using nodal-attached geometry and a frictional pressure vector generated by an auxiliary kernel through a user-friendly action. Results are diffed against non-action output'
  rel_err = 1.0e-5
  abs_zero = 1.0e-5
  max_parallel = 8
[]

[frictional-mortar-3d-interp-geometry]
  type = 'Exodiff'
  input = 'frictional-mortar-3d-interp-geometry.i'
  min_ad_size = 100
  exodiff = 'frictional-mortar-3d-interp-geometry_out.e'
  requirement = 'The system shall solve a 3D frictional bouncing block problem with mortar constraint using weakly interpolated mesh geometry, such as normal and tangent vectors'
  rel_err = 1.0e-5
  abs_zero = 1.0e-5
  max_parallel = 8
[]

[./contact_initial_adaptivity]
  type = 'Exodiff'
  input = 'contact_initial_adaptivity.i'
  exodiff = 'contact_initial_adaptivity_out.e'
  abs_zero = 1e-10
  superlu = true
  issues = '#6799'
  requirement = "Contact shall be enforced on new nodes created due to mesh refinement"
  design = 'Constraints/index.md'
  max_threads = 1 # See https://github.com/idaholab/moose/issues/18685
[../]

[grid_sequencing]
  type = Exodiff
  input = grid-sequencing.i
  exodiff = grid-sequencing_exo.e
  requirement = 'The system shall use grid sequencing in order to improve the performance of the nonlinear solve in a frictional contact problem'
  cli_args = 'Postprocessors/active=contact'
  ad_indexing_type = global
  # Fix below when I understand duplicate element centroids
  map = False
  partial = True
  mesh_mode = replicated # Fix when mortar segment mesh is distributed
[]

[ranfs-ping-pong]
  type = RunException
  expect_err = 'Solve failed and timestep already at or below dtmin, cannot continue!'
  expect_out = 'Time Step 12'
  input = ranfs-ping-pong.i
  requirement = 'The node-face discretization with a RANFS formulation for frictionless mechanical contact shall be susceptible to ping-ponging, specifically in this case to a secondary node oscillating back and forth between different primary faces'
  heavy = True
[]

[kinematic-ping-pong]
  type = RunException
  expect_err = 'Solve failed and timestep already at or below dtmin, cannot continue!'
  expect_out = 'Time Step 12'
  input = kinematic-ping-pong.i
  requirement = 'The node-face discretization with a kinematic formulation for frictionless mechanical contact shall be susceptible to ping-ponging, specifically in this case to a secondary node oscillating back and forth between different primary faces'
  heavy = True
[]

[mortar-no-ping-pong_weighted]
  type = 'Exodiff'
  input = 'mortar-no-ping-pong_weighted.i'
  exodiff = 'mortar-no-ping-pong_weighted_exo.e'
  requirement = 'A variational consistent mortar formulation with dual bases for frictionless mechanical contact shall not show any ping-ponging behavior'
  heavy = True
  ad_indexing_type = 'global'
  rel_err = 6e-5
  abs_zero = 5e-9
[]

[frictionless-nodal-min-lm-mortar-disp-coarse]
  type = 'Exodiff'
  input = 'frictionless-nodal-min-lm-mortar-disp.i'
  exodiff = 'frictionless-nodal-min-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictionless bouncing block problem with nodal constraint enforcement, mortar application of forces, and min NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
[]

[frictionless-nodal-fb-lm-mortar-disp-coarse]
  type = 'Exodiff'
  input = 'frictionless-nodal-fb-lm-mortar-disp.i'
  exodiff = 'frictionless-nodal-fb-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictionless bouncing block problem with nodal constraint enforcement, mortar application of forces, and fb NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
[]

[frictionless-nodal-min-lm-nodal-disp-coarse]
  type = 'Exodiff'
  input = 'frictionless-nodal-min-lm-nodal-disp.i'
  exodiff = 'frictionless-nodal-min-lm-nodal-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictionless bouncing block problem with nodal constraint enforcement, nodal application of forces, and min NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
[]

[frictionless-nodal-fb-lm-nodal-disp-coarse]
  type = 'Exodiff'
  input = 'frictionless-nodal-fb-lm-nodal-disp.i'
  exodiff = 'frictionless-nodal-fb-lm-nodal-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictionless bouncing block problem with nodal constraint enforcement, nodal application of forces, and fb NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
[]

[frictionless-mortar-min-lm-mortar-disp-coarse]
  allow_test_objects = true
  type = 'Exodiff'
  input = 'frictionless-mortar-min-lm-mortar-disp.i'
  exodiff = 'frictionless-mortar-min-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictionless bouncing block problem with mortar constraint enforcement, mortar application of forces, and min NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
[]

[frictionless-mortar-fb-lm-mortar-disp-coarse]
  type = 'Exodiff'
  input = 'frictionless-mortar-fb-lm-mortar-disp.i'
  exodiff = 'frictionless-mortar-fb-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictionless bouncing block problem with mortar constraint enforcement, mortar application of forces, and fb NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
[]

[frictional-mortar-min-lm-mortar-disp-coarse]
  type = 'Exodiff'
  input = 'frictional-mortar-min-lm-mortar-disp.i'
  exodiff = 'frictional-mortar-min-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictional bouncing block problem with mortar constraint enforcement, mortar application of forces, and min NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
  max_parallel = 1
  abs_zero = 1e-9
[]

[frictional-mortar-fb-lm-mortar-disp-coarse]
  type = 'Exodiff'
  input = 'frictional-mortar-fb-lm-mortar-disp.i'
  exodiff = 'frictional-mortar-fb-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictional bouncing block problem with mortar constraint enforcement, mortar application of forces, and fb NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
  allow_test_objects = true
[]

[frictional-nodal-min-lm-mortar-disp-coarse]
  type = 'Exodiff'
  input = 'frictional-nodal-min-lm-mortar-disp.i'
  exodiff = 'frictional-nodal-min-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictional bouncing block problem with nodal constraint enforcement, mortar application of forces, and min NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active='' Outputs/exodus/sync_only=true --redirect-output output.log"
  allow_warnings = True
  max_parallel = 1
  method = 'OPT'
  allow_test_objects = true
  abs_zero = 1e-8
[]

[frictional-nodal-min-normal-lm-mortar-fb-tangential-lm-mortar-disp-coarse]
  type = 'Exodiff'
  input = 'frictional-nodal-min-normal-lm-mortar-fb-tangential-lm-mortar-disp.i'
  exodiff = 'frictional-nodal-min-normal-lm-mortar-fb-tangential-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictional bouncing block problem with nodal constraint enforcment for the normal LM using min NCP, mortar constraint enforcement of the tangential LM with fb, and mortar application of forces'
  cli_args = "Executioner/end_time=50 Postprocessors/active=''"
[]

[ranfs-project-svd]
  type = 'Exodiff'
  input = 'frictionless-nodal-reduced-active-set.i'
  exodiff = 'ranfs-project-svd.e'
  absent_out = '[1-9]+[0-9]* of 16 singular values'
  expect_out = '0 of 16 singular values'
  detail = 'nonsingular both with bounds projection and'
  cli_args = 'Outputs/file_base=ranfs-project-svd -pc_type svd -pc_svd_monitor'
[]

[ranfs-dont-project-svd]
  type = 'Exodiff'
  input = 'frictionless-nodal-reduced-active-set.i'
  exodiff = 'ranfs-dont-project-svd.e'
  detail = 'without bounds projection and be'
  cli_args = 'Outputs/file_base=ranfs-dont-project-svd Executioner/line_search=none -pc_type svd -pc_svd_monitor'
[]

[ranfs-project-amg]
  type = 'Exodiff'
  input = 'frictionless-nodal-reduced-active-set.i'
  exodiff = 'ranfs-project-amg.e'
  detail = 'solvable with amg both with bounds projection'
  cli_args = 'Outputs/file_base=ranfs-project-amg'
  expect_out = '[bB][oO][oO][mM][eE][rR][aA][mM][gG]'
[]

[ranfs-dont-project-amg]
  type = 'Exodiff'
  input = 'frictionless-nodal-reduced-active-set.i'
  exodiff = 'ranfs-dont-project-amg.e'
  detail = 'and without bounds projection.'
  cli_args = 'Outputs/file_base=ranfs-dont-project-amg Executioner/line_search=none'
  expect_out = '[bB][oO][oO][mM][eE][rR][aA][mM][gG]'
[]

[ranfs-jacobian]
  type = 'PetscJacobianTester'
  input = 'frictionless-nodal-reduced-active-set.i'
  difference_tol = 1e-7
  cli_args = "Outputs/active='' Executioner/line_search=none"
  expect_out = '[bB][oO][oO][mM][eE][rR][aA][mM][gG]'
  detail = "The system's RANFS scheme shall have a perfect Jacobian for mechanical contact that only has one non-zero normal component"
[]

[full-bouncing-block]
  type = 'Exodiff'
  input = 'bouncing-block-ranfs.i'
  exodiff = 'bouncing-block-ranfs_exo.e'
  detail = 'The system shall be able to detect when a secondary node is ping-ponging back and forth between different primary faces and consequently tie the locations of the secondary and corresponding primary node using Lagrange Multipliers corresponding to equality constraints, e.g. more RANFS'
  cli_args = "Postprocessors/active=''"
  custom_cmp = 'bouncing-block-ranfs.cmp'
  heavy = true
[]

[bouncing-block]
  type = 'Exodiff'
  input = 'bouncing-block-ranfs.i'
  exodiff = 'bouncing-block-coarse.e'
  abs_zero = 1e-9
  detail = 'The system shall be able to solve a smaller model of the full ping-ponging problem'
  cli_args = "Postprocessors/active='' Mesh/file=long-bottom-block-no-lower-d-coarse.e Executioner/end_time=100 Outputs/file_base=bouncing-block-coarse"
[]

[physics]
  type = 'Exodiff'
  input = tied-nodes.i
  exodiff = 'tied-nodes_exo.e'
  absent_out = '[1-9]+[0-9]* of 16 singular values'
  expect_out = '0 of 16 singular values'
  detail = 'tie nodes together and'
  allow_test_objects = True
  cli_args = '-pc_type svd -pc_svd_monitor'
[]

[jac]
  type = 'PetscJacobianTester'
  input = tied-nodes.i
  cli_args = "Outputs/active=''"
  detail = 'have a perfect Jacobian'
  allow_test_objects = True
[]

[equal_order]
  type = 'Exodiff'
  input = frictionless-weighted-gap.i
  exodiff = 'frictionless-weighted-gap_out.e'
  detail = 'using equal, first order bases for displacements and the lagrange multiplier'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
  ad_indexing_type = 'global'
[]

[mixed_order]
  type = 'Exodiff'
  input = frictionless-weighted-gap-mixed-basis.i
  exodiff = 'frictionless-weighted-gap-mixed-basis_out.e'
  detail = 'using a second order basis for displacements and a first order basis for the lagrange multiplier'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
  ad_indexing_type = 'global'
  min_ad_size = '73'
[]

[equal_order_edge_dropping]
  type = 'Exodiff'
  input = frictionless-weighted-gap.i
  exodiff = 'frictionless-weighted-gap-edge-dropping_out.e'
  detail = 'using equal, first order bases for displacements and the lagrange multiplier with correct edge dropping'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Constraints/weighted_gap_lm/correct_edge_dropping=true Executioner/end_time=50 Outputs/exodus=true Outputs/file_base=frictionless-weighted-gap-edge-dropping_out Postprocessors/active=''"
  ad_indexing_type = 'global'
[]

[verbose]
  type = 'Exodiff'
  input = frictional-nodal-min-normal-lm-mortar-pdass-tangential-lm-mortar-disp.i
  exodiff = 'frictional-nodal-min-normal-lm-mortar-pdass-tangential-lm-mortar-disp_out.e'
  detail = 'using verbose input file'
  cli_args = "Postprocessors/active=''"
  ad_indexing_type = 'global'
  abs_zero = 5e-9
[]

[action]
  type = 'Exodiff'
  input = frictional-nodal-min-normal-lm-mortar-pdass-tangential-lm-mortar-action.i
  exodiff = 'frictional-nodal-min-normal-lm-mortar-pdass-tangential-lm-mortar-action_out.e'
  detail = 'using the contact action'
  cli_args = "Postprocessors/active=''"
  ad_indexing_type = 'global'
  abs_zero = 5e-9
[]

[split_secondary]
  type = 'Exodiff'
  input = 'variational-frictional.i'
  exodiff = 'variational-frictional_out.e'
  ad_indexing_type = 'global'
  cli_args = 'Outputs/hide=procid'
  requirement = 'The system shall be able to solve a frictional, variationally consistent, mortar mechanical contact problem in which the secondary side of the contact interface is split between processes when run in parallel.'
  method = "!dbg"
  valgrind = 'heavy'
[]

[mixed_inactive_nodes]
  type = 'Exodiff'
  input = 'mixed-weighted-gap-swapped.i'
  exodiff = 'mixed-weighted-gap-swapped_out.e'
  requirement = 'The system shall not attempt to zero Lagrange multipliers that do not exist on inactive nodes.'
  method = "!dbg"
  valgrind = 'heavy'
  ad_indexing_type = 'global'
  min_ad_size = '100'
[]

[./catch_release]
  type = 'Exodiff'
  input = 'catch_release.i'
  exodiff = 'catch_release_out.e'
  custom_cmp = 'catch_release.exodiff'
  # This test fails with a negative Jacobian for older versions of
  # PETSc, possibly due to the changes in 3061bbd5d. It had already
  # been restricted for many years to PETSc >= 3.1, most likely due
  # to numerical instability.
  petsc_version = '>=3.5.0'
  requirement = 'The contact system shall enforce three-dimensional block to block interaction using a penalty approach.'
[../]

[./contact_displacements]
  type = 'RunException'
  input = 'contact_displacements.i'
  expect_err = "Contact requires updated coordinates.  Use the"
[../]

[dm_contact]
  type = 'CSVDiff'
  input = 'dm_mechanical_contact.i'
  csvdiff = 'dual_mortar_frictionless_contact_out.csv'
  rel_err = 1e-4
  cli_args = 'Contact/leftright/use_dual=true Outputs/file_base=dual_mortar_frictionless_contact_out'
  requirement = 'The system shall converge and match the solution produced by standard mortar contact.'
  allow_warnings = true
  mesh_mode = REPLICATED
  ad_indexing_type = 'global'
[]

[std_contact]
  type = 'CSVDiff'
  input = 'dm_mechanical_contact.i'
  csvdiff = 'std_mortar_frictionless_contact_out.csv'
  rel_err = 1e-4
  cli_args = 'Contact/leftright/use_dual=false Outputs/file_base=std_mortar_frictionless_contact_out'
  requirement = 'The system shall converge and match the solution produced by dual mortar contact.'
  allow_warnings = true
  mesh_mode = REPLICATED
  ad_indexing_type = 'global'
[]

[dm_contact_precon]
  type = 'CSVDiff'
  input = 'dm_mechanical_contact_precon.i'
  csvdiff = 'dual_mortar_frictionless_contact_out.csv'
  rel_err = 1e-4
  cli_args = 'Outputs/file_base=dual_mortar_frictionless_contact_out Preconditioning/vcp/adaptive_condensation=true'
  requirement = 'The system shall converge and match the solution with the standard methods using variable condenstation with AMG.'
  allow_warnings = true
  mesh_mode = REPLICATED
  ad_indexing_type = 'global'
[]

[dm_contact_precon_no_adaptivity]
  type = 'CSVDiff'
  input = 'dm_mechanical_contact_precon.i'
  csvdiff = 'dual_mortar_frictionless_contact_out.csv'
  rel_err = 1e-4
  cli_args = 'Outputs/file_base=dual_mortar_frictionless_contact_out Preconditioning/vcp/adaptive_condensation=false'
  requirement = 'The system shall converge and match the solution with the standard methods using variable condenstation with AMG, by always condensing out the LMs.'
  allow_warnings = true
  mesh_mode = REPLICATED
  ad_indexing_type = 'global'
[]

[dm_contact_precon_full_dinv]
  type = 'CSVDiff'
  input = 'dm_mechanical_contact_precon.i'
  csvdiff = 'dual_mortar_frictionless_contact_out.csv'
  rel_err = 1e-4
  cli_args = 'Outputs/file_base=dual_mortar_frictionless_contact_out Preconditioning/vcp/is_lm_coupling_diagonal=false'
  requirement = 'The system shall converge and match the solution with the standard methods using variable condenstation with AMG, by using LU to solve for the LM variable (not assuming diagonal coupling with the primal variable).'
  allow_warnings = true
  mesh_mode = REPLICATED
  ad_indexing_type = 'global'
[]

[2block-3D]
  type = 'Exodiff'
  input = '2blocks3d.i'
  exodiff = '2blocks3d_out.e'
  max_parallel = 4
  requirement = 'The system shall allow for split preconditioning based on contact surfaces.'
[]

[./single_point_2d]
  type = 'Exodiff'
  input = 'single_point_2d.i'
  exodiff = 'single_point_2d_out.e'
  custom_cmp = 'single_point_2d.cmp'
  superlu = true
  requirement = 'The contact system shall enforce 2D single-point contact with significant accumulated slip.'
[../]

[./single_point_2d_predictor]
  type = 'Exodiff'
  input = 'single_point_2d_predictor.i'
  exodiff = 'single_point_2d_predictor_out.e'
  custom_cmp = 'single_point_2d.cmp'
  superlu = true
  requirement = 'The contact system shall enforce 2D single-point contact with significant accumulated slip. With predictor solver options.'
[../]

[./single_point_2d_tp]
  type = 'Exodiff'
  input = 'single_point_2d_tp.i'
  exodiff = 'single_point_2d_tp_out.e'
  custom_cmp = 'single_point_2d.cmp'
  superlu = true
  requirement = 'The contact system shall enforce 2D single-point contact with significant accumulated slip when formulation selected is tangential_penalty contact.'
[../]

[./sliding_elastic_blocks_2d]
  type = 'Exodiff'
  input = 'sliding_elastic_blocks_2d.i'
  exodiff = 'sliding_elastic_blocks_2d_out.e'
  custom_cmp = 'sliding_elastic_blocks_2d.cmp'
  superlu = true
  requirement = 'The contact system shall enforce 2D line contact between quads with significant accumulated slip.'
[../]

[./sliding_elastic_blocks_2d_tp]
  type = 'Exodiff'
  input = 'sliding_elastic_blocks_2d_tp.i'
  exodiff = 'sliding_elastic_blocks_2d_tp_out.e'
  custom_cmp = 'sliding_elastic_blocks_2d.cmp'
  superlu = true
  requirement = 'The contact system shall enforce 2D line contact between quads with significant accumulated slip, when formulation selected is tangential_penalty.'
[../]

[./glued_contact_mechanical_constraint]
  type = 'Exodiff'
  input = 'glued_contact_mechanical_constraint_test.i'
  exodiff = 'mechanical_constraint_out.e'
  abs_zero = 5e-05
  max_parallel = 1
  requirement = 'The contact system shall enforce a glued contact constraint that ties together two blocks that are separated by an initial gap when the come in contact with each other so that the blocks move together.'
[../]

[./hertz_contact_rz]
  input = 'hertz_contact_rz.i'
  type = 'RunApp'
  requirement = 'The Contact system shall simulate Hertz contact between sphere and plane as a 2D axisymmetric problem with Quad4 elements.'
[../]

[./hertz_contact_rz_quad8]
  input = 'hertz_contact_rz_quad8.i'
  type = 'RunApp'
  requirement = 'The Contact system shall simulate Hertz contact between sphere and plane as a 2D axisymmetric problem with Quad8 elements.'
[../]

[./hertz_contact]
  input = 'hertz_contact.i'
  type = 'RunApp'
  heavy = true
  requirement = 'The Contact system shall simulate Hertz contact between sphere and plane as a 3D problem with Hex8 elements.'
[../]

[./hertz_contact_hex27]
  input = 'hertz_contact_hex27.i'
  type = 'RunApp'
  heavy = true
  requirement = 'The Contact system shall simulate Hertz contact between sphere and plane as a 3D problem with Hex27 elements.'
[../]

[./slip]
  type = Exodiff
  input = 'incremental_slip.i'
  exodiff = 'incremental_slip_out.e'
  abs_zero = 5e-09
  superlu = true
[../]

[auto-scaling]
  type = Exodiff
  input = bouncing-block-kinematic.i
  exodiff = bouncing-block-kinematic_exo.e
  requirement = 'The system shall be able to apply automatic scaling in conjunection with kinematic contact constraint enforcement and show no penetration and exhibit good nonlinear convergence'
  cli_args = "Postprocessors/active=''"
[]

[no-scaling]
  type = Exodiff
  input = bouncing-block-kinematic.i
  exodiff = no-scaling.e
  requirement = 'The system shall yield the same physical results when solving a kinematic contact problem with and without automatic scaling'
  cli_args = "Postprocessors/active='' Executioner/automatic_scaling=false Outputs/exo/file_base=no-scaling"
[]

[./frictionless_kinematic]
  type = 'Exodiff'
  input = 'frictionless_kinematic.i'
  exodiff = 'frictionless_kinematic_out.e'
  max_parallel = 1 # -pc_type lu
  requirement = "The contact system shall enforce a frictionless mechanical contact condition
                   between two blocks with a combination of normal and tangential motion using
                   a kinematic enforcement with the Constraint system."
[../]

[./frictionless_kinematic_gap_offset]
  type = 'Exodiff'
  input = 'frictionless_kinematic_gap_offsets.i'
  exodiff = 'frictionless_kinematic_gap_offsets_out.e'
  max_parallel = 1 # -pc_type lu
  requirement = "The contact system shall enforce a frictionless mechanical contact condition
                   between two blocks with gap offsets on both primary and secondary blocks using
                   a kinematic enforcement with the DiracKernel system."
  issues = '#14089'
[../]

[./frictionless_penalty]
  type = 'Exodiff'
  input = 'frictionless_penalty.i'
  exodiff = 'frictionless_penalty_out.e'
  max_parallel = 1 # -pc_type lu
  requirement = "The contact system shall enforce a frictionless mechanical contact condition
                   between two blocks with a combination of normal and tangential motion using
                   a penalty enforcement with the Constraint system."
[../]

[./glued_kinematic]
  type = 'Exodiff'
  input = 'glued_kinematic.i'
  exodiff = 'glued_kinematic_out.e'
  max_parallel = 1 # -pc_type lu
  requirement = "The contact system shall enforce a glued mechanical contact condition
                   between two blocks with a combination of normal and tangential motion using
                   a kinematic enforcement with the Constraint system."
[../]

[./glued_penalty]
  type = 'Exodiff'
  input = 'glued_penalty.i'
  exodiff = 'glued_penalty_out.e'
  max_parallel = 1 # -pc_type lu
  requirement = "The contact system shall enforce a glued mechanical contact condition
                   between two blocks with a combination of normal and tangential motion using
                   a penalty enforcement with the Constraint system."
[../]

[block-dynamics-aux-vel]
  type = 'Exodiff'
  exodiff = 'block-dynamics-aux-vel_out.e'
  input = 'block-dynamics-aux-vel.i'
  rel_err = 5.0e-4
  abs_zero = 1.0e-5
  min_parallel = 1
  max_parallel = 4
  min_ad_size = '50'
  ad_indexing_type = 'global'
  requirement = 'Contact module shall compute nodal weighted gap distance velocity values via a mortar auxiliary kernel.'
[]

[block-dynamics-aux-wear]
  type = 'Exodiff'
  exodiff = 'block-dynamics-aux-wear_out.e'
  input = 'block-dynamics-aux-wear.i'
  rel_err = 5.0e-4
  abs_zero = 1.0e-5
  min_parallel = 1
  max_parallel = 4
  min_ad_size = '50'
  ad_indexing_type = 'global'
  valgrind = 'none'
  requirement = 'Contact module shall compute nodal wear depth values in accordance with Archard equation via a mortar auxiliary kernel.'

[]

[block-dynamics-aux-wear-vel]
  type = 'Exodiff'
  exodiff = 'block-dynamics-aux-wear-vel_out.e'
  input = 'block-dynamics-aux-wear-vel.i'
  rel_err = 1.0e-3
  abs_zero = 1.0e-4
  min_parallel = 1
  max_parallel = 4
  min_ad_size = '50'
  ad_indexing_type = 'global'
  heavy = true
  method = 'opt'
  valgrind = 'none'
  requirement = 'Contact module shall compute nodal wear depth values in accordance with Archard equation and the mortar gap velocity via mortar auxiliary kernels in the same input file.'

[]

[frictional-mortar-3d-status]
  type = 'Exodiff'
  exodiff = 'frictional-mortar-3d-status_out.e'
  input = 'frictional-mortar-3d-status.i'
  rel_err = 1.0e-3
  abs_zero = 1.0e-4
  min_ad_size = '100'
  ad_indexing_type = 'global'
  cli_args = 'AuxKernels/frictional_state/tangent_two=mortar_tangential_3d_lm'
  requirement = 'Contact module shall compute nodal frictional status of mortar surfaces for a simple problem in which nodes are in stick and slip states.'
[]

[frictional-mortar-3d-status-error]
  type = RunException
  input = 'frictional-mortar-3d-status.i'
  expect_err = 'tangent_two: MortarFrictionalStateAux auxiliary kernel requires a second tangent Lagrange multiplier for three-dimensional problems'
  min_parallel = 1
  max_parallel = 2
  min_ad_size = '100'
  ad_indexing_type = 'global'
  requirement = 'Contact module shall generate an informative error if the nodal frictional status auxiliary kernel is not provided a second frictional Lagrange multiplier for a three-dimensional problem.'
[]

[block-dynamics-aux-fretting-wear-test-quick]
  type = 'Exodiff'
  exodiff = 'block-dynamics-aux-fretting-wear-test-quick_out.e'
  input = 'block-dynamics-aux-fretting-wear-test.i'
  rel_err = 1.0e-3
  abs_zero = 1.0e-4
  min_parallel = 1
  max_parallel = 4
  min_ad_size = '50'
  ad_indexing_type = 'global'
  cli_args = 'Executioner/end_time=0.15 Outputs/file_base=block-dynamics-aux-fretting-wear-test-quick_out'
  requirement = 'Contact module shall compute nodal wear depth values in accordance with Archard equation and the mortar gap velocity via mortar auxiliary kernels while including these computations in the definition of mortar normal contact constraints in an asymmetric problem for a short simulation.'
[]

[block-dynamics]
  type = 'Exodiff'
  exodiff = 'block-dynamics_out.e'
  input = 'block-dynamics.i'
  rel_err = 5.0e-4
  abs_zero = 1.0e-5
  min_parallel = 1
  max_parallel = 4
  min_ad_size = '50'
  ad_indexing_type = 'global'
  cli_args = 'Executioner/dt=0.05 Executioner/end_time=0.3'
  requirement = 'The system shall solve mortar frictionless contact between two blocks with weighted gap time stabilization using mortar nodal geometry.'
[]

[block-dynamics-action]
  type = 'Exodiff'
  exodiff = 'block-dynamics-action_out.e'
  input = 'block-dynamics-action.i'
  rel_err = 5.0e-4
  abs_zero = 1.0e-5
  min_parallel = 1
  max_parallel = 4
  min_ad_size = '50'
  ad_indexing_type = 'global'
  cli_args = 'Contact/mechanical/mortar_dynamics=true Contact/mechanical/newmark_beta=0.25  Contact/mechanical/newmark_gamma=0.5 Executioner/dt=0.05 Executioner/end_time=0.3'
  requirement = 'The system shall solve mortar frictionless contact between two blocks with weighted gap time stabilization using mortar nodal geometry via the contact mechanics action.'
[]

[block-dynamics-action-beta-error]
  type = RunException
  input = 'block-dynamics-action.i'
  expect_err = 'newmark_beta can only be used with the mortar_dynamics option'
  cli_args = 'Contact/mechanical/mortar_dynamics=false Contact/mechanical/newmark_beta=0.25 Executioner/end_time=0.3'
  min_ad_size = '50'
  ad_indexing_type = 'global'
  requirement = 'The system shall generate an error if dynamics is not specifically requested and Newmark-beta integration parameter beta is provided.'
[]

[block-dynamics-action-gamma-error]
  type = RunException
  input = 'block-dynamics-action.i'
  expect_err = 'newmark_gamma can only be used with the mortar_dynamics option'
  cli_args = 'Contact/mechanical/mortar_dynamics=false Contact/mechanical/newmark_gamma=0.5 Executioner/end_time=0.3'
  min_ad_size = '50'
  ad_indexing_type = 'global'
  requirement = 'The system shall generate an error if dynamics is not specifically requested and Newmark-beta integration parameter gamma is provided.'

[]

[block-dynamics-friction]
  type = 'Exodiff'
  exodiff = 'block-dynamics-friction_out.e'
  input = 'block-dynamics-friction.i'
  rel_err = 5.0e-4
  abs_zero = 1.0e-5
  min_parallel = 1
  max_parallel = 4
  min_ad_size = '50'
  ad_indexing_type = 'global'
  cli_args = 'Executioner/dt=0.05 Executioner/end_time=0.3'
  requirement = 'The system shall solve mortar frictional contact between two blocks with weighted gap time stabilization using mortar nodal geometry.'
[]

[block-dynamics-friction-action]
  type = 'Exodiff'
  exodiff = 'block-dynamics-friction-action_out.e'
  input = 'block-dynamics-friction-action.i'
  rel_err = 5.0e-4
  abs_zero = 1.0e-5
  min_ad_size = '50'
  ad_indexing_type = 'global'
  cli_args = 'Executioner/dt=0.05 Executioner/end_time=0.3'
  requirement = 'The system shall solve mortar frictional contact between two blocks with weighted gap time stabilization using mortar nodal geometry via the contact action.'
[]

[block-dynamics-friction-creep]
  type = 'Exodiff'
  exodiff = 'block-dynamics-friction-creep_out.e'
  input = 'block-dynamics-friction-creep.i'
  rel_err = 5.0e-4
  abs_zero = 1.0e-5
  min_parallel = 1
  min_ad_size = '50'
  ad_indexing_type = 'global'
  requirement = 'The system shall simulate mortar frictional contact between two blocks with weighted gap time stabilization using mortar nodal geometry with a creep material model.'
[]

[frictional-mortar-3d-dynamics]
  type = 'CSVDiff'
  input = 'frictional-mortar-3d-dynamics.i'
  csvdiff = 'frictional-mortar-3d-dynamics_out_frictional-pressure_0001.csv frictional-mortar-3d-dynamics_out_frictional-pressure-3d_0001.csv frictional-mortar-3d-dynamics_out_tangent_x_0001.csv frictional-mortar-3d-dynamics_out_tangent_y_0001.csv'
  requirement = 'The system shall solve a dynamic 3D frictional bouncing block problem with mortar constraint using nodal-attached geometry'
  rel_err = 1.0e-5
  abs_zero = 1.0e-5
  min_ad_size = 100
  max_parallel = 8
  method = 'opt'
[]

[frictional-mortar-3d-dynamics-light]
  type = 'Exodiff'
  exodiff = 'frictional-mortar-3d-dynamics-light_out.e'
  input = 'frictional-mortar-3d-dynamics-light.i'
  rel_err = 5.0e-4
  abs_zero = 1.0e-5
  min_ad_size = 100
  ad_indexing_type = 'global'
  requirement = 'The system shall solve a dynamic 3D frictional one-element bouncing block problem with mortar constraint using nodal-attached geometry and the correct edge dropping setting'
[]

[frictional-mortar-3d-dynamics-light-miss-dropping]
  type = 'Exodiff'
  exodiff = 'frictional-mortar-3d-dynamics-light_out.e'
  input = 'frictional-mortar-3d-dynamics-light.i'
  rel_err = 5.0e-4
  abs_zero = 1.0e-5
  min_ad_size = 100
  ad_indexing_type = 'global'
  cli_args = 'Constraints/friction/correct_edge_dropping=false Constraints/normal_x/correct_edge_dropping=false Constraints/normal_y/correct_edge_dropping=false Constraints/normal_z/correct_edge_dropping=false Constraints/tangential_x/correct_edge_dropping=false Constraints/tangential_y/correct_edge_dropping=false Constraints/tangential_z/correct_edge_dropping=false Constraints/tangential_dir_x/correct_edge_dropping=false Constraints/tangential_dir_y/correct_edge_dropping=false Constraints/tangential_dir_z/correct_edge_dropping=false'
  requirement = 'The system shall solve a dynamic 3D frictional one-element bouncing block problem with mortar constraint using nodal-attached geometry and the incorrect edge dropping setting'
[]

[frictional-mortar-3d-dynamics-light-function]
  type = 'Exodiff'
  exodiff = 'frictional-mortar-3d-dynamics-light-function_out.e'
  input = 'frictional-mortar-3d-dynamics-light-function.i'
  rel_err = 5.0e-4
  abs_zero = 1.0e-5
  min_ad_size = 100
  ad_indexing_type = 'global'
  max_parallel = '2'
  requirement = 'The system shall solve a dynamic 3D frictional one-element bouncing block problem with mortar constraint using nodal-attached geometry and a friction coefficient that depends on the normal contact pressure and the relative tangential velocity'
[]

[frictional_bouncing_block_action_restart_1]
  type = 'Exodiff'
  input = 'frictional_bouncing_block_action_restart_1.i'
  exodiff = 'frictional_bouncing_block_action_restart_1_out.e'
  abs_zero = 1.0e-4
  rel_err = 1.0e-4
  ad_indexing_type = 'global'
  requirement = 'The system shall be able to run a two-dimensional frictional model using the contact action for mortar applications'
[]

[frictional_bouncing_block_action_restart_2]
  type = 'Exodiff'
  input = 'frictional_bouncing_block_action_restart_2.i'
  exodiff = 'frictional_bouncing_block_action_restart_2_out.e'
  recover = false
  abs_zero = 1.0e-4
  rel_err = 1.0e-4
  ad_indexing_type = 'global'
  requirement = 'The system shall be able to restart a mortar mechanical contact simulation via the action without generating additional lower dimensional subdomains which may be unused'
[]

[small]
  detail = 'using the small strain formulation.'
  type = 'CSVDiff'
  input = 'small.i'
  csvdiff = 'small.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite]
  detail = 'using the finite strain formulation.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'finite.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite_scaling]
  detail = 'using the finite strain formulation with automatic scaling.'
  type = 'CSVDiff'
  input = 'finite.i'
  cli_args = 'Executioner/automatic_scaling=true Executioner/compute_scaling_once=false name=finite_scaling'
  csvdiff = 'finite_scaling.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite_rr]
  detail = 'using the finite strain formulation and reference residual.'
  type = 'CSVDiff'
  input = 'finite_rr.i'
  csvdiff = 'finite_rr.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[small-jac]
  type = 'PetscJacobianTester'
  prereq = soft/small
  input = 'small.i'
  ratio_tol = 3e-3
  difference_tol = 1
  cli_args = 'Executioner/num_steps=3 -snes_test_err 1e-9'
  run_sim = true
  heavy = true
  detail = 'using the small strain formulation and calculate a perfect Jacobian.'
  ad_indexing_type = 'global'
[]

[finite-jac]
  type = 'PetscJacobianTester'
  prereq = soft/finite
  input = 'finite.i'
  ratio_tol = 3e-3
  difference_tol = 1
  cli_args = 'Executioner/num_steps=3 -snes_test_err 1e-9'
  run_sim = true
  heavy = true
  detail = 'using the finite strain formulation and calculate a perfect Jacobian.'
  ad_indexing_type = 'global'
[]

[stiff_stiff]
  detail = 'using with a stiff block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'stiff_stiff.csv'
  cli_args = 'E_block=1e9 E_plank=1e9 name=stiff_stiff'
  heavy = true
  ad_indexing_type = 'global'
[]

[soft_stiff]
  detail = 'using with a soft block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'soft_stiff.csv'
  cli_args = 'E_block=1e7 E_plank=1e9 name=soft_stiff'
  heavy = true
  ad_indexing_type = 'global'
[]

[small]
  detail = 'using the small strain formulation.'
  type = 'CSVDiff'
  input = 'small.i'
  csvdiff = 'small.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite]
  detail = 'using the finite strain formulation.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'finite.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite_scaling]
  detail = 'using the finite strain formulation with automatic scaling.'
  type = 'CSVDiff'
  input = 'finite.i'
  cli_args = 'Executioner/automatic_scaling=true Executioner/compute_scaling_once=false name=finite_scaling'
  csvdiff = 'finite_scaling.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite_rr]
  detail = 'using the finite strain formulation and reference residual.'
  type = 'CSVDiff'
  input = 'finite_rr.i'
  csvdiff = 'finite_rr.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[small-jac]
  type = 'PetscJacobianTester'
  prereq = soft/small
  input = 'small.i'
  ratio_tol = 1e-3
  difference_tol = 1
  cli_args = 'Executioner/num_steps=3 -snes_test_err 1e-9'
  run_sim = true
  heavy = true
  detail = 'using the small strain formulation and calculate a perfect Jacobian.'
  ad_indexing_type = 'global'
[]

[finite-jac]
  type = 'PetscJacobianTester'
  prereq = soft/finite
  input = 'finite.i'
  ratio_tol = 1e-3
  difference_tol = 1
  cli_args = 'Executioner/num_steps=2 -snes_test_err 1e-9'
  run_sim = true
  heavy = true
  detail = 'using the finite strain formulation and calculate a perfect Jacobian.'
  ad_indexing_type = 'global'
[]

[stiff_stiff]
  detail = 'using with a stiff block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'stiff_stiff.csv'
  cli_args = 'E_block=1e9 E_plank=1e9 name=stiff_stiff'
  heavy = true
  ad_indexing_type = 'global'
[]

[soft_stiff]
  detail = 'using with a soft block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'soft_stiff.csv'
  cli_args = 'E_block=1e7 E_plank=1e9 name=soft_stiff'
  heavy = true
  ad_indexing_type = 'global'
[]

[small]
  detail = 'using the small strain formulation.'
  type = 'CSVDiff'
  input = 'small.i'
  csvdiff = 'small.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite]
  detail = 'using the finite strain formulation.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'finite.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite_scaling]
  detail = 'using the finite strain formulation with automatic scaling.'
  type = 'CSVDiff'
  input = 'finite.i'
  cli_args = 'Executioner/automatic_scaling=true Executioner/compute_scaling_once=false name=finite_scaling'
  csvdiff = 'finite_scaling.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite_rr]
  detail = 'using the finite strain formulation and reference residual.'
  type = 'CSVDiff'
  input = 'finite_rr.i'
  csvdiff = 'finite_rr.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[stiff_stiff]
  detail = 'using with a stiff block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'stiff_stiff.csv'
  cli_args = 'E_block=1e9 E_plank=1e9 name=stiff_stiff'
  heavy = true
  ad_indexing_type = 'global'
[]

[soft_stiff]
  detail = 'using with a soft block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'soft_stiff.csv'
  cli_args = 'E_block=1e7 E_plank=1e9 name=soft_stiff'
  heavy = true
  ad_indexing_type = 'global'
[]

[small]
  detail = 'using the small strain formulation.'
  type = 'CSVDiff'
  input = 'small.i'
  csvdiff = 'small.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite]
  detail = 'using the finite strain formulation.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'finite.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite_scaling]
  detail = 'using the finite strain formulation with automatic scaling.'
  type = 'CSVDiff'
  input = 'finite.i'
  cli_args = 'Executioner/automatic_scaling=true Executioner/compute_scaling_once=false name=finite_scaling'
  csvdiff = 'finite_scaling.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite_rr]
  detail = 'using the finite strain formulation and reference residual.'
  type = 'CSVDiff'
  input = 'finite_rr.i'
  csvdiff = 'finite_rr.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[stiff_stiff]
  detail = 'using with a stiff block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'stiff_stiff.csv'
  cli_args = 'E_block=1e9 E_plank=1e9 name=stiff_stiff'
  heavy = true
  ad_indexing_type = 'global'
[]

[soft_stiff]
  detail = 'using with a soft block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'soft_stiff.csv'
  cli_args = 'E_block=1e7 E_plank=1e9 name=soft_stiff'
  heavy = true
  ad_indexing_type = 'global'
[]

[small]
  detail = 'using the small strain formulation.'
  type = 'CSVDiff'
  input = 'small.i'
  csvdiff = 'small.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite]
  detail = 'using the finite strain formulation.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'finite.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite_scaling]
  detail = 'using the finite strain formulation with automatic scaling.'
  type = 'CSVDiff'
  input = 'finite.i'
  cli_args = 'Executioner/automatic_scaling=true Executioner/compute_scaling_once=false name=finite_scaling'
  csvdiff = 'finite_scaling.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite_rr]
  detail = 'using the finite strain formulation and reference residual.'
  type = 'CSVDiff'
  input = 'finite_rr.i'
  csvdiff = 'finite_rr.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[small-jac]
  type = 'PetscJacobianTester'
  prereq = soft/small
  input = 'small.i'
  ratio_tol = 4e-3
  difference_tol = 1
  cli_args = 'Executioner/num_steps=3 -snes_test_err 1e-9'
  run_sim = true
  heavy = true
  ad_indexing_type = 'global'
  detail = 'using the small strain formulation and calculate a perfect Jacobian.'
[]

[finite-jac]
  type = 'PetscJacobianTester'
  prereq = soft/finite
  input = 'finite.i'
  ratio_tol = 5e-3
  difference_tol = 1
  cli_args = 'Executioner/num_steps=3 -snes_test_err 1e-9'
  run_sim = true
  heavy = true
  ad_indexing_type = 'global'
  detail = 'using the finite strain formulation and calculate a perfect Jacobian.'
[]

[stiff_stiff]
  detail = 'using with a stiff block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'stiff_stiff.csv'
  cli_args = 'E_block=1e9 E_plank=1e9 name=stiff_stiff'
  heavy = true
  ad_indexing_type = 'global'
[]

[soft_stiff]
  detail = 'using with a soft block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'soft_stiff.csv'
  cli_args = 'E_block=1e7 E_plank=1e9 name=soft_stiff'
  heavy = true
  ad_indexing_type = 'global'
[]

[small]
  detail = 'using the small strain formulation.'
  type = 'CSVDiff'
  input = 'small.i'
  csvdiff = 'small.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite]
  detail = 'using the finite strain formulation.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'finite.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite_scaling]
  detail = 'using the finite strain formulation with automatic scaling.'
  type = 'CSVDiff'
  input = 'finite.i'
  cli_args = 'Executioner/automatic_scaling=true Executioner/compute_scaling_once=false name=finite_scaling'
  csvdiff = 'finite_scaling.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite_rr]
  detail = 'using the finite strain formulation and reference residual.'
  type = 'CSVDiff'
  input = 'finite_rr.i'
  csvdiff = 'finite_rr.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[small-jac]
  type = 'PetscJacobianTester'
  prereq = soft/small
  input = 'small.i'
  ratio_tol = 2e-3
  difference_tol = 1
  cli_args = 'Executioner/num_steps=3 -snes_test_err 1e-9'
  run_sim = true
  heavy = true
  detail = 'using the small strain formulation and calculate a perfect Jacobian.'
  ad_indexing_type = 'global'
[]

[finite-jac]
  type = 'PetscJacobianTester'
  prereq = soft/finite
  input = 'finite.i'
  ratio_tol = 2e-3
  difference_tol = 1
  cli_args = 'Executioner/num_steps=3 -snes_test_err 1e-9'
  run_sim = true
  heavy = true
  detail = 'using the finite strain formulation and calculate a perfect Jacobian.'
  ad_indexing_type = 'global'
[]

[stiff_stiff]
  detail = 'using with a stiff block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'stiff_stiff.csv'
  cli_args = 'E_block=1e9 E_plank=1e9 name=stiff_stiff'
  heavy = true
  ad_indexing_type = 'global'
[]

[soft_stiff]
  detail = 'using with a soft block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'soft_stiff.csv'
  cli_args = 'E_block=1e7 E_plank=1e9 name=soft_stiff'
  heavy = true
  ad_indexing_type = 'global'
[]

[small]
  detail = 'using the small strain formulation.'
  type = 'CSVDiff'
  input = 'small.i'
  csvdiff = 'small.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite]
  detail = 'using the finite strain formulation.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'finite.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite_scaling]
  detail = 'using the finite strain formulation with automatic scaling.'
  type = 'CSVDiff'
  input = 'finite.i'
  cli_args = 'Executioner/automatic_scaling=true Executioner/compute_scaling_once=false name=finite_scaling'
  csvdiff = 'finite_scaling.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite_rr]
  detail = 'using the finite strain formulation and reference residual.'
  type = 'CSVDiff'
  input = 'finite_rr.i'
  csvdiff = 'finite_rr.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[stiff_stiff]
  detail = 'using with a stiff block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'stiff_stiff.csv'
  cli_args = 'E_block=1e9 E_plank=1e9 name=stiff_stiff'
  heavy = true
  ad_indexing_type = 'global'
[]

[soft_stiff]
  detail = 'using with a soft block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'soft_stiff.csv'
  cli_args = 'E_block=1e7 E_plank=1e9 name=soft_stiff'
  heavy = true
  ad_indexing_type = 'global'
[]

[small]
  detail = 'using the small strain formulation.'
  type = 'CSVDiff'
  input = 'small.i'
  csvdiff = 'small.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite]
  detail = 'using the finite strain formulation.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'finite.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite_scaling]
  detail = 'using the finite strain formulation with automatic scaling.'
  type = 'CSVDiff'
  input = 'finite.i'
  cli_args = 'Executioner/automatic_scaling=true Executioner/compute_scaling_once=false name=finite_scaling'
  csvdiff = 'finite_scaling.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[finite_rr]
  detail = 'using the finite strain formulation and reference residual.'
  type = 'CSVDiff'
  input = 'finite_rr.i'
  csvdiff = 'finite_rr.csv'
  heavy = true
  ad_indexing_type = 'global'
[]

[stiff_stiff]
  detail = 'using with a stiff block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'stiff_stiff.csv'
  cli_args = 'E_block=1e9 E_plank=1e9 name=stiff_stiff'
  heavy = true
  ad_indexing_type = 'global'
[]

[soft_stiff]
  detail = 'using with a soft block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'soft_stiff.csv'
  cli_args = 'E_block=1e7 E_plank=1e9 name=soft_stiff'
  heavy = true
  ad_indexing_type = 'global'
[]

[./horizontal_blocks_mortar_TM]
  type = 'RunApp'
  input = horizontal_blocks_mortar_TM.i
  ad_indexing_type = 'global'
[../]

[multiple_pairs]
  type = Exodiff
  input = 'multiple_pairs.i'
  exodiff = 'multiple_pairs_out.e'
  abs_zero = 1.0e-8
  rel_err = 1.0e-5
  cli_args = "Postprocessors/active=''"
  requirement = 'Contact module action shall allow for multiple contact pairs when selecting node-face formulation.'
[]

[./3D]
  type = 'Exodiff'
  input = 'nodal_area_3D.i'
  exodiff = 'nodal_area_3D_out.e'
  requirement = 'The system shall compute the nodal area for use with contact calculations in 3D.'
[../]

[./3D2Proc]
  type = 'Exodiff'
  input = 'nodal_area_3D.i'
  exodiff = 'nodal_area_3D_out.e'
  prereq = 3D
  min_parallel = 2
  requirement = 'The system shall compute the nodal area in parallel for use with contact calculations in 3D.'
[../]

[./2D]
  type = 'Exodiff'
  input = 'nodal_area_2D.i'
  exodiff = 'nodal_area_2D_out.e'
  requirement = 'The system shall compute the nodal area for use with contact calculations in 2D.'
[../]

[./2D2Proc]
  type = 'Exodiff'
  input = 'nodal_area_2D.i'
  exodiff = 'nodal_area_2D_out.e'
  prereq = 2D
  min_parallel = 2
  requirement = 'The system shall compute the nodal area in parallel for use with contact calculations in 2D.'
[../]

[./Hex20]
  type = 'Exodiff'
  input = 'nodal_area_Hex20.i'
  exodiff = 'nodal_area_Hex20_out.e'
  custom_cmp = 'nodal_area_Hex20.exodiff'
  requirement = 'The system shall compute the nodal area for Hex20 elements for use with contact calculations.'
[../]

[./Hex20_2]
  type = 'Exodiff'
  input = 'nodal_area_Hex20.i'
  cli_args = 'Contact/dummy_name/model=frictionless Outputs/file_base=nodal_area_Hex20_out'
  exodiff = 'nodal_area_Hex20_out.e'
  custom_cmp = 'nodal_area_Hex20.exodiff'
  prereq = Hex20
  requirement = 'The system shall compute the nodal area for Hex20 elements for use with frictionless contact calculations.'
[../]

[./Hex20_3]
  type = 'Exodiff'
  input = 'nodal_area_Hex20_3.i'
  exodiff = 'nodal_area_Hex20_3_out.e'
  custom_cmp = 'nodal_area_Hex20_3.exodiff'
  # This test fails with a negative Jacobian for older versions of
  # PETSc, possibly due to the changes in 3061bbd5d.
  petsc_version = '>=3.5.0'
  requirement = 'The system shall compute the nodal area for Hex20 elements for use with penalty contact calculations.'
[../]

[./Hex202Proc]
  type = 'Exodiff'
  input = 'nodal_area_Hex20.i'
  exodiff = 'nodal_area_Hex20_out.e'
  custom_cmp = 'nodal_area_Hex20.exodiff'
  prereq = Hex20_2
  min_parallel = 2
  requirement = 'The system shall compute the nodal area in parallel for Hex20 elements for use with contact calculations.'
[../]

[./Hex27]
  type = 'Exodiff'
  input = 'nodal_area_Hex27.i'
  exodiff = 'nodal_area_Hex27_out.e'
  custom_cmp = 'nodal_area_Hex27.exodiff'
  requirement = 'The system shall compute the nodal area for Hex27 elements for use with contact calculations.'
[../]

[./Hex272Proc]
  type = 'Exodiff'
  input = 'nodal_area_Hex27.i'
  exodiff = 'nodal_area_Hex27_out.e'
  custom_cmp = 'nodal_area_Hex27.exodiff'
  prereq = Hex27
  min_parallel = 2
  requirement = 'The system shall compute the nodal area in parallel for Hex27 elements for use with contact calculations.'
[../]

[nonsingular]
  type = Exodiff
  input = frictional-mortar.i
  exodiff = frictional-mortar_exodus.e
  design = 'ComputeFrictionalForceLMMechanicalContact.md'
  issues = '#17495'
  requirement = 'The system shall not generate singular Jacobians in frictional mortar contact.'
  expect_out = '0 of 42 singular values'
  absent_out = '[1-9]+[0-9]* of 42 singular values'
  ad_indexing_type = 'global'
  mesh_mode = 'REPLICATED'
[]

[exo]
  type = Exodiff
  input = normal-nodal-lm-tan-tolerance.i
  exodiff = normal-nodal-lm-tan-tolerance_out.e
  requirement = 'The system shall provide an option for extending primary surfaces when doing mortar contact to ensure that secondary nodes which project near the edge of the primary surface are registered with the near-surface.'
  method = '!DBG' # In debug mode, you get a ton of warnings out of libMesh when we Assembly::reinitNeighborAtPhysical for the secondary node that doesn't actually project onto the primary face; but that's the point of fudging with the tangential_tolerance parameter
  allow_warnings = true # Warning about slowness of standard finite differencing
  rel_err = 2e-3 # Even with nl_rel_tol of 1e-12 significant differences in parallel
  mesh_mode = REPLICATED
  valgrind = HEAVY
[]

[./normalized_penalty]
  type = Exodiff
  input = 'normalized_penalty.i'
  exodiff = 'normalized_penalty_out.e'
  abs_zero = 1e-7
  max_parallel = 1 # -pc_type lu
  requirement = 'The contact system shall yield repeatable results for 2D contact with elements of various aspect ratios. Penalty contact.'
[../]

[./normalized_penalty_Q8]
  type = Exodiff
  input = 'normalized_penalty_Q8.i'
  exodiff = 'normalized_penalty_Q8_out.e'
  abs_zero = 1e-7
  max_parallel = 1 # -pc_type lu
  requirement = 'The contact system shall yield repeatable results for 2D contact with Q8 elements of various aspect ratios. Penalty contact.'
[../]

[./normalized_penalty_kin]
  type = Exodiff
  input = 'normalized_penalty_kin.i'
  exodiff = 'normalized_penalty_kin_out.e'
  abs_zero = 1e-7
  max_parallel = 1 # -pc_type lu
  # This test no longer works on older versions of PETSc due to the
  # changes in 3061bbd5d.
  petsc_version = '>=3.5.0'
  requirement = 'The contact system shall yield repeatable results for 2D contact with elements of various aspect ratios. Kinematic contact.'
[../]

[./normalized_penalty_kin_Q8]
  type = Exodiff
  input = 'normalized_penalty_kin_Q8.i'
  exodiff = 'normalized_penalty_kin_Q8_out.e'
  abs_zero = 1e-7
  max_parallel = 1 # -pc_type lu
  # This test no longer works on older versions of PETSc due to the
  # changes in 3061bbd5d.
  petsc_version = '>=3.5.0'
  requirement = 'The contact system shall yield repeatable results for 2D contact with Q8 elements of various aspect ratios. Kinematic contact.'
[../]

[cylinder_friction]
  type = 'CSVDiff'
  csvdiff = 'cylinder_friction_check_friction_FINAL.csv cylinder_friction_check_x_disp_FINAL.csv cylinder_friction_check_y_disp_FINAL.csv cylinder_friction_check_cont_press_FINAL.csv'
  input = 'cylinder_friction.i'
  rel_err = 5.0e-4
  min_parallel = 4
  max_parallel = 4
  heavy = true
  min_ad_size = '50'
  ad_indexing_type = 'global'
  requirement = 'Contact module shall solve frictional contact between a cylinder and a plane.'
[]

[ironing]
  type = 'CSVDiff'
  csvdiff = 'ironing_chkfile_cont_press_0009.csv ironing_chkfile_friction_0009.csv'
  input = 'ironing.i'
  rel_err = 1.0e-4
  min_parallel = 4
  heavy = true
  min_ad_size = '200'
  ad_indexing_type = 'global'
  requirement = 'Contact module shall solve frictional contact between a semicircular tool and flexible base material.'
[]

[frictional_bouncing_block]
  type = 'CSVDiff'
  csvdiff = 'frictional_bouncing_block_checkfile_cont_press_0029.csv frictional_bouncing_block_checkfile_friction_0029.csv'
  input = 'frictional_bouncing_block.i'
  rel_err = 1.0e-4
  min_ad_size = '50'
  heavy = true
  ad_indexing_type = 'global'
  requirement = 'Contact module shall solve frictional contact between a bouncing block and flexible base material.'
[]

[frictional_bouncing_block_action]
  type = 'CSVDiff'
  csvdiff = 'frictional_bouncing_block_action_checkfile_cont_press_0029.csv frictional_bouncing_block_action_checkfile_friction_0029.csv'
  input = 'frictional_bouncing_block_action.i'
  rel_err = 1.0e-4
  min_ad_size = '50'
  heavy = true
  ad_indexing_type = 'global'
  requirement = 'Contact module shall solve frictional contact between a bouncing block and flexible base material verifying setup in the contact action.'
[]

[frictional_bouncing_block_edge_dropping]
  type = 'CSVDiff'
  csvdiff = 'frictional_bouncing_block_edge_dropping_cont_press_0023.csv frictional_bouncing_block_edge_dropping_friction_0023.csv'
  input = 'frictional_bouncing_block.i'
  rel_err = 1.0e-4
  min_ad_size = '50'
  method = '!dbg'
  valgrind = 'none'
  ad_indexing_type = 'global'
  cli_args = "Constraints/weighted_gap_lm/correct_edge_dropping=true Executioner/end_time=5.5 Outputs/file_base=frictional_bouncing_block_edge_dropping"
  requirement = 'Contact module shall solve frictional contact between a bouncing block and a flexible substrate when correct edge dropping is enabled. Additional requirement is that the correct edge dropping treatment must yield same results as incorrect edge dropping treatment when there is not actual edge dropping, e.g.: All secondary surface projects to the primary surface.'
[]

[./pressureAugLag]
  type = 'Exodiff'
  input = 'pressureAugLag.i'
  exodiff = 'pressureAugLag_out.e'
  custom_cmp = 'pressure.exodiff'
  max_parallel = 1
  requirement = 'The contact system shall reproduce contact pressure results among various formulation types. Augmented Lagrangian formulation.'
[../]

[./pressurePenalty]
  type = 'Exodiff'
  input = 'pressurePenalty.i'
  exodiff = 'pressurePenalty_out.e'
  custom_cmp = 'pressure.exodiff'
  max_parallel = 1
  requirement = 'The contact system shall reproduce contact pressure results among various formulation types. Penalty.'
[../]

[./pressurePenalty_mechanical_constraint]
  type = 'Exodiff'
  input = 'pressurePenalty_mechanical_constraint.i'
  exodiff = 'pressurePenalty_mechanical_constraint_out.e'
  custom_cmp = 'pressure.exodiff'
  max_parallel = 1
  requirement = 'The contact system shall reproduce contact pressure results among various formulation types. Mechanical constraint.'
[../]

[auto-scaling]
  type = Exodiff
  input = bouncing-block-ranfs.i
  exodiff = bouncing-block-ranfs_exo.e
  requirement = 'The system shall be able to apply automatic scaling in conjunction with ranfs contact'
  cli_args = "Postprocessors/active=''"
  abs_zero = 1e-8
[]

[no-scaling]
  type = Exodiff
  input = bouncing-block-ranfs.i
  exodiff = no-scaling.e
  requirement = 'The system shall be able to solve ranfs contact with no scaling'
  cli_args = "Executioner/automatic_scaling=false Outputs/file_base=no-scaling Postprocessors/active=''"
  abs_zero = 1e-8
[]

[./ring_contact]
  type = 'Exodiff'
  input = 'ring_contact.i'
  exodiff = 'ring_contact_out.e'
  abs_zero = 5e-07
  use_old_floor = True
  max_parallel = 1
  requirement = 'The contact system shall enforce contact between three-dimensional non-conformal surfaces with Hex20 elements.'
[../]

[./simple_contact_test]
  type = Exodiff
  input = 'simple_contact_test.i'
  exodiff = 'simple_contact_test_out.e'
[../]

[./simple_contact_test2]
  type = Exodiff
  input = 'simple_contact_test2.i'
  exodiff = 'simple_contact_test2_out.e'
[../]

[./simple_contact_rz_test]
  type = Exodiff
  input = 'simple_contact_rz_test.i'
  exodiff = 'simple_contact_rz_test_out.e'
  abs_zero = 1e-7
  scale_refine = 1
  max_parallel = 1
[../]

[./simple_contact_rspherical]
  type = Exodiff
  input = 'simple_contact_rspherical.i'
  exodiff = 'simple_contact_rspherical_out.e'
  max_parallel = 1
[../]

[frictionless-compress-2d]
  type = 'CSVDiff'
  input = 'two_equal_blocks_compress_2d.i'
  csvdiff = 'two_equal_blocks_compress_2d_out.csv'
  mesh_mode = 'REPLICATED'
  requirement = 'The system shall simulate correct contact behavior in 2D when two blocks with the same height come into contact using the dual basis'
  min_ad_size = 100
  abs_zero = 1.0e-8
[]

[frictionless-compress-2d-nondual]
  type = 'CSVDiff'
  input = 'two_equal_blocks_compress_2d.i'
  csvdiff = 'two_equal_blocks_compress_2d_nondual_out.csv'
  mesh_mode = 'REPLICATED'
  requirement = 'The system shall simulate correct contact behavior in 2D when two blocks with the same height come into contact using the standard (non-dual) basis'
  cli_args = 'Variables/normal_lm/use_dual=false Outputs/file_base="two_equal_blocks_compress_2d_nondual_out"'
  allow_test_objects = true
  min_ad_size = 100
  abs_zero = 1.0e-8
[]

[frictionless-compress-3d]
  type = 'CSVDiff'
  input = 'two_equal_blocks_compress_3d.i'
  csvdiff = 'two_equal_blocks_compress_3d_out.csv'
  mesh_mode = 'REPLICATED'
  requirement = 'The system shall simulate correct contact behavior in 3D when two blocks with the same height come into contact using the dual basis'
  allow_test_objects = true
  min_ad_size = 100
  abs_zero = 1.0e-8
[]

[frictionless-compress-3d-nondual]
  type = 'CSVDiff'
  input = 'two_equal_blocks_compress_3d.i'
  csvdiff = 'two_equal_blocks_compress_3d_nondual_out.csv'
  mesh_mode = 'REPLICATED'
  requirement = 'The system shall simulate correct contact behavior in 3D when two blocks with the same height come into contact using the standard (non-dual) basis'
  cli_args = 'Variables/normal_lm/use_dual=false Outputs/file_base="two_equal_blocks_compress_3d_nondual_out"'
  allow_test_objects = true
  min_ad_size = 100
  abs_zero = 1.0e-8
[]

[frictionless-slide-2d]
  type = 'CSVDiff'
  input = 'two_equal_blocks_slide_2d.i'
  csvdiff = 'two_equal_blocks_slide_2d_out.csv'
  mesh_mode = 'REPLICATED'
  requirement = 'The system shall simulate correct contact behavior in 2D when the node from a secondary mortar element does not project to the primary surface using the dual basis'
  allow_test_objects = true
  min_ad_size = 100
[]

[frictionless-slide-2d-nondual]
  type = 'CSVDiff'
  input = 'two_equal_blocks_slide_2d.i'
  csvdiff = 'two_equal_blocks_slide_2d_nondual_out.csv'
  mesh_mode = 'REPLICATED'
  requirement = 'The system shall simulate correct contact behavior in 2D when the node from a secondary mortar element does not project to the primary surface using the standard (non-dual) basis'
  cli_args = 'Variables/normal_lm/use_dual=false Outputs/file_base="two_equal_blocks_slide_2d_nondual_out"'
  allow_test_objects = true
  min_ad_size = 100
[]

[frictionless-slide-3d]
  type = 'CSVDiff'
  input = 'two_equal_blocks_slide_3d.i'
  csvdiff = 'two_equal_blocks_slide_3d_out.csv'
  mesh_mode = 'REPLICATED'
  requirement = 'The system shall simulate correct contact behavior in 3D when the node from a secondary mortar element does not project to the primary surface'
  allow_test_objects = true
  min_ad_size = 100
[]

[frictionless-slide-3d-nondual]
  type = 'CSVDiff'
  input = 'two_equal_blocks_slide_3d.i'
  csvdiff = 'two_equal_blocks_slide_3d_nondual_out.csv'
  mesh_mode = 'REPLICATED'
  requirement = 'The system shall simulate correct contact behavior in 3D when the node from a secondary mortar element does not project to the primary surface using the standard (non-dual) basis'
  cli_args = 'Variables/normal_lm/use_dual=false Outputs/file_base="two_equal_blocks_slide_3d_nondual_out"'
  allow_test_objects = true
  min_ad_size = 100
[]

[frictionless_lm]
  type = 'Exodiff'
  input = 'frictionless_lm.i'
  exodiff = 'frictionless_lm_out.e'
  heavy = true
  abs_zero = 1e-6
  design = 'Constraints/index.md'
  issues = '#13080 #20773'
  requirement = 'We shall be able to run our canonical frictional sliding block problem with lagrange multipliers and the mortar method'
  mesh_mode = 'replicated'
  ad_indexing_type = 'global'
  cli_args = "Outputs/active='out'"
[]

[./frictional_02_penalty]
  type = 'Exodiff'
  input = 'frictional_02_penalty.i'
  exodiff = 'frictional_02_penalty_out.e'
  heavy = true
  min_parallel = 4
  abs_zero = 1e-6
  max_time = 800
  issues = '#2816'
  design = 'MechanicalContactConstraint.md'
  requirement = 'We shall be able to solve the Coulomb friction sliding block problem using the penalty method and a friction coefficient of .2'
[../]

[./frictional_04_penalty]
  type = 'Exodiff'
  input = 'frictional_04_penalty.i'
  exodiff = 'frictional_04_penalty_out.e'
  heavy = true
  min_parallel = 4
  abs_zero = 1e-6
  max_time = 800
  issues = '#2816'
  design = 'MechanicalContactConstraint.md'
  requirement = 'We shall be able to solve the Coulomb friction sliding block problem using the penalty method and a friction coefficient of .4'
[../]

[./frictionless_kinematic]
  type = 'Exodiff'
  input = 'frictionless_kinematic.i'
  exodiff = 'frictionless_kinematic_out.e'
  heavy = true
  min_parallel = 4
  abs_zero = 1e-6
  max_time = 800
  custom_cmp = 'frictionless_kinematic_out.cmp'
  issues = '#2816'
  design = 'MechanicalContactConstraint.md'
  requirement = 'We shall be able to solve the frictionless sliding block problem using a kinematic constraint formulation.'
[../]

[./frictionless_kinematic_scaling]
  type = 'Exodiff'
  input = 'frictionless_kinematic.i'
  exodiff = 'frictionless_kinematic_out.e'
  heavy = true
  min_parallel = 4
  abs_zero = 1e-6
  max_time = 800
  cli_args = 'Variables/disp_x/scaling=1e-6 Variables/disp_y/scaling=1e-6'
  prereq = 'frictionless_kinematic'
  custom_cmp = 'frictionless_kinematic_out.cmp'
  requirement = "Kinematic contact shall produce the same results regardless of whether
                   variable scaling is used or not"
  design = 'MechanicalContactConstraint.md'
  issues = '#11601'
[../]

[./frictionless_penalty]
  type = 'Exodiff'
  input = 'frictionless_penalty.i'
  exodiff = 'frictionless_penalty_out.e'
  heavy = true
  superlu = true
  min_parallel = 4
  abs_zero = 1e-6
  max_time = 800
  issues = '#2816'
  design = 'MechanicalContactConstraint.md'
  requirement = 'We shall be able to solve the frictionless sliding block problem using a penalty constraint formulation.'
[../]

[./frictionless_penalty_contact_line_search]
  type = 'Exodiff'
  input = 'frictionless_penalty_contact_line_search.i'
  exodiff = 'frictionless_penalty_contact_line_search_out.e'
  heavy = true
  superlu = true
  min_parallel = 4
  abs_zero = 1e-6
  max_time = 800
  # The linesearch used by this test requires PETSc >= 3.6.0.
  petsc_version = '>=3.6.0'
  issues = '#10950'
  design = 'ContactLineSearch.md'
  requirement = 'We shall be able to solve the frictionless sliding block problem with a line serach customized for mechanical contact.'
[../]

[./frictional_mu02_penalty]
  type = 'Exodiff'
  input = 'frictional_02_penalty.i'
  exodiff = 'frictional_02_penalty_out.e'
  heavy = true
  superlu = true
  min_parallel = 4
  abs_zero = 1e-7
  max_time = 800
  allow_warnings = true
  rel_err = 1e-4
  requirement = 'The system shall support mechanics frictional contact problems'
[../]

[./frictional_mu04_penalty]
  type = 'Exodiff'
  input = 'frictional_04_penalty.i'
  exodiff = 'frictional_04_penalty_out.e'
  heavy = true
  min_parallel = 4
  abs_zero = 1e-7
  rel_err = 1e-4
  max_time = 800
  #    deleted = '#6139'
  requirement = 'The system shall support mechanics frictional contact problems'
[../]

[./frictionless_kinematic]
  type = 'Exodiff'
  input = 'frictionless_kinematic.i'
  exodiff = 'frictionless_kinematic_out.e'
  heavy = true
  min_parallel = 4
  abs_zero = 1e-7
  rel_err = 1e-4
  max_time = 800
  #    deleted = '#6139'
  requirement = 'The system shall support mechanics frictionless contact problems'
[../]

[./frictionless_penalty]
  type = 'Exodiff'
  input = 'frictionless_penalty.i'
  exodiff = 'frictionless_penalty_out.e'
  heavy = true
  superlu = true
  min_parallel = 4
  abs_zero = 1e-7
  rel_err = 1e-4
  max_time = 800
  requirement = 'The system shall support mechanics frictionless contact problems'
[../]

[./frictionless_aug]
  type = 'Exodiff'
  input = 'frictionless_aug.i'
  exodiff = 'frictionless_aug_out.e'
  heavy = true
  superlu = true
  min_parallel = 4
  abs_zero = 1e-4
  rel_err = 1e-4
  max_time = 800
  requirement = 'The system shall support mechanics frictionless contact problems'
[../]

[./frictional_mu02_aug]
  type = 'Exodiff'
  input = 'frictional_02_aug.i'
  exodiff = 'frictional_02_aug_out.e'
  heavy = true
  superlu = true
  min_parallel = 4
  rel_err = 1e-4
  abs_zero = 1e-4
  max_time = 800
  requirement = 'The system shall support mechanics frictional contact problems'
[../]

[auto-scaling]
  type = Exodiff
  input = bouncing-block-tan-pen.i
  exodiff = bouncing-block-tan-pen_exo.e
  requirement = 'The system shall be able to apply automatic scaling in conjunection with tangential penalty contact constraint enforcement and show no penetration and exhibit good nonlinear convergence'
  cli_args = "Postprocessors/active=''"
  abs_zero = 1e-8
[]

[no-scaling]
  type = Exodiff
  input = bouncing-block-tan-pen.i
  exodiff = no-scaling.e
  requirement = 'The system shall yield the same physical results when solving a tangential penalty contact problem with and without automatic scaling'
  cli_args = "Postprocessors/active='' Executioner/automatic_scaling=false Outputs/exo/file_base=no-scaling"
  abs_zero = 1e-8
[]

[./4ElemTensionRelease]
  type = 'Exodiff'
  input = 4ElemTensionRelease.i
  exodiff = 4ElemTensionRelease_out.e
  custom_cmp = '4ElemTensionRelease.exodiff'
  min_parallel = 4
  requirement = 'The contact system shall enforce and release contact conditions. 4 elements.'
  # This test has substantial diffs in older versions of PETSc due to the changes
  # in the way residuals are computed in contact problems in 3061bbd5d.
  petsc_version = '>=3.5.0'
[../]

[./4ElemTensionRelease_mechanical_constraint]
  type = 'Exodiff'
  input = 4ElemTensionRelease_mechanical_constraint.i
  exodiff = 4ElemTensionRelease_mechanical_constraint_out.e
  custom_cmp = '4ElemTensionRelease.exodiff'
  min_parallel = 4
  requirement = 'The contact system shall enforce and release contact conditions. 4 elements and mechanical constraints.'
  # This test has substantial diffs in older versions of PETSc due to the changes
  # in the way residuals are computed in contact problems in 3061bbd5d.
  petsc_version = '>=3.5.0'
[../]

[./4ElemTensionRelease_no_new_nonzeros]
  type = 'Exodiff'
  input = 4ElemTensionRelease_mechanical_constraint.i
  exodiff = 4ElemTensionRelease_no_new_nonzeros_out.e
  custom_cmp = '4ElemTensionRelease.exodiff'
  cli_args = 'Problem/error_on_jacobian_nonzero_reallocation=true Outputs/file_base=4ElemTensionRelease_no_new_nonzeros_out'
  max_parallel = 1
  # This test has substantial diffs in older versions of PETSc due to the changes
  # in the way residuals are computed in contact problems in 3061bbd5d.
  petsc_version = '>=3.5.0'
  requirement = 'The contact system shall enforce and release contact conditions. 4 elements and ensure no new Jacobian allocations.'

[../]

[./8ElemTensionRelease]
  type = 'Exodiff'
  input = 8ElemTensionRelease.i
  exodiff = 8ElemTensionRelease_out.e
  custom_cmp = '8ElemTensionRelease.exodiff'
  min_parallel = 4
  max_parallel = 4
  mesh_mode = REPLICATED
  # This test has substantial diffs in older versions of PETSc due to the changes
  # in the way residuals are computed in contact problems in 3061bbd5d.
  petsc_version = '>=3.5.0'
  requirement = 'The contact system shall enforce and release contact conditions. 8 elements.'

[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template1.i'
  csvdiff = 'hertz_cyl_half_1deg_glued_kin_check2.csv hertz_cyl_half_1deg_glued_kin_check_cont_press_0010.csv hertz_cyl_half_1deg_glued_kin_check_x_disp_0010.csv hertz_cyl_half_1deg_glued_kin_check_y_disp_0010.csv'
  cli_args = 'Executioner/end_time=1.0 Contact/interface/model=glued Contact/interface/formulation=kinematic
    Outputs/file_base=hertz_cyl_half_1deg_glued_kin_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_glued_kin_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_glued_kin_check2'
  rel_err = 2e-4
  abs_zero = 2e-6
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D Hertz half-symmetry cylindrical contact problem."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template1.i'
  csvdiff = 'hertz_cyl_half_1deg_glued_pen_check2.csv hertz_cyl_half_1deg_glued_pen_check_cont_press_0010.csv hertz_cyl_half_1deg_glued_pen_check_x_disp_0010.csv hertz_cyl_half_1deg_glued_pen_check_y_disp_0010.csv'
  cli_args = 'Executioner/end_time=1.0 Contact/interface/model=glued Contact/interface/formulation=penalty
    Outputs/file_base=hertz_cyl_half_1deg_glued_pen_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_glued_pen_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_glued_pen_check2'
  rel_err = 2e-4
  abs_zero = 1e-6
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 2D Hertz half-symmetry cylindrical contact problem."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template1.i'
  csvdiff = 'hertz_cyl_half_1deg_frictionless_kin_check2.csv
    hertz_cyl_half_1deg_frictionless_kin_check_cont_press_0035.csv
    hertz_cyl_half_1deg_frictionless_kin_check_x_disp_0035.csv
    hertz_cyl_half_1deg_frictionless_kin_check_y_disp_0035.csv'
  cli_args = 'Contact/interface/model=frictionless Contact/interface/formulation=kinematic
    Outputs/file_base=hertz_cyl_half_1deg_frictionless_kin_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_frictionless_kin_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_frictionless_kin_check2
    Outputs/chkfile/start_time=3.49'
  rel_err = 1e-5
  abs_zero = 3e-6
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 2D Hertz half-symmetry cylindrical contact problem."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template1.i'
  csvdiff = 'hertz_cyl_half_1deg_frictionless_pen_check2.csv
    hertz_cyl_half_1deg_frictionless_pen_check_cont_press_0035.csv
    hertz_cyl_half_1deg_frictionless_pen_check_x_disp_0035.csv
    hertz_cyl_half_1deg_frictionless_pen_check_y_disp_0035.csv'
  cli_args = 'Contact/interface/model=frictionless Contact/interface/formulation=penalty
    Outputs/file_base=hertz_cyl_half_1deg_frictionless_pen_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_frictionless_pen_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_frictionless_pen_check2
    Outputs/chkfile/start_time=3.49'
  rel_err = 1e-5
  abs_zero = 3e-6
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 2D Hertz half-symmetry cylindrical contact problem."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template1.i'
  csvdiff = 'hertz_cyl_half_1deg_frictionless_aug_check2.csv
    hertz_cyl_half_1deg_frictionless_aug_check_cont_press_0035.csv
    hertz_cyl_half_1deg_frictionless_aug_check_x_disp_0035.csv
    hertz_cyl_half_1deg_frictionless_aug_check_y_disp_0035.csv'
  cli_args = 'Problem/type=AugmentedLagrangianContactProblem Problem/maximum_lagrangian_update_iterations=200
    Contact/interface/model=frictionless
    Contact/interface/formulation=augmented_lagrange
    Contact/interface/al_penetration_tolerance=1e-12
    Outputs/file_base=hertz_cyl_half_1deg_frictionless_aug_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_frictionless_aug_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_frictionless_aug_check2
    Outputs/chkfile/start_time=3.49'
  rel_err = 1e-5
  abs_zero = 3e-6
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 2D Hertz half-symmetry cylindrical contact problem."
[../]

[./mu_0_pen]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template3.i'
  csvdiff = 'hertz_cyl_half_1deg_mu_0_pen_check2.csv hertz_cyl_half_1deg_mu_0_pen_check_cont_press_0020.csv hertz_cyl_half_1deg_mu_0_pen_check_x_disp_0020.csv hertz_cyl_half_1deg_mu_0_pen_check_y_disp_0020.csv'
  cli_args = 'Contact/interface/friction_coefficient=0.0
    Outputs/file_base=hertz_cyl_half_1deg_mu_0_pen_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_mu_0_pen_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_mu_0_pen_check2
    Outputs/chkfile/start_time=1.99'
  rel_err = 1e-5
  abs_zero = 1e-4
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D Hertz half-symmetry cylindrical contact problem with friction coefficient
    of 0."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template3.i'
  csvdiff = 'hertz_cyl_half_1deg_mu_0_2_pen_check2.csv hertz_cyl_half_1deg_mu_0_2_pen_check_cont_press_0060.csv hertz_cyl_half_1deg_mu_0_2_pen_check_x_disp_0060.csv hertz_cyl_half_1deg_mu_0_2_pen_check_y_disp_0060.csv'
  cli_args = 'Executioner/end_time=6.0 Contact/interface/friction_coefficient=0.2
    Outputs/file_base=hertz_cyl_half_1deg_mu_0_2_pen_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_mu_0_2_pen_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_mu_0_2_pen_check2
    Outputs/chkfile/start_time=5.99'
  max_time = 1500
  rel_err = 1e-5
  abs_zero = 3e-5
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D Hertz half-symmetry cylindrical contact problem with friction coefficient
    of 0.2."
[../]

[./mu_1_0_pen]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template3.i'
  csvdiff = 'hertz_cyl_half_1deg_mu_1_0_pen_check2.csv hertz_cyl_half_1deg_mu_1_0_pen_check_cont_press_0010.csv hertz_cyl_half_1deg_mu_1_0_pen_check_x_disp_0010.csv hertz_cyl_half_1deg_mu_1_0_pen_check_y_disp_0010.csv'
  cli_args = 'Executioner/end_time=1.0 Contact/interface/friction_coefficient=1.0 Contact/interface/penalty=5e+8
    Outputs/file_base=hertz_cyl_half_1deg_mu_1_0_pen_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_mu_1_0_pen_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_mu_1_0_pen_check2
    Outputs/chkfile/start_time=0.99'
  rel_err = 1e-5
  abs_zero = 1e-6
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D Hertz half-symmetry cylindrical contact problem with friction coefficient
    of 1.0."
[../]

[./mu_0_2_aug]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template1.i'
  csvdiff = 'hertz_cyl_half_1deg_mu_0_2_aug_check2.csv hertz_cyl_half_1deg_mu_0_2_aug_check_cont_press_0030.csv hertz_cyl_half_1deg_mu_0_2_aug_check_x_disp_0030.csv hertz_cyl_half_1deg_mu_0_2_aug_check_y_disp_0030.csv'
  cli_args = 'Problem/type=AugmentedLagrangianContactProblem Problem/maximum_lagrangian_update_iterations=200
    Contact/interface/formulation=augmented_lagrange Contact/interface/al_penetration_tolerance=1e-12
    Executioner/end_time=3.0 Contact/interface/friction_coefficient=0.2
    Outputs/file_base=hertz_cyl_half_1deg_mu_0_2_aug_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_mu_0_2_aug_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_mu_0_2_aug_check2
    Outputs/chkfile/start_time=2.99'
  max_time = 1500
  rel_err = 1e-5
  abs_zero = 3e-5
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictional, Augmented
    Lagrange contact for 2D Hertz half-symmetry cylindrical contact problem."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template1.i'
  csvdiff = 'hertz_cyl_half_1deg_glued_kin_check2.csv hertz_cyl_half_1deg_glued_kin_check_cont_press_0010.csv hertz_cyl_half_1deg_glued_kin_check_x_disp_0010.csv hertz_cyl_half_1deg_glued_kin_check_y_disp_0010.csv'
  cli_args = 'Executioner/end_time=1.0 Contact/interface/model=glued Contact/interface/formulation=kinematic
    Outputs/file_base=hertz_cyl_half_1deg_glued_kin_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_glued_kin_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_glued_kin_check2'
  rel_err = 1e-5
  abs_zero = 1e-6
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D Hertz half-symmetry cylindrical contact problem using higher order QUAD8 elements."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template1.i'
  csvdiff = 'hertz_cyl_half_1deg_glued_pen_check2.csv hertz_cyl_half_1deg_glued_pen_check_cont_press_0010.csv hertz_cyl_half_1deg_glued_pen_check_x_disp_0010.csv hertz_cyl_half_1deg_glued_pen_check_y_disp_0010.csv'
  cli_args = 'Executioner/end_time=1.0 Contact/interface/model=glued Contact/interface/formulation=penalty
    Outputs/file_base=hertz_cyl_half_1deg_glued_pen_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_glued_pen_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_glued_pen_check2'
  rel_err = 1e-5
  abs_zero = 1e-5
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 2D Hertz half-symmetry cylindrical contact problem using higher order QUAD8 elements."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template1.i'
  csvdiff = 'hertz_cyl_half_1deg_frictionless_kin_check2.csv
    hertz_cyl_half_1deg_frictionless_kin_check_cont_press_0035.csv
    hertz_cyl_half_1deg_frictionless_kin_check_x_disp_0035.csv
    hertz_cyl_half_1deg_frictionless_kin_check_y_disp_0035.csv'
  cli_args = 'Contact/interface/model=frictionless Contact/interface/formulation=kinematic
    Outputs/file_base=hertz_cyl_half_1deg_frictionless_kin_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_frictionless_kin_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_frictionless_kin_check2
    Outputs/chkfile/start_time=3.49'
  max_time = 1000
  rel_err = 1e-5
  abs_zero = 2e-6
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 2D Hertz half-symmetry cylindrical contact problem using higher order QUAD8 elements."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template1.i'
  csvdiff = 'hertz_cyl_half_1deg_frictionless_pen_check2.csv
    hertz_cyl_half_1deg_frictionless_pen_check_cont_press_0035.csv
    hertz_cyl_half_1deg_frictionless_pen_check_x_disp_0035.csv
    hertz_cyl_half_1deg_frictionless_pen_check_y_disp_0035.csv'
  cli_args = 'Contact/interface/model=frictionless Contact/interface/formulation=penalty
    Outputs/file_base=hertz_cyl_half_1deg_frictionless_pen_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_frictionless_pen_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_frictionless_pen_check2
    Outputs/chkfile/start_time=3.49'
  max_time = 1000
  rel_err = 3e-4
  abs_zero = 2e-6
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 2D Hertz half-symmetry cylindrical contact problem using higher order QUAD8 elements."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template1.i'
  csvdiff = 'hertz_cyl_half_1deg_frictionless_aug_check2.csv
    hertz_cyl_half_1deg_frictionless_aug_check_cont_press_0035.csv
    hertz_cyl_half_1deg_frictionless_aug_check_x_disp_0035.csv
    hertz_cyl_half_1deg_frictionless_aug_check_y_disp_0035.csv'
  cli_args = 'Problem/type=AugmentedLagrangianContactProblem Problem/maximum_lagrangian_update_iterations=200
    Contact/interface/model=frictionless Contact/interface/formulation=augmented_lagrange
    Contact/interface/al_penetration_tolerance=1e-12
    Outputs/file_base=hertz_cyl_half_1deg_frictionless_aug_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_frictionless_aug_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_frictionless_aug_check2
    Outputs/chkfile/start_time=3.49'
  max_time = 1000
  rel_err = 5e-4
  abs_zero = 2e-6
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 2D Hertz half-symmetry cylindrical contact problem using higher order QUAD8 elements."
[../]

[./mu_0_pen]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template3.i'
  csvdiff = 'hertz_cyl_half_1deg_mu_0_pen_check2.csv hertz_cyl_half_1deg_mu_0_pen_check_cont_press_0020.csv hertz_cyl_half_1deg_mu_0_pen_check_x_disp_0020.csv hertz_cyl_half_1deg_mu_0_pen_check_y_disp_0020.csv'
  cli_args = 'Contact/interface/friction_coefficient=0.0
    Outputs/file_base=hertz_cyl_half_1deg_mu_0_pen_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_mu_0_pen_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_mu_0_pen_check2
    Outputs/chkfile/start_time=1.99'
  max_time = 1000
  rel_err = 3e-5
  abs_zero = 2e-6
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D Hertz half-symmetry cylindrical contact problem using higher order QUAD8
    elements and with a friction coefficient of 0."
[../]

[./mu_1_0_pen]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template3.i'
  csvdiff = 'hertz_cyl_half_1deg_mu_1_0_pen_check2.csv hertz_cyl_half_1deg_mu_1_0_pen_check_cont_press_0010.csv hertz_cyl_half_1deg_mu_1_0_pen_check_x_disp_0010.csv hertz_cyl_half_1deg_mu_1_0_pen_check_y_disp_0010.csv'
  cli_args = 'Executioner/end_time=1.0 Contact/interface/friction_coefficient=1.0 Contact/interface/penalty=2e+8
    Outputs/file_base=hertz_cyl_half_1deg_mu_1_0_pen_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_mu_1_0_pen_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_mu_1_0_pen_check2
    Outputs/chkfile/start_time=0.99'
  max_time = 1000
  rel_err = 1e-5
  abs_zero = 1e-6
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D Hertz half-symmetry cylindrical contact problem using higher order QUAD8
    elements and with a friction coefficient of 1.0."
[../]

[./mu_0_2_aug]
  type = 'CSVDiff'
  input = 'hertz_cyl_half_1deg_template1.i'
  csvdiff = 'hertz_cyl_half_1deg_mu_0_2_aug_check2.csv hertz_cyl_half_1deg_mu_0_2_aug_check_cont_press_0030.csv hertz_cyl_half_1deg_mu_0_2_aug_check_x_disp_0030.csv hertz_cyl_half_1deg_mu_0_2_aug_check_y_disp_0030.csv'
  cli_args = 'Problem/type=AugmentedLagrangianContactProblem Problem/maximum_lagrangian_update_iterations=200
    Contact/interface/formulation=augmented_lagrange Contact/interface/al_penetration_tolerance=1e-12
    Executioner/end_time=3.0 Contact/interface/friction_coefficient=0.2
    Outputs/file_base=hertz_cyl_half_1deg_mu_0_2_aug_out Outputs/chkfile/file_base=hertz_cyl_half_1deg_mu_0_2_aug_check
    Outputs/chkfile2/file_base=hertz_cyl_half_1deg_mu_0_2_aug_check2
    Outputs/chkfile/start_time=2.99'
  max_time = 1500
  rel_err = 2e-5
  abs_zero = 1e-5
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictional, Augmented
    Lagrange contact for 2D Hertz half-symmetry cylindrical contact problem using higher order QUAD8 elements."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'hertz_cyl_qsym_1deg_template1.i'
  csvdiff = 'hertz_cyl_qsym_1deg_glued_kin_check.csv hertz_cyl_qsym_1deg_glued_kin_check_cont_press_0010.csv hertz_cyl_qsym_1deg_glued_kin_check_x_disp_0010.csv hertz_cyl_qsym_1deg_glued_kin_check_y_disp_0010.csv'
  cli_args = 'Contact/interface/model=glued Contact/interface/formulation=kinematic Outputs/file_base=hertz_cyl_qsym_1deg_glued_kin_out Outputs/chkfile/file_base=hertz_cyl_qsym_1deg_glued_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D Hertz quarter-symmetry cylindrical contact problem."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'hertz_cyl_qsym_1deg_template1.i'
  csvdiff = 'hertz_cyl_qsym_1deg_glued_pen_check.csv hertz_cyl_qsym_1deg_glued_pen_check_cont_press_0010.csv hertz_cyl_qsym_1deg_glued_pen_check_x_disp_0010.csv hertz_cyl_qsym_1deg_glued_pen_check_y_disp_0010.csv'
  cli_args = 'Contact/interface/model=glued Contact/interface/formulation=penalty Outputs/file_base=hertz_cyl_qsym_1deg_glued_pen_out Outputs/chkfile/file_base=hertz_cyl_qsym_1deg_glued_pen_check'
  rel_err = 2e-5
  abs_zero = 1e-8
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 2D Hertz quarter-symmetry cylindrical contact problem."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'hertz_cyl_qsym_1deg_template1.i'
  csvdiff = 'hertz_cyl_qsym_1deg_frictionless_kin_check.csv hertz_cyl_qsym_1deg_frictionless_kin_check_cont_press_0010.csv hertz_cyl_qsym_1deg_frictionless_kin_check_x_disp_0010.csv hertz_cyl_qsym_1deg_frictionless_kin_check_y_disp_0010.csv'
  cli_args = 'Contact/interface/model=frictionless Contact/interface/formulation=kinematic Outputs/file_base=hertz_cyl_qsym_1deg_frictionless_kin_out Outputs/chkfile/file_base=hertz_cyl_qsym_1deg_frictionless_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 2D Hertz quarter-symmetry cylindrical contact problem."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'hertz_cyl_qsym_1deg_template1.i'
  csvdiff = 'hertz_cyl_qsym_1deg_frictionless_pen_check.csv hertz_cyl_qsym_1deg_frictionless_pen_check_cont_press_0010.csv hertz_cyl_qsym_1deg_frictionless_pen_check_x_disp_0010.csv hertz_cyl_qsym_1deg_frictionless_pen_check_y_disp_0010.csv'
  cli_args = 'Contact/interface/model=frictionless Contact/interface/formulation=penalty Outputs/file_base=hertz_cyl_qsym_1deg_frictionless_pen_out Outputs/chkfile/file_base=hertz_cyl_qsym_1deg_frictionless_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 2D Hertz quarter-symmetry cylindrical contact problem."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'hertz_cyl_qsym_1deg_template1.i'
  csvdiff = 'hertz_cyl_qsym_1deg_frictionless_aug_check.csv hertz_cyl_qsym_1deg_frictionless_aug_check_cont_press_0010.csv hertz_cyl_qsym_1deg_frictionless_aug_check_x_disp_0010.csv hertz_cyl_qsym_1deg_frictionless_aug_check_y_disp_0010.csv'
  cli_args = 'Problem/type=AugmentedLagrangianContactProblem Problem/maximum_lagrangian_update_iterations=200
    Contact/interface/model=frictionless Contact/interface/formulation=augmented_lagrange Contact/interface/al_penetration_tolerance=1e-12
    Outputs/file_base=hertz_cyl_qsym_1deg_frictionless_aug_out Outputs/chkfile/file_base=hertz_cyl_qsym_1deg_frictionless_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 2D Hertz quarter-symmetry cylindrical contact problem."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'hertz_cyl_qsym_1deg_template1.i'
  csvdiff = 'hertz_cyl_qsym_1deg_glued_kin_check.csv hertz_cyl_qsym_1deg_glued_kin_check_cont_press_0010.csv hertz_cyl_qsym_1deg_glued_kin_check_x_disp_0010.csv hertz_cyl_qsym_1deg_glued_kin_check_y_disp_0010.csv'
  cli_args = 'Contact/interface/model=glued Contact/interface/formulation=kinematic Outputs/file_base=hertz_cyl_qsym_1deg_glued_kin_out Outputs/chkfile/file_base=hertz_cyl_qsym_1deg_glued_kin_check'
  max_time = 1000
  rel_err = 3e-5
  abs_zero = 1e-6
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D Hertz quarter-symmetry cylindrical contact problem using higher order QUAD8 elements."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'hertz_cyl_qsym_1deg_template1.i'
  csvdiff = 'hertz_cyl_qsym_1deg_glued_pen_check.csv hertz_cyl_qsym_1deg_glued_pen_check_cont_press_0010.csv hertz_cyl_qsym_1deg_glued_pen_check_x_disp_0010.csv hertz_cyl_qsym_1deg_glued_pen_check_y_disp_0010.csv'
  cli_args = 'Contact/interface/model=glued Contact/interface/formulation=penalty Outputs/file_base=hertz_cyl_qsym_1deg_glued_pen_out Outputs/chkfile/file_base=hertz_cyl_qsym_1deg_glued_pen_check'
  max_time = 1000
  rel_err = 3e-5
  abs_zero = 1e-6
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 2D Hertz quarter-symmetry cylindrical contact problem using higher order QUAD8 elements."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'hertz_cyl_qsym_1deg_template1.i'
  csvdiff = 'hertz_cyl_qsym_1deg_frictionless_kin_check.csv hertz_cyl_qsym_1deg_frictionless_kin_check_cont_press_0010.csv hertz_cyl_qsym_1deg_frictionless_kin_check_x_disp_0010.csv hertz_cyl_qsym_1deg_frictionless_kin_check_y_disp_0010.csv'
  cli_args = 'Contact/interface/model=frictionless Contact/interface/formulation=kinematic Outputs/file_base=hertz_cyl_qsym_1deg_frictionless_kin_out Outputs/chkfile/file_base=hertz_cyl_qsym_1deg_frictionless_kin_check'
  max_time = 1000
  rel_err = 4e-4
  abs_zero = 1e-6
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 2D Hertz quarter-symmetry cylindrical contact problem using higher order QUAD8 elements."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'hertz_cyl_qsym_1deg_template1.i'
  csvdiff = 'hertz_cyl_qsym_1deg_frictionless_pen_check.csv hertz_cyl_qsym_1deg_frictionless_pen_check_cont_press_0010.csv hertz_cyl_qsym_1deg_frictionless_pen_check_x_disp_0010.csv hertz_cyl_qsym_1deg_frictionless_pen_check_y_disp_0010.csv'
  cli_args = 'Contact/interface/model=frictionless Contact/interface/formulation=penalty Outputs/file_base=hertz_cyl_qsym_1deg_frictionless_pen_out Outputs/chkfile/file_base=hertz_cyl_qsym_1deg_frictionless_pen_check'
  max_time = 1000
  rel_err = 4e-4
  abs_zero = 5e-6
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 2D Hertz quarter-symmetry cylindrical contact problem using higher order QUAD8 elements."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'hertz_cyl_qsym_1deg_template1.i'
  csvdiff = 'hertz_cyl_qsym_1deg_frictionless_aug_check.csv hertz_cyl_qsym_1deg_frictionless_aug_check_cont_press_0010.csv hertz_cyl_qsym_1deg_frictionless_aug_check_x_disp_0010.csv hertz_cyl_qsym_1deg_frictionless_aug_check_y_disp_0010.csv'
  cli_args = 'Problem/type=AugmentedLagrangianContactProblem Problem/maximum_lagrangian_update_iterations=200
    Contact/interface/model=frictionless Contact/interface/formulation=augmented_lagrange Contact/interface/al_penetration_tolerance=1e-12
    Outputs/file_base=hertz_cyl_qsym_1deg_frictionless_aug_out Outputs/chkfile/file_base=hertz_cyl_qsym_1deg_frictionless_aug_check'
  rel_err = 4e-4
  abs_zero = 1e-8
  max_parallel = 1
  heavy = true
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 2D Hertz quarter-symmetry cylindrical contact problem using higher order QUAD8 elements."
[../]

[./test]
  type = 'Exodiff'
  input = 'overclosure.i'
  exodiff = 'overclosure_exodus.e'
  rel_err = 1e-5
  abs_zero = 1e-6
  superlu = true
  requirement = "The Contact system shall enforce eliminating initial
                   overclosure between the primary and secondary surfaces."
[../]

[./sliding_blocks]
  type = 'Exodiff'
  input = 'sliding_update.i'
  exodiff = 'sliding_update_out.e'
  abs_zero = 1e-7
  superlu = true
  requirement = "The Contact system shall enforce the automatic patch update
                   using the iteration option."
[../]

[./error_test]
  type = RunException
  input = 'sliding_update.i'
  cli_args = 'Mesh/ghosting_patch_size=2
                Outputs/file_base=sliding_error_out'
  expect_err = "Error in NearestNodeLocator : The nearest neighbor lies outside the ghosted set of elements. Increase the ghosting_patch_size parameter in the mesh block and try again."
  prereq = sliding_blocks
  min_parallel = 2
  mesh_mode = REPLICATED
  requirement = "The Contact system shall enforce that the nearest neighbor
                   node is inside the ghosted set of elements."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'brick1_template1.i'
  csvdiff = 'brick1_glued_kin_check.csv brick1_glued_kin_check_cont_press_0001.csv brick1_glued_kin_check_x_disp_0001.csv brick1_glued_kin_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=kinematic Outputs/file_base=brick1_glued_kin_out Outputs/chkfile/file_base=brick1_glued_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 3D brick geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'brick1_template1.i'
  csvdiff = 'brick1_glued_pen_check.csv brick1_glued_pen_check_cont_press_0001.csv brick1_glued_pen_check_x_disp_0001.csv brick1_glued_pen_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=penalty Outputs/file_base=brick1_glued_pen_out Outputs/chkfile/file_base=brick1_glued_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 3D brick geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'brick1_template1.i'
  csvdiff = 'brick1_frictionless_kin_check.csv brick1_frictionless_kin_check_cont_press_0001.csv brick1_frictionless_kin_check_x_disp_0001.csv brick1_frictionless_kin_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic Outputs/file_base=brick1_frictionless_kin_out Outputs/chkfile/file_base=brick1_frictionless_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 3D brick geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'brick1_template1.i'
  csvdiff = 'brick1_frictionless_pen_check.csv brick1_frictionless_pen_check_cont_press_0001.csv brick1_frictionless_pen_check_x_disp_0001.csv brick1_frictionless_pen_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Outputs/file_base=brick1_frictionless_pen_out Outputs/chkfile/file_base=brick1_frictionless_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 3D brick geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'brick1_template2.i'
  csvdiff = 'brick1_frictionless_aug_check.csv brick1_frictionless_aug_check_cont_press_0001.csv brick1_frictionless_aug_check_x_disp_0001.csv brick1_frictionless_aug_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=augmented_lagrange Outputs/file_base=brick1_frictionless_aug_out Outputs/chkfile/file_base=brick1_frictionless_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 3D brick geometry (NAFEMS CGS1 contact patch test
    with matched nodes)."
[../]

[./mu_0_2_aug]
  type = 'CSVDiff'
  input = 'brick1_template2.i'
  csvdiff = 'brick1_mu_0_2_aug_check.csv brick1_mu_0_2_aug_check_cont_press_0001.csv brick1_mu_0_2_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=coulomb Contact/leftright/formulation=augmented_lagrange Contact/leftright/friction_coefficient=0.2 Contact/leftright/al_incremental_slip_tolerance=1e-6 Contact/leftright/al_frictional_force_tolerance=5e-2 Outputs/file_base=brick1_mu_0_2_aug_out Outputs/chkfile/file_base=brick1_mu_0_2_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, Augmented
    Lagrange contact for 3D brick geometry (NAFEMS CGS1 contact patch test
    with matched nodes)."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'brick1_mu_0_2_pen.i'
  csvdiff = 'brick1_mu_0_2_pen_check.csv brick1_mu_0_2_pen_check_cont_press_0001.csv brick1_mu_0_2_pen_check_x_disp_0001.csv'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 3D brick geometry (NAFEMS CGS1 contact patch test
    with matched nodes)."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'brick2_template1.i'
  csvdiff = 'brick2_glued_kin_check.csv brick2_glued_kin_check_cont_press_0001.csv brick2_glued_kin_check_x_disp_0001.csv brick2_glued_kin_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=kinematic Outputs/file_base=brick2_glued_kin_out Outputs/chkfile/file_base=brick2_glued_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 3D brick geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'brick2_template1.i'
  csvdiff = 'brick2_glued_pen_check.csv brick2_glued_pen_check_cont_press_0001.csv brick2_glued_pen_check_x_disp_0001.csv brick2_glued_pen_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=penalty Outputs/file_base=brick2_glued_pen_out Outputs/chkfile/file_base=brick2_glued_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 3D brick geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'brick2_template1.i'
  csvdiff = 'brick2_frictionless_kin_check.csv brick2_frictionless_kin_check_cont_press_0001.csv brick2_frictionless_kin_check_x_disp_0001.csv brick2_frictionless_kin_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic Outputs/file_base=brick2_frictionless_kin_out Outputs/chkfile/file_base=brick2_frictionless_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 3D brick geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'brick2_template1.i'
  csvdiff = 'brick2_frictionless_pen_check.csv brick2_frictionless_pen_check_cont_press_0001.csv brick2_frictionless_pen_check_x_disp_0001.csv brick2_frictionless_pen_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Outputs/file_base=brick2_frictionless_pen_out Outputs/chkfile/file_base=brick2_frictionless_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 3D brick geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'brick2_template2.i'
  csvdiff = 'brick2_frictionless_aug_check.csv brick2_frictionless_aug_check_cont_press_0001.csv brick2_frictionless_aug_check_x_disp_0001.csv brick2_frictionless_aug_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=augmented_lagrange Outputs/file_base=brick2_frictionless_aug_out Outputs/chkfile/file_base=brick2_frictionless_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 3D brick geometry (NAFEMS CGS1 contact patch test
    with mismatched nodes)."
[../]

[./mu_0_2_aug]
  type = 'CSVDiff'
  input = 'brick2_template2.i'
  csvdiff = 'brick2_mu_0_2_aug_check.csv brick2_mu_0_2_aug_check_cont_press_0001.csv brick2_mu_0_2_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=coulomb Contact/leftright/formulation=augmented_lagrange Contact/leftright/friction_coefficient=0.2 Contact/leftright/al_incremental_slip_tolerance=1e-6 Contact/leftright/al_frictional_force_tolerance=5e-2 Outputs/file_base=brick2_mu_0_2_aug_out Outputs/chkfile/file_base=brick2_mu_0_2_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, Augmented
    Lagrange contact for 3D brick geometry (NAFEMS CGS1 contact patch test
    with mismatched nodes)."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'brick2_mu_0_2_pen.i'
  csvdiff = 'brick2_mu_0_2_pen_check.csv brick2_mu_0_2_pen_check_cont_press_0001.csv brick2_mu_0_2_pen_check_x_disp_0001.csv'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 3D brick geometry (NAFEMS CGS1 contact patch test
    with mismatched nodes)."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'brick3_template1.i'
  csvdiff = 'brick3_frictionless_kin_check.csv brick3_frictionless_kin_check_cont_press_0001.csv brick3_frictionless_kin_check_x_disp_0001.csv brick3_frictionless_kin_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic Outputs/file_base=brick3_frictionless_kin_out Outputs/chkfile/file_base=brick3_frictionless_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 3D HEX20 brick geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'brick3_template1.i'
  csvdiff = 'brick3_frictionless_pen_check.csv brick3_frictionless_pen_check_cont_press_0001.csv brick3_frictionless_pen_check_x_disp_0001.csv brick3_frictionless_pen_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Outputs/file_base=brick3_frictionless_pen_out Outputs/chkfile/file_base=brick3_frictionless_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 3D HEX20 brick geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'brick3_template2.i'
  csvdiff = 'brick3_frictionless_aug_check.csv brick3_frictionless_aug_check_cont_press_0001.csv brick3_frictionless_aug_check_x_disp_0001.csv brick3_frictionless_aug_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=augmented_lagrange Outputs/file_base=brick3_frictionless_aug_out Outputs/chkfile/file_base=brick3_frictionless_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 3D HEX20 brick geometry (NAFEMS CGS1 contact patch
    test with matched nodes)."
[../]

[./mu_0_2_aug]
  type = 'CSVDiff'
  input = 'brick3_template2.i'
  csvdiff = 'brick3_mu_0_2_aug_check.csv brick3_mu_0_2_aug_check_cont_press_0001.csv brick3_mu_0_2_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=coulomb Contact/leftright/formulation=augmented_lagrange Contact/leftright/friction_coefficient=0.2 Contact/leftright/al_incremental_slip_tolerance=1e-6 Contact/leftright/al_frictional_force_tolerance=5e-2 Outputs/file_base=brick3_mu_0_2_aug_out Outputs/chkfile/file_base=brick3_mu_0_2_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, Augmented
    Lagrange contact for 3D HEX20 brick geometry (NAFEMS CGS1 contact patch
    test with matched nodes)."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'brick3_mu_0_2_pen.i'
  csvdiff = 'brick3_mu_0_2_pen_check.csv brick3_mu_0_2_pen_check_cont_press_0001.csv brick3_mu_0_2_pen_check_x_disp_0001.csv'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 3D HEX20 brick geometry (NAFEMS CGS1 contact patch
    test with matched nodes)."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'brick4_template1.i'
  csvdiff = 'brick4_glued_kin_check.csv brick4_glued_kin_check_cont_press_0001.csv brick4_glued_kin_check_x_disp_0001.csv brick4_glued_kin_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=kinematic Outputs/file_base=brick4_glued_kin_out Outputs/chkfile/file_base=brick4_glued_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 3D HEX20 brick geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'brick4_template1.i'
  csvdiff = 'brick4_glued_pen_check.csv brick4_glued_pen_check_cont_press_0001.csv brick4_glued_pen_check_x_disp_0001.csv brick4_glued_pen_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=penalty Outputs/file_base=brick4_glued_pen_out Outputs/chkfile/file_base=brick4_glued_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 3D HEX20 brick geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'brick4_template1.i'
  csvdiff = 'brick4_frictionless_kin_check.csv brick4_frictionless_kin_check_cont_press_0001.csv brick4_frictionless_kin_check_x_disp_0001.csv brick4_frictionless_kin_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic Outputs/file_base=brick4_frictionless_kin_out Outputs/chkfile/file_base=brick4_frictionless_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 3D HEX20 brick geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'brick4_template1.i'
  csvdiff = 'brick4_frictionless_pen_check.csv brick4_frictionless_pen_check_cont_press_0001.csv brick4_frictionless_pen_check_x_disp_0001.csv brick4_frictionless_pen_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Outputs/file_base=brick4_frictionless_pen_out Outputs/chkfile/file_base=brick4_frictionless_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 3D HEX20 brick geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'brick4_template2.i'
  csvdiff = 'brick4_frictionless_aug_check.csv brick4_frictionless_aug_check_cont_press_0001.csv brick4_frictionless_aug_check_x_disp_0001.csv brick4_frictionless_aug_check_y_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=augmented_lagrange Outputs/file_base=brick4_frictionless_aug_out Outputs/chkfile/file_base=brick4_frictionless_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-7
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 3D HEX20 brick geometry (NAFEMS CGS1 contact patch
    test with mismatched nodes)."
[../]

[./mu_0_2_aug]
  type = 'CSVDiff'
  input = 'brick4_template2.i'
  csvdiff = 'brick4_mu_0_2_aug_check.csv brick4_mu_0_2_aug_check_cont_press_0001.csv brick4_mu_0_2_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=coulomb Contact/leftright/formulation=augmented_lagrange Contact/leftright/friction_coefficient=0.2 Contact/leftright/al_incremental_slip_tolerance=1e-6 Contact/leftright/al_frictional_force_tolerance=5e-2 Outputs/file_base=brick4_mu_0_2_aug_out Outputs/chkfile/file_base=brick4_mu_0_2_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, Augmented
    Lagrange contact for 3D HEX20 brick geometry (NAFEMS CGS1 contact patch
    test with mismatched nodes)."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'brick4_mu_0_2_pen.i'
  csvdiff = 'brick4_mu_0_2_pen_check.csv brick4_mu_0_2_pen_check_cont_press_0001.csv brick4_mu_0_2_pen_check_x_disp_0001.csv'
  rel_err = 1e-5
  abs_zero = 2e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 3D HEX20 brick geometry (NAFEMS CGS1 contact patch
    test with mismatched nodes)."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'cyl1_template1.i'
  csvdiff = 'cyl1_glued_kin_check.csv cyl1_glued_kin_check_cont_press_0001.csv cyl1_glued_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=kinematic Outputs/file_base=cyl1_glued_kin_out Outputs/chkfile/file_base=cyl1_glued_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D axisymmetric geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'cyl1_template1.i'
  csvdiff = 'cyl1_glued_pen_check.csv cyl1_glued_pen_check_cont_press_0001.csv cyl1_glued_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=penalty Outputs/file_base=cyl1_glued_pen_out Outputs/chkfile/file_base=cyl1_glued_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 2D axisymmetric geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'cyl1_template1.i'
  csvdiff = 'cyl1_frictionless_kin_check.csv cyl1_frictionless_kin_check_cont_press_0001.csv cyl1_frictionless_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic Outputs/file_base=cyl1_frictionless_kin_out Outputs/chkfile/file_base=cyl1_frictionless_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 2D axisymmetric geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'cyl1_template1.i'
  csvdiff = 'cyl1_frictionless_pen_check.csv cyl1_frictionless_pen_check_cont_press_0001.csv cyl1_frictionless_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Outputs/file_base=cyl1_frictionless_pen_out Outputs/chkfile/file_base=cyl1_frictionless_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 2D axisymmetric geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'cyl1_template2.i'
  csvdiff = 'cyl1_frictionless_aug_check.csv cyl1_frictionless_aug_check_cont_press_0001.csv cyl1_frictionless_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=augmented_lagrange Outputs/file_base=cyl1_frictionless_aug_out Outputs/chkfile/file_base=cyl1_frictionless_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty
    Augmented Lagrange contact for 2D axisymmetric geometry (NAFEMS CGS1
    contact patch test with matched nodes)."
[../]

[./mu_0_2_aug]
  type = 'CSVDiff'
  input = 'cyl1_template2.i'
  csvdiff = 'cyl1_mu_0_2_aug_check.csv cyl1_mu_0_2_aug_check_cont_press_0001.csv cyl1_mu_0_2_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=coulomb Contact/leftright/formulation=augmented_lagrange Contact/leftright/friction_coefficient=0.2 Contact/leftright/al_incremental_slip_tolerance=1e-6 Contact/leftright/al_frictional_force_tolerance=5e-2 Outputs/file_base=cyl1_mu_0_2_aug_out Outputs/chkfile/file_base=cyl1_mu_0_2_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, Augmented
    Lagrange contact for 2D axisymmetric geometry (NAFEMS CGS1 contact patch test
    with matched nodes)."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'cyl1_mu_0_2_pen.i'
  csvdiff = 'cyl1_mu_0_2_pen_check.csv cyl1_mu_0_2_pen_check_cont_press_0001.csv cyl1_mu_0_2_pen_check_x_disp_0001.csv'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D axisymmetric geometry (NAFEMS CGS1 contact patch test
    with matched nodes)."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'cyl2_template1.i'
  csvdiff = 'cyl2_glued_kin_check.csv cyl2_glued_kin_check_cont_press_0001.csv cyl2_glued_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=kinematic Outputs/file_base=cyl2_glued_kin_out Outputs/chkfile/file_base=cyl2_glued_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D axisymmetric geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'cyl2_template1.i'
  csvdiff = 'cyl2_glued_pen_check.csv cyl2_glued_pen_check_cont_press_0001.csv cyl2_glued_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=penalty Outputs/file_base=cyl2_glued_pen_out Outputs/chkfile/file_base=cyl2_glued_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 2D axisymmetric geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'cyl2_template1.i'
  csvdiff = 'cyl2_frictionless_kin_check.csv cyl2_frictionless_kin_check_cont_press_0001.csv cyl2_frictionless_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic Outputs/file_base=cyl2_frictionless_kin_out Outputs/chkfile/file_base=cyl2_frictionless_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 2D axisymmetric geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'cyl2_template1.i'
  csvdiff = 'cyl2_frictionless_pen_check.csv cyl2_frictionless_pen_check_cont_press_0001.csv cyl2_frictionless_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Outputs/file_base=cyl2_frictionless_pen_out Outputs/chkfile/file_base=cyl2_frictionless_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 2D axisymmetric geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'cyl2_template2.i'
  csvdiff = 'cyl2_frictionless_aug_check.csv cyl2_frictionless_aug_check_cont_press_0001.csv cyl2_frictionless_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=augmented_lagrange Outputs/file_base=cyl2_frictionless_aug_out Outputs/chkfile/file_base=cyl2_frictionless_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 2D axisymmetric geometry (NAFEMS CGS1 contact patch
    test with mismatched nodes)."
[../]

[./mu_0_2_aug]
  type = 'CSVDiff'
  input = 'cyl2_template2.i'
  csvdiff = 'cyl2_mu_0_2_aug_check.csv cyl2_mu_0_2_aug_check_cont_press_0001.csv cyl2_mu_0_2_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=coulomb Contact/leftright/formulation=augmented_lagrange Contact/leftright/friction_coefficient=0.2 Contact/leftright/al_incremental_slip_tolerance=1e-6 Contact/leftright/al_frictional_force_tolerance=5e-2 Outputs/file_base=cyl2_mu_0_2_aug_out Outputs/chkfile/file_base=cyl2_mu_0_2_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, Augmented
    Lagrange contact for 2D axisymmetric geometry (NAFEMS CGS1 contact patch
    test with mismatched nodes)."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'cyl2_mu_0_2_pen.i'
  csvdiff = 'cyl2_mu_0_2_pen_check.csv cyl2_mu_0_2_pen_check_cont_press_0001.csv cyl2_mu_0_2_pen_check_x_disp_0001.csv'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D axisymmetric geometry (NAFEMS CGS1 contact patch
    test with mismatched nodes)."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'cyl3_template1.i'
  csvdiff = 'cyl3_glued_kin_check.csv cyl3_glued_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=kinematic
                Outputs/file_base=cyl3_glued_kin_out Outputs/chkfile/file_base=cyl3_glued_kin_check
                Executioner/petsc_options_iname="-pc_type -pc_factor_mat_solver_package" Executioner/petsc_options_value="lu superlu_dist"'
  rel_err = 1e-5
  abs_zero = 2e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D QUAD8 axisymmetric geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'cyl3_template1.i'
  csvdiff = 'cyl3_glued_pen_check.csv cyl3_glued_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=penalty
                Contact/leftright/normalize_penalty=true
                Outputs/chkfile/file_base=cyl3_glued_pen_check
                Executioner/petsc_options_iname="-pc_type -pc_factor_mat_solver_package"
                Executioner/petsc_options_value="lu superlu_dist"'
  rel_err = 1e-5
  abs_zero = 2e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 2D QUAD8 axisymmetric geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'cyl3_template1.i'
  csvdiff = 'cyl3_frictionless_kin_check.csv cyl3_frictionless_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic
                Outputs/file_base=cyl3_frictionless_kin_out Outputs/chkfile/file_base=cyl3_frictionless_kin_check
                Executioner/petsc_options_iname="-pc_type -pc_factor_mat_solver_package"
                Executioner/petsc_options_value="lu superlu_dist"'
  rel_err = 2e-5
  abs_zero = 1e-4
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  # This test is fails reliably on INTEL compilers with PETSc 3.7.4 and 3.7.5, so restrict the test to not INTEL
  compiler = 'GCC CLANG'
  petsc_version = '>=3.7.6'
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 2D QUAD8 axisymmetric geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'cyl3_template1.i'
  csvdiff = 'cyl3_frictionless_pen_check.csv cyl3_frictionless_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Contact/leftright/normalize_penalty=true
                Outputs/file_base=cyl3_frictionless_pen_out Outputs/chkfile/file_base=cyl3_frictionless_pen_check
                Executioner/petsc_options_iname="-pc_type -pc_factor_mat_solver_package"
                Executioner/petsc_options_value="lu superlu_dist"'
  rel_err = 1e-5
  abs_zero = 2e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 2D QUAD8 axisymmetric geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'cyl3_template2.i'
  csvdiff = 'cyl3_frictionless_aug_check.csv cyl3_frictionless_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=augmented_lagrange
                Contact/leftright/normalize_penalty=true Outputs/file_base=cyl3_frictionless_aug_out
                Outputs/chkfile/file_base=cyl3_frictionless_aug_check'
  rel_err = 2e-5
  abs_zero = 2e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 2D QUAD8 axisymmetric geometry (NAFEMS CGS1 contact
    patch test with matched nodes)."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'cyl3_mu_0_2_pen.i'
  csvdiff = 'cyl3_mu_0_2_pen_check.csv cyl3_mu_0_2_pen_check_cont_press_0001.csv cyl3_mu_0_2_pen_check_x_disp_0001.csv'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D QUAD8 axisymmetric geometry (NAFEMS CGS1 contact
    patch test with matched nodes)."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'cyl4_template1.i'
  csvdiff = 'cyl4_glued_kin_check.csv cyl4_glued_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=kinematic Outputs/file_base=cyl4_glued_kin_out Outputs/chkfile/file_base=cyl4_glued_kin_check'
  rel_err = 1e-5
  abs_zero = 2e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D QUAD8 axisymmetric geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'cyl4_template1.i'
  csvdiff = 'cyl4_glued_pen_check.csv cyl4_glued_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=penalty Contact/leftright/normalize_penalty=true Outputs/file_base=cyl4_glued_pen_out Outputs/chkfile/file_base=cyl4_glued_pen_check'
  rel_err = 1e-5
  abs_zero = 2e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 2D QUAD8 axisymmetric geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'cyl4_template1.i'
  csvdiff = 'cyl4_frictionless_kin_check.csv cyl4_frictionless_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic Outputs/file_base=cyl4_frictionless_kin_out Outputs/chkfile/file_base=cyl4_frictionless_kin_check'
  rel_err = 1e-5
  abs_zero = 2e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 2D QUAD8 axisymmetric geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'cyl4_template1.i'
  csvdiff = 'cyl4_frictionless_pen_check.csv cyl4_frictionless_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Contact/leftright/normalize_penalty=true Outputs/file_base=cyl4_frictionless_pen_out Outputs/chkfile/file_base=cyl4_frictionless_pen_check'
  rel_err = 1e-5
  abs_zero = 2e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 2D QUAD8 axisymmetric geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'cyl4_template2.i'
  csvdiff = 'cyl4_frictionless_aug_check.csv cyl4_frictionless_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=augmented_lagrange Contact/leftright/normalize_penalty=true Outputs/file_base=cyl4_frictionless_aug_out Outputs/chkfile/file_base=cyl4_frictionless_aug_check'
  rel_err = 1e-5
  abs_zero = 2e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 2D QUAD8 axisymmetric geometry (NAFEMS CGS1 contact
    patch test with mismatched nodes)."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'cyl4_mu_0_2_pen.i'
  csvdiff = 'cyl4_mu_0_2_pen_check.csv cyl4_mu_0_2_pen_check_cont_press_0001.csv cyl4_mu_0_2_pen_check_x_disp_0001.csv'
  rel_err = 1e-5
  abs_zero = 1e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D QUAD8 axisymmetric geometry (NAFEMS CGS1 contact
    patch test with mismatched nodes)."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'plane1_template1.i'
  csvdiff = 'plane1_glued_kin_check.csv plane1_glued_kin_check_cont_press_0001.csv plane1_glued_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=kinematic Outputs/file_base=plane1_glued_kin_out Outputs/chkfile/file_base=plane1_glued_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D plane geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'plane1_template1.i'
  csvdiff = 'plane1_glued_pen_check.csv plane1_glued_pen_check_cont_press_0001.csv plane1_glued_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=penalty Outputs/file_base=plane1_glued_pen_out Outputs/chkfile/file_base=plane1_glued_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact for
    2D plane geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'plane1_template1.i'
  csvdiff = 'plane1_frictionless_kin_check.csv plane1_frictionless_kin_check_cont_press_0001.csv plane1_frictionless_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic Outputs/file_base=plane1_frictionless_kin_out Outputs/chkfile/file_base=plane1_frictionless_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, kinematic
    contact for 2D plane geometry (NAFEMS CGS1 contact patch test with matched
    nodes)."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'plane1_template1.i'
  csvdiff = 'plane1_frictionless_pen_check.csv plane1_frictionless_pen_check_cont_press_0001.csv plane1_frictionless_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Outputs/file_base=plane1_frictionless_pen_out Outputs/chkfile/file_base=plane1_frictionless_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty
    contact for 2D plane geometry (NAFEMS CGS1 contact patch test with matched
    nodes)."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'plane1_template2.i'
  csvdiff = 'plane1_frictionless_aug_check.csv plane1_frictionless_aug_check_cont_press_0001.csv plane1_frictionless_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=augmented_lagrange Outputs/file_base=plane1_frictionless_aug_out Outputs/chkfile/file_base=plane1_frictionless_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 2D plane geometry (NAFEMS CGS1 contact patch test
    with matched nodes)."
[../]

[./mu_0_2_aug]
  type = 'CSVDiff'
  input = 'plane1_template2.i'
  csvdiff = 'plane1_mu_0_2_aug_check.csv plane1_mu_0_2_aug_check_cont_press_0001.csv plane1_mu_0_2_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=coulomb Contact/leftright/formulation=augmented_lagrange Contact/leftright/friction_coefficient=0.2 Contact/leftright/al_incremental_slip_tolerance=1e-6 Contact/leftright/al_frictional_force_tolerance=5e-2 Outputs/file_base=plane1_mu_0_2_aug_out Outputs/chkfile/file_base=plane1_mu_0_2_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, Augmented
    Lagrange contact for 2D plane geometry (NAFEMS CGS1 contact patch test
    with matched nodes)."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'plane1_mu_0_2_pen.i'
  csvdiff = 'plane1_mu_0_2_pen_check.csv plane1_mu_0_2_pen_check_cont_press_0001.csv plane1_mu_0_2_pen_check_x_disp_0001.csv'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D plane geometry (NAFEMS CGS1 contact patch test with matched
    nodes)."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'plane2_template1.i'
  csvdiff = 'plane2_glued_kin_check.csv plane2_glued_kin_check_cont_press_0001.csv plane2_glued_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=kinematic Outputs/file_base=plane2_glued_kin_out Outputs/chkfile/file_base=plane2_glued_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D plane geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'plane2_template1.i'
  csvdiff = 'plane2_glued_pen_check.csv plane2_glued_pen_check_cont_press_0001.csv plane2_glued_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=penalty Outputs/file_base=plane2_glued_pen_out Outputs/chkfile/file_base=plane2_glued_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact for
    2D plane geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'plane2_template1.i'
  csvdiff = 'plane2_frictionless_kin_check.csv plane2_frictionless_kin_check_cont_press_0001.csv plane2_frictionless_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic Outputs/file_base=plane2_frictionless_kin_out Outputs/chkfile/file_base=plane2_frictionless_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, kinematic
    contact for 2D plane geometry (NAFEMS CGS1 contact patch test with mismatched
    nodes)."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'plane2_template1.i'
  csvdiff = 'plane2_frictionless_pen_check.csv plane2_frictionless_pen_check_cont_press_0001.csv plane2_frictionless_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Outputs/file_base=plane2_frictionless_pen_out Outputs/chkfile/file_base=plane2_frictionless_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty
    contact for 2D plane geometry (NAFEMS CGS1 contact patch test with mismatched
    nodes)."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'plane2_template2.i'
  csvdiff = 'plane2_frictionless_aug_check.csv plane2_frictionless_aug_check_cont_press_0001.csv plane2_frictionless_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=augmented_lagrange Outputs/file_base=plane2_frictionless_aug_out Outputs/chkfile/file_base=plane2_frictionless_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 2D plane geometry (NAFEMS CGS1 contact patch test
    with mismatched nodes)."
[../]

[./mu_0_2_aug]
  type = 'CSVDiff'
  input = 'plane2_template2.i'
  csvdiff = 'plane2_mu_0_2_aug_check.csv plane2_mu_0_2_aug_check_cont_press_0001.csv plane2_mu_0_2_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=coulomb Contact/leftright/formulation=augmented_lagrange Contact/leftright/friction_coefficient=0.2 Contact/leftright/al_incremental_slip_tolerance=1e-6 Contact/leftright/al_frictional_force_tolerance=5e-2 Outputs/file_base=plane2_mu_0_2_aug_out Outputs/chkfile/file_base=plane2_mu_0_2_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, Augmented
    Lagrange contact for 2D plane geometry (NAFEMS CGS1 contact patch test
    with mismatched nodes)."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'plane2_mu_0_2_pen.i'
  csvdiff = 'plane2_mu_0_2_pen_check.csv plane2_mu_0_2_pen_check_cont_press_0001.csv plane2_mu_0_2_pen_check_x_disp_0001.csv'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D plane geometry (NAFEMS CGS1 contact patch test with mismatched
    nodes)."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'plane3_template1.i'
  csvdiff = 'plane3_glued_kin_check.csv plane3_glued_kin_check_cont_press_0001.csv plane3_glued_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=kinematic Outputs/file_base=plane3_glued_kin_out Outputs/chkfile/file_base=plane3_glued_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  petsc_version = '>=3.7.6'
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D QUAD8 plane geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'plane3_template1.i'
  csvdiff = 'plane3_glued_pen_check.csv plane3_glued_pen_check_cont_press_0001.csv plane3_glued_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=penalty Outputs/file_base=plane3_glued_pen_out Outputs/chkfile/file_base=plane3_glued_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact for
    2D QUAD8 plane geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'plane3_template1.i'
  csvdiff = 'plane3_frictionless_kin_check.csv plane3_frictionless_kin_check_cont_press_0001.csv plane3_frictionless_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic Outputs/file_base=plane3_frictionless_kin_out Outputs/chkfile/file_base=plane3_frictionless_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  petsc_version = '>=3.7.6'
  requirement = "The Contact system shall enforce frictionless, kinematic contact for
    2D QUAD8 plane geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'plane3_template1.i'
  csvdiff = 'plane3_frictionless_pen_check.csv plane3_frictionless_pen_check_cont_press_0001.csv plane3_frictionless_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Outputs/file_base=plane3_frictionless_pen_out Outputs/chkfile/file_base=plane3_frictionless_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact for
    2D QUAD8 plane geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'plane3_template2.i'
  csvdiff = 'plane3_frictionless_aug_check.csv plane3_frictionless_aug_check_cont_press_0001.csv plane3_frictionless_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=augmented_lagrange Outputs/file_base=plane3_frictionless_aug_out Outputs/chkfile/file_base=plane3_frictionless_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 2D QUAD8 plane geometry (NAFEMS CGS1 contact patch test
    with matched nodes)."
[../]

[./mu_0_2_aug]
  type = 'CSVDiff'
  input = 'plane3_template2.i'
  csvdiff = 'plane3_mu_0_2_aug_check.csv plane3_mu_0_2_aug_check_cont_press_0001.csv plane3_mu_0_2_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=coulomb Contact/leftright/formulation=augmented_lagrange Contact/leftright/friction_coefficient=0.2 Contact/leftright/al_incremental_slip_tolerance=1e-6 Contact/leftright/al_frictional_force_tolerance=5e-2 Outputs/file_base=plane3_mu_0_2_aug_out Outputs/chkfile/file_base=plane3_mu_0_2_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, Augmented
    Lagrange contact for 2D QUAD8 plane geometry (NAFEMS CGS1 contact patch test
    with matched nodes)."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'plane3_mu_0_2_pen.i'
  csvdiff = 'plane3_mu_0_2_pen_check.csv plane3_mu_0_2_pen_check_cont_press_0001.csv plane3_mu_0_2_pen_check_x_disp_0001.csv'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D QUAD8 plane geometry (NAFEMS CGS1 contact patch test with matched
    nodes)."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'plane4_template1.i'
  csvdiff = 'plane4_glued_kin_check.csv plane4_glued_kin_check_cont_press_0001.csv plane4_glued_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=kinematic Outputs/file_base=plane4_glued_kin_out Outputs/chkfile/file_base=plane4_glued_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D QUAD8 plane geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'plane4_template1.i'
  csvdiff = 'plane4_glued_pen_check.csv plane4_glued_pen_check_cont_press_0001.csv plane4_glued_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=penalty Outputs/file_base=plane4_glued_pen_out Outputs/chkfile/file_base=plane4_glued_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 2D QUAD8 plane geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'plane4_template1.i'
  csvdiff = 'plane4_frictionless_kin_check.csv plane4_frictionless_kin_check_cont_press_0001.csv plane4_frictionless_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic Outputs/file_base=plane4_frictionless_kin_out Outputs/chkfile/file_base=plane4_frictionless_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 2D QUAD8 plane geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'plane4_template1.i'
  csvdiff = 'plane4_frictionless_pen_check.csv plane4_frictionless_pen_check_cont_press_0001.csv plane4_frictionless_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Outputs/file_base=plane4_frictionless_pen_out Outputs/chkfile/file_base=plane4_frictionless_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 2D QUAD8 plane geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'plane4_template2.i'
  csvdiff = 'plane4_frictionless_aug_check.csv plane4_frictionless_aug_check_cont_press_0001.csv plane4_frictionless_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=augmented_lagrange Outputs/file_base=plane4_frictionless_aug_out Outputs/chkfile/file_base=plane4_frictionless_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 2D QUAD8 plane geometry (NAFEMS CGS1 contact patch test
    with mismatched nodes)."
[../]

[./mu_0_2_aug]
  type = 'CSVDiff'
  input = 'plane4_template2.i'
  csvdiff = 'plane4_mu_0_2_aug_check.csv plane4_mu_0_2_aug_check_cont_press_0001.csv plane4_mu_0_2_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=coulomb Contact/leftright/formulation=augmented_lagrange Contact/leftright/friction_coefficient=0.2 Contact/leftright/al_incremental_slip_tolerance=1e-6 Contact/leftright/al_frictional_force_tolerance=5e-2 Outputs/file_base=plane4_mu_0_2_aug_out Outputs/chkfile/file_base=plane4_mu_0_2_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, Augmented
    Lagrange contact for 2D QUAD8 plane geometry (NAFEMS CGS1 contact patch test
    with matched nodes)."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'plane4_mu_0_2_pen.i'
  csvdiff = 'plane4_mu_0_2_pen_check.csv plane4_mu_0_2_pen_check_cont_press_0001.csv plane4_mu_0_2_pen_check_x_disp_0001.csv'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D QUAD8 plane geometry (NAFEMS CGS1 contact patch test with
    mismatched nodes)."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'ring1_template1.i'
  csvdiff = 'ring1_glued_kin_check.csv ring1_glued_kin_check_cont_press_0001.csv ring1_glued_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=kinematic Outputs/file_base=ring1_glued_kin_out Outputs/chkfile/file_base=ring1_glued_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D axisymmetric ring geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'ring1_template1.i'
  csvdiff = 'ring1_glued_pen_check.csv ring1_glued_pen_check_cont_press_0001.csv ring1_glued_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=penalty Outputs/file_base=ring1_glued_pen_out Outputs/chkfile/file_base=ring1_glued_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 2D axisymmetric ring geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'ring1_template1.i'
  csvdiff = 'ring1_frictionless_kin_check.csv ring1_frictionless_kin_check_cont_press_0001.csv ring1_frictionless_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic Outputs/file_base=ring1_frictionless_kin_out Outputs/chkfile/file_base=ring1_frictionless_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frctionless, kinematic contact
    for 2D axisymmetric ring geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'ring1_template1.i'
  csvdiff = 'ring1_frictionless_pen_check.csv ring1_frictionless_pen_check_cont_press_0001.csv ring1_frictionless_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Outputs/file_base=ring1_frictionless_pen_out Outputs/chkfile/file_base=ring1_frictionless_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frctionless, penalty contact
    for 2D axisymmetric ring geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'ring1_template2.i'
  csvdiff = 'ring1_frictionless_aug_check.csv ring1_frictionless_aug_check_cont_press_0001.csv ring1_frictionless_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=augmented_lagrange Outputs/file_base=ring1_frictionless_aug_out Outputs/chkfile/file_base=ring1_frictionless_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 2D axisymmetric ring geometry (NAFEMS CGS1 contact
    patch test with matched nodes)."
[../]

[./mu_0_2_aug]
  type = 'CSVDiff'
  input = 'ring1_template2.i'
  csvdiff = 'ring1_mu_0_2_aug_check.csv ring1_mu_0_2_aug_check_cont_press_0001.csv ring1_mu_0_2_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=coulomb Contact/leftright/formulation=augmented_lagrange Contact/leftright/friction_coefficient=0.2 Contact/leftright/al_incremental_slip_tolerance=1e-6 Contact/leftright/al_frictional_force_tolerance=5e-2 Outputs/file_base=ring1_mu_0_2_aug_out Outputs/chkfile/file_base=ring1_mu_0_2_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, Augmented
    Lagrange contact for 2D axisymmetric ring geometry (NAFEMS CGS1 contact
    patch test with matched nodes)."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'ring1_mu_0_2_pen.i'
  csvdiff = 'ring1_mu_0_2_pen_check.csv ring1_mu_0_2_pen_check_cont_press_0001.csv ring1_mu_0_2_pen_check_x_disp_0001.csv'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D axisymmetric ring geometry (NAFEMS CGS1 contact
    patch test with matched nodes)."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'ring2_template1.i'
  csvdiff = 'ring2_glued_kin_check.csv ring2_glued_kin_check_cont_press_0001.csv ring2_glued_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=kinematic Outputs/file_base=ring2_glued_kin_out Outputs/chkfile/file_base=ring2_glued_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D axisymmetric ring geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'ring2_template1.i'
  csvdiff = 'ring2_glued_pen_check.csv ring2_glued_pen_check_cont_press_0001.csv ring2_glued_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=penalty Outputs/file_base=ring2_glued_pen_out Outputs/chkfile/file_base=ring2_glued_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 2D axisymmetric ring geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'ring2_template1.i'
  csvdiff = 'ring2_frictionless_kin_check.csv ring2_frictionless_kin_check_cont_press_0001.csv ring2_frictionless_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic Outputs/file_base=ring2_frictionless_kin_out Outputs/chkfile/file_base=ring2_frictionless_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 2D axisymmetric ring geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'ring2_template1.i'
  csvdiff = 'ring2_frictionless_pen_check.csv ring2_frictionless_pen_check_cont_press_0001.csv ring2_frictionless_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Outputs/file_base=ring2_frictionless_pen_out Outputs/chkfile/file_base=ring2_frictionless_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 2D axisymmetric ring geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'ring2_template2.i'
  csvdiff = 'ring2_frictionless_aug_check.csv ring2_frictionless_aug_check_cont_press_0001.csv ring2_frictionless_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=augmented_lagrange Outputs/file_base=ring2_frictionless_aug_out Outputs/chkfile/file_base=ring2_frictionless_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 2D axisymmetric ring geometry (NAFEMS CGS1 contact
    patch test with mismatched nodes)."
[../]

[./mu_0_2_aug]
  type = 'CSVDiff'
  input = 'ring2_template2.i'
  csvdiff = 'ring2_mu_0_2_aug_check.csv ring2_mu_0_2_aug_check_cont_press_0001.csv ring2_mu_0_2_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=coulomb Contact/leftright/formulation=augmented_lagrange Contact/leftright/friction_coefficient=0.2 Contact/leftright/al_incremental_slip_tolerance=1e-6 Contact/leftright/al_frictional_force_tolerance=5e-2 Outputs/file_base=ring2_mu_0_2_aug_out Outputs/chkfile/file_base=ring2_mu_0_2_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, Augmented
    Lagrange contact for 2D axisymmetric ring geometry (NAFEMS CGS1 contact
    patch test with mismatched nodes)."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'ring2_mu_0_2_pen.i'
  csvdiff = 'ring2_mu_0_2_pen_check.csv ring2_mu_0_2_pen_check_cont_press_0001.csv ring2_mu_0_2_pen_check_x_disp_0001.csv'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D axisymmetric ring geometry (NAFEMS CGS1 contact
    patch test with mismatched nodes)."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'ring3_template1.i'
  csvdiff = 'ring3_glued_kin_check.csv ring3_glued_kin_check_cont_press_0001.csv ring3_glued_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=kinematic Outputs/file_base=ring3_glued_kin_out Outputs/chkfile/file_base=ring3_glued_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D QUAD8 axisymmetric ring geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'ring3_template1.i'
  csvdiff = 'ring3_glued_pen_check.csv ring3_glued_pen_check_cont_press_0001.csv ring3_glued_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=penalty Outputs/file_base=ring3_glued_pen_out Outputs/chkfile/file_base=ring3_glued_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 2D QUAD8 axisymmetric ring geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'ring3_template1.i'
  csvdiff = 'ring3_frictionless_kin_check.csv ring3_frictionless_kin_check_cont_press_0001.csv ring3_frictionless_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic Outputs/file_base=ring3_frictionless_kin_out Outputs/chkfile/file_base=ring3_frictionless_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  petsc_version = '>=3.7.6'
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 2D QUAD8 axisymmetric ring geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'ring3_template1.i'
  csvdiff = 'ring3_frictionless_pen_check.csv ring3_frictionless_pen_check_cont_press_0001.csv ring3_frictionless_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Outputs/file_base=ring3_frictionless_pen_out Outputs/chkfile/file_base=ring3_frictionless_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 2D QUAD8 axisymmetric ring geometry (NAFEMS CGS1 contact patch test with matched nodes)."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'ring3_template2.i'
  csvdiff = 'ring3_frictionless_aug_check.csv ring3_frictionless_aug_check_cont_press_0001.csv ring3_frictionless_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=augmented_lagrange Outputs/file_base=ring3_frictionless_aug_out Outputs/chkfile/file_base=ring3_frictionless_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 2D QUAD8 axisymmetric ring geometry (NAFEMS CGS1
    contact patch test with matched nodes)."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'ring3_mu_0_2_pen.i'
  csvdiff = 'ring3_mu_0_2_pen_check.csv ring3_mu_0_2_pen_check_cont_press_0001.csv ring3_mu_0_2_pen_check_x_disp_0001.csv'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D QUAD8 axisymmetric ring geometry (NAFEMS CGS1
    contact patch test with matched nodes)."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'ring4_template1.i'
  csvdiff = 'ring4_glued_kin_check.csv ring4_glued_kin_check_cont_press_0001.csv ring4_glued_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=kinematic Outputs/file_base=ring4_glued_kin_out Outputs/chkfile/file_base=ring4_glued_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D QUAD8 axisymmetric ring geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'ring4_template1.i'
  csvdiff = 'ring4_glued_pen_check.csv ring4_glued_pen_check_cont_press_0001.csv ring4_glued_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=glued Contact/leftright/formulation=penalty Outputs/file_base=ring4_glued_pen_out Outputs/chkfile/file_base=ring4_glued_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 2D QUAD8 axisymmetric ring geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'ring4_template1.i'
  csvdiff = 'ring4_frictionless_kin_check.csv ring4_frictionless_kin_check_cont_press_0001.csv ring4_frictionless_kin_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic Outputs/file_base=ring4_frictionless_kin_out Outputs/chkfile/file_base=ring4_frictionless_kin_check'
  rel_err = 1e-5
  abs_zero = 1e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 2D QUAD8 axisymmetric ring geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'ring4_template1.i'
  csvdiff = 'ring4_frictionless_pen_check.csv ring4_frictionless_pen_check_cont_press_0001.csv ring4_frictionless_pen_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Outputs/file_base=ring4_frictionless_pen_out Outputs/chkfile/file_base=ring4_frictionless_pen_check'
  rel_err = 1e-5
  abs_zero = 1e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 2D QUAD8 axisymmetric ring geometry (NAFEMS CGS1 contact patch test with mismatched nodes)."
[../]

[./frictionless_aug]
  type = 'CSVDiff'
  input = 'ring4_template2.i'
  csvdiff = 'ring4_frictionless_aug_check.csv ring4_frictionless_aug_check_cont_press_0001.csv ring4_frictionless_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=frictionless Contact/leftright/formulation=augmented_lagrange Outputs/file_base=ring4_frictionless_aug_out Outputs/chkfile/file_base=ring4_frictionless_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictionless, Augmented
    Lagrange contact for 2D QUAD8 axisymmetric ring geometry (NAFEMS CGS1
    contact patch test with mismatched nodes)."
[../]

[./mu_0_2_aug]
  type = 'CSVDiff'
  input = 'ring4_template2.i'
  csvdiff = 'ring4_mu_0_2_aug_check.csv ring4_mu_0_2_aug_check_cont_press_0001.csv ring4_mu_0_2_aug_check_x_disp_0001.csv'
  cli_args = 'Contact/leftright/model=coulomb Contact/leftright/formulation=augmented_lagrange Contact/leftright/friction_coefficient=0.2 Contact/leftright/al_incremental_slip_tolerance=1e-6 Contact/leftright/al_frictional_force_tolerance=5e-2 Outputs/file_base=ring4_mu_0_2_aug_out Outputs/chkfile/file_base=ring4_mu_0_2_aug_check'
  rel_err = 1e-5
  abs_zero = 1e-8
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, Augmented
    Lagrange contact for 2D QUAD8 axisymmetric ring geometry (NAFEMS CGS1
    contact patch test with mismatched nodes)."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'ring4_mu_0_2_pen.i'
  csvdiff = 'ring4_mu_0_2_pen_check.csv ring4_mu_0_2_pen_check_cont_press_0001.csv ring4_mu_0_2_pen_check_x_disp_0001.csv'
  rel_err = 1e-5
  abs_zero = 1e-6
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  requirement = "The Contact system shall enforce frictional, penalty
    contact for 2D QUAD8 axisymmetric ring geometry (NAFEMS CGS1
    contact patch test with mismatched nodes)."
[../]

[./glued_kin]
  type = 'CSVDiff'
  input = 'single_point_2d.i'
  csvdiff = 'single_point_2d_out_glued_kin.csv'
  cli_args = 'Executioner/end_time=0.002 Contact/leftright/model=glued Contact/leftright/formulation=kinematic'
  rel_err = 5e-4
  abs_zero = 3e-5
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  petsc_version = '>=3.7.6'
  compiler = '!INTEL'
  allow_warnings = true
  requirement = "The Contact system shall enforce glued, kinematic contact
    for 2D plane strain single point contact model."
[../]

[./glued_pen]
  type = 'CSVDiff'
  input = 'single_point_2d.i'
  csvdiff = 'single_point_2d_out_glued_pen.csv'
  cli_args = 'Executioner/end_time=0.002 Executioner/l_tol=1e-4 Executioner/nl_rel_tol=1e-8 Executioner/nl_abs_tol=1e-10
    Contact/leftright/formulation=penalty Contact/leftright/penalty=2.0e12 Outputs/file_base=single_point_2d_out_glued_pen'
  rel_err = 5e-4
  abs_zero = 3e-5
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  allow_warnings = true
  requirement = "The Contact system shall enforce glued, penalty contact
    for 2D plane strain single point contact model."
[../]

[./frictionless_kin]
  type = 'CSVDiff'
  input = 'single_point_2d.i'
  csvdiff = 'single_point_2d_out_frictionless_kin.csv'
  cli_args = 'Executioner/end_time=0.002 Executioner/nl_rel_tol=1e-8 Executioner/nl_abs_tol=1e-10 Contact/leftright/model=frictionless Contact/leftright/formulation=kinematic Outputs/file_base=single_point_2d_out_frictionless_kin'
  rel_err = 1e-5
  abs_zero = 1e-5
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  allow_warnings = true
  requirement = "The Contact system shall enforce frictionless, kinematic contact
    for 2D plane strain single point contact model."
[../]

[./frictionless_pen]
  type = 'CSVDiff'
  input = 'single_point_2d.i'
  csvdiff = 'single_point_2d_out_frictionless_pen.csv'
  cli_args = 'Executioner/end_time=0.002 Executioner/nl_rel_tol=1e-8 Executioner/nl_abs_tol=1e-10 Contact/leftright/model=frictionless Contact/leftright/formulation=penalty Contact/leftright/penalty=5.0e12 Outputs/file_base=single_point_2d_out_frictionless_pen'
  rel_err = 1e-4
  abs_zero = 1e-5
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  allow_warnings = true
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 2D plane strain single point contact model."
[../]

[./frictionless_pen_contact_line_search]
  type = 'CSVDiff'
  input = 'single_point_2d_contact_line_search.i'
  csvdiff = 'single_point_2d_contact_line_search_out.csv'
  allow_warnings = true
  abs_zero = 1e-8
  rel_err = 5e-4
  mesh_mode = REPLICATED
  # The linesearch used by this test requires PETSc >= 3.6.0.
  petsc_version = '>=3.6.0'
  requirement = "The Contact system shall enforce frictionless, penalty contact
    for 2D plane strain single point contact model using the contact line
    search solver options."
  compiler = 'GCC CLANG'
[../]

[./mu_0_0_kin]
  type = 'CSVDiff'
  input = 'single_point_2d.i'
  csvdiff = 'single_point_2d_out_fric_0_0_kin.csv'
  cli_args = 'Executioner/end_time=0.002 Executioner/nl_rel_tol=1e-8 Executioner/nl_abs_tol=1e-10 Contact/leftright/model=coulomb Contact/leftright/formulation=kinematic Contact/leftright/friction_coefficient=0.0 Outputs/file_base=single_point_2d_out_fric_0_0_kin'
  rel_err = 1e-5
  abs_zero = 1e-5
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  allow_warnings = true
  requirement = "The Contact system shall enforce frictional, kinematic contact
    for 2D plane strain single point contact model."
[../]

[./mu_0_0_pen]
  type = 'CSVDiff'
  input = 'single_point_2d.i'
  csvdiff = 'single_point_2d_out_fric_0_0_pen.csv'
  cli_args = 'Executioner/end_time=0.002 Executioner/nl_rel_tol=1e-8 Executioner/nl_abs_tol=1e-10 Contact/leftright/model=coulomb Contact/leftright/formulation=penalty Contact/leftright/friction_coefficient=0.0 Contact/leftright/penalty=8.0e12 Outputs/file_base=single_point_2d_out_fric_0_0_pen'
  rel_err = 1e-4
  abs_zero = 1e-5
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  allow_warnings = true
  requirement = "The Contact system shall enforce frictional, penalty contact
    for 2D plane strain single point contact model."
[../]

[./mu_0_2_kin]
  type = 'CSVDiff'
  input = 'single_point_2d_frictional.i'
  csvdiff = 'single_point_2d_out_frictional_0_2_kin.csv'
  cli_args = 'Executioner/end_time=0.002 Executioner/nl_rel_tol=1e-8 Executioner/nl_abs_tol=1e-10 Contact/leftright/friction_coefficient=0.2 Outputs/file_base=single_point_2d_out_frictional_0_2_kin'
  rel_err = 1e-4
  abs_zero = 4e-5
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  allow_warnings = true
  requirement = "The Contact system shall enforce frictional, kinematic contact
    for 2D plane strain single point contact model with a non-zero friction
    coefficient."
[../]

[./mu_0_2_pen]
  type = 'CSVDiff'
  input = 'single_point_2d_frictional.i'
  csvdiff = 'single_point_2d_out_frictional_0_2_pen.csv'
  cli_args = 'Executioner/end_time=0.002 Executioner/nl_rel_tol=1e-8 Executioner/nl_abs_tol=1e-10 Contact/leftright/formulation=penalty
    Contact/leftright/model=coulomb Contact/leftright/friction_coefficient=0.2 Contact/leftright/penalty=1.0e12 Outputs/file_base=single_point_2d_out_frictional_0_2_pen'
  rel_err = 1e-4
  abs_zero = 3e-5
  max_parallel = 1
  mesh_mode = REPLICATED
  superlu = true
  allow_warnings = true
  requirement = "The Contact system shall enforce frictional, penalty contact
    for 2D plane strain single point contact model with a non-zero friction
    coefficient."
[../]