Contact Requirements Traceability Matrix

[frictionless-mortar-3d]
  type = 'Exodiff'
  input = 'frictionless-mortar-3d.i'
  min_ad_size = 150
  exodiff = 'frictionless-mortar-3d_out.e'
  requirement = 'The system shall solve a 3D frictionless bouncing block problem with mortar constraint'
  method = '!dbg'
  valgrind = 'none'
[]

[frictionless-mortar-3d_pg]
  type = 'CSVDiff'
  input = 'frictionless-mortar-3d_pg.i'
  min_ad_size = 150
  csvdiff = 'frictionless-mortar-3d_pg_out.csv'
  requirement = 'The system shall solve a 3D frictionless bouncing block problem with mortar constraint using Petrov-Galerkin approach.'
  method = '!dbg'
  valgrind = 'none'
[]

[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.'
[]

[frictionless-mortar-3d-action]
  type = 'CSVDiff'
  input = 'frictionless-mortar-3d-action.i'
  min_ad_size = 150
  csvdiff = 'frictionless-mortar-3d-action_out_contact-pressure_0002.csv'
  mesh_mode = 'REPLICATED' # see issue 24255
  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
  method = '!dbg'
  valgrind = 'none'
[]

[frictionless-mortar-3d-pg-action]
  prereq = 'frictionless-mortar-3d_pg'
  type = 'CSVDiff'
  input = 'frictionless-mortar-3d-action.i'
  cli_args = 'Contact/mortar/use_petrov_galerkin=true Outputs/file_base="frictionless-mortar-3d_pg_out" Postprocessors/active="contact lambda"'
  min_ad_size = 100
  csvdiff = 'frictionless-mortar-3d_pg_out.csv'
  requirement = 'The system shall solve a 3D frictionless bouncing block problem with mortar constraint using Petrov-Galerkin approach in contact action.'
[]

[frictionless-mortar-3d-action-correct-edge-dropping]
  type = 'CSVDiff'
  input = 'frictionless-mortar-3d-action.i'
  min_ad_size = 150
  csvdiff = 'frictionless-mortar-3d-action-correct-edge-dropping_out_contact-pressure_0002.csv'
  mesh_mode = 'REPLICATED' # see issue 24255
  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
  method = '!dbg'
  valgrind = 'none'
[]

[frictionless-mortar-3d-debug-mesh]
  type = 'Exodiff'
  input = 'frictionless-mortar-3d.i'
  min_ad_size = 150
  exodiff = 'mortar_segment_mesh.e'
  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 UserObjects/weighted_gap_uo/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
[]

[frictional-mortar-3d]
  type = 'Exodiff'
  input = 'frictional-mortar-3d.i'
  min_ad_size = 150
  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
  method = '!dbg'
  valgrind = 'none'
[]

[frictional-mortar-3d-pg]
  type = 'CSVDiff'
  input = 'frictional-mortar-3d_pg.i'
  min_ad_size = 150
  csvdiff = 'frictional-mortar-3d_pg_out_contact-pressure_0001.csv frictional-mortar-3d_pg_out_frictional-pressure-3d_0001.csv frictional-mortar-3d_pg_out_frictional-pressure_0001.csv'
  requirement = 'The system shall solve a 3D frictional bouncing block problem with mortar constraint using the Petrov-Galerkin approach.'
  rel_err = 1.0e-5
  abs_zero = 1.0e-5
  method = '!dbg'
  valgrind = 'none'
[]

[frictional-mortar-3d-function]
  type = 'Exodiff'
  input = 'frictional-mortar-3d-function.i'
  min_ad_size = 150
  exodiff = 'frictional-mortar-3d-function_out.e'
  requirement = 'The system shall solve a 3D frictional block problem with mortar constraint using a simple frictional model dependent on relative velocity and contact pressure.'
  rel_err = 1.0e-5
  abs_zero = 1.0e-5
  method = '!dbg'
  valgrind = 'none'
[]

[frictional-mortar-3d-action]
  type = 'CSVDiff'
  input = 'frictional-mortar-3d-action.i'
  min_ad_size = 150
  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
  method = '!dbg'
  valgrind = 'none'
[]

[frictional-mortar-3d-pg-action]
  prereq = 'frictional-mortar-3d-pg'
  type = 'CSVDiff'
  input = 'frictional-mortar-3d-action.i'
  cli_args = 'Contact/mortar/use_petrov_galerkin=true Outputs/file_base=frictional-mortar-3d_pg_out'
  min_ad_size = 150
  csvdiff = 'frictional-mortar-3d_pg_out_contact-pressure_0001.csv frictional-mortar-3d_pg_out_frictional-pressure-3d_0001.csv frictional-mortar-3d_pg_out_frictional-pressure_0001.csv'
  requirement = 'The system shall solve a 3D frictional bouncing block problem with mortar constraint using the Petrov-Galerkin approach in contact action.'
  rel_err = 1.0e-5
  abs_zero = 1.0e-5
[]

[half_sphere_nodal_geometry]
  type = 'Exodiff'
  exodiff = 'nodal_geometry_only.e'
  input = 'half_sphere_nodal_geometry.i'
  min_ad_size = 100
  valgrind = 'heavy'
  requirement = 'The system shall generate consistent mortar nodal geometry (normal and tangent vectors) on a spherical surface.'
  rel_err = 1.0e-5
  abs_zero = 1.0e-5
  allow_warnings = true
  max_parallel = 1
[]

[frictional-mortar-3d-penalty]
  type = 'Exodiff'
  input = 'frictional-mortar-3d-penalty.i'
  min_ad_size = 150
  exodiff = 'frictional-mortar-3d-penalty_out.e'
  requirement = 'The system shall solve a 3D frictional problem with mortar constraint using a penalty approach'
  rel_err = 1.0e-5
  abs_zero = 1.0e-5
[]

[frictionless-mortar-3d-penalty]
  type = 'Exodiff'
  input = 'frictionless-mortar-3d-penalty.i'
  min_ad_size = 200
  exodiff = 'frictionless-mortar-3d-penalty_out.e'
  requirement = 'The system shall solve a 3D frictionless problem with mortar constraint using a penalty approach'
  rel_err = 1.0e-5
  abs_zero = 1.0e-5
[]

[frictional-mortar-3d-al]
  type = 'Exodiff'
  input = 'frictional-mortar-3d-al.i'
  exodiff = 'frictional-mortar-3d-al_out.e'
  requirement = 'The system shall be able to solve a 3D frictional problem with mortar constraint using a augmented Lagrange approach and converging to prescribed tolerance.'
  rel_err = 1.0e-5
  abs_zero = 1.0e-5
  valgrind = none
  method = 'opt'
  heavy = true
  min_ad_size = 200
  max_time = 1000
[]

[./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'
  # 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
  rel_err = 6e-5
  abs_zero = 5e-9
[]

[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=''"
[]

[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=''"
  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=''"
[]

[penalty]
  type = 'Exodiff'
  input = frictionless-penalty-weighted-gap.i
  exodiff = 'frictionless-penalty-weighted-gap_out.e'
  detail = 'using a first order basis for displacements and penalty multiplication times the negative gap distance to form the contact force.'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
[]

[penalty_action]
  type = 'Exodiff'
  input = frictionless-penalty-weighted-gap-action.i
  exodiff = 'frictionless-penalty-weighted-gap-action_out.e'
  detail = 'using a first order basis for displacements and penalty multiplication times the negative gap distance to form the contact force via the contact action.'
  map = true
  cli_args = "Mesh/mesh_file/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active='' Outputs/hide='nodal_area contact_pressure penetration'"
[]

[very_parallel_penalty]
  requirement = 'The system shall be able to communicate semilocal weighted gaps back to any process that contributed to computing said weighted gaps.'
  issues = '#21329'
  type = 'Exodiff'
  input = frictionless-penalty-weighted-gap.i
  exodiff = 'parallel-frictionless-penalty-weighted-gap_out.e'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active='' Outputs/file_base=parallel-frictionless-penalty-weighted-gap_out"
  min_parallel = 13
[]

[very_parallel_frictional_penalty]
  requirement = 'The system shall be able to communicate semilocal weighted gaps and velocities back to any process that contributed to computing said weighted quantities. We check against the serial numerical results.'
  issues = '#21329'
  type = 'Exodiff'
  input = frictional-penalty-weighted-vel.i
  exodiff = 'parallel-frictional-penalty-weighted-vel_out.e'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active='' Outputs/file_base=parallel-frictional-penalty-weighted-vel_out"
  min_parallel = 13
[]

[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=''"
  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=''"
  abs_zero = 5e-9
[]

[split_secondary]
  type = 'Exodiff'
  input = 'variational-frictional.i'
  exodiff = 'variational-frictional_out.e'
  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'
[]

[split_secondary_action]
  type = 'Exodiff'
  input = 'variational-frictional-action.i'
  exodiff = 'variational-frictional-action_exodus.e'
  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 while using the contact action to build the set of constraints, user objects, and application of generalized forces.'
  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'
  min_ad_size = '100'
[]

[./catch_release]
  type = 'Exodiff'
  input = 'catch_release.i'
  exodiff = 'catch_release_out.e'
  custom_cmp = 'catch_release.exodiff'
  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"
[../]

[only_czm]
  type = Exodiff
  input = bilinear_mixed_mortar_only_czm.i
  exodiff = bilinear_mixed_mortar_only_czm_out.e
  requirement = 'The system shall be able to apply cohesive zone modeling with bilinear mixed mode traction to a simple model.'
  abs_zero = 1.0e-6
  map = false
[]

[czm_contact]
  type = Exodiff
  input = bilinear_mixed_mortar_only_czm.i
  exodiff = bilinear_mixed_mortar_only_czm_contact_out.e
  requirement = 'The system shall be able to apply cohesive zone modeling with bilinear mixed mode traction and mechanical contact to a simple model.'
  cli_args = "UserObjects/czm_uo/penalty=1e6 UserObjects/czm_uo/penalty_friction=0.0 Outputs/file_base=bilinear_mixed_mortar_only_czm_contact_out"
  abs_zero = 1.0e-6
  map = false
[]

[mortar_czm]
  type = Exodiff
  input = mortar_czm.i
  exodiff = mortar_czm_out.e
  requirement = 'The system shall be able to apply reproduce the setup of the interface cohesive zone modeling with a mortar formulation.'
  map = false
[]

[mortar_czm_analysis]
  type = Exodiff
  input = mortar_czm_analysis.i
  exodiff = mortar_czm_analysis_out.e
  requirement = 'The system shall be able to apply reproduce the setup of the interface cohesive zone modeling with a mortar formulation and output analysis information to the Exodus mesh.'
  map = false
[]

[bilinear_mixed]
  type = CSVDiff
  input = bilinear_mixed.i
  csvdiff = bilinear_mixed_out.csv
  requirement = 'The system shall be able to apply a cohesive zone model based on interface objects.'
  rel_err = 1.0e-5
  abs_zero = 0.001
  heavy = true
  method = 'opt'
[]

[bilinear_mixed_compare]
  type = CSVDiff
  input = bilinear_mixed_compare.i
  csvdiff = bilinear_mixed_compare_out.csv
  requirement = 'The system shall be able to apply a mortar cohesive zone model and yield results similar to those of interface objects.'
  # This test compares different formulations
  rel_err = 0.035
  abs_zero = 0.001
  heavy = true
  method = 'opt'
[]

[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
[]

[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
[]

[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
[]

[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
[]

[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
[]

[block_penalty]
  type = 'CSVDiff'
  input = 'first_test.i'
  csvdiff = 'first_test_out.csv'
  abs_zero = 1.0e-4
  rel_err = 5.0e-5
  allow_warnings = true
  requirement = 'The system shall be able to solve a simple few-element normal contact problem using explicit dynamics.'
[]

[test_balance]
  type = 'Exodiff'
  input = 'test_balance.i'
  exodiff = 'test_balance_out.e'
  abs_zero = 1.0e-4
  allow_warnings = true
  method = '!dbg'
  requirement = 'The system shall be able to solve a simple few-element normal contact problem using explicit dynamics solving uncoupled, local equations of momentum balance.'
[]

[test_balance_overwrite]
  type = 'CSVDiff'
  input = 'test_balance.i'
  csvdiff = 'test_balance_overwrite_out.csv'
  cli_args = 'ExplicitDynamicsContact/my_contact/overwrite_current_solution=true Outputs/file_base=test_balance_overwrite_out'
  abs_zero = 1.0e-4
  allow_warnings = true
  method = '!dbg'
  requirement = 'The system shall be able to solve a simple few-element normal contact problem using explicit dynamics solving uncoupled, local equations of momentum balance  and overwrite boundary variables after applying the time stepper scheme.'
[]

[test_balance_short]
  type = 'Exodiff'
  input = 'test_balance.i'
  exodiff = 'test_balance_short_out.e'
  abs_zero = 1.0e-4
  allow_warnings = true
  cli_args = 'Outputs/file_base=test_balance_short_out Executioner/end_time=-0.009'
  requirement = 'The system shall be able to solve a simple few-element normal contact problem using explicit dynamics solving uncoupled, local equations of momentum balance in debug mode.'
[]

[test_balance_overwrite_short]
  type = 'CSVDiff'
  input = 'test_balance.i'
  csvdiff = 'test_balance_overwrite_short_out.csv'
  cli_args = 'ExplicitDynamicsContact/my_contact/overwrite_current_solution=true Outputs/file_base=test_balance_overwrite_short_out Executioner/end_time=-0.009'
  abs_zero = 1.0e-4
  allow_warnings = true
  requirement = 'The system shall be able to solve a simple few-element normal contact problem using explicit dynamics solving uncoupled, local equations of momentum balance  and overwrite boundary variables after applying the time stepper scheme in debug mode.'
[]

[settlement]
  type = 'CSVDiff'
  input = 'settlement.i'
  csvdiff = 'settlement_out.csv'
  abs_zero = 1.0e-4
  allow_warnings = true
  requirement = 'The system shall be able to solve a simple few-element normal contact problem using explicit dynamics solving uncoupled, local equations of momentum balance during an impact-settling under gravity acceleration.'
  heavy = true
[]

[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.'
[]

[mortar-FS]
  type = 'Exodiff'
  input = 'frictionless_mortar_FS.i'
  exodiff = 'frictionless_mortar_FS_out.e'
  abs_zero = 1e-8
  requirement = 'The system shall allow split preconditioning of frictionless mortar contact.'
  petsc_version = '>= 3.16.6' # early versions of petsc can't handle change in sparsity pattern in field split
[]

[mortar-FS-frictional]
  type = 'Exodiff'
  input = 'frictional_mortar_FS.i'
  exodiff = 'frictional_mortar_FS_out.e'
  abs_zero = 1e-8
  requirement = 'The system shall allow split preconditioning of frictional mortar contact.'
  petsc_version = '>= 3.16.6' # early versions of petsc can't handle change in sparsity pattern in field split
[]

[./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."
[../]

[normal_pressure]
  type = RunException
  input = aux_errors.i
  cli_args = 'AuxKernels/gap/contact_quantity=normal_pressure'
  detail = 'normal_pressure'
  expect_err = "The 'NORMAL_PRESSURE' quantity is only provided by a 'PenaltyWeightedGapUserObject' or derived object\."
[]

[accumulated_slip_one]
  type = RunException
  input = aux_errors.i
  cli_args = 'AuxKernels/gap/contact_quantity=accumulated_slip_one'
  detail = 'accumulated_slip_one'
  expect_err = "The 'ACCUMULATED_SLIP_ONE' quantity is only provided by a 'PenaltyFrictionUserObject' or derived object\."
[]

[tangential_pressure_one]
  type = RunException
  input = aux_errors.i
  cli_args = 'AuxKernels/gap/contact_quantity=tangential_pressure_one'
  detail = 'tangential_pressure_one'
  expect_err = "The 'TANGENTIAL_PRESSURE_ONE' quantity is only provided by a 'PenaltyFrictionUserObject' or derived object\."
[]

[tangential_velocity_one]
  type = RunException
  input = aux_errors.i
  cli_args = 'AuxKernels/gap/contact_quantity=tangential_velocity_one'
  detail = 'tangential_velocity_one'
  expect_err = "The 'TANGENTIAL_VELOCITY_ONE' quantity is only provided by a 'WeightedVelocitiesUserObject' or derived object\."
[]

[accumulated_slip_two]
  type = RunException
  input = aux_errors.i
  cli_args = 'AuxKernels/gap/contact_quantity=accumulated_slip_two'
  detail = 'accumulated_slip_two'
  expect_err = "The 'ACCUMULATED_SLIP_TWO' quantity is only provided by a 'PenaltyFrictionUserObject' or derived object\."
[]

[tangential_pressure_two]
  type = RunException
  input = aux_errors.i
  cli_args = 'AuxKernels/gap/contact_quantity=tangential_pressure_two'
  detail = 'tangential_pressure_two'
  expect_err = "The 'TANGENTIAL_PRESSURE_TWO' quantity is only provided by a 'PenaltyFrictionUserObject' or derived object\."
[]

[tangential_velocity_two]
  type = RunException
  input = aux_errors.i
  cli_args = 'AuxKernels/gap/contact_quantity=tangential_velocity_two'
  detail = 'tangential_velocity_two'
  expect_err = "The 'TANGENTIAL_VELOCITY_TWO' quantity is only provided by a 'WeightedVelocitiesUserObject' or derived object\."
[]

[normal_gap]
  type = RunException
  input = aux_errors.i
  cli_args = 'AuxKernels/gap/contact_quantity=normal_gap'
  detail = 'normal_gap'
  expect_err = "The 'NORMAL_GAP' quantity is only provided by a 'WeightedGapUserObject' or derived object\."
[]

[normal_lm]
  type = RunException
  input = aux_errors.i
  cli_args = 'AuxKernels/gap/contact_quantity=normal_lm'
  detail = 'normal_lm'
  expect_err = "The 'NORMAL_LM' quantity is only provided by a 'PenaltyWeightedGapUserObject' or derived object\."
[]

[delta_tangential_lm_one]
  type = RunException
  input = aux_errors.i
  cli_args = 'AuxKernels/gap/contact_quantity=delta_tangential_lm_one'
  detail = 'delta_tangential_lm_one'
  expect_err = "The 'DELTA_TANGENTIAL_LM_ONE' quantity is only provided by a 'PenaltyFrictionUserObject' or derived object\."
[]

[delta_tangential_lm_two]
  type = RunException
  input = aux_errors.i
  cli_args = 'AuxKernels/gap/contact_quantity=delta_tangential_lm_two'
  detail = 'delta_tangential_lm_two'
  expect_err = "The 'DELTA_TANGENTIAL_LM_TWO' quantity is only provided by a 'PenaltyFrictionUserObject' or derived object\."
[]

[active_set]
  type = RunException
  input = aux_errors.i
  cli_args = 'AuxKernels/gap/contact_quantity=active_set'
  detail = 'active_set'
  expect_err = "The 'ACTIVE_SET' quantity is only provided by a 'PenaltyWeightedGapUserObject' or derived object\."
[]

[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'
  requirement = 'The contact module shall be able to 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'
  valgrind = 'none'
  requirement = 'The contact module shall be able to compute nodal wear depth values in accordance with Archard equation via a mortar auxiliary kernel.'
  method = '!dbg'
[]

[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 = 5.0e-3
  min_parallel = 1
  max_parallel = 4
  min_ad_size = '50'
  heavy = true
  method = 'opt'
  valgrind = 'none'
  requirement = 'The contact module shall be able to 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'
  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'
  requirement = 'The contact module shall be able to 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'
  cli_args = 'Executioner/end_time=0.2 Outputs/file_base=block-dynamics-aux-fretting-wear-test-quick_out'
  requirement = 'The contact module shall be able to 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-aux-fretting-wear-test-quick-multiple-projections]
  type = 'Exodiff'
  exodiff = 'block-dynamics-aux-fretting-wear-test-quick-multiple-projections_out.e'
  input = 'block-dynamics-aux-fretting-wear-test.i'
  rel_err = 1.0e-3
  abs_zero = 2.0e-4
  min_ad_size = '50'
  map = False
  partial = True
  cli_args = 'Mesh/file=long-bottom-block-1elem-blocks-multiple-projections.e Executioner/end_time=0.2 Outputs/file_base=block-dynamics-aux-fretting-wear-test-quick-multiple-projections_out'
  requirement = 'The contact module shall be able to 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 where secondary elements/nodes may project in an oblique manner and where sidesets only span the possible contact patch area. '
[]

[block-dynamics-aux-fretting-wear-test-quick-multiple-projections-lowerd]
  type = 'Exodiff'
  exodiff = 'block-dynamics-aux-fretting-wear-test-quick-multiple-projections-lowerd_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'
  map = False
  partial = True
  cli_args = 'Mesh/file=long-bottom-block-1elem-blocks-multiple-projections-lowerd.e Executioner/end_time=0.2 Constraints/weighted_gap_lm/correct_edge_dropping=true Constraints/weighted_gap_lm/debug_mesh=true UserObjects/weighted_vel_uo/correct_edge_dropping=false UserObjects/weighted_vel_uo/debug_mesh=true Constraints/normal_x/correct_edge_dropping=true Constraints/normal_x/debug_mesh=true Constraints/normal_y/correct_edge_dropping=true Constraints/normal_y/debug_mesh=true Constraints/tangential_x/correct_edge_dropping=true Constraints/tangential_x/debug_mesh=true Constraints/tangential_y/correct_edge_dropping=true Constraints/tangential_y/debug_mesh=true AuxKernels/gap_vel/correct_edge_dropping=true AuxKernels/gap_vel/debug_mesh=true AuxKernels/worn_depth/correct_edge_dropping=true AuxKernels/worn_depth/debug_mesh=true Outputs/file_base=block-dynamics-aux-fretting-wear-test-quick-multiple-projections-lowerd_out'
  requirement = 'The contact module shall be able to 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 where secondary elements/nodes may project in an oblique manner due to wrap-around lower-dimensional domains and sidesets. This test checks for the ability to discard oblique segments with a default minimum projection angle and the fact that wrapping around the lowerd mesh modifies contact geometry (normals, e.g.) and this has an impact on the contact force direction at the geometry discontinuity.'
[]

[block-dynamics-aux-fretting-wear-test-action-quick]
  type = 'Exodiff'
  exodiff = 'block-dynamics-aux-fretting-wear-test-action-quick_out.e'
  input = 'block-dynamics-aux-fretting-wear-test-action.i'
  rel_err = 1.0e-3
  abs_zero = 5.0e-4
  min_parallel = 1
  max_parallel = 4
  min_ad_size = '50'
  map = False
  partial = True
  cli_args = 'Executioner/end_time=0.15 Outputs/file_base=block-dynamics-aux-fretting-wear-test-action-quick_out'
  requirement = 'The contact module shall be able to 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 using the contact action.'
[]

[block-dynamics-aux-fretting-wear-test-projection_angle]
  type = 'Exodiff'
  exodiff = 'mortar_segment_mesh.e'
  input = 'block-dynamics-aux-fretting-wear-test-projection_angle.i'
  rel_err = 1.0e-3
  abs_zero = 1.0e-4
  max_parallel = 1
  min_ad_size = '50'
  map = False
  partial = True
  requirement = 'The system shall allow for a user controllable input such that the mortar segment mesh does not contain segments built from oblique projections. Therefore, the system shall generate many segments from oblique projections if a very small minimum projection angle is chosen. '
[]

[pressure-aux-frictionless]
  type = CSVDiff
  input = 'pressure-aux-frictionless.i'
  csvdiff = 'pressure-aux-frictionless_out_normal_lm_0011.csv'
  abs_zero = 1.0e-4
  min_ad_size = '50'
  mesh_mode = 'REPLICATED'
  requirement = 'The system shall produce the correct normal contact pressure value for frictionless mortar contact when using an auxiliary kernel.'
[]

[pressure-aux-frictionless-3d]
  type = CSVDiff
  input = 'pressure-aux-frictionless-3d.i'
  csvdiff = 'pressure-aux-frictionless-3d_out_normal_lm_0003.csv'
  abs_zero = 1.0e-4
  min_ad_size = '150'
  mesh_mode = 'REPLICATED'
  requirement = 'The system shall be able to  produce the correct normal contact pressure value in 3D for frictionless mortar contact when using an auxiliary kernel.'
[]

[pressure-aux-friction]
  type = CSVDiff
  input = 'pressure-aux-friction.i'
  csvdiff = 'pressure-aux-friction_out_tangent_lm_0011.csv'
  abs_zero = 1.0e-4
  min_ad_size = '50'
  mesh_mode = 'REPLICATED'
  requirement = 'The system shall be able to produce the correct tangential contact pressure value for frictional mortar contact when using an auxiliary kernel.'
[]

[pressure-aux-friction-3d]
  type = CSVDiff
  input = 'pressure-aux-friction-3d.i'
  csvdiff = 'pressure-aux-friction-3d_out_tangent2_lm_0003.csv'
  abs_zero = 1.0e-4
  min_ad_size = '150'
  mesh_mode = 'REPLICATED'
  requirement = 'The system shall be able to produce the correct tangential contact pressure value in 3D for frictional mortar contact when using an auxiliary kernel.'
[]

[frictionless-weighted-gap-lm]
  type = 'Exodiff'
  exodiff = 'frictionless-weighted-gap-lm_out.e'
  input = 'frictionless-weighted-gap-lm.i'
  rel_err = 5.0e-4
  min_parallel = 4
  max_parallel = 4
  min_ad_size = '50'
  valgrind = 'heavy'
  installation_type = in_tree
  requirement = 'The contact module shall be able to solve a two-dimensional frictionless problem with Cartesian Lagrange multipliers.'
[]

[frictionless-mortar-3d]
  type = 'Exodiff'
  exodiff = 'frictionless-mortar-3d_out.e'
  input = 'frictionless-mortar-3d.i'
  rel_err = 5.0e-4
  max_parallel = 4
  min_ad_size = '150'
  mesh_mode = 'REPLICATED'
  requirement = 'The contact module shall be able to solve a three-dimensional frictionless problem with Cartesian Lagrange multipliers.'
[]

[friction-mortar-3d]
  type = 'Exodiff'
  exodiff = 'frictionless-mortar-3d-friction_out.e'
  input = 'frictionless-mortar-3d-friction.i'
  rel_err = 5.0e-4
  max_parallel = 4
  min_ad_size = '150'
  mesh_mode = 'REPLICATED'
  requirement = 'Contact module shall solve a three-dimensional frictional problem with Cartesian Lagrange multipliers.'
[]

[two_block_1st_order_constraint_lm_xy_friction]
  type = 'Exodiff'
  exodiff = 'two_block_1st_order_constraint_lm_xy_friction_out.e'
  input = 'two_block_1st_order_constraint_lm_xy_friction.i'
  rel_err = 5.0e-4
  max_parallel = 4
  min_ad_size = '50'
  requirement = 'The contact module shall be able to solve a flat surface-flat surface frictional problem with Cartesian Lagrange multipliers.'
[]

[cylinder_friction_cartesian]
  type = 'Exodiff'
  exodiff = 'cylinder_friction_cartesian_out.e'
  input = 'cylinder_friction_cartesian.i'
  rel_err = 5.0e-4
  abs_zero = 1.0e-4
  max_parallel = 4
  min_ad_size = '50'
  method = 'opt'
  requirement = 'The contact module shall be able to solve a cylinder-on-plane plane strain frictional problem with Cartesian Lagrange multipliers.'
[]

[two_block_1st_order_constraint_lm_xy_friction_vcp]
  type = 'Exodiff'
  exodiff = 'two_block_1st_order_constraint_lm_xy_friction_vcp_out.e'
  input = 'two_block_1st_order_constraint_lm_xy_friction_vcp.i'
  rel_err = 5.0e-4
  min_ad_size = '50'
  mesh_mode = 'REPLICATED'
  requirement = 'The contact module shall be able to solve a flat surface-flat surface frictional problem with Cartesian Lagrange multipliers using VCP.'
  # This test is causing sporadic test failures (see issue #24986)
  # Limit the parallel test as a temporary solution here
  max_parallel = 1
[]

[cylinder_friction_cartesian_vcp]
  type = 'Exodiff'
  exodiff = 'cylinder_friction_cartesian_vcp_out.e'
  input = 'cylinder_friction_cartesian_vcp.i'
  rel_err = 5.0e-4
  abs_zero = 1.0e-4
  max_parallel = 4
  min_ad_size = '50'
  method = 'opt'
  mesh_mode = 'REPLICATED'
  requirement = 'The contact module shall be able to solve a cylinder-on-plane plane strain frictional problem with Cartesian Lagrange multipliers using VCP.'
[]

[block-dynamics]
  type = 'Exodiff'
  exodiff = 'block-dynamics_out.e'
  input = 'block-dynamics.i'
  rel_err = 5.0e-4
  abs_zero = 1.0e-5
  min_ad_size = '50'
  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-reference]
  type = 'CSVDiff'
  input = 'block-dynamics-reference.i'
  csvdiff = 'block-dynamics-reference_out.csv'
  rel_err = 5.0e-4
  abs_zero = 1.0e-5
  min_ad_size = '50'
  requirement = 'The system shall solve dynamics mortar contact between two blocks with weighted gap time stabilization and match reference kinetic and elastic energy results.'
[]

[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_ad_size = '50'
  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'
  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'
  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_ad_size = '50'
  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'
  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'
  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 = 150
  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
  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
  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
  requirement = 'The system shall solve a dynamic 3D frictional one-element bouncing block problem with mortar constraints using nodal 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
  requirement = 'The system shall be able to run a two-dimensional frictional model using the contact action for mortar applications'
  valgrind = 'heavy'
  method = '!dbg'
[]

[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
  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'
  prereq = 'frictional_bouncing_block_action_restart_1'
  valgrind = 'heavy'
  method = '!dbg'
  # skip test if test is being run out-of-tree. Issue Ref: #25279
  installation_type = in_tree
[]

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

[finite_short]
  detail = 'using the finite strain formulation for a limited time simulation.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'finite_short.csv'
  cli_args = 'name=finite_short Executioner/num_steps=10 Executioner/nl_abs_tol=1.0e-5'
  valgrind = 'heavy'
  method = '!dbg'
[]

[finite_rr]
  detail = 'using the finite strain formulation and reference residual.'
  type = 'CSVDiff'
  input = 'finite_rr.i'
  cli_args = 'Executioner/nl_abs_tol=1.0e-6 '
  csvdiff = 'finite_rr.csv'
  heavy = true
[]

[finite_rr_full]
  detail = 'using the finite strain formulation and reference residual with extra_vector_tags passed to the constraints.'
  type = 'CSVDiff'
  input = 'finite_rr.i'
  csvdiff = 'finite_rr_full.csv'
  cli_args = 'Contact/frictional/extra_vector_tags=ref name=finite_rr_full Executioner/nl_abs_tol=1.0e-2 Executioner/num_steps=10'
  valgrind = 'heavy'
  method = '!dbg'
[]

[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 Contact/frictional/c_tangential=1e0 Executioner/nl_abs_tol=1.0e-4'
  heavy = true
[]

[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
[]

[finite]
  detail = 'using the finite strain formulation.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'finite_second.csv'
  cli_args = 'name=finite_second order=SECOND elem=QUAD8'
  heavy = true
[]

[finite_rr]
  detail = 'using the finite strain formulation and reference residual.'
  type = 'CSVDiff'
  input = 'finite_rr.i'
  csvdiff = 'finite_rr_second.csv'
  heavy = true
  cli_args = 'name=finite_rr_second order=SECOND elem=QUAD8'
[]

[stiff_stiff]
  detail = 'using with a stiff block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite_stiff.i'
  csvdiff = 'stiff_stiff_second.csv'
  cli_args = 'E_block=1e9 E_plank=1e9 name=stiff_stiff_second Contact/frictional/c_normal=1e7 Contact/frictional/c_tangential=1e-2 order=SECOND elem=QUAD8 Executioner/nl_abs_tol=1.0e-2'
  heavy = true
[]

[soft_stiff]
  detail = 'using with a soft block and a stiff plank.'
  type = 'CSVDiff'
  input = 'finite_stiff.i'
  csvdiff = 'soft_stiff_second.csv'
  cli_args = 'E_block=1e7 E_plank=1e9 name=soft_stiff_second order=SECOND elem=QUAD8 Contact/frictional/c_normal=1e6 Contact/frictional/c_tangential=1e-2 Executioner/nl_abs_tol=1.0e-4'
  heavy = true
[]

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

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

[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
[]

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

[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.'
[]

[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.'
[]

[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
[]

[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
[]

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

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

[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
  valgrind = 'heavy'
[]

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

[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
  valgrind = 'heavy'
  detail = 'using the small strain formulation and calculate a perfect Jacobian.'
[]

[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
  valgrind = 'heavy'
  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
[]

[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
[]

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

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

[finite_short]
  detail = 'using the finite strain formulation for a limited time simulation.'
  type = 'CSVDiff'
  input = 'finite.i'
  csvdiff = 'finite_short.csv'
  cli_args = 'name=finite_short Executioner/num_steps=10'
[]

[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
[]

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

[finite_rr_full]
  detail = 'using the finite strain formulation and reference residual with extra_vector_tags passed to the constraints.'
  type = 'CSVDiff'
  input = 'finite_rr.i'
  csvdiff = 'finite_rr_full.csv'
  cli_args = 'Contact/frictionless/extra_vector_tags=ref name=finite_rr_full Executioner/num_steps=10'
[]

[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
[]

[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
[]

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

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

[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
[]

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

[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
[]

[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
[]

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

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

[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
[]

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

[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
  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
  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
[]

[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
[]

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

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

[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
[]

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

[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.'
[]

[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.'
[]

[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
[]

[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
[]

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

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

[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
[]

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

[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
[]

[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
[]

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

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

[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
[]

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

[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
[]

[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
[]

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

[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 = 'The system shall set up node-face interactions between multiple surface pairs using a contact action syntax.'
[]

[multiple_pairs_cleanup]
  type = Exodiff
  input = 'multiple_pairs.i'
  exodiff = 'multiple_pairs_cleanup_out.e'
  abs_zero = 1.0e-8
  rel_err = 1.0e-5
  cli_args = "Postprocessors/active='' Contact/action_name/primary='20 20 101 20' Contact/action_name/secondary='10 101 20 101' Outputs/file_base=multiple_pairs_cleanup_out"
  requirement = 'The system shall set up node-face interactions between multiple surface pairs using a contact action syntax and remove repeated contact pairs.'
[]

[multiple_pairs_mortar]
  type = Exodiff
  input = 'multiple_pairs_mortar.i'
  exodiff = 'multiple_pairs_mortar_out.e'
  abs_zero = 1.0e-8
  rel_err = 1.0e-5
  cli_args = "Postprocessors/active=''"
  map = False
  partial = True
  requirement = 'The system shall set up mortar interactions between multiple surface pairs using a contact action syntax with distinct input blocks.'
[]

[multiple_pairs_mortar_friction]
  type = Exodiff
  input = 'multiple_pairs_mortar_friction.i'
  exodiff = 'multiple_pairs_mortar_friction_out.e'
  abs_zero = 1.0e-8
  rel_err = 1.0e-5
  cli_args = "Postprocessors/active=''"
  map = False
  partial = True
  requirement = 'The system shall set up frictional mortar interactions between multiple surface pairs using a contact action syntax with distinct input blocks.'
[]

[three_hexagons_coarse]
  type = 'Exodiff'
  input = 'three_hexagons_coarse.i'
  exodiff = 'three_hexagons_coarse_out.e'
  requirement = 'The system shall assign and accumulate a single contact pressure auxiliary variable for all contact pairs supplied to a single contact action object.'
  abs_zero = 1.0e-4
  rel_err = 1.0e-4
  allow_warnings = true
[]

[three_hexagons_coarse_automatic_pair_reference]
  type = 'Exodiff'
  issues = '#24122'
  input = 'three_hexagons_coarse.i'
  exodiff = 'three_hexagons_coarse_automatic_pair_reference_out.e'
  cli_args = "Postprocessors/active='' Contact/contact_pressure/secondary='102 102 201' Contact/contact_pressure/primary='201 301 301' Outputs/hide=penetration Outputs/file_base=three_hexagons_coarse_automatic_pair_reference_out"
  requirement = 'The system shall assign and accumulate a single contact pressure auxiliary variable for all contact pairs supplied to a single contact action object with a different combination of primary-secondary pairs.'
  abs_zero = 1.0e-4
  rel_err = 1.0e-4
  allow_warnings = true
[]

[three_hexagons_coarse_automatic_pair]
  type = 'Exodiff'
  issues = '#24122'
  input = 'three_hexagons_coarse_automatic_pair.i'
  exodiff = 'three_hexagons_coarse_automatic_pair_out.e'
  cli_args = "Contact/contact_pressure/automatic_pairing_boundaries='102 201 301' Contact/contact_pressure/automatic_pairing_method='CENTROID' Outputs/hide=penetration Postprocessors/active=''"
  requirement = 'The system shall be able to, if the user provides a reference distance, assign primary-secondary contact pairs automatically.'
  abs_zero = 1.0e-4
  rel_err = 1.0e-4
  allow_warnings = true
[]

[three_hexagons_coarse_automatic_pair_node_proximity]
  type = 'Exodiff'
  issues = '#24609'
  input = 'three_hexagons_coarse_automatic_pair.i'
  exodiff = 'three_hexagons_coarse_automatic_pair_node_proximity_out.e'
  cli_args = "Contact/contact_pressure/automatic_pairing_boundaries='102 201 301' Contact/contact_pressure/automatic_pairing_method=NODE Contact/contact_pressure/automatic_pairing_distance=0.002001 Outputs/file_base=three_hexagons_coarse_automatic_pair_node_proximity_out Outputs/hide=penetration"
  requirement = 'The system shall be able to, if the user provides a reference distance, assign primary-secondary contact pairs automatically using proximity between nodes that belong to arbitrary boundaries.'
  abs_zero = 1.0e-4
  rel_err = 1.0e-4
  allow_warnings = true
[]

[three_hexagons_coarse_automatic_pair_node_proximity_small_distance]
  type = 'Exodiff'
  issues = '#24609'
  input = 'three_hexagons_coarse_automatic_pair.i'
  exodiff = 'three_hexagons_coarse_automatic_pair_node_proximity_small_distance_out.e'
  cli_args = "Contact/contact_pressure/automatic_pairing_boundaries='102 201 301' Contact/contact_pressure/automatic_pairing_method=NODE Contact/contact_pressure/automatic_pairing_distance=0.0013 Outputs/file_base=three_hexagons_coarse_automatic_pair_node_proximity_small_distance_out Outputs/hide=penetration"
  requirement = 'The system shall be able to, if the user provides a too-small reference distance, avoid assigning primary-secondary contact pairs using proximity between nodes.'
  abs_zero = 1.0e-4
  rel_err = 1.0e-4
  allow_warnings = true
[]

[three_hexagons_coarse_automatic_pair_node_proximity_large_distance]
  type = 'Exodiff'
  issues = '#24609'
  input = 'three_hexagons_coarse_automatic_pair.i'
  exodiff = 'three_hexagons_coarse_automatic_pair_node_proximity_large_distance_out.e'
  cli_args = "Contact/contact_pressure/automatic_pairing_boundaries='102 201 301' Contact/contact_pressure/automatic_pairing_method=NODE Contact/contact_pressure/automatic_pairing_distance=4.0 Outputs/file_base=three_hexagons_coarse_automatic_pair_node_proximity_large_distance_out Outputs/hide=penetration"
  requirement = 'The system shall be able to, if the user provides a very large reference distance, assign primary-secondary contact pairs automatically using proximity between nodes that belong to arbitrary boundaries.'
  abs_zero = 1.0e-4
  rel_err = 1.0e-4
  allow_warnings = true
[]

[three_hexagons_coarse_automatic_pair_error]
  type = 'RunException'
  issues = '#24122'
  input = 'three_hexagons_coarse_automatic_pair.i'
  expect_err = "'primary' and 'secondary' surfaces or a list of boundaries"
  requirement = 'The system shall generate an error if insufficient information to automatically define contact surface pairs is provided.'
  allow_warnings = true
[]

[three_hexagons_coarse_various_actions]
  type = 'Exodiff'
  input = 'three_hexagons_coarse_various_actions.i'
  exodiff = 'three_hexagons_coarse_various_actions_out.e'
  requirement = 'The system shall assign and accumulate a single contact pressure auxiliary variable for all contact pairs supplied to a distinct contact action object.'
  abs_zero = 1.0e-4
  rel_err = 1.0e-4
  allow_warnings = true
[]

[split_sidesets]
  type = 'Exodiff'
  input = 'split_sidesets.i'
  exodiff = 'split_sidesets_out.e'
  requirement = 'The system shall compute mechanical contact with multiple pairs, which are defined with separate contact sidesets.'
  abs_zero = 1.0e-4
  rel_err = 1.0e-4
  allow_warnings = true
[]

[continuous_sidesets]
  type = 'Exodiff'
  input = 'continuous_sidesets.i'
  exodiff = 'continuous_sidesets_out.e'
  requirement = 'The system shall compute mechanical contact allowing a secondary surface to enforce constraints with multiple primary surfaces.'
  abs_zero = 1.0e-4
  rel_err = 1.0e-4
  allow_warnings = true
[]

[./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'
  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'
  mesh_mode = 'REPLICATED'
[]

[./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
  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
  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_x_disp_FINAL.csv cylinder_friction_check_cont_press_FINAL.csv cylinder_friction_check_y_disp_FINAL.csv cylinder_friction_check_friction_FINAL.csv'
  input = 'cylinder_friction.i'
  rel_err = 5.0e-4
  min_parallel = 4
  heavy = true
  requirement = 'The contact module shall be able to solve frictional contact between a cylinder and a plane.'
[]

[cylinder_friction_penalty]
  type = 'CSVDiff'
  csvdiff = 'cylinder_friction_penalty_check_x_disp_FINAL.csv cylinder_friction_penalty_check_cont_press_FINAL.csv cylinder_friction_penalty_check_y_disp_FINAL.csv cylinder_friction_penalty_check_friction_FINAL.csv'
  input = 'cylinder_friction_penalty.i'
  rel_err = 7.5e-4
  abs_zero = 1.0e-4
  heavy = true
  requirement = "The contact module shall be able to approximate frictional contact solution between a cylinder and a plane using penalty mortar contact."
[]

[cylinder_friction_penalty_adaptivity]
  type = 'CSVDiff'
  csvdiff = cylinder_friction_penalty_adaptivity_surface_FINAL.csv
  input = 'cylinder_friction_penalty_adaptivity.i'
  rel_err = 5.0e-4
  abs_zero = 1.0e-2
  valgrind = none
  heavy = true
  requirement = "The contact module shall be able to approximate frictional contact solution between a cylinder and a plane using penalty mortar contact after performing initial adaptivity around the contact area."
[]

[cylinder_friction_penalty_al_normal]
  type = 'CSVDiff'
  csvdiff = 'cylinder_friction_penalty_normal_al_vectorpp_output_surface_FINAL.csv'
  input = 'cylinder_friction_penalty_normal_al.i'
  rel_err = 5.0e-4
  abs_zero = 1.0e-2
  heavy = true
  method = opt
  valgrind = none
  requirement = "The contact module shall be able to approximate normal contact in a two-dimensional cylinder via an augmented Lagrange approach to constraint enforcement."
[]

[action_error]
  type = 'RunException'
  input = 'cylinder_friction_penalty_frictional_al_action.i'
  expect_err = "Augmented Lagrange parameter was specified, but the selected problem type does not support Augmented Lagrange iterations."
  cli_args = 'Problem/type=FEProblem'
  requirement = "The contact action shall be alert the user if Augmented Lagrage parameters are supplied, but the selected Problem class does not support AL."
[]

[cylinder_friction_penalty_al_frictional]
  type = 'CSVDiff'
  csvdiff = 'cylinder_friction_penalty_frictional_al_vectorpp_output_surface_FINAL.csv'
  input = 'cylinder_friction_penalty_frictional_al.i'
  rel_err = 1.0e-4
  abs_zero = 1.0e-2
  heavy = true
  valgrind = none
  method = opt
  requirement = "The contact module shall be able to approximate frictional contact in a two-dimensional cylinder via an augmented Lagrange approach to constraint enforcement."
[]

[cylinder_friction_penalty_frictional_al_tight_slip]
  type = 'CSVDiff'
  csvdiff = 'cylinder_friction_penalty_frictional_al_tight_slip_vectorpp_output_surface_FINAL.csv'
  input = 'cylinder_friction_penalty_frictional_al_tight_slip.i'
  rel_err = 1.0e-3
  abs_zero = 1.0e-4
  heavy = true
  valgrind = none
  method = opt
  requirement = "The contact module shall be able to approximate frictional contact in a two-dimensional cylinder via an augmented Lagrange approach to constraint enforcement with a tight slip distance tolerance."
[]

[cylinder_friction_penalty_al_frictional_action]
  type = 'CSVDiff'
  csvdiff = 'cylinder_friction_penalty_frictional_al_action_vectorpp_output_surface_FINAL.csv'
  input = 'cylinder_friction_penalty_frictional_al_action.i'
  rel_err = 1.0e-4
  abs_zero = 1.0e-2
  heavy = true
  valgrind = none
  method = opt
  requirement = "The contact module shall be able to approximate frictional contact in a two-dimensional cylinder via an augmented Lagrange approach to constraint enforcement using the contact action."
[]

[cylinder_friction_penalty_al_frictional_action_amg]
  type = 'CSVDiff'
  csvdiff = 'cylinder_friction_penalty_frictional_al_action_amg_vectorpp_output_surface_FINAL.csv'
  input = 'cylinder_friction_penalty_frictional_al_action_amg.i'
  rel_err = 1.0e-4
  abs_zero = 1.0e-2
  requirement = "The contact module shall be able to approximate frictional contact in a two-dimensional cylinder via an augmented Lagrange approach to constraint enforcement using the contact action and use the hypre algebraic multigrid preconditioner since the augmented Lagrange (Uzawa) constraint enfrocement approach allows for maintaining reasonable system condition numbers."
[]

[cylinder_friction_penalty_frictional_al_action_amg_tight]
  type = 'CSVDiff'
  csvdiff = 'cylinder_friction_penalty_frictional_al_action_amg_tight_vectorpp_output_surface_FINAL.csv'
  input = 'cylinder_friction_penalty_frictional_al_action_amg_tight.i'
  rel_err = 1.0e-4
  abs_zero = 1.0e-2
  method = opt
  valgrind = none
  requirement = "The contact module shall be able to approximate frictional contact in a two-dimensional cylinder via an augmented Lagrange approach to constraint enforcement using the contact action and use the hypre algebraic multigrid preconditioner since the augmented Lagrange (Uzawa) constraint enfrocement approach allows for maintaining reasonable system condition numbers and selecting the 'Bussetta' and 'Simple' adaptivity strategies for the normal and frictional penalty values, respectively."
[]

[cylinder_friction_penalty_al_frictional_action_amg_bussetta_simple]
  type = 'CSVDiff'
  csvdiff = 'cylinder_friction_penalty_frictional_al_action_amg_bussetta_simple_vectorpp_output_surface_FINAL.csv'
  input = 'cylinder_friction_penalty_frictional_al_action_amg_bussetta_simple.i'
  rel_err = 1.0e-4
  abs_zero = 1.0e-2
  requirement = "The contact module shall be able to approximate frictional contact in a two-dimensional cylinder via an augmented Lagrange approach to constraint enforcement using the contact action and use the hypre algebraic multigrid preconditioner since the augmented Lagrange (Uzawa) constraint enfrocement approach allows for maintaining reasonable system condition numbers and selecting the 'Bussetta' and 'FrictionLimit' adaptivity strategies for the normal and frictional penalty values, respectively."
  method = '!dbg'
[]

[cylinder_friction_penalty_al_separation]
  type = 'CSVDiff'
  csvdiff = 'cylinder_friction_penalty_normal_al_backup_vectorpp_output_surface_FINAL.csv'
  input = 'cylinder_friction_penalty_normal_al_backup.i'
  rel_err = 5.0e-4
  abs_zero = 1.0e-2
  heavy = true
  valgrind = none
  method = opt
  requirement = "The contact module shall be able to approximate normal contact in a two-dimensional cylinder via an augmented Lagrange approach to constraint enforcement and allow for separation."
[]

[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'
  requirement = 'The contact module shall be able to solve frictional contact between a semicircular tool and flexible base material.'
[]

[ironing_penalty]
  type = 'CSVDiff'
  csvdiff = 'ironing_penalty_out_penalty_normal_pressure_0012.csv'
  input = 'ironing_penalty.i'
  rel_err = 1.0e-4
  abs_zero = 1.0e-4
  min_ad_size = '200'
  cli_args = 'Executioner/end_time=0.15'
  requirement = "The contact module shall be able to solve frictional contact between a semicircular tool and flexible base material using a penalty approach to mortar contact."
[]

[ironing_penalty_action]
  type = 'CSVDiff'
  csvdiff = 'ironing_penalty_action_out_penalty_normal_pressure_0012.csv'
  input = 'ironing_penalty_action.i'
  rel_err = 1.0e-4
  abs_zero = 1.0e-4
  min_ad_size = '200'
  cli_args = 'Executioner/end_time=0.15'
  requirement = "The contact module shall be able to solve frictional contact between a semicircular tool and flexible base material using a penalty approach to mortar contact through the contact action."
[]

[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
  heavy = true
  requirement = 'The contact module shall be able to 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
  heavy = true
  requirement = "The contact module shall be able to 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
  method = '!dbg'
  valgrind = 'none'
  cli_args = "Constraints/weighted_gap_lm/correct_edge_dropping=true UserObjects/weighted_vel_uo/correct_edge_dropping=true Executioner/end_time=5.5 Outputs/file_base=frictional_bouncing_block_edge_dropping"
  requirement = "The contact module shall be able to solve frictional contact between a bouncing block and a flexible substrate when correct edge dropping is enabled. An additional requirement is that the correct edge dropping treatment must yield same results as an incorrect edge dropping treatment when there is not edge dropping, e.g. when the entire 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'
  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_pg]
  type = 'CSVDiff'
  input = 'two_equal_blocks_compress_2d_pg.i'
  csvdiff = 'two_equal_blocks_compress_2d_pg_out.csv'
  requirement = 'The system shall simulate correct contact behavior in 2D when two blocks with the same height come into contact using the Petrov-Galerkin approach.'
  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'
  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'
  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_pg]
  type = 'CSVDiff'
  input = 'two_equal_blocks_compress_3d_pg.i'
  csvdiff = 'two_equal_blocks_compress_3d_pg_out.csv'
  requirement = 'The system shall simulate correct contact behavior in 3D when two blocks with the same height come into contact using the Petrov-Galerkin approach.'
  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'
  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'
  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
  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
  allow_warnings = true
  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
  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.'
[../]

[./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.'
[../]

[./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
  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
  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 = 2e-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_always]
  type = 'Exodiff'
  input = 'sliding_update.i'
  exodiff = 'sliding_update_out.e'
  cli_args = "Mesh/patch_update_strategy='always'"
  abs_zero = 1e-7
  requirement = "The system shall enforce the automatic patch update using the 'always' option."
[]

[sliding_blocks]
  type = 'Exodiff'
  input = 'sliding_update.i'
  exodiff = 'sliding_blocks_iteration.e'
  cli_args = "Mesh/patch_update_strategy='iteration' Outputs/file_base=sliding_blocks_iteration"
  abs_zero = 1e-7
  requirement = "The system shall enforce the automatic patch update using the iteration option and give results equal to the 'always' option."
[]

[error_test]
  type = RunException
  input = 'sliding_update.i'
  cli_args = "Mesh/ghosting_patch_size=2 Mesh/patch_update_strategy='iteration' 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
  # max parallel set to avoid not running into error for more procs
  max_parallel = 2
  mesh_mode = REPLICATED
  requirement = "The system shall enforce that the nearest neighbor
                   node is inside the ghosted set of elements."
[]

[iteration_adaptivity_parallel]
  type = 'Exodiff'
  input = 'iteration_adaptivity_parallel.i'
  exodiff = 'iteration_adaptivity_parallel_out.e'
  # Issue #24824: Unable to make it distributed at the moment.
  mesh_mode = replicated
  map = false
  requirement = "The system shall be able to couple iteration patch update strategies with initial adaptivity when using geometric searches with a replicated mesh."
[]

[iteration_adaptivity_parallel_node_face]
  type = 'Exodiff'
  input = 'iteration_adaptivity_parallel_node_face.i'
  exodiff = 'iteration_adaptivity_parallel_node_face_out.e'
  abs_zero = 5.0e-8
  rel_err = 5.0e-5
  requirement = "The system shall be able to couple iteration patch update strategies with initial adaptivity when using geometric searches using node face contact with a replicated or distributed mesh."
[]

[./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)."
[../]