Ray Tracing Requirements Traceability Matrix

[missing]
  type = RunException
  input = 'add_raybc_action.i'
  cli_args = 'RayBCs/active=missing_study_by_name'
  expect_err = "NullRayBC 'missing_study_by_name': Could not find the requested study 'dummy'."
  allow_test_objects = true

  detail = 'could not be found '
[]

[not_a_study]
  type = RunException
  input = 'add_raybc_action.i'
  cli_args = 'RayBCs/active=not_a_study UserObjects/active=not_a_study'
  expect_err = "NullRayBC 'not_a_study' requested the study not_a_study but the provided study is not a RayTracingStudy-derived object."
  allow_test_objects = true

  detail = 'or is not a RayTracingStudy.'
[]

[multiple]
  type = RunException
  input = 'add_raybc_action.i'
  cli_args = 'RayBCs/active=multiple_studies UserObjects/active="study another_study"'
  expect_err = "While constructing the NullRayBC 'multiple_studies', multiple RayTracingStudy objects were found."
  allow_test_objects = true

  detail = 'and multiple studies were found '
[]

[missing]
  type = RunException
  input = 'add_raybc_action.i'
  cli_args = 'RayBCs/active=missing_study'
  expect_err = "NullRayBC 'missing_study' did not provide a RayTracingStudy to associate with via the 'study' parameter and a study was not found."
  allow_test_objects = true

  detail = 'or when no studies were found.'
[]

[missing]
  type = RunException
  input = 'add_raykernel_action.i'
  cli_args = 'RayKernels/active=missing_study_by_name'
  expect_err = "NullRayKernel 'missing_study_by_name': Could not find the requested study 'dummy'."
  allow_test_objects = true

  detail = 'could not be found '
[]

[not_a_study]
  type = RunException
  input = 'add_raykernel_action.i'
  cli_args = 'RayKernels/active=not_a_study UserObjects/active=not_a_study'
  expect_err = "NullRayKernel 'not_a_study' requested the study not_a_study but the provided study is not a RayTracingStudy-derived object."
  allow_test_objects = true

  detail = 'or is not a RayTracingStudy.'
[]

[multiple]
  type = RunException
  input = 'add_raykernel_action.i'
  cli_args = 'RayKernels/active=multiple_studies UserObjects/active="study another_study"'
  expect_err = "While constructing the NullRayKernel 'multiple_studies', multiple RayTracingStudy objects were found."
  allow_test_objects = true

  detail = 'and multiple studies were found '
[]

[missing]
  type = RunException
  input = 'add_raykernel_action.i'
  cli_args = 'RayKernels/active=missing_study'
  expect_err = "NullRayKernel 'missing_study' did not provide a RayTracingStudy to associate with via the 'study' parameter and a study was not found."
  allow_test_objects = true

  detail = 'or when no studies were found.'
[]

[get_study_bad_param]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'RayKernels/test/get_study_bad=true RayKernels/test/study=study'
  allow_test_objects = true
  expect_err = 'Supplied study of type RayTracingStudyTest is not the required study type RepeatableRayStudy'
  detail = 'the study is provided via param and is of the wrong type'
[]

[get_study_bad_default]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'RayKernels/test/get_study_bad=true'
  allow_test_objects = true
  expect_err = 'Supplied study of type RayTracingStudyTest is not the required study type RepeatableRayStudy'
  detail = 'the study is provided by default and is of the wrong type'
[]

[rz]
  type = 'CSVDiff'
  input = 'rz_line_integral.i'
  csvdiff = 'rz_line_integral_out.csv'

  detail = 'rz coordinate system '
[]

[rspherical]
  type = 'CSVDiff'
  input = 'rspherical_line_integral.i'
  csvdiff = 'rspherical_line_integral_out.csv'

  detail = 'and in the rspherical coordinate system.'
[]

[quad4]
  type = 'Exodiff'
  input = 'ray_mesh_output.i'
  exodiff = 'quad4_rays.e'
  cli_args = 'Outputs/rays/file_base=quad4_rays'

  detail = 'QUAD4,'
[]

[tri3]
  type = 'Exodiff'
  input = 'ray_mesh_output.i'
  exodiff = 'tri3_rays.e'
  cli_args = 'Mesh/gmg/elem_type=TRI3
                  Outputs/rays/file_base=tri3_rays'
  detail = 'TRI3,'
[]

[hex8]
  type = 'Exodiff'
  input = 'ray_mesh_output_3d.i'
  exodiff = 'hex8_rays.e'
  cli_args = 'Outputs/rays/file_base=hex8_rays'
  detail = 'HEX8, and'
[]

[tet4]
  type = 'Exodiff'
  input = 'ray_mesh_output_3d.i'
  exodiff = 'tet4_rays.e'
  cli_args = 'Mesh/gmg/elem_type=tet4
                  Outputs/rays/file_base=tet4_rays'
  detail = 'TET4.'
[]

[transient]
  type = Exodiff
  input = 'ray_mesh_output_transient.i'
  exodiff = 'ray_mesh_output_transient_rays.e
               ray_mesh_output_transient_rays.e-s0001
               ray_mesh_output_transient_rays.e-s0002
               ray_mesh_output_transient_rays.e-s0003'
  allow_test_objects = true
  requirement = 'The system shall support the mesh output of traced rays into the exodus format for a transient simulation'
  design = 'RayTracingExodus.md RayTracingMeshOutput.md'
[]

[nemesis]
  type = Exodiff
  input = 'ray_mesh_output.i'
  cli_args = 'Mesh/Partitioner/type=GridPartitioner
                Mesh/Partitioner/nx=2
                Mesh/Partitioner/ny=2
                Outputs/rays/type=RayTracingNemesis
                Outputs/rays/output_properties="ray_id intersections pid processor_crossings trajectory_changes"
                Outputs/rays/file_base=nemesis_rays'
  exodiff = 'nemesis_rays.e.4.0 nemesis_rays.e.4.1 nemesis_rays.e.4.2 nemesis_rays.e.4.3'
  min_parallel = 4
  max_parallel = 4
  recover = false
  requirement = 'The system shall support the mesh output of traced rays using the Nemesis format.'
  design = 'RayTracingNemesis.md RayTracingMeshOutput.md'
[]

[nemesis_transient]
  type = Exodiff
  input = 'ray_mesh_output_transient.i'
  cli_args = 'Mesh/Partitioner/type=GridPartitioner
                Mesh/Partitioner/nx=2
                Mesh/Partitioner/ny=2
                Outputs/rays/type=RayTracingNemesis
                Outputs/rays/output_properties="intersections pid processor_crossings trajectory_changes"
                Outputs/rays/file_base=ray_mesh_output_transient_nemesis_rays'
  # Missing some files here beacuse exodiff doesn't like diffing empty output,
  # which is the case for the early transient when not all procs are hit
  exodiff = 'ray_mesh_output_transient_nemesis_rays.e.4.0
               ray_mesh_output_transient_nemesis_rays.e.4.1
               ray_mesh_output_transient_nemesis_rays.e-s0001.4.0
               ray_mesh_output_transient_nemesis_rays.e-s0001.4.1
               ray_mesh_output_transient_nemesis_rays.e-s0002.4.0
               ray_mesh_output_transient_nemesis_rays.e-s0002.4.1
               ray_mesh_output_transient_nemesis_rays.e-s0002.4.3
               ray_mesh_output_transient_nemesis_rays.e-s0003.4.0
               ray_mesh_output_transient_nemesis_rays.e-s0003.4.1
               ray_mesh_output_transient_nemesis_rays.e-s0003.4.2
               ray_mesh_output_transient_nemesis_rays.e-s0003.4.3'
  allow_test_objects = true
  min_parallel = 4
  max_parallel = 4
  requirement = 'The system shall support the mesh output of traced rays into the nemesis format for a transient simulation'
  design = 'RayTracingNemesis.md RayTracingMeshOutput.md'
  # cannot do recover because we don't have a single exodus file, the first two timesteps won't exist
  recover = false
[]

[exodus]
  type = Exodiff
  input = 'ray_mesh_output.i'
  cli_args = 'Mesh/Partitioner/type=GridPartitioner
                  Mesh/Partitioner/ny=3
                  UserObjects/study/segments_on_cache_traces=false
                  Outputs/rays/output_properties="ray_id intersections pid processor_crossings trajectory_changes"
                  Outputs/rays/file_base=no_segments_rays'
  exodiff = 'no_segments_rays.e'
  min_parallel = 3
  max_parallel = 3
  recover = false
  detail = 'Exodus format and'
[]

[nemesis]
  type = Exodiff
  input = 'ray_mesh_output.i'
  cli_args = 'Mesh/Partitioner/type=GridPartitioner
                  Mesh/Partitioner/ny=2
                  Outputs/rays/type=RayTracingNemesis
                  UserObjects/study/segments_on_cache_traces=false
                  Outputs/rays/output_properties="ray_id intersections pid processor_crossings trajectory_changes"
                  Outputs/rays/file_base=no_segments_nemesis_rays'
  exodiff = 'no_segments_nemesis_rays.e.2.0 no_segments_nemesis_rays.e.2.1'
  min_parallel = 2
  max_parallel = 2
  recover = false
  detail = 'in Nemesis format.'
[]

[exodus]
  type = Exodiff
  input = 'ray_mesh_output_data.i'
  exodiff = 'ray_mesh_output_data_rays.e ray_mesh_output_data_rays_nodal.e'
  detail = 'exodus and'
[]

[nemesis]
  type = Exodiff
  input = 'ray_mesh_output_data.i'
  cli_args = 'Mesh/Partitioner/type=GridPartitioner
                  Mesh/Partitioner/nx=2
                  Outputs/rays/type=RayTracingNemesis
                  Outputs/rays_nodal/type=RayTracingNemesis'
  exodiff = 'ray_mesh_output_data_rays.e.2.0 ray_mesh_output_data_rays.e.2.1 ray_mesh_output_data_rays_nodal.e.2.0 ray_mesh_output_data_rays_nodal.e.2.1'
  min_parallel = 2
  max_parallel = 2
  recover = false
  detail = 'nemesis formats.'
[]

[data]
  type = Exodiff
  input = 'ray_mesh_output_data.i'
  exodiff = 'ray_mesh_output_data_specific.e ray_mesh_output_data_specific_nodal.e'
  cli_args = "UserObjects/study/ray_data_names='foo'
                  UserObjects/study/initial_ray_data='5; 6; 7; 8'
                  Outputs/rays/output_data=false
                  Outputs/rays/output_data_names='foo'
                  Outputs/rays/file_base=ray_mesh_output_data_specific
                  Outputs/rays_nodal/output_data=false
                  Outputs/rays_nodal/output_data_names='foo'
                  Outputs/rays_nodal/file_base=ray_mesh_output_data_specific_nodal"
  detail = 'ray data'
[]

[aux_data]
  type = Exodiff
  input = 'ray_mesh_output_data.i'
  exodiff = 'ray_mesh_output_aux_data_specific.e'
  cli_args = "UserObjects/study/ray_aux_data_names='test_aux foo'
                  UserObjects/study/initial_ray_aux_data='1 5; 2 6; 3 7; 4 8'
                  Outputs/active=rays
                  Outputs/rays/output_aux_data=false
                  Outputs/rays/output_aux_data_names='foo'
                  Outputs/rays/file_base=ray_mesh_output_aux_data_specific"
  detail = 'ray aux data'
[]

[no_cached_data]
  type = RunException
  input = 'ray_mesh_output_errors.i'
  cli_args = 'Outputs/rays/output_data=true'
  expect_err = "In order to output Ray data in output 'rays', the RayTracingStudy 'study' must set data_on_cache_traces = true"

  detail = 'ray data is requested to be output but it is not cached,'
[]

[no_cached_aux_data]
  type = RunException
  input = 'ray_mesh_output_errors.i'
  cli_args = 'Outputs/rays/output_aux_data=true'
  expect_err = "In order to output Ray aux data in output 'rays', the RayTracingStudy 'study' must set aux_data_on_cache_traces = true"

  detail = 'aux data is requested to be output but it is not cached,'
[]

[data_nodal_no_data]
  type = RunException
  input = 'ray_mesh_output_errors.i'
  cli_args = 'UserObjects/study/data_on_cache_traces=true
                  Outputs/rays/output_data_nodal=true'
  expect_err = "Cannot be used unless there is data to output"

  detail = 'nodal data is requested but data is not requested,'
[]

[nodal_no_segments]
  type = RunException
  input = 'ray_mesh_output_errors.i'
  cli_args = 'UserObjects/study/data_on_cache_traces=true
                  UserObjects/study/segments_on_cache_traces=false
                  Outputs/rays/output_data=true
                  Outputs/rays/output_data_nodal=true'
  expect_err = "Not supported when study segments_on_cache_traces = false"

  detail = 'nodal data is requested to be output but segment-wise information is not cached,'
[]

[data_all_and_specific]
  type = RunException
  input = 'ray_mesh_output_errors.i'
  cli_args = 'UserObjects/study/data_on_cache_traces=true
                  Outputs/rays/output_data=true
                  Outputs/rays/output_data_names=unused'
  expect_err = "Cannot be used in addition to 'output_data_names'"

  detail = 'all data and specific data are both requested,'
[]

[aux_data_all_and_specific]
  type = RunException
  input = 'ray_mesh_output_errors.i'
  cli_args = 'UserObjects/study/aux_data_on_cache_traces=true
                  Outputs/rays/output_aux_data=true
                  Outputs/rays/output_aux_data_names=unused'
  expect_err = "Cannot be used in addition to 'output_aux_data_names'"

  detail = 'all aux data and specific aux data are both requested'
[]

[missing_data]
  type = RunException
  input = 'ray_mesh_output_errors.i'
  cli_args = 'UserObjects/study/data_on_cache_traces=true
                  UserObjects/study/always_cache_traces=true
                  Outputs/rays/output_data_names=no_exist'
  expect_err = "The ray data 'no_exist' is not registered"

  detail = 'requesting to output data and the data does not exist'
[]

[missing_aux_data]
  type = RunException
  input = 'ray_mesh_output_errors.i'
  cli_args = 'UserObjects/study/aux_data_on_cache_traces=true
                  UserObjects/study/always_cache_traces=true
                  Outputs/rays/output_aux_data_names=no_exist'
  expect_err = "The ray aux data 'no_exist' is not registered"

  detail = 'requesting to output aux data and the aux data does not exist'
[]

[no_segments_warning]
  type = RunException
  input = 'ray_mesh_output_errors.i'
  expect_err = "No cached trace segments were found in the study 'study'."

  design = 'RayTracingMeshOutput.md'
  requirement = 'The system shall report a warning when rays are requested to be output in a mesh format but cached ray inforamtion is not found.'
[]

[stationary]
  type = Exodiff
  input = 'stationary.i'
  exodiff = 'stationary_rays.e stationary_rays.e-s0001 stationary_rays.e-s0002'
  allow_test_objects = true
  design = 'RayTracingMeshOutput.md'
  requirement = 'The system shall support the mesh output of stationary rays'
[]

[test]
  type = 'CSVDiff'
  input = 'ray_data_value.i'
  csvdiff = 'ray_data_value_out.csv'
  requirement = 'The system shall support the ability to output the data stored on a ray.'
[]

[no_registration]
  type = RunException
  input = 'ray_data_value_errors.i'
  expect_err = 'LotsOfRaysRayStudy "lots" does not support Ray registration'
  cli_args = 'UserObjects/active=lots Postprocessors/active=no_registration'
  allow_test_objects = true

  detail = 'the RayTracingStudy does not support accessing rays by name,'
[]

[no_banking]
  type = RunException
  input = 'ray_data_value_errors.i'
  expect_err = "DisableRayBankingStudy \"no_banking_study\" does not bank Rays on completion."
  cli_args = 'UserObjects/active=no_banking_study Postprocessors/active=no_banking'
  allow_test_objects = true

  detail = 'the RayTracingStudy does not support the banking of rays on completion,'
[]

[ray_name_not_found]
  type = RunException
  input = 'ray_data_value_errors.i'
  expect_err = "Could not find a registered Ray with the name 'dummy'"
  cli_args = 'UserObjects/active=repeatable Postprocessors/active=ray_name_not_found'
  allow_test_objects = true

  detail = 'the named ray was not found,'
[]

[data_name_not_found]
  type = RunException
  input = 'ray_data_value_errors.i'
  expect_err = "The RepeatableRayStudy \"repeatable\" does not have Ray data associated with the name 'dummy'"
  cli_args = 'UserObjects/active=repeatable Postprocessors/active=data_name_not_found'
  allow_test_objects = true

  detail = 'the named ray data was not found,'
[]

[aux_data_name_not_found]
  type = RunException
  input = 'ray_data_value_errors.i'
  expect_err = "The RepeatableRayStudy \"repeatable\" does not have Ray aux data associated with the name 'dummy'"
  cli_args = 'UserObjects/active=repeatable Postprocessors/active=aux_data_name_not_found'
  allow_test_objects = true

  detail = 'the named ray aux data was not found,'
[]

[id_not_found]
  type = RunException
  input = 'ray_data_value_errors.i'
  expect_err = "Could not find a Ray with the ID 1 in the Ray banks."
  cli_args = 'UserObjects/active=repeatable Postprocessors/active=id_not_found RayBCs/active=kill'
  allow_test_objects = true

  detail = 'the ray with the given id was not found,'
[]

[neither_provided]
  type = RunException
  input = 'ray_data_value_errors.i'
  expect_err = "Must have either the 'ray_id' or the 'ray_name' param set."
  cli_args = 'UserObjects/active=repeatable Postprocessors/active=neither_provided RayBCs/active=kill'
  allow_test_objects = true

  detail = 'neither an ID or name was provided, and'
[]

[both_provided]
  type = RunException
  input = 'ray_data_value_errors.i'
  expect_err = "Either 'ray_id' or 'ray_name' must be set, but not both"
  cli_args = 'UserObjects/active=repeatable Postprocessors/active=both_provided RayBCs/active=kill'
  allow_test_objects = true

  detail = 'both an ID and name was provided.'
[]

[not_integral_ray_kernel]
  type = RunException
  input = 'ray_integral_value_errors.i'
  cli_args = 'RayKernels/active=null Postprocessors/active=not_integral_ray_kernel'
  expect_err = "NullRayKernel is not derived from an IntegralRayKernel."

  detail = 'the provided IntegralRayKernel is not an IntegralRayKernel,'
[]

[kernel_not_found]
  type = RunException
  input = 'ray_integral_value_errors.i'
  cli_args = 'RayKernels/active=null Postprocessors/active=kernel_not_found'
  expect_err = "The RayKernel by the name 'dummy' was not found."

  detail = 'the provided IntegralRayKernel was not found,'
[]

[ray_not_found]
  type = RunException
  input = 'ray_integral_value_errors.i'
  cli_args = 'RayKernels/active=variable_integral Postprocessors/active=ray_not_found'
  expect_err = "Could not find a Ray named 'dummy'"

  detail = 'a Ray with the given name was not found,'
[]

[no_registration]
  type = RunException
  input = 'ray_integral_value_errors.i'
  cli_args = 'UserObjects/active=lots RayKernels/active=variable_integral Postprocessors/active=no_registration'
  expect_err = "Cannot be used because the supplied LotsOfRaysRayStudy does not have Ray registration enabled."
  allow_test_objects = true

  detail = 'the provided study does not support Ray registration,'
[]

[no_banking]
  type = RunException
  input = 'ray_integral_value_errors.i'
  cli_args = 'UserObjects/active=no_banking_study RayKernels/active=variable_integral Postprocessors/active=no_banking'
  expect_err = "Cannot be used because the supplied DisableRayBankingStudy does not bank Rays on completion."
  allow_test_objects = true

  detail = 'or the provided study does not support Ray banking.'
[]

[test]
  type = 'CSVDiff'
  input = 'ray_tracing_study_result.i'
  csvdiff = 'ray_tracing_study_result_out.csv'
  allow_test_objects = true

  # ignore the partitioning-specific columns
  override_columns = 'total_processor_crossings max_processor_crossings'
  override_rel_err = '1e12 1e12'
  override_abs_zero = '1e-10 1e-10'

  requirement = 'The system shall support the output of the total number of rays started, the total distance traveled by rays, and the number of processor crossings the rays have encountered.'
[]

[test]
  type = CSVDiff
  input = 'ray_bc_dependencies.i'
  csvdiff = 'ray_bc_dependencies_out.csv'
  allow_test_objects = true

  requirement = 'The system shall support the dependency resolution of RayBCs.'
[]

[missing]
  type = RunException
  input = 'ray_bc_dependencies.i'
  cli_args = 'RayBCs/add_1/depends_on=dummy'
  allow_test_objects = true
  expect_err = 'The RayBoundaryCondition \'dummy\' does not exist'

  requirement = 'The system shall throw a reasonable error when adding a dependency for a RayBC that does not exist.'
[]

[should_not_continue]
  type = RunException
  input = 'raybc_errors.i'
  cli_args = "RayBCs/active='kill change'
                  RayBCs/change/depends_on=kill"
  expect_err = "Cannot changeRayDirection\(\) for a Ray that should not continue"
  allow_test_objects = true

  detail = 'the ray was set to not continue by another RayBC,'
[]

[again]
  type = RunException
  input = 'raybc_errors.i'
  cli_args = "RayBCs/active='change change_again'"
  expect_err = "Cannot change direction for a ray whose direction has already been changed"
  allow_test_objects = true

  detail = 'the ray already had its direction changed,'
[]

[end_set]
  type = RunException
  input = 'raybc_errors.i'
  cli_args = "Mesh/active='gmg internal'
                  UserObjects/active=set_end_study
                  RayBCs/active=change_internal"
  expect_err = "Cannot change the direction of a Ray whose end point is set upon generation"
  allow_test_objects = true

  detail = 'the end point of the ray was set upon generation,'
[]

[zero]
  type = RunException
  input = 'raybc_errors.i'
  cli_args = "RayBCs/active=change_zero"
  expect_err = "Cannot set zero vector direction for a Ray"
  allow_test_objects = true

  detail = 'and when the direction is changed to the zero vector.'
[]

[test]
  type = 'CSVDiff'
  input = 'reflect_ray_bc.i'
  csvdiff = 'reflect_ray_bc_out.csv'
  allow_test_objects = true

  requirement = 'The system shall support the specular reflection of rays on a boundary.'
[]

[non_planar_warning]
  type = RunException
  input = 'reflect_ray_bc_nonplanar.i'
  expect_err = 'A Ray is being reflected on a non-planar side.'
  allow_test_objects = true

  requirement = 'The system shall report a reasonable warning when using an approximation to reflect a ray on a side that is not planar.'
[]

[xyz_only]
  type = RunException
  input = 'errors.i'
  cli_args = 'Problem/coord_type=RZ'
  expect_err = 'Not valid on coordinate systems other than XYZ'
  detail = 'is used with a coordinate system other than xyz or'
[]

[implicit_only]
  type = RunException
  input = 'errors.i'
  cli_args = 'RayKernels/line_source/implicit=false'
  expect_err = 'ADRayKernels do not currently support explicit solves.'
  detail = 'is used in an explicit manner.'
[]

[const_monomial_only]
  type = RunException
  input = 'const_monomial_only.i'
  expect_err = 'Only CONSTANT MONOMIAL variables are supported'
  requirement = 'The system shall throw a reasonable error if a non-supported auxiliary variable is used with a ray aux kernel.'
[]

[test]
  type = 'Exodiff'
  input = 'coupled_line_source_ray_kernel.i'
  exodiff = 'coupled_line_source_ray_kernel_out.e'
  allow_test_objects = true

  requirement = 'The system shall support the coupling of variables on rays that contribute to the residual and Jacobian.'
[]

[test_jac]
  type = 'PetscJacobianTester'
  input = 'coupled_line_source_ray_kernel.i'
  run_sim = True
  ratio_tol = 1e-7
  difference_tol = 1e-6
  allow_test_objects = true

  requirement = 'The system shall be able to produce the exact Jacobian with coupled variables on rays.'
[]

[test_ad]
  type = 'Exodiff'
  input = 'coupled_line_source_ray_kernel.i'
  exodiff = 'coupled_line_source_ray_kernel_out.e'
  cli_args = 'RayKernels/active=source_ad'
  allow_test_objects = true

  requirement = 'The system shall support the coupling of variables using automatic differentiation on rays that contribute to the residual and Jacobian.'
[]

[test_ad_jac]
  type = 'PetscJacobianTester'
  input = 'coupled_line_source_ray_kernel.i'
  run_sim = True
  ratio_tol = 1e-7
  difference_tol = 1e-6
  cli_args = 'RayKernels/active=source_ad'
  allow_test_objects = true

  requirement = 'The system shall be able to produce the exact Jacobian using automatic differentiation on rays.'
[]

[test]
  type = CSVDiff
  input = 'ray_kernel_dependencies.i'
  csvdiff = 'ray_kernel_dependencies_out.csv'
  allow_test_objects = true

  requirement = 'The system shall support the dependency resolution of RayKernels.'
[]

[missing]
  type = RunException
  input = 'ray_kernel_dependencies.i'
  cli_args = 'RayKernels/add_1/depends_on=dummy'
  allow_test_objects = true
  expect_err = 'The RayKernel \'dummy\' does not exist'

  requirement = 'The system shall throw a reasonable error when adding a dependency for a RayKernel that does not exist.'
[]

[should_not_continue]
  type = RunException
  input = 'raykernel_errors.i'
  cli_args = "RayKernels/active='kill change_after_kill'"
  expect_err = "Cannot changeRayStartDirection\(\) for a Ray that should not continue."
  allow_test_objects = true

  detail = 'the ray was set to not continue by another RayKernel,'
[]

[should_not_continue_at_end]
  type = RunException
  input = 'raykernel_errors.i'
  cli_args = "UserObjects/active=end_study
                  RayKernels/active='kill change_after_kill'"
  expect_err = "It has also hit its user-set end point."
  allow_test_objects = true

  detail = 'the ray was set to not continue by another RayKernel and it is at its end point,'
[]

[changed_again]
  type = RunException
  input = 'raykernel_errors.i'
  cli_args = "RayKernels/active='change change_again'"
  expect_err = "Cannot change a Ray\'s trajectory when its trajectory has already been changed"
  allow_test_objects = true

  detail = 'the ray\'s trajectory was already changed,'
[]

[end_set]
  type = RunException
  input = 'raykernel_errors.i'
  cli_args = "UserObjects/active=end_study
                  UserObjects/end_study/end_points='1 0 0'
                  RayKernels/active='change'"
  expect_err = "Cannot change the direction of a Ray whose end point is set upon generation \(via setStartingEndPoint\(\)\)."
  allow_test_objects = true

  detail = 'the ray\'s end point has been set,'
[]

[out_of_elem]
  type = RunException
  input = 'raykernel_errors.i'
  cli_args = "RayKernels/active='change_out_of_elem'"
  expect_err = "is not within the element that the RayKernel was executed on"
  allow_test_objects = true

  detail = 'the ray\'s start point is not within the current element, and'
[]

[zero_direction]
  type = RunException
  input = 'raykernel_errors.i'
  cli_args = "RayKernels/active='change_zero'"
  expect_err = "Cannot set zero vector direction for a Ray"
  allow_test_objects = true

  detail = 'when the ray\'s new direction is set to the zero vector.'
[]

[test]
  type = 'CSVDiff'
  input = 'function_integral_ray_kernel.i'
  csvdiff = 'function_integral_ray_kernel_out.csv'
  requirement = 'The system shall provide the ability to compute the integral of a Function along a line.'
[]

[average]
  type = 'CSVDiff'
  input = 'function_integral_ray_kernel.i'
  csvdiff = 'function_integral_ray_kernel_average.csv'
  cli_args = 'Outputs/file_base=function_integral_ray_kernel_average
                RayKernels/function_integral/average=true'
  requirement = 'The system shall provide the ability to compute the average value of a Function along a line.'
[]

[test]
  type = 'Exodiff'
  input = 'line_source_ray_kernel.i'
  exodiff = 'line_source_ray_kernel_out.e'
  requirement = 'The system shall provide an ability to produce a source term along a line.'
[]

[fe]
  type = 'Exodiff'
  input = 'simple_diffusion_line_source.i'
  exodiff = 'simple_diffusion_line_source_out.e'
  detail = 'finite elements '
[]

[fv]
  type = 'Exodiff'
  input = 'fv_simple_diffusion_line_source.i'
  exodiff = 'fv_simple_diffusion_line_source_out.e'
  detail = 'and finite volumes.'
[]

[test]
  type = 'CSVDiff'
  input = 'material_integral_ray_kernel.i'
  csvdiff = 'material_integral_ray_kernel_out.csv'

  requirement = 'The system shall provide the ability to compute the integral of a material property along a line.'
[]

[average]
  type = 'CSVDiff'
  input = 'material_integral_ray_kernel.i'
  csvdiff = 'material_integral_ray_kernel_average.csv'
  cli_args = 'Outputs/file_base=material_integral_ray_kernel_average
                RayKernels/material_integral/average=true'

  requirement = 'The system shall provide the ability to compute the average value of a material property along a line.'
[]

[test]
  type = 'Exodiff'
  input = 'ray_distance_aux.i'
  exodiff = 'ray_distance_aux_out.e'
  allow_test_objects = true
  requirement = 'The system shall be able to store the distance traversed by a ray in each element in an auxiliary field.'
[]

[xyz_only]
  type = RunException
  input = 'errors.i'
  cli_args = 'Problem/coord_type=RZ'
  expect_err = 'Not valid on coordinate systems other than XYZ'
  detail = 'is utilized with a coordinate system other than xyz or'
[]

[bad_execute_on]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/execute_on=initial'
  expect_err = 'In this case, the study must use the execute_on = PRE_KERNELS'
  detail = 'if the execute flag for the associated ray study is not set correctly to use the RayKernel.'
[]

[test]
  type = 'CSVDiff'
  input = 'variable_integral_ray_kernel.i'
  csvdiff = 'variable_integral_ray_kernel_out.csv'
  design = 'VariableIntegralRayKernel.md'
  requirement = 'The system shall provide the ability to compute the integral of a variable along a line.'
[]

[average]
  type = 'CSVDiff'
  input = 'variable_integral_ray_kernel.i'
  csvdiff = 'variable_integral_ray_kernel_average.csv'
  cli_args = 'Outputs/file_base=variable_integral_ray_kernel_average
                RayKernels/variable_integral/average=true
                RayKernels/aux_variable_integral/average=true'
  design = 'VariableIntegralRayKernel.md'
  requirement = 'The system shall provide the ability to compute the average value of a variable along a line.'
[]

[simple_diffusion_rays]
  type = 'Exodiff'
  input = 'simple_diffusion_line_integral.i'
  exodiff = 'simple_diffusion_line_integral_rays.e'
  design = 'VariableIntegralRayKernel.md RayTracingMeshOutput.md'
  cli_args = 'UserObjects/study/always_cache_traces=true
                UserObjects/study/data_on_cache_traces=true
                Outputs/rays/output_data=true
                Outputs/rays/execute_on=TIMESTEP_END'
  requirement = 'The system shall provide the ability to output the segment wise accumulated integral of a variable along a line with a 2D, finite-element diffusion problem.'
[]

[simple_diffusion_rays_nodal]
  type = 'Exodiff'
  input = 'simple_diffusion_line_integral.i'
  exodiff = 'simple_diffusion_line_integral_rays_nodal.e'
  design = 'VariableIntegralRayKernel.md RayTracingMeshOutput.md'
  cli_args = 'UserObjects/study/always_cache_traces=true
                UserObjects/study/data_on_cache_traces=true
                Outputs/rays/execute_on=TIMESTEP_END
                Outputs/rays/output_data=true
                Outputs/rays/output_data_nodal=true
                Outputs/rays/file_base=simple_diffusion_line_integral_rays_nodal'
  requirement = 'The system shall provide the ability to output the segment wise, linearly approximated, accumulated integral of a variable along a line with a 2D, finite-element diffusion problem.'
[]

[fe]
  type = 'CSVDiff'
  input = 'simple_diffusion_line_integral.i'
  csvdiff = 'simple_diffusion_line_integral_out.csv'
  detail = 'finite elements '
[]

[fv]
  type = 'CSVDiff'
  input = 'fv_simple_diffusion_line_integral.i'
  csvdiff = 'fv_simple_diffusion_line_integral_out.csv'
  detail = 'and finite volumes.'
[]

[get_info]
  type = 'RunApp'
  input = 'ray_lots.i'
  allow_test_objects = true
  cli_args = 'UserObjects/study/get_info=true'
  requirement = 'The system shall support the console output of information pertaining to a Ray.'
[]

[equality]
  type = RunApp
  input = 'ray_lots.i'
  allow_test_objects = true
  cli_args = "UserObjects/study/equality=true
                UserObjects/other_study/type=RepeatableRayStudy
                UserObjects/other_study/start_points='0 0 0'
                UserObjects/other_study/end_points='0.1 0.1 0'
                UserObjects/other_study/names=ray
                UserObjects/other_study/ray_kernel_coverage_check=false"
  requirement = 'The system shall provide a method for determining the equality of rays'
  expect_out = 'Ray equality test: same_passes = 408, duplicate_passes = 7344'
[]

[at_end_without_set]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/at_end_without_set=true'
  expect_err = 'Cannot use Ray::atEnd\(\) for a Ray that does not have an end set'
  allow_test_objects = true
  detail = 'checking if a ray is at its end point without an end point being set'
[]

[end_point_without_set]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/end_point_without_set=true'
  expect_err = 'Cannot use Ray::endPoint\(\) for a Ray that does not have an end set'
  allow_test_objects = true
  detail = 'retrieving a ray\'s end point when its end point has not been set'
[]

[set_start_again]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_start_again=true'
  expect_err = 'Starting point was already set via Ray::setStart\(\) and is being changed.'
  allow_test_objects = true
  detail = 'setting a ray\'s start point that has already been set'
[]

[set_direction_again]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_direction_again=true'
  expect_err = 'Cannot change a Ray\'s starting direction using Ray::setStartingDirection\(\)'
  allow_test_objects = true
  detail = 'setting a ray\'s starting direction multiple times without clearing it'
[]

[set_start_fail_bbox]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_start_fail_bbox=true'
  expect_err = 'Mesh does not contain starting point.'
  allow_test_objects = true
  detail = 'setting a ray\'s start point that is not contained within the mesh bounding box'
[]

[set_side_without_elem]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_side_without_elem=true'
  expect_err = 'Starting incoming side is set but starting element is not set.'
  allow_test_objects = true
  detail = 'setting a ray\'s starting incoming side without setting its starting element'
[]

[set_invalid_side]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_invalid_side=true'
  expect_err = 'Starting incoming side is not valid for its starting element.'
  allow_test_objects = true
  detail = 'setting a ray\'s starting incoming side to an invalid side'
[]

[set_bad_side]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_bad_side=true'
  expect_err = 'Starting incoming side does not contain the starting point.'
  allow_test_objects = true
  detail = 'setting a ray\'s starting incoming side to one that does not contain the starting point'
[]

[set_bad_start]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_bad_start=true'
  expect_err = 'Starting element does not contain the starting point.'
  allow_test_objects = true
  detail = 'setting a ray\'s starting point to a point that is not contained within a given starting element'
[]

[set_direction_before_start]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_direction_before_start=true'
  expect_err = 'Cannot use Ray::setStartingDirection\(\) before Ray::setStart\(\).'
  allow_test_objects = true
  detail = 'setting a ray\'s starting direction before setting its starting information,'
[]

[set_zero_direction]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_zero_direction=true'
  expect_err = 'Starting direction in Ray::setStartingDirection\(\) is the zero vector.'
  allow_test_objects = true
  detail = 'setting a ray\'s starting direction to the zero vector'
[]

[set_end_before_start]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_end_before_start=true'
  expect_err = 'Cannot use Ray::setStartingEndPoint\(\) before Ray::setStart\(\).'
  allow_test_objects = true
  detail = 'setting a ray\'s end point before setting its start point'
[]

[set_end_equal_start]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_end_equal_start=true'
  expect_err = 'End point is equal to the start point in Ray::setStartingEndPoint\(\).'
  allow_test_objects = true
  detail = 'setting a ray\'s end point to the same as its start point'
[]

[set_end_with_direction]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_end_with_direction=true'
  expect_err = 'Cannot use Ray::setStartingEndPoint\(\) after Ray::setStartingDirection\(\).'
  allow_test_objects = true
  detail = 'setting a ray\'s end point and starting direction'
[]

[set_distance_with_end]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_distance_with_end=true'
  expect_err = 'Cannot use Ray::setStartingMaxDistance\(\) after Ray::setStartingEndPoint\(\).'
  allow_test_objects = true
  detail = 'setting a ray\'s max distance after setting its end point'
[]

[set_end_with_distance]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_end_with_distance=true'
  expect_err = 'Cannot use Ray::setStartingEndPoint\(\) after Ray::setStartingMaxDistance\(\)'
  allow_test_objects = true
  detail = 'setting a ray\'s end point after setting its max distance'
[]

[set_end_fail_bbox]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_end_fail_bbox=true'
  expect_err = 'End point is not within the mesh for Ray::setStartingEndPoint\(\).'
  allow_test_objects = true
  detail = 'setting a ray\'s end point to a point that is not within the mesh bounding box'
[]

[set_distance_before_start]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_distance_before_start=true'
  expect_err = 'Cannot use Ray::setStartingMaxDistance\(\) before Ray::setStart\(\).'
  allow_test_objects = true
  detail = 'setting a ray\'s starting maximum distance before setting its start information'
[]

[set_distance_negative]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_distance_negative=true'
  expect_err = 'Starting max distance is <= 0'
  allow_test_objects = true
  detail = 'setting a ray\'s starting maximum distance to a negative value'
[]

[set_start_inactive]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_start_inactive=true
                  Adaptivity/steps=1
                  Adaptivity/marker=uniform
                  Adaptivity/max_h_level=2
                  Adaptivity/initial_steps=2
                  Adaptivity/Markers/uniform/type=UniformMarker
                  Adaptivity/Markers/uniform/mark=refine'
  expect_err = 'Starting element is not active'
  allow_test_objects = true
  detail = 'setting the starting element for a ray to an inactive element'
[]

[test]
  type = 'Exodiff'
  input = 'stationary.i'
  exodiff = 'stationary_out.e'
  allow_test_objects = true
  requirement = 'The system shall support the tracing of stationary rays'
[]

[stationary_before_start]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_stationary_before_start=true'
  expect_err = 'Cannot use Ray::setStationary\(\) before Ray::setStart\(\)'
  allow_test_objects = true
  detail = 'its starting point has not been set yet'
[]

[stationary_with_direction]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_stationary_with_direction=true'
  expect_err = 'Cannot use Ray::setStationary\(\) with Ray::setStartingDirection\(\)'
  allow_test_objects = true
  detail = 'its starting direction has already been set'
[]

[stationary_with_end]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/set_stationary_with_end=true'
  expect_err = 'Cannot use Ray::setStationary\(\) after Ray::setStartingEndPoint\(\)'
  allow_test_objects = true
  detail = 'its starting end point has already been set'
[]

[quad4]
  type = 'CSVDiff'
  input = 'adaptivity_2d.i'
  cli_args = 'Mesh/gmg/elem_type=QUAD4 Outputs/file_base=adaptivity_quad4_out'
  csvdiff = 'adaptivity_quad4_out.csv'
  allow_test_objects = true
  detail = 'QUAD4, '
[]

[tri3]
  type = 'CSVDiff'
  input = 'adaptivity_2d.i'
  csvdiff = 'adaptivity_tri3_out.csv'
  cli_args = 'Mesh/gmg/elem_type=TRI3 Outputs/file_base=adaptivity_tri3_out'
  allow_test_objects = true
  detail = 'TRI3, '
[]

[hex8]
  type = 'CSVDiff'
  input = 'adaptivity_3d.i'
  csvdiff = 'adaptivity_hex8_out.csv'
  cli_args = 'Mesh/gmg/elem_type=HEX8 Outputs/file_base=adaptivity_hex8_out'
  allow_test_objects = true
  detail = 'HEX8, '
[]

[edge2]
  type = 'CSVDiff'
  input = 'adaptivity_1d.i'
  csvdiff = 'adaptivity_edge2_out.csv'
  cli_args = 'Mesh/gmg/elem_type=EDGE2 Outputs/file_base=adaptivity_edge2_out'
  allow_test_objects = true
  detail = 'and EDGE2.'
[]

[2d]
  type = 'CSVDiff'
  input = 'backface_culling.i'
  csvdiff = 'backface_culling_2d_out.csv'
  abs_zero = 1e-8
  cli_args = 'Mesh/gmg/dim=2
                  RayBCs/active=kill_2d
                  Outputs/file_base=backface_culling_2d_out'
  allow_test_objects = true
  detail = 'two-dimensional meshes '
[]

[3d]
  type = 'CSVDiff'
  input = 'backface_culling.i'
  csvdiff = 'backface_culling_3d_out.csv'
  abs_zero = 1e-8
  cli_args = 'Mesh/gmg/dim=3
                  RayBCs/active=kill_3d
                  Outputs/file_base=backface_culling_3d_out
                  UserObjects/study/edge_to_edge=true'
  allow_test_objects = true
  detail = 'and for three-dimensional meshes.'
[]

[test]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'lots_out.csv'
  cli_args = 'Mesh/gmg/elem_type=edge2
                Mesh/gmg/dim=1
                RayBCs/active=kill_1d
                UserObjects/lots/vertex_to_vertex=true
                UserObjects/lots/centroid_to_vertex=true
                UserObjects/lots/centroid_to_centroid=true
                Outputs/file_base=lots_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'
  requirement = 'The system shall be able to trace rays within EDGE2 elements.'
[]

[test]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'lots_out.csv'
  cli_args = 'Mesh/gmg/elem_type=edge3
                Mesh/gmg/dim=1
                RayBCs/active=kill_1d
                UserObjects/lots/vertex_to_vertex=true
                UserObjects/lots/centroid_to_vertex=true
                UserObjects/lots/centroid_to_centroid=true
                Outputs/file_base=lots_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'
  requirement = 'The system shall be able to trace rays within EDGE3 elements.'
[]

[test]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'lots_out.csv'
  cli_args = 'Mesh/gmg/elem_type=edge4
                Mesh/gmg/dim=1
                RayBCs/active=kill_1d
                UserObjects/lots/vertex_to_vertex=true
                UserObjects/lots/centroid_to_vertex=true
                UserObjects/lots/centroid_to_centroid=true
                Outputs/file_base=lots_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'
  requirement = 'The system shall be able to trace rays within EDGE4 elements.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ntn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex20
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=ntn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex20
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex20
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[ete]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ete_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex20
                  UserObjects/lots/edge_to_edge=true
                  Outputs/file_base=ete_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary edges in the direction of all other edge centroids in the same element, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex20
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex20
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ntn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex27
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=ntn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex27
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex27
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[ete]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ete_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex27
                  UserObjects/lots/edge_to_edge=true
                  Outputs/file_base=ete_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary edges in the direction of all other edge centroids in the same element, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex27
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex27
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ntn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex8
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=ntn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex8
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex8
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[ete]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ete_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex8
                  UserObjects/lots/edge_to_edge=true
                  Outputs/file_base=ete_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary edges in the direction of all other edge centroids in the same element, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex8
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=hex8
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[1d]
  type = 'CSVDiff'
  input = 'internal_sidesets_1d.i'
  csvdiff = 'internal_sidesets_1d_kill_out.csv'
  cli_args = 'Outputs/file_base=internal_sidesets_1d_kill_out
                  RayBCs/active=kill_internal'
  allow_test_objects = true

  detail = 'one-dimensional meshes, '
[]

[2d]
  type = 'CSVDiff'
  input = 'internal_sidesets_2d.i'
  csvdiff = 'internal_sidesets_2d_kill_out.csv'
  cli_args = 'Outputs/file_base=internal_sidesets_2d_kill_out
                  RayBCs/active=kill_internal'
  allow_test_objects = true

  detail = 'two-dimensional meshes, '
[]

[3d]
  type = 'CSVDiff'
  input = 'internal_sidesets_3d.i'
  csvdiff = 'internal_sidesets_3d_kill_out.csv'
  cli_args = 'Outputs/file_base=internal_sidesets_3d_kill_out
                  RayBCs/active=kill_internal'
  allow_test_objects = true

  detail = 'and three-dimensional meshes.'
[]

[1d]
  type = 'CSVDiff'
  input = 'internal_sidesets_1d.i'
  csvdiff = 'internal_sidesets_1d_reflect_out.csv'
  cli_args = 'Outputs/file_base=internal_sidesets_1d_reflect_out
                  RayBCs/active="kill_external reflect_internal"'
  allow_test_objects = true

  detail = 'one-dimensional meshes, '
[]

[2d]
  type = 'CSVDiff'
  input = 'internal_sidesets_2d.i'
  csvdiff = 'internal_sidesets_2d_reflect_out.csv'
  cli_args = 'Outputs/file_base=internal_sidesets_2d_reflect_out
                  RayBCs/active="kill_external reflect_internal"'
  allow_test_objects = true

  detail = 'two-dimensional meshes, '
[]

[3d]
  type = 'CSVDiff'
  input = 'internal_sidesets_3d.i'
  csvdiff = 'internal_sidesets_3d_reflect_out.csv'
  cli_args = 'Outputs/file_base=internal_sidesets_3d_reflect_out
                  RayBCs/active="kill_external reflect_internal"'
  allow_test_objects = true

  detail = 'and three-dimensional meshes.'
[]

[internal_sidesets_disabled]
  type = RunException
  input = 'internal_sidesets_1d.i'
  cli_args = 'UserObjects/study/use_internal_sidesets=false
                RayBCs/active=kill_internal'
  expect_err = 'RayBCs are defined on internal sidesets, but the study is not set to use internal sidesets during tracing.'
  allow_test_objects = true

  requirement = 'The system shall throw an error if RayBCs are defined on internal sidesets, but the study does not have internal sideset capability enabled.'
[]

[not_subdomain_bounded]
  type = RunException
  input = 'internal_sidesets_2d.i'
  cli_args = 'Mesh/change_subdomain/type=RenameBlockGenerator
                Mesh/change_subdomain/input=interior_top
                Mesh/change_subdomain/old_block=5
                Mesh/change_subdomain/new_block=1'
  expect_err = 'RayBCs exist on internal sidesets that are not bounded by a different'
  allow_test_objects = true

  requirement = 'The system shall throw an error if RayBCs are defined on internal sidesets and said internal sidesets are not bounded by different subdomains.'
[]

[test]
  type = 'CSVDiff'
  input = 'nonplanar.i'
  csvdiff = 'nonplanar_out.csv'
  allow_test_objects = true
  abs_zero = 1e-8
  requirement = 'The system shall be able to trace rays on 3D meshes that have nonplanar faces.'
[]

[warning]
  type = RunException
  input = 'nonplanar.i'
  cli_args = 'UserObjects/lots/warn_non_planar=true'
  expect_err = 'Ray tracing on non-planar faces is an approximation and may fail.'
  allow_test_objects = true
  requirement = 'The system shall by default warn that tracing on non-planar faces is an approximation.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ntn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=prism15
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=ntn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=prism15
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=4
                  Mesh/gmg/ny=4
                  Mesh/gmg/nz=4
                  Mesh/gmg/elem_type=prism15
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[ete]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ete_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=prism15
                  UserObjects/lots/edge_to_edge=true
                  Outputs/file_base=ete_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary edges in the direction of all other edge centroids in the same element, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=prism15
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=prism15
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ntn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=prism18
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=ntn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=prism18
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=4
                  Mesh/gmg/ny=4
                  Mesh/gmg/nz=4
                  Mesh/gmg/elem_type=prism18
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[ete]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ete_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=prism18
                  UserObjects/lots/edge_to_edge=true
                  Outputs/file_base=ete_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary edges in the direction of all other edge centroids in the same element, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=prism18
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=prism18
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ntn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=prism6
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=ntn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=prism6
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=4
                  Mesh/gmg/ny=4
                  Mesh/gmg/nz=4
                  Mesh/gmg/elem_type=prism6
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[ete]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ete_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=prism6
                  UserObjects/lots/edge_to_edge=true
                  Outputs/file_base=ete_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary edges in the direction of all other edge centroids in the same element, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=prism6
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=prism6
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ntn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=pyramid13
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=ntn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=pyramid13
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=4
                  Mesh/gmg/ny=4
                  Mesh/gmg/ny=4
                  Mesh/gmg/elem_type=pyramid13
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[ete]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ete_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=pyramid13
                  UserObjects/lots/edge_to_edge=true
                  Outputs/file_base=ete_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary edges in the direction of all other edge centroids in the same element, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=pyramid13
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=pyramid13
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ntn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=pyramid14
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=ntn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=pyramid14
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=4
                  Mesh/gmg/ny=4
                  Mesh/gmg/ny=4
                  Mesh/gmg/elem_type=pyramid14
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[ete]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ete_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=pyramid14
                  UserObjects/lots/edge_to_edge=true
                  Outputs/file_base=ete_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary edges in the direction of all other edge centroids in the same element, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=pyramid14
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=pyramid14
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ntn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=pyramid5
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=ntn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctn_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=pyramid5
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctn_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=4
                  Mesh/gmg/ny=4
                  Mesh/gmg/ny=4
                  Mesh/gmg/elem_type=pyramid5
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[ete]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ete_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=pyramid5
                  UserObjects/lots/edge_to_edge=true
                  Outputs/file_base=ete_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary edges in the direction of all other edge centroids in the same element, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=pyramid5
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/elem_type=pyramid5
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'vtv_out.csv'
  cli_args = 'Mesh/gmg/elem_type=quad4
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=vtv_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctv_out.csv'
  cli_args = 'Mesh/gmg/elem_type=quad4
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctv_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/elem_type=quad4
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/elem_type=quad4
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/elem_type=quad4
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'vtv_out.csv'
  cli_args = 'Mesh/gmg/elem_type=quad8
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=vtv_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctv_out.csv'
  cli_args = 'Mesh/gmg/elem_type=quad8
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctv_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/elem_type=quad8
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/elem_type=quad8
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/elem_type=quad8
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'vtv_out.csv'
  cli_args = 'Mesh/gmg/elem_type=quad9
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=vtv_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctv_out.csv'
  cli_args = 'Mesh/gmg/elem_type=quad9
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctv_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/elem_type=quad9
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/elem_type=quad9
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/elem_type=quad9
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[both]
  type = RunException
  input = 'raybc_check.i'
  expect_err = 'Boundaries that did not have any RayBCs\:.*1 \(right\).*2 \(top\)'
  detail = 'the boundary does not have any RayBCs '
[]

[one]
  type = RunException
  input = 'raybc_check.i'
  expect_err = 'Boundaries that did not have any RayBCs\:.*1 \(right\)'
  cli_args = 'RayBCs/active=top'
  detail = 'and then the boundary does not have any RayBCs that change the trajectory of the ray or end the ray.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'vtv_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet10
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=vtv_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctv_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet10
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctv_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet10
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[ete]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ete_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet10
                  UserObjects/lots/edge_to_edge=true
                  Outputs/file_base=ete_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary edges in the direction of all other edge centroids in the same element, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet10
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet10
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'vtv_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet14
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=vtv_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctv_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet14
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctv_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet14
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[ete]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ete_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet14
                  UserObjects/lots/edge_to_edge=true
                  Outputs/file_base=ete_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary edges in the direction of all other edge centroids in the same element, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet14
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet14
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'vtv_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet4
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=vtv_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctv_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet4
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctv_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet4
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[ete]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ete_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet4
                  UserObjects/lots/edge_to_edge=true
                  Outputs/file_base=ete_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary edges in the direction of all other edge centroids in the same element, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet4
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/dim=3
                  Mesh/gmg/nx=3
                  Mesh/gmg/ny=3
                  Mesh/gmg/nz=3
                  Mesh/gmg/elem_type=tet4
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out
                  RayBCs/active=kill_3d'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'vtv_out.csv'
  cli_args = 'Mesh/gmg/elem_type=tri3
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=vtv_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctv_out.csv'
  cli_args = 'Mesh/gmg/elem_type=tri3
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctv_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/elem_type=tri3
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/elem_type=tri3
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/elem_type=tri3
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'vtv_out.csv'
  cli_args = 'Mesh/gmg/elem_type=tri6
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=vtv_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctv_out.csv'
  cli_args = 'Mesh/gmg/elem_type=tri6
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctv_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/elem_type=tri6
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/elem_type=tri6
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/elem_type=tri6
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[vtv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'vtv_out.csv'
  cli_args = 'Mesh/gmg/elem_type=tri7
                  UserObjects/lots/vertex_to_vertex=true
                  Outputs/file_base=vtv_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary vertices on boundary sides in the direction of the vertices on the other side of the same boundary element, '
[]

[ctv]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctv_out.csv'
  cli_args = 'Mesh/gmg/elem_type=tri7
                  UserObjects/lots/centroid_to_vertex=true
                  Outputs/file_base=ctv_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of boundary sides in the direction of all other vertices in the same boundary element, '
[]

[ctc]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'ctc_out.csv'
  cli_args = 'Mesh/gmg/elem_type=tri7
                  UserObjects/lots/centroid_to_centroid=true
                  Outputs/file_base=ctc_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from the centroids of all boundary sides to the centroids of all other boundary elements, '
[]

[side_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'side_aq_out.csv'
  cli_args = 'Mesh/gmg/elem_type=tri7
                  UserObjects/lots/side_aq=true
                  Outputs/file_base=side_aq_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'from boundary sides centroids using an angular quadrature, '
[]

[centroid_aq]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'centroid_aq_out.csv'
  cli_args = 'Mesh/gmg/elem_type=tri7
                  UserObjects/lots/centroid_aq=true
                  Outputs/file_base=centroid_aq_out'
  allow_test_objects = true
  abs_zero = 1e-8
  method = '!dbg'

  detail = 'and from boundary element centroids using an angular quadrature.'
[]

[test]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'unset_incoming_side_out.csv'
  cli_args = 'Mesh/gmg/elem_type=quad4
                UserObjects/lots/vertex_to_vertex=true
                UserObjects/lots/centroid_to_vertex=true
                UserObjects/lots/centroid_to_centroid=true
                UserObjects/lots/side_aq=true
                UserObjects/lots/set_incoming_side=false
                Outputs/file_base=unset_incoming_side_out'
  allow_test_objects = true
  abs_zero = 1e-8
  requirement = 'The system shall support not setting the incoming side at the beginning of a trace.'
[]

[2d]
  type = 'Exodiff'
  input = 'cone_ray_study.i'
  exodiff = 'cone_ray_study_out.e cone_ray_study_rays.e'
  detail = 'two-dimensional meshes '
[]

[3d]
  type = 'Exodiff'
  input = 'cone_ray_study_3d.i'
  exodiff = 'cone_ray_study_3d_out.e cone_ray_study_3d_rays.e'
  detail = 'and in three-dimensional meshes.'
[]

[direction_size]
  type = 'RunException'
  input = 'cone_ray_study_errors.i'
  cli_args = "UserObjects/study/directions='2 1 0 1 0 0'"
  expect_err = 'Not the same size as start_points'
  detail = 'the provided directions are inconsistent in size'
[]

[scaling_factors_size]
  type = 'RunException'
  input = 'cone_ray_study_errors.i'
  cli_args = "UserObjects/study/scaling_factors='1 1'"
  expect_err = 'Not the same size as start_points'
  detail = 'the provided scaling factors are inconsistent in size'
[]

[half_cone_angles_size]
  type = 'RunException'
  input = 'cone_ray_study_errors.i'
  cli_args = "UserObjects/study/half_cone_angles='1 1'"
  expect_err = 'Not the same size as start_points'
  detail = 'the provided half cone angles are inconsistent in size'
[]

[half_cone_angles_valid]
  type = 'RunException'
  input = 'cone_ray_study_errors.i'
  cli_args = "UserObjects/study/half_cone_angles='-10'"
  expect_err = 'Must be > 0 and <= 90 degrees'
  detail = 'the provided half cone angles are not valid angles'
[]

[polar_quad_orders_size]
  type = 'RunException'
  input = 'cone_ray_study_errors.i'
  cli_args = "UserObjects/study/polar_quad_orders='1 2'"
  expect_err = 'Not the same size as start_points'
  detail = 'the provided polar quadrature orders are inconsistent in size'
[]

[azimuthal_quad_orders_size]
  type = 'RunException'
  input = 'cone_ray_study_errors.i'
  cli_args = "UserObjects/study/azimuthal_quad_orders='1 2'"
  expect_err = 'Not the same size as start_points'
  detail = 'the provided azimuthal quadrature orders are inconsistent in size'
[]

[1d]
  type = 'RunException'
  input = 'cone_ray_study_errors.i'
  cli_args = 'Mesh/active=1d Mesh/1d/type=GeneratedMeshGenerator Mesh/1d/dim=1'
  expect_err = 'Does not support 1D'
  detail = 'it is run with an invalid mesh dimension'
[]

[2d]
  type = 'CSVDiff'
  input = 'bc_create_ray.i'
  csvdiff = 'bc_create_ray_2d_out.csv'
  allow_test_objects = true
  cli_args = 'Outputs/file_base=bc_create_ray_2d_out'
  detail = 'two-dimensional meshes '
[]

[3d]
  type = 'CSVDiff'
  input = 'bc_create_ray.i'
  csvdiff = 'bc_create_ray_3d_out.csv'
  allow_test_objects = true
  cli_args = 'Outputs/file_base=bc_create_ray_3d_out
                  Mesh/active=gmg_3d
                  RayBCs/active="kill_3d create_3d"'
  detail = 'and in three-dimensional meshes.'
[]

[lifo]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'lots_out.csv'
  cli_args = 'UserObjects/lots/work_buffer_type=lifo
                  Outputs/file_base=lots_out'
  allow_test_objects = true
  detail = 'last in, first out buffer '
[]

[circular]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'lots_out.csv'
  cli_args = 'UserObjects/lots/work_buffer_type=circular
                  Outputs/file_base=lots_out'
  allow_test_objects = true
  detail = 'and with a circular buffer.'
[]

[duplicate_ray]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/study/add_duplicate_ray=true'
  allow_test_objects = true
  expect_err = 'Multiple shared_ptrs were found that point to the same Ray'
  detail = 'the same ray is added to be traced multiple times,'
[]

[local_non_unique_id]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/study/add_local_non_unique_id_ray=true
                  Mesh/gmg/nx=50'
  allow_test_objects = true
  expect_err = 'Multiple Rays exist with ID 0 on processor'
  detail = 'when local rays that have the same ID are added to be traced,'
[]

[global_non_unique_id]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/study/add_global_non_unique_id_ray=true
                  Mesh/gmg/nx=50'
  allow_test_objects = true
  min_parallel = 2
  expect_err = 'Ray with ID 0 exists on ranks'
  detail = 'when rays that have the same ID on different processors are added to be traced,'
[]

[ray_kernel_coverage_check]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/ray_kernel_coverage_check=true'
  expect_err = 'Subdomains { 0 } do not have RayKernels defined!'
  allow_test_objects = true
  detail = 'when the domain is not covered with ray kernels,'
[]

[unneeded_rays]
  type = RunException
  input = 'errors.i'
  cli_args = 'RayKernels/null/type=NullRayKernel RayKernels/null/rays=dummy'
  expect_err = 'Rays cannot be supplied when the study does not require Ray registration'
  allow_test_objects = true
  detail = 'when ray names are provided to a ray tracing object and the study does not register ray names,'
[]

[ray_not_found]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/active=repeatable
                  RayKernels/null/type=NullRayKernel
                  RayKernels/null/rays=dummy'
  expect_err = 'Supplied ray \'dummy\' is not a registered Ray'
  allow_test_objects = true
  detail = 'when a ray name is provided to a ray tracing object and said name is not registered to a ray,'
[]

[mixed_pre_kernel_execute]
  type = RunException
  input = 'errors.i'
  cli_args = "UserObjects/study/execute_on='INITIAL PRE_KERNELS'"
  expect_err = 'PRE_KERNELS cannot be mixed with any other execution flag.'
  allow_test_objects = true
  detail = 'when the study is set to execute with residual kernels and non-residual kernels,'
[]

[eigenvalue]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/execute_on=PRE_KERNELS
                  Executioner/type=Eigenvalue'
  expect_err = 'is not supported for an eigenvalue solve.'
  allow_test_objects = true
  slepc = true
  detail = 'when the study is set to execute with residual kernels and an eigenvalue executioner,'
[]

[error_if_tracing]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/ray_error_if_tracing=true'
  expect_err = 'Cannot use Ray::clearStartingInfo\(\) after it has started tracing'
  allow_test_objects = true
  detail = 'when the trajectory of a ray is changed after it has started tracing,'
[]

[reset_counters]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/ray_reset_counters=true'
  expect_err = 'Ray::resetCounters\(\) can only be used during generateRays()'
  allow_test_objects = true
  detail = 'when a ray\'s counters are reset when not generating rays'
[]

[adapt_bad_type]
  type = RunException
  input = 'errors.i'
  cli_args = 'Mesh/gmg/elem_type=EDGE3
                  Adaptivity/steps=1
                  Adaptivity/marker=uniform
                  Adaptivity/max_h_level=2
                  Adaptivity/initial_steps=1
                  Adaptivity/Markers/uniform/type=UniformMarker
                  Adaptivity/Markers/uniform/mark=refine'
  allow_test_objects = true
  expect_err = 'Element type EDGE3 is not supported in ray tracing with adaptivity'
  detail = 'when tracing in a mesh with adaptivity and an unsupported element type'
[]

[hmax_vary]
  type = RunException
  input = 'errors.i'
  cli_args = 'Mesh/active=2d
                  Mesh/2d/type=GeneratedMeshGenerator
                  Mesh/2d/dim=2
                  Mesh/2d/xmax=1000'
  allow_test_objects = true
  expect_err = 'Element hmin varies significantly from subdomain hmax.'
  detail = 'when tracing in a mesh when the average subdomain hmax is a poor approximation due to stretched elements'
[]

[register_ray_data_late]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/register_ray_data_late=true'
  expect_err = 'Cannot register Ray data after initialSetup()'
  allow_test_objects = true
  detail = 'when ray data is registered too late'
[]

[register_ray_data_same_name]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/register_ray_data_same_name=true'
  expect_err = 'Cannot register Ray aux data with name foo because Ray \(non-aux\) data already exists with said name.'
  allow_test_objects = true
  detail = 'when ray data is registered with the same name across systems'
[]

[ray_data_index_bad_name]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/ray_data_index_bad_name=true'
  expect_err = 'Unknown Ray data with name foo'
  allow_test_objects = true
  detail = 'when requesting a ray data index for an invalid name'
[]

[ray_data_name_bad_index]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/ray_data_name_bad_index=true'
  expect_err = 'Unknown Ray data with index 123'
  allow_test_objects = true
  detail = 'when requesting a ray data name for an invalid index'
[]

[get_ray_kernels_early]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/get_ray_kernels_early=true'
  expect_err = 'Should not call getRayKernels\(\) before initialSetup\(\)'
  allow_test_objects = true
  detail = 'when requesting ray kernels before they are setup'
[]

[get_ray_bcs_early]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/get_ray_bcs_early=true'
  expect_err = 'Should not call getRayBCs\(\) before initialSetup\(\)'
  allow_test_objects = true
  detail = 'when requesting ray bcs before they are setup'
[]

[get_bank_disabled]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/type=RayTracingStudyNoBankingTest'
  expect_err = 'The Ray bank is not available because the private parameter \'_bank_rays_on_completion\' is set to false.'
  allow_test_objects = true
  detail = 'when obtaining the ray bank when banking is disabled'
[]

[get_bank_bad]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/get_ray_bank_generating=true'
  expect_err = 'Cannot get the Ray bank during generation or propagation.'
  allow_test_objects = true
  detail = 'when obtaining the ray bank when it is not available'
[]

[register_ray_no_registration]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/register_ray_no_registration=true'
  expect_err = 'Cannot use registerRay\(\) with Ray registration disabled'
  allow_test_objects = true
  detail = 'when registering a ray when ray registration is disabled'
[]

[registered_ray_id_no_registration]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/registered_ray_id_no_registration=true'
  expect_err = 'Should not use registeredRayID\(\) with Ray registration disabled'
  allow_test_objects = true
  detail = 'when getting a registered ray ID when ray registration is disabled'
[]

[registered_ray_name_no_registration]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/registered_ray_name_no_registration=true'
  expect_err = 'Should not use registeredRayName\(\) with Ray registration disabled'
  allow_test_objects = true
  detail = 'when getting a registered ray name when ray registration is disabled'
[]

[registered_ray_id_missing]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/type=RayTracingStudyWithRegistrationTest
                  UserObjects/study/registered_ray_id_missing=true'
  expect_err = 'Attempted to obtain ID of registered Ray foo, but a Ray with said name is not registered.'
  allow_test_objects = true
  detail = 'when requesting a registered ray ID for a ray that is not registered'
[]

[registered_ray_name_missing]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/type=RayTracingStudyWithRegistrationTest
                  UserObjects/study/registered_ray_name_missing=true'
  expect_err = 'Attempted to obtain name of registered Ray with ID 0, but a Ray with said ID is not registered.'
  allow_test_objects = true
  detail = 'when requesting a registered ray name for a ray that is not registered'
[]

[reserve_bad]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/reserve_bad=true'
  expect_err = 'Can only reserve in Ray buffer during generateRays\(\)'
  allow_test_objects = true
  detail = 'when reserving space in the ray buffer outside of registration, and'
[]

[subdomain_hmax_missing]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/subdomain_hmax_missing=true'
  expect_err = 'Subdomain 1337 not found in subdomain hmax map'
  allow_test_objects = true
  detail = 'when requesting a subdomain hmax for an invalid subdomain'
[]

[elem_extrema_build_edge]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/elem_extrema_build_edge=true'
  expect_err = 'Element does not contain vertices in ElemExtrema'
  allow_test_objects = true
  # will hit assertions before the error outside of opt
  method = opt
  detail = 'building an element extrema edge for invalid vertices'
[]

[get_elem_normals_unimplemented]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/get_elem_normals_unimplemented=true'
  expect_err = 'Unimplemented element normal caching in RayTracingStudyTest::getElemNormals\(\)'
  allow_test_objects = true
  detail = 'requesting cached element normals when element normal caching has not been setup'
[]

[ray_data_index_other_exists]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/ray_data_index_other_exists=true'
  expect_err = 'However, Ray aux data with said name was found.'
  allow_test_objects = true
  detail = 'getting the index of ray data when aux data with the same name exists'
[]

[ray_data_aux_index_other_exists]
  type = RunException
  input = 'errors.i'
  cli_args = 'UserObjects/study/ray_data_aux_index_other_exists=true'
  expect_err = 'However, Ray non-aux data with said name was found.'
  allow_test_objects = true
  detail = 'getting the index of ray aux data when data with the same name exists'
[]

[test]
  type = 'CSVDiff'
  input = 'kernel_change_ray.i'
  csvdiff = 'kernel_change_ray_out.csv'
  allow_test_objects = true
  issues = '#16028'
  design = 'RayKernels/index.md RayTracingStudy.md'
  requirement = 'The system shall be able to change the trajectory of rays during tracing within RayKernels.'
[]

[2d]
  type = 'CSVDiff'
  input = 'kernel_create_ray.i'
  csvdiff = 'kernel_create_ray_2d_out.csv'
  allow_test_objects = true
  cli_args = 'Outputs/file_base=kernel_create_ray_2d_out'
  detail = 'two-dimensional meshes '
[]

[3d]
  type = 'CSVDiff'
  input = 'kernel_create_ray.i'
  csvdiff = 'kernel_create_ray_3d_out.csv'
  allow_test_objects = true
  cli_args = 'Outputs/file_base=kernel_create_ray_3d_out
                  Mesh/active=gmg_3d
                  RayBCs/active=kill_3d'
  detail = 'and in three-dimensional meshes.'
[]

[smart]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'lots_out.csv'
  cli_args = 'UserObjects/lots/method=smart
                  Outputs/file_base=lots_out'
  allow_test_objects = true
  detail = 'an asynchronous method, '
[]

[bs]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'lots_out.csv'
  cli_args = 'UserObjects/lots/method=bs
                  UserObjects/lots/allow_new_work_during_execution=false
                  Outputs/file_base=lots_out'
  allow_test_objects = true
  detail = 'a bulk-synchronous method, '
[]

[harm]
  type = 'CSVDiff'
  input = '../lots.i'
  csvdiff = 'lots_out.csv'
  cli_args = 'UserObjects/lots/method=harm
                  UserObjects/lots/allow_new_work_during_execution=false
                  Outputs/file_base=lots_out'
  allow_test_objects = true
  detail = 'and the harm method.'
[]

[test]
  type = 'Exodiff'
  input = 'multiple_subdomains.i'
  exodiff = 'multiple_subdomains_out.e'
  allow_test_objects = true
  issues = '#16028'
  design = 'RayTracingStudy.md'
  requirement = 'The system shall support the tracing of rays across multiple subdomains.'
[]

[sized]
  type = 'CSVDiff'
  input = 'ray_data.i'
  csvdiff = 'ray_data_out.csv'
  allow_test_objects = true
  detail = 'sizing the Ray data on acquire '
[]

[unsized]
  type = 'CSVDiff'
  input = 'ray_data.i'
  # needs to be a separate file because recover tests will fail
  # when appending to an output from another test in the same directory
  csvdiff = 'ray_data_unsized_out.csv'
  allow_test_objects = true
  detail = 'and when not sizing the Ray data on acquire.'
  cli_args = 'UserObjects/test1/use_unsized_rays=true
                  UserObjects/test2/use_unsized_rays=true
                  UserObjects/test3/use_unsized_rays=true
                  UserObjects/test4/use_unsized_rays=true
                  Outputs/file_base=ray_data_unsized_out'
[]

[initial]
  type = CSVDiff
  input = registration_restart.i
  csvdiff = initial.csv
  cli_args = 'Outputs/file_base=initial'
  recover = false
  requirement = 'The system shall support the storage of registered rays via the restart and recover system'
[]

[restart]
  type = CSVDiff
  input = registration_restart.i
  csvdiff = restart.csv
  prereq = initial
  cli_args = 'Outputs/file_base=restart Executioner/start_time=1 Problem/force_restart=true Problem/restart_file_base=initial_cp/LATEST'
  recover = false
  requirement = 'The system shall support the loading of registered rays via the restart and recover system'
[]

[test]
  type = 'CSVDiff'
  input = 'reuse_rays.i'
  csvdiff = 'reuse_rays_out.csv'
  allow_test_objects = true
  issues = '#16028'
  design = 'RayTracingStudy.md'
  requirement = 'The system shall support the re-use of traced rays.'
[]

[no_bc]
  type = 'RunException'
  input = 'tolerate_failure.i'
  expect_err = 'Ray on processor .* and thread .* failed to trace'
  requirement = 'The system shall support the failure of a trace gracefully via a warning.'
  allow_warnings = true
  should_crash = false
[]

[max_distance]
  type = 'CSVDiff'
  input = 'max_distance.i'
  csvdiff = 'max_distance_out.csv'
  issues = '#16028'
  design = 'RayTracingStudy.md'
  requirement = 'The system shall support ending rays by a user-set maximum distance.'
[]

[directions_or_end_points]
  type = RunException
  input = 'errors.i'
  cli_args = "UserObjects/study/end_points='1 0 0'
                  UserObjects/study/directions='1 0 0'"
  expect_err = "Can only use 'directions' or 'end_points', but not both"

  detail = 'both directions and points are provided,'
[]

[neither_directions_or_end_points]
  type = RunException
  input = 'errors.i'
  expect_err = "Must set 'end_points' or 'directions'"

  detail = 'neither directions and points are provided,'
[]

[start_size]
  type = RunException
  input = 'errors.i'
  cli_args = "UserObjects/study/end_points='1 0 0'
                  UserObjects/study/start_points='1 0 0 1 0 0'"
  expect_err = "Not the same size as names"

  detail = 'the number of start points provided is inconsistent,'
[]

[directions_size]
  type = RunException
  input = 'errors.i'
  cli_args = "UserObjects/study/directions='1 0 0 1 0 0'"
  expect_err = "Not the same size as names"

  detail = 'the number of directions provided is inconsistent,'
[]

[max_distance_with_end]
  type = RunException
  input = 'errors.i'
  cli_args = "UserObjects/study/end_points='1 0 0'
                  UserObjects/study/max_distances=1"
  expect_err = "Can only be used when trajectories are set with the 'directions' parameter"

  detail = 'both end points and maximum distances are provided,'
[]

[max_distance_size]
  type = RunException
  input = 'errors.i'
  cli_args = "UserObjects/study/directions='1 0 0'
                  UserObjects/study/max_distances='1 2'"
  expect_err = "Must be the same size as 'start_points'"

  detail = 'the number of maximum distances provided is inconsistent,'
[]

[max_distance_positive]
  type = RunException
  input = 'errors.i'
  cli_args = "UserObjects/study/directions='1 0 0'
                  UserObjects/study/max_distances='-1'"
  expect_err = "Values must be positive"

  detail = 'the provided maximum distance is non-positive,'
[]

[end_points_size]
  type = RunException
  input = 'errors.i'
  cli_args = "UserObjects/study/end_points='1 0 0 1 0 0'"
  expect_err = "Not the same size as names"

  detail = 'the number of end points provided is inconsistent,'
[]

[initial_data_size]
  type = RunException
  input = 'errors.i'
  cli_args = "UserObjects/study/end_points='1 0 0'
                  UserObjects/study/ray_data_names=dummy
                  UserObjects/study/initial_ray_data='0; 0'"
  expect_err = "Data for 2 ray\(s\) was provided, but 1 ray\(s\) were defined"

  detail = 'the initial data size is inconsistent with the number of rays,'
[]

[initial_data_individual_size]
  type = RunException
  input = 'errors.i'
  cli_args = "UserObjects/study/end_points='1 0 0'
                  UserObjects/study/ray_data_names=dummy
                  UserObjects/study/initial_ray_data='0 0'"
  expect_err = "Data for index 0 \(ray 'ray'\) is not the size of 'ray_data_names'"

  detail = 'the initial data size is inconsistent with the provided data names,'
[]

[data_names_not_set]
  type = RunException
  input = 'errors.i'
  cli_args = "UserObjects/study/end_points='1 0 0'
                  UserObjects/study/initial_ray_data='1'"
  expect_err = "Can only be used if 'ray_data_names' is set"

  detail = 'initial data is provided but data names are not provided,'
[]

[initial_aux_data_size]
  type = RunException
  input = 'errors.i'
  cli_args = "UserObjects/study/end_points='1 0 0'
                  UserObjects/study/ray_aux_data_names=dummy
                  UserObjects/study/initial_ray_aux_data='0; 0'"
  expect_err = "Aux data for 2 ray\(s\) was provided, but 1 ray\(s\) were defined"

  detail = 'the initial aux data size is inconsistent with the number of rays,'
[]

[initial_aux_data_individual_size]
  type = RunException
  input = 'errors.i'
  cli_args = "UserObjects/study/end_points='1 0 0'
                  UserObjects/study/ray_aux_data_names=dummy
                  UserObjects/study/initial_ray_aux_data='0 0'"
  expect_err = "Aux data for index 0 \(ray 'ray'\) is not the size of 'ray_aux_data_names'"

  detail = 'the initial aux data size is inconsistent with the provided aux data names,'
[]

[aux_data_names_not_set]
  type = RunException
  input = 'errors.i'
  cli_args = "UserObjects/study/end_points='1 0 0'
                  UserObjects/study/initial_ray_aux_data='1'"
  expect_err = "Can only be used if 'ray_aux_data_names' is set"

  detail = 'and aux initial data is provided but aux data names are not provided.'
[]

[checkpoint]
  type = 'RunApp'
  input = 'recover.i'
  allow_test_objects = true
  cli_args = 'Outputs/checkpoint=true Executioner/num_steps=1'
  recover = false
  detail = 'outputting of checkpoint files and'
[]

[recover]
  type = 'CSVDiff'
  input = 'recover.i'
  csvdiff = 'recover_out.csv'
  cli_args = '--recover'
  allow_test_objects = true
  prereq = 'test/checkpoint'
  delete_output_before_running = false
  recover = false
  detail = 'shall be capable of restarting from the output data.'
[]

[replicated_no_claim]
  requirement = 'The RepeatableRayStudyBase shall report a useful warning when rays are to be replicated but claiming is disabled.'
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/study/type=RepeatableRayStudyDefineNoClaimTest'
  allow_test_objects = True
  expect_err = '_define_rays_replicated is being set to false.'
[]

[non_unique_id_error]
  requirement = 'The RepeatableRayStudyBase shall return a useful error when a ray is added with a non-unique ID.'
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/study/create_non_unique_id_rays=true'
  allow_test_objects = true
  expect_err = 'Multiple Rays exist with ID 0 on processor'
[]

[additional]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/study/create_additional_ray_pid=0'
  min_parallel = 2
  max_parallel = 2
  allow_test_objects = true
  expect_err = 'The size of _rays on rank 1 does not match the size of rays on rank 0.'
  detail = 'additional rays are created in a non-replicated manner on a subset of processors,'
[]

[non_replicated_id]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/study/create_non_replicated_id_ray_pid=0'
  min_parallel = 2
  max_parallel = 2
  allow_test_objects = true
  expect_err = 'A Ray was found on rank 1 with an ID that does not exist on rank 0.'
  detail = 'a ray exists with a non-replicated ID on a processor,'
[]

[non_replicated]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/study/create_non_replicated_ray_pid=0'
  min_parallel = 2
  max_parallel = 2
  allow_test_objects = true
  expect_err = 'A Ray was found on rank 1 that does not exist on rank 0.'
  detail = 'a non-replicated ray exists on a processor,'
[]

[starting_elem_set]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/study/define_ray_with_starting_elem=true'
  allow_test_objects = true
  expect_err = 'A Ray was found in _rays after defineRays\(\) that has a starting element or incoming side set.'
  detail = 'and a ray is created with starting element information when it should not.'
[]

[no_rays]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/study/define_no_rays=true'
  allow_test_objects = true
  expect_err = 'No Rays were moved to _rays in defineRays\(\)'

  detail = 'no rays were defined or'
[]

[null_ray]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/study/define_nullptr_ray=true'
  allow_test_objects = true
  expect_err = 'A nullptr Ray was found in _rays after defineRays\(\).'

  detail = 'a null ray was defined.'
[]

[pic]
  type = 'Exodiff'
  input = 'pic_ray_study.i'
  exodiff = 'pic_ray_study_rays.e
                 pic_ray_study_rays.e-s0001
                 pic_ray_study_rays.e-s0002
                 pic_ray_study_rays.e-s0003
                 pic_ray_study_rays.e-s0004
                 pic_ray_study_rays.e-s0005
                 pic_ray_study_rays.e-s0006
                 pic_ray_study_rays.e-s0007
                 pic_ray_study_rays.e-s0008
                 pic_ray_study_rays.e-s0009'
  allow_test_objects = tru
  requirement = 'The system shall support the propagation of rays for particle-in-cell capability by tracing rays via a temporally and spatially varying velocity in time'
[]

[non_positive_polar_order]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/test/non_positive_polar_order=true'
  allow_test_objects = true
  expect_err = 'polar_order must be positive in RayTracingAngularQuadrature'
  detail = 'the provided polar order is non-positive'
[]

[non_positive_azimuthal_order]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/test/non_positive_azimuthal_order=true'
  allow_test_objects = true
  expect_err = 'azimuthal_order must be positive in RayTracingAngularQuadrature'
  detail = 'the provided azimuthal order is non-positive'
[]

[mu_min_larger]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/test/mu_min_larger=true'
  allow_test_objects = true
  expect_err = 'mu_min must be < mu_max in RayTracingAngularQuadrature'
  detail = 'the minimum polar cosine is greater than the maximum polar cosine'
[]

[mu_min_too_small]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/test/mu_min_too_small=true'
  allow_test_objects = true
  expect_err = 'mu_min must be >= -1 in RayTracingAngularQuadrature'
  detail = 'the minimum polar cosine is less than negative one'
[]

[mu_max_too_big]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/test/mu_max_too_big=true'
  allow_test_objects = true
  expect_err = 'mu_max must be <= 1 in RayTracingAngularQuadrature'
  detail = 'the maximum polar cosine is greater than one'
[]

[dim1]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/test/dim1=true'
  allow_test_objects = true
  expect_err = 'RayTracingAngularQuadrature only supports dimensions 2 and 3'
  detail = 'the dimension is less than two'
[]

[non_positive_gauss_legendre_order]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/test/non_positive_gauss_legendre_order=true'
  allow_test_objects = true
  expect_err = 'Order must be positive in gaussLegendre\(\)'
  detail = 'the Gauss-Legendre order is non-positive'
[]

[check_direction]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/test/check_direction=true'
  allow_test_objects = true
  expect_err = 'RayTracingAngularQuadrature does not have direction 1337'
  detail = 'an invalid direction index is requested'
[]

[orthonormal_vector_zero]
  type = 'RunException'
  input = 'errors.i'
  cli_args = 'UserObjects/test/orthonormal_vector_zero=true'
  allow_test_objects = true
  expect_err = 'Vector v has norm close to 0 in orthonormalVector\(\)'
  detail = 'an orthonormal vector is requested from a vector that has a zero norm'
[]

[test]
  type = 'RunApp' # exodiff
  input = 'per_processor_ray_tracing_results_vector_postprocessor.i'
  # TODO: find out a test for this that is not platform dependent so we can do exodiff
  # exodiff = 'per_processor_ray_tracing_results_vector_postprocessor_out.e'
  # min_parallel = 4
  # max_parallel = 4
  # min_threads = 1
  # max_threads = 1
  allow_test_objects = true
  issues = '#16028'
  design = 'PerProcessorRayTracingResultsVectorPostprocessor.md'
  requirement = 'The system shall support the output of information pretaining to the tracing of rays on a per-processor basis.'
[]