From ee7e7e9815350d5537acd1a52c9e1b8b5bd7f9ce Mon Sep 17 00:00:00 2001 From: Tyagi Date: Mon, 5 Oct 2026 14:41:55 +1100 Subject: [PATCH] test: provide a hierarchy for the constrained traction benchmark Keep the SolCx fine grid at 32x32 elements by refining a 16x16 base once. This supplies geometric multigrid to the high-contrast augmented saddle-point solve instead of falling back to single-level GAMG. The unchanged upstream test times out after 180 seconds on Gadi. A guarded diagnostic spends 21,666 velocity iterations in 90 seconds; the hierarchy completes and satisfies all existing traction, negative-control and topography bounds in 35-40 seconds with PETSc 3.25.4 int32 and int64. No tolerances, physical parameters, accuracy bounds or production defaults are changed. Underworld development team with AI support from Claude Code --- tests/test_1063_constrained_traction.py | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/tests/test_1063_constrained_traction.py b/tests/test_1063_constrained_traction.py index 003e747bb..c7a641b00 100644 --- a/tests/test_1063_constrained_traction.py +++ b/tests/test_1063_constrained_traction.py @@ -30,8 +30,14 @@ def _solve(): + # Keep the 32x32 fine mesh while providing geometric multigrid for the + # high-contrast augmented system; a single-level GAMG solve stalls on Gadi. mesh = uw.meshing.StructuredQuadBox( - elementRes=(RES, RES), minCoords=(0.0, 0.0), maxCoords=(1.0, 1.0), qdegree=3 + elementRes=(RES // 2, RES // 2), + refinement=1, + minCoords=(0.0, 0.0), + maxCoords=(1.0, 1.0), + qdegree=3, ) sol = A.SolCx(mesh, eta_A=1.0, eta_B=ETA_B, x_c=0.5, n=1)