Restore vertex default and pin legacy mortar baselines

Author: Andreas Henrik Frederiksen

framework/src/constraints/MortarInterfaceWarehouse.C

@@ -19,8 +19,8 @@ namespace
 {
 struct ParsedMortarTriangulationOptions
 {
-  MortarSegmentTriangulationMode mode = MortarSegmentTriangulationMode::centroid;
-  std::string canonical_name = "centroid";
+  MortarSegmentTriangulationMode mode = MortarSegmentTriangulationMode::vertex;
+  std::string canonical_name = "vertex";
 };
 
 ParsedMortarTriangulationOptions

framework/src/constraints/MortarSegmentHelper.C

@@ -10,9 +10,11 @@
 #include "MooseError.h"
 
 #include "libmesh/int_range.h"
+#if defined(LIBMESH_HAVE_TRIANGLE) || defined(LIBMESH_HAVE_POLY2TRI)
 #include "libmesh/replicated_mesh.h"
 #include "libmesh/mesh_triangle_interface.h"
 #include "libmesh/poly2tri_triangulator.h"
+#endif
 
 #include <algorithm>
 #include <array>

framework/src/interfaces/MortarConsumerInterface.C

@@ -97,7 +97,7 @@ MortarConsumerInterface::validParams()
 #else
       "vertex centroid ear_clipping",
 #endif
-      "centroid");
+      "vertex");
   triangulation.addDocumentation(
       "vertex",
       "Triangulate clipped 3D mortar polygons by forming a fan from an existing polygon vertex.");
@@ -115,8 +115,7 @@ MortarConsumerInterface::validParams()
   params.addParam<MooseEnum>(
       "triangulation",
       triangulation,
-      "Strategy used to triangulate clipped 3D mortar polygons into mortar segments. The default "
-      "is 'centroid' to preserve the legacy 3D mortar segmentation behavior.");
+      "Strategy used to triangulate clipped 3D mortar polygons into mortar segments.");
   params.addParam<bool>(
       "triangulate_triangles",
       false,

modules/contact/doc/content/syntax/Contact/index.md

@@ -28,8 +28,7 @@ In addition, face/face contact using a mortar method can also be specified using
 
 For 3D mortar contact, the triangulation of clipped overlap polygons can be selected with the
 `triangulation` parameter. The `triangulate_triangles` parameter controls whether polygons that
-are already triangles are left intact or centroid-subdivided again. By default, clipped polygons
-use centroid triangulation and already-triangular polygons are left intact.
+are already triangles are left intact or centroid-subdivided again.
 
 ## Constructed Objects
 

modules/contact/src/actions/ContactAction.C

@@ -241,7 +241,7 @@ ContactAction::validParams()
 #else
       "vertex centroid ear_clipping",
 #endif
-      "centroid");
+      "vertex");
   triangulation.addDocumentation(
       "vertex",
       "Triangulate clipped 3D mortar polygons by forming a fan from an existing polygon vertex.");
@@ -259,8 +259,7 @@ ContactAction::validParams()
   params.addParam<MooseEnum>(
       "triangulation",
       triangulation,
-      "Strategy used to triangulate clipped 3D mortar polygons into mortar segments. The default "
-      "is 'centroid' to preserve the legacy 3D mortar segmentation behavior.");
+      "Strategy used to triangulate clipped 3D mortar polygons into mortar segments.");
   params.addParam<bool>(
       "triangulate_triangles",
       false,

test/tests/mortar/3d-periodic/tests

@@ -5,6 +5,8 @@
     type = 'Exodiff'
     input = 'periodic.i'
     exodiff = 'periodic_out.e'
+    # Preserve the historical solution baseline while the default triangulation is vertex.
+    cli_args = 'Constraints/ev/triangulation=centroid'
     requirement = 'We shall be able to enforce periodic boundary conditions on 3d domains using mortar constraints.'
     abs_zero = 1e-08
     map = False

test/tests/mortar/continuity-3d-non-conforming/tests

@@ -5,7 +5,9 @@
     type = 'Exodiff'
     input = 'continuity_sphere_hex.i'
     exodiff = 'continuity_sphere_hex8_out.e'
-    cli_args = 'Outputs/file_base=continuity_sphere_hex8_out'
+    # Preserve the historical solution baseline while the default triangulation is vertex.
+    cli_args = 'Outputs/file_base=continuity_sphere_hex8_out '
+               'Constraints/mortar/triangulation=centroid'
     requirement = 'The system shall be able to produce the expected result for a solution continuity '
                   'test case using the mortar method on a HEX8 mesh with curved geometry, e.g. the '
                   'primal variable values across the mortar interface shall be the same.'
@@ -152,7 +154,9 @@
     type = 'Exodiff'
     input = 'continuity_penalty_sphere_hex.i'
     exodiff = 'continuity_penalty_sphere_hex8_out.e'
-    cli_args = 'Outputs/file_base=continuity_penalty_sphere_hex8_out'
+    # Preserve the historical solution baseline while the default triangulation is vertex.
+    cli_args = 'Outputs/file_base=continuity_penalty_sphere_hex8_out '
+               'Constraints/mortar/triangulation=centroid'
     requirement = 'The system shall be able to produce the expected result for a solution continuity '
                   'test case using the penalty-based mortar method on a HEX8 mesh with curved '
                   'geometry, e.g. the primal variable values across the mortar interface shall be '