test: add scaled N256 Taylor-Green diagnostic

tests/validation/taylor_green/run_case.py

@@ -18,6 +18,7 @@ SOURCE_SHA256 = "a5dbfd804befce868067d79a23e5fd47dce7f505b3f9b1bb2145fb57cb5f217
 BUILD_CAPTURE_SHA256 = "3dbaa6547cb0a4cb1541054f3b64ad8b6e18cc70527ac6166e56b944e9bd464a"
 REFERENCE_HARNESS = "11901a7f77557a901e017a171eeae1758d584fa5"
 INITIALIZATION = "physical FP64 equilibrium G at actual cell centers; rho=1+3p"
+SCALED_DIAGNOSTIC = "scaled_n256"
 
 
 def sample_steps(n, early_pressure=False):
@@ -145,10 +146,12 @@ def main():
     p.add_argument("--run-gpu", action="store_true")
     p.add_argument("--dry-run", action="store_true")
     p.add_argument("--re", type=int, choices=[100, 1000], required=True)
-    p.add_argument("--n", type=int, choices=[32, 64, 128, 256], required=True)
+    p.add_argument("--n", type=int, choices=[32, 64, 128, 256, 512], required=True)
     p.add_argument("--depth", type=int, choices=[1, 2, 3], required=True)
     p.add_argument("--output", type=Path, required=True)
     p.add_argument("--explosion", choices=["constant", "projected_linear"], default="constant")
+    p.add_argument("--diagnostic", choices=[SCALED_DIAGNOSTIC],
+                   help="explicit four-case N256/fine512 same-placement diagnostic")
     p.add_argument("--early-pressure", action="store_true",
                    help="also save 65 scalar pressure samples at t*=k*.001875, k=0..64")
     args = p.parse_args()
@@ -156,6 +159,13 @@ def main():
         p.error("Output already exists; choose a new path to retain prior evidence")
     if args.explosion != "constant" and args.depth != 2:
         p.error("projected_linear requires depth 2")
+    if args.diagnostic:
+        cases = {(256, 1, "constant"), (512, 1, "constant"),
+                 (256, 2, "constant"), (256, 2, "projected_linear")}
+        if args.re != 100 or (args.n, args.depth, args.explosion) not in cases or not args.early_pressure:
+            p.error("scaled_n256 requires one of its four Re100 cases and --early-pressure")
+    elif args.n == 512:
+        p.error("N512 requires --diagnostic scaled_n256")
     selected = os.environ.get("CUDA_VISIBLE_DEVICES", "")
     if not args.dry_run and (not args.run_gpu or not selected or "," in selected or selected == "-1"):
         p.error("--run-gpu and exactly one explicit CUDA_VISIBLE_DEVICES are required")
@@ -164,7 +174,8 @@ def main():
     source = verify_source(repo, args.expected_revision)
     runtime = runtime_source(repo)
     if args.dry_run:
-        print(json.dumps(dict(source=source, runtime_source=runtime, dry_run=True)))
+        print(json.dumps(dict(source=source, runtime_source=runtime, dry_run=True,
+                              **({"diagnostic": args.diagnostic} if args.diagnostic else {}))))
         return 0
     from CelerisLab.cuda.context import CudaContext
     from CelerisLab.lbm.blocks import BlockLayout2D
@@ -212,6 +223,8 @@ def main():
                             sampled_max_pressure_l2=None),
         completed=False,
     )
+    if args.diagnostic:
+        report["diagnostic"] = args.diagnostic
     args.output.parent.mkdir(parents=True, exist_ok=True)
     def save():
         temporary = args.output.with_suffix(args.output.suffix + ".tmp")