#!/usr/bin/env python3

import contextlib
import base64
from collections import Counter
import hashlib
import io
import json
import math
import os
from pathlib import Path
import re
import sys
import tempfile
import unittest
from unittest import mock


sys.path.insert(0, str(Path(__file__).resolve().parents[1]))

from postgis_exampletest import ExampleTester, QueryRows, parse_source_example


def xyz_vertices(points, ring=False):
    return [
        {"path": [1, index] if ring else [index], "point": point}
        for index, point in enumerate(points, 1)
    ]


class ExampleTestParserTest(unittest.TestCase):
    def setUp(self):
        self.tester = ExampleTester.__new__(ExampleTester)

    def test_empty_psql_table_is_empty_rows(self):
        rows = self.tester.expected_rows_from_psql_lines(
            [
                " name | default_version | installed_version",
                "------+-----------------+-------------------",
                "(0 rows)",
            ]
        )
        self.assertEqual([], rows)

    def test_non_empty_psql_table_rows_are_parsed(self):
        rows = self.tester.expected_rows_from_psql_lines(
            [
                " name    | default_version | installed_version",
                "---------+-----------------+-------------------",
                " postgis | 3.7.0dev        | null",
                "(1 row)",
            ]
        )
        self.assertEqual([["postgis", "3.7.0dev", "null"]], rows)

    def test_unicode_psql_table_rows_are_parsed(self):
        rows = self.tester.expected_rows_from_psql_lines(
            [
                "cube_surface_area │ solid_surface_area",
                "──────────────────┼────────────────────",
                "6                 │ 0",
            ]
        )
        self.assertEqual([["6", "0"]], rows)

    def test_psql_headers_are_preserved_for_visual_labels(self):
        lines = [
            " center | nearest ",
            "--------+---------",
            " POINT(0 0) | POINT(1 1)",
            "(1 row)",
        ]
        self.assertEqual(["center", "nearest"], self.tester.psql_headers_from_lines(lines))

    def test_psql_expanded_record_is_transposed_into_one_row(self):
        lines = [
            "-[ RECORD 1 ]-------------------",
            "auth_name | EPSG",
            "auth_srid | 3005",
            "srname    | NAD83 / BC Albers",
        ]
        self.assertEqual(
            [["EPSG", "3005", "NAD83 / BC Albers"]],
            self.tester.expected_rows_from_psql_lines(lines),
        )
        self.assertEqual(
            ["auth_name", "auth_srid", "srname"],
            self.tester.psql_headers_from_lines(lines),
        )

    def test_multiple_psql_expanded_records_become_multiple_rows(self):
        lines = [
            "-[ RECORD 1 ]---",
            "name | first",
            "srid | 1",
            "-[ RECORD 2 ]---",
            "name | second",
            "srid | 2",
        ]
        self.assertEqual(
            [["first", "1"], ["second", "2"]],
            self.tester.expected_rows_from_psql_lines(lines),
        )
        self.assertEqual(["name", "srid"], self.tester.psql_headers_from_lines(lines))

    def test_multiline_programlisting_is_trimmed(self):
        query, expected = parse_source_example(
            """
SELECT ST_AsText(
    ST_Point(1, 2));
""",
            """
POINT(1 2)
""",
        )
        self.assertEqual("SELECT ST_AsText(\n    ST_Point(1, 2))", query)
        self.assertEqual([["POINT(1 2)"]], expected)

    def test_function_whitespace_is_allowed(self):
        self.assertTrue(self.tester.query_is_auto_safe("SELECT ST_AsText ('POINT(1 2)'::geometry);"))

    def test_geometry_candidates_are_lexical_and_comment_aware(self):
        sql = """-- 'POINT(9 9)'
SELECT 'POINT(1 2)', $$LINESTRING(0 0,1 1)$$,
       'not geometry', /* 'POINT(8 8)' */ 'CURVEPOLYGON(CIRCULARSTRING(0 0,1 1,2 0,1 -1,0 0))';
"""
        self.assertEqual(
            [
                "POINT(1 2)",
                "LINESTRING(0 0,1 1)",
                "CURVEPOLYGON(CIRCULARSTRING(0 0,1 1,2 0,1 -1,0 0))",
            ],
            self.tester.geometry_candidates(sql, quoted_only=True),
        )

    def test_geometry_candidates_reject_unclosed_wkt(self):
        with self.assertRaisesRegex(RuntimeError, "Unclosed geometry candidate"):
            self.tester.geometry_candidates("POLYGON((0 0,1 0,0 0)")

    def test_hexwkb_geometry_recognizes_endian_flags_offsets_and_extended_types(self):
        def payload(type_word, byte_order=1):
            endian = "little" if byte_order == 1 else "big"
            encoded_type = type_word & 0x0FFFFFFF
            dimension_offset = encoded_type // 1000
            base_type = encoded_type % 1000 if dimension_offset in (1, 2, 3) else encoded_type
            dimensions = 2
            dimensions += int(bool(type_word & 0x80000000) or dimension_offset in (1, 3))
            dimensions += int(bool(type_word & 0x40000000) or dimension_offset in (2, 3))
            body = b""
            if type_word & 0x20000000:
                body += (4326).to_bytes(4, endian)
            if base_type == 1:
                body += b"\x00" * (dimensions * 8)
            elif base_type == 21:
                body += b"\x00" * 12
            else:
                body += b"\x00" * 4
            return (bytes([byte_order]) + type_word.to_bytes(4, endian) + body).hex()

        type_names = {
            1: "POINT", 2: "LINESTRING", 3: "POLYGON", 4: "MULTIPOINT",
            5: "MULTILINESTRING", 6: "MULTIPOLYGON", 7: "GEOMETRYCOLLECTION",
            8: "CIRCULARSTRING", 9: "COMPOUNDCURVE", 10: "CURVEPOLYGON",
            11: "MULTICURVE", 12: "MULTISURFACE", 15: "POLYHEDRALSURFACE",
            16: "TIN", 17: "TRIANGLE", 21: "NURBSCURVE",
        }
        cases = [
            (payload(type_code, byte_order), expected_type)
            for type_code, expected_type in type_names.items()
            for byte_order in (0, 1)
        ]
        cases.extend((
            (payload(0xE0000001), "POINT"),
            (payload(3003, byte_order=0), "POLYGON"),
            (payload(0x20000000 | 2021), "NURBSCURVE"),
            (payload(0x20000000 | 3021, byte_order=0), "NURBSCURVE"),
        ))
        for value, expected_type in cases:
            with self.subTest(value=value):
                self.assertEqual(
                    expected_type,
                    self.tester.hexwkb_geometry(value)["type"],
                )

    def test_hexwkb_geometry_rejects_false_positives_and_bbox_flag(self):
        self.assertIsNone(self.tester.hexwkb_geometry("0103"))
        self.assertIsNone(self.tester.hexwkb_geometry("0103000000"))
        self.assertIsNone(self.tester.hexwkb_geometry("010300000"))
        self.assertIsNone(self.tester.hexwkb_geometry("0203000000"))
        self.assertIsNone(self.tester.hexwkb_geometry("0100000000"))
        self.assertIsNone(self.tester.hexwkb_geometry("0112000000"))
        self.assertIsNone(self.tester.hexwkb_geometry("0103000010"))
        self.assertIsNone(self.tester.hexwkb_geometry("01E9030080" + "00" * 32))
        self.assertIsNone(self.tester.hexwkb_geometry("0115000080" + "00" * 24))
        self.assertIsNone(self.tester.hexwkb_geometry("01E9030020" + "00" * 32))
        self.assertIsNone(self.tester.hexwkb_geometry("0115000000" + "00" * 4))
        self.assertIsNone(self.tester.hexwkb_geometry("89504E470D0A1A0A"))
        self.assertIsNone(self.tester.hexwkb_geometry("not hex"))

    def test_hexwkb_geometry_accepts_bytea_prefix_only_when_requested(self):
        self.assertEqual(
            {"hexwkb": "010300000000000000", "type": "POLYGON"},
            self.tester.hexwkb_geometry("  \\x010300000000000000  "),
        )
        self.assertIsNone(
            self.tester.hexwkb_geometry("\\x010300000000000000", allow_bytea_prefix=False)
        )

    def test_geometry_output_hex_requires_geometry_type_metadata(self):
        value = "010300000000000000"
        self.assertEqual("POLYGON", self.tester.geometry_output_hex(value, "geometry")["type"])
        self.assertEqual("POLYGON", self.tester.geometry_output_hex(value, "geography")["type"])
        self.assertIsNone(self.tester.geometry_output_hex(value, "text"))
        self.assertIsNone(self.tester.geometry_output_hex("\\x" + value, "bytea"))

    def test_textual_surface_wkt_compares_canonically(self):
        actual = QueryRows(
            [[
                "POLYHEDRALSURFACE Z ("
                "((0 0 0,0 0 2,0 2 0,0 0 0)),"
                "((0 0 0,0 2 0,2 0 0,0 0 0))"
                ")"
            ]],
            types=["text"],
        )
        expected = [[
            "POLYHEDRALSURFACE Z ("
            "((0 0 0,0 2 0,2 0 0,0 0 0)),"
            "((0 0 0,0 0 2,0 2 0,0 0 0))"
            ")"
        ]]
        self.assertTrue(self.tester.typed_rows_equal("unused", actual, expected))

    def test_multiline_single_column_wkt_rows_are_coalesced(self):
        rows = QueryRows([
            ["POLYGON((234776.9 899563.7,234896.5 899456.7,"],
            ["234914 899436.4,234946.6 899356.9,"],
            ["234776.9 899563.7))"],
        ])
        self.assertEqual(
            [[
                "POLYGON((234776.9 899563.7,234896.5 899456.7,"
                "234914 899436.4,234946.6 899356.9,"
                "234776.9 899563.7))"
            ]],
            self.tester.coalesce_single_column_wkt_rows(rows),
        )

    def test_function_type_substring_does_not_swallow_argument_geometry(self):
        self.assertEqual(
            ["POLYGON((0 0,2 0,0 2,0 0))", "POINT(1 1)"],
            self.tester.geometry_candidates(
                "SELECT ST_ClosestPoint('POLYGON((0 0,2 0,0 2,0 0))', 'POINT(1 1)')",
                quoted_only=True,
            ),
        )

    def test_visual_skip_keeps_example_test_without_rendering(self):
        xml = """<book xmlns="http://docbook.org/ns/docbook">
  <refentry xml:id="example">
    <refsection>
      <programlisting role="visual-skip">SELECT ST_AsText('LINESTRING(0 0,1 1)'::geometry);</programlisting>
      <screen>LINESTRING(0 0,1 1)</screen>
    </refsection>
  </refentry>
</book>"""
        with tempfile.NamedTemporaryFile("w", suffix=".xml", encoding="utf-8") as source:
            source.write(xml)
            source.flush()
            example = ExampleTester(source.name).examples()[0]
        self.assertTrue(example["valid"])
        self.assertIsNone(example["visual_id"])

    def test_explicit_visual_can_be_selected_from_binary_image_output(self):
        xml = """<book xmlns="http://docbook.org/ns/docbook">
  <refentry xml:id="image-example">
    <refsection>
      <programlisting>SELECT ST_AsPNG(ST_AsRaster(ST_Buffer(ST_Point(1, 5), 10), 150, 150)) AS image;</programlisting>
      <screen role="visual-primary"> image
-------
 PNG image, 150 x 150 pixels
      </screen>
    </refsection>
  </refentry>
</book>"""
        with tempfile.NamedTemporaryFile("w", suffix=".xml", encoding="utf-8") as source:
            source.write(xml)
            source.flush()
            example = ExampleTester(source.name).examples()[0]
        self.assertEqual("explicit", example["visual_kind"])
        self.assertEqual("visual-image-example-01", example["visual_id"])

    def test_visual_separate_output_role_is_preserved(self):
        xml = """<book xmlns="http://docbook.org/ns/docbook">
  <refentry xml:id="comparison">
    <refsection>
      <programlisting>SELECT ST_AsText('LINESTRING(0 0,1 1)'), ST_AsText('POINT(0 0)');</programlisting>
      <screen role="visual-separate-output">LINESTRING(0 0,1 1) | POINT(0 0)</screen>
    </refsection>
  </refentry>
</book>"""
        with tempfile.NamedTemporaryFile("w", suffix=".xml", encoding="utf-8") as source:
            source.write(xml)
            source.flush()
            example = ExampleTester(source.name).examples()[0]
        self.assertTrue(example["visual_separate_output"])

    def test_visual_overlay_role_is_preserved(self):
        xml = """<book xmlns="http://docbook.org/ns/docbook">
  <refentry xml:id="overlay">
    <refsection>
      <programlisting>SELECT ST_AsText(ST_Buffer('POINT(0 0)'::geometry, 1));</programlisting>
      <screen role="visual-overlay">POLYGON((0 0,1 0,0 1,0 0))</screen>
    </refsection>
  </refentry>
</book>"""
        with tempfile.NamedTemporaryFile("w", suffix=".xml", encoding="utf-8") as source:
            source.write(xml)
            source.flush()
            example = ExampleTester(source.name).examples()[0]
        self.assertTrue(example["visual_overlay"])

    def test_visual_layout_roles_do_not_disable_self_contained_queries(self):
        for role, flag in (
            ("visual-overlay", "visual_overlay"),
            ("visual-direction", "visual_direction"),
            ("visual-separate-output", "visual_separate_output"),
            ("visual-tessellation", "visual_tessellation"),
        ):
            with self.subTest(role=role):
                xml = f"""<book xmlns="http://docbook.org/ns/docbook">
  <refentry xml:id="alpha-shape">
    <refsection>
      <programlisting>SELECT ST_AsText(CG_AlphaShape('MULTIPOINT((0 0),(1 0),(0 1))'::geometry, 1));</programlisting>
      <screen role="{role}">POLYGON((0 0,1 0,0 1,0 0))</screen>
    </refsection>
  </refentry>
</book>"""
                with tempfile.NamedTemporaryFile("w", suffix=".xml", encoding="utf-8") as source:
                    source.write(xml)
                    source.flush()
                    example = ExampleTester(source.name).examples()[0]
                self.assertFalse(example["documented_only"])
                self.assertEqual("explicit", example["visual_kind"])
                self.assertTrue(example[flag])

    def test_documented_output_role_keeps_explicit_visual_without_execution(self):
        xml = """<book xmlns="http://docbook.org/ns/docbook">
  <refentry xml:id="preview">
    <refsection>
      <programlisting role="documented-output">SELECT ExpensiveVersionSensitiveGeometry() AS geom;</programlisting>
      <screen role="visual-primary">POLYGON((0 0,1 0,0 1,0 0))</screen>
    </refsection>
  </refentry>
</book>"""
        with tempfile.NamedTemporaryFile("w", suffix=".xml", encoding="utf-8") as source:
            source.write(xml)
            source.flush()
            example = ExampleTester(source.name).examples()[0]
        self.assertTrue(example["documented_only"])
        self.assertEqual("explicit", example["visual_kind"])
        self.assertEqual("visual-preview-01", example["visual_id"])

    def test_documented_output_role_keeps_auto_visual_without_execution(self):
        xml = """<book xmlns="http://docbook.org/ns/docbook">
  <refentry xml:id="preview">
    <refsection>
      <programlisting role="documented-output">SELECT ST_Buffer('POINT(0 0)'::geometry, 1);</programlisting>
      <screen>POLYGON((1 0,0 1,-1 0,0 -1,1 0))</screen>
    </refsection>
  </refentry>
</book>"""
        with tempfile.NamedTemporaryFile("w", suffix=".xml", encoding="utf-8") as source:
            source.write(xml)
            source.flush()
            example = ExampleTester(source.name).examples()[0]
        self.assertTrue(example["documented_only"])
        self.assertEqual("geometry-output", example["visual_kind"])
        self.assertEqual("visual-preview-01", example["visual_id"])

    def test_self_contained_functions_are_not_skipped_by_name(self):
        tester = ExampleTester.__new__(ExampleTester)
        for function in (
            "ST_MinimumSpanningTree",
            "ST_ClusterKMeans",
            "CG_3DConvexHull",
            "CG_AlphaShape",
            "CG_OptimalConvexPartition",
        ):
            with self.subTest(function=function):
                self.assertEqual("", tester.obvious_skip_reason(f"SELECT {function}(geom)"))

    def test_capability_roles_are_parsed_and_compared_generically(self):
        xml = """<book xmlns="http://docbook.org/ns/docbook">
  <refentry xml:id="capability"><refsection>
    <programlisting role="requires-geos-3.15 requires-proj-9.2 requires-cgal-6.1 requires-protobuf-1.1.0 requires-raster-1.0.0">SELECT 1;</programlisting>
    <screen>1</screen>
  </refsection></refentry>
</book>"""
        with tempfile.NamedTemporaryFile("w", suffix=".xml", encoding="utf-8") as source:
            source.write(xml)
            source.flush()
            tester = ExampleTester(source.name)
            example = tester.examples()[0]
        self.assertEqual(
            ["geos", "proj", "cgal", "protobuf", "raster"],
            [item["name"] for item in example["requirements"]],
        )
        self.assertTrue(tester.requirements_satisfied(
            example["requirements"], {
                "geos": (3, 15, 1),
                "proj": (9, 2),
                "cgal": (6, 1),
                "protobuf": (1, 5, 2),
                "raster": (1, 0, 0),
            }
        ))
        self.assertFalse(tester.requirements_satisfied(
            example["requirements"], {
                "geos": (3, 14, 9),
                "proj": (9, 5),
                "cgal": (6, 2),
                "protobuf": (1, 5, 2),
                "raster": (1, 0, 0),
            }
        ))
        self.assertEqual(
            (6, 1, 0),
            tester.capability_version(
                "cgal",
                'SFCGAL="2.3.0" CGAL="6.1.0" Boost="1.83.0"',
            ),
        )
        self.assertEqual(
            (),
            tester.capability_version("cgal", 'SFCGAL="2.3.0" Boost="1.83.0"'),
        )

    def test_capability_requirements_are_inferred_from_known_functions(self):
        xml = """<book xmlns="http://docbook.org/ns/docbook">
  <refentry xml:id="capability"><refsection>
    <programlisting>SELECT ST_MakeEmptyRaster(1, 1, 0, 0, 1, -1, 0, 0, 0),
                           ST_CoverageEdges('MULTIPOLYGON EMPTY'::geometry),
                           CG_PolygonRepair('POLYGON EMPTY'::geometry);</programlisting>
    <screen>dummy</screen>
  </refsection></refentry>
</book>"""
        with tempfile.NamedTemporaryFile("w", suffix=".xml", encoding="utf-8") as source:
            source.write(xml)
            source.flush()
            example = ExampleTester(source.name).examples()[0]

        requirements = {item["name"]: item["minimum"] for item in example["requirements"]}
        self.assertEqual((1, 0, 0), requirements["raster"])
        self.assertEqual((3, 15, 0), requirements["geos"])
        self.assertEqual((2, 3, 0), requirements["sfcgal"])
        self.assertEqual((6, 0, 0), requirements["cgal"])

    def test_explicit_capability_requirement_overrides_inferred_minimum(self):
        xml = """<book xmlns="http://docbook.org/ns/docbook">
  <refentry xml:id="capability"><refsection>
    <programlisting role="requires-geos-3.12">SELECT ST_LargestEmptyCircle(
      'MULTIPOINT((0 0),(1 1))'::geometry
    );</programlisting>
    <screen>dummy</screen>
  </refsection></refentry>
</book>"""
        with tempfile.NamedTemporaryFile("w", suffix=".xml", encoding="utf-8") as source:
            source.write(xml)
            source.flush()
            example = ExampleTester(source.name).examples()[0]

        self.assertEqual(
            {"geos": (3, 12)},
            {
                requirement["name"]: requirement["minimum"]
                for requirement in example["requirements"]
            },
        )

    def test_cg_3dbuffer_requires_sfcgal_2_0(self):
        tester = ExampleTester.__new__(ExampleTester)

        requirements = tester.inferred_capability_requirements(
            "SELECT CG_3DBuffer('POINT EMPTY'::geometry, 1, 8, 0)"
        )

        self.assertEqual([{
            "name": "sfcgal",
            "minimum": (2, 0, 0),
            "role": "auto-requires-sfcgal-2.0.0",
        }], requirements)

    def test_legacy_2d_wkt_output_compares_against_2d_projection(self):
        tester = ExampleTester.__new__(ExampleTester)
        query = tester.geometry_comparison_query(
            "0101000040000000000000F03F00000000000000400000000000000840",
            expected_wkt="POINT(1 2)",
            actual_type="POINT",
        )
        self.assertIn("ST_Force2D(ST_GeomFromEWKB", query)
        self.assertIn("ST_Force2D(ST_GeomFromEWKT", query)

    def test_implicit_3d_wkt_output_keeps_z_for_historical_postgis_spelling(self):
        tester = ExampleTester.__new__(ExampleTester)
        query = tester.geometry_comparison_query(
            "0101000080000000000000F03F00000000000000400000000000000840",
            expected_wkt="POINT(1 2 3)",
            actual_type="POINT",
        )
        self.assertNotIn("ST_Force2D", query)

    def test_explicit_geometry_output_precision_allows_hausdorff_tolerance(self):
        tester = ExampleTester.__new__(ExampleTester)
        query = tester.geometry_comparison_query(
            "0101000000000000000000F03F0000000000000040",
            expected_wkt="POINT(1.000000 2.000000)",
            actual_type="POINT",
            expected_wkt_digits=6,
        )
        self.assertIn("ST_HausdorffDistance(actual, expected) <= 9.9999999999999995e-07", query)

    def test_surface_output_precision_preserves_canonical_face_comparison(self):
        tester = ExampleTester.__new__(ExampleTester)
        query = tester.geometry_comparison_query(
            "00",
            expected_wkt="POLYHEDRALSURFACE Z EMPTY",
            actual_type="POLYHEDRALSURFACE",
            expected_wkt_digits=0,
        )
        self.assertEqual(2, query.count("jsonb_agg(face_key ORDER BY face_key)"))
        self.assertEqual(2, query.count("ST_AsEWKT(point.geom, 0)"))
        self.assertNotIn("ST_HausdorffDistance", query)

    def test_psql_environment_sets_default_query_timeout(self):
        tester = ExampleTester.__new__(ExampleTester)
        with mock.patch.dict(os.environ, {}, clear=True):
            env = tester.psql_environment()
        self.assertIn("postgis.gdal_enabled_drivers=PNG", env["PGOPTIONS"])
        self.assertIn("statement_timeout=60s", env["PGOPTIONS"])

    def test_psql_environment_allows_timeout_override(self):
        tester = ExampleTester.__new__(ExampleTester)
        with mock.patch.dict(
            os.environ,
            {"POSTGIS_EXAMPLETEST_STATEMENT_TIMEOUT": "5s"},
            clear=True,
        ):
            env = tester.psql_environment()
        self.assertIn("statement_timeout=5s", env["PGOPTIONS"])

    def test_psql_subprocess_timeout_follows_statement_timeout(self):
        tester = ExampleTester.__new__(ExampleTester)
        with mock.patch.dict(
            os.environ,
            {"POSTGIS_EXAMPLETEST_STATEMENT_TIMEOUT": "5s"},
            clear=True,
        ):
            self.assertEqual(20.0, tester.psql_subprocess_timeout())

    def test_psql_subprocess_timeout_can_be_disabled(self):
        tester = ExampleTester.__new__(ExampleTester)
        with mock.patch.dict(
            os.environ,
            {"POSTGIS_EXAMPLETEST_SUBPROCESS_TIMEOUT": "off"},
            clear=True,
        ):
            self.assertIsNone(tester.psql_subprocess_timeout())


class ExampleTestComparisonTest(unittest.TestCase):
    def setUp(self):
        self.tester = ExampleTester.__new__(ExampleTester)

    def test_polygon_ring_start_is_canonicalized(self):
        self.assertTrue(
            self.tester.values_equal(
                "POLYGON((0 0,1 0,1 1,0 0))",
                "POLYGON((1 0,1 1,0 0,1 0))",
            )
        )
        self.assertFalse(
            self.tester.values_equal(
                "POLYGON((0 0,1 0,1 1,0 0))",
                "POLYGON((0 0,2 0,1 1,0 0))",
            )
        )

    def test_multipoint_wrapped_point_spelling_is_canonicalized(self):
        self.assertTrue(
            self.tester.values_equal(
                "MULTIPOINT(0.5 0.5,2.5 2.5)",
                "MULTIPOINT((0.5 0.5),(2.5 2.5))",
            )
        )
        self.assertFalse(
            self.tester.values_equal(
                "MULTIPOINT(0.5 0.5,2.5 2.5)",
                "MULTIPOINT(2.5 2.5,0.5 0.5)",
            )
        )

    def test_png_bytea_is_summarized_without_losing_visual_payload(self):
        raw = b"\x89PNG\r\n\x1a\n" + b"\x00" * 8 + (150).to_bytes(4, "big") + (90).to_bytes(4, "big")
        actual = QueryRows([["\\x" + raw.hex()]], ["image"])
        prepared = self.tester.prepare_visual_rows(actual, ["bytea"])
        self.assertEqual([["PNG image, 150 x 90 pixels"]], prepared)
        self.assertEqual(1, len(prepared.visuals))
        self.assertEqual("Output", prepared.visuals[0]["label"])
        self.assertEqual("image/png", prepared.visuals[0]["media_type"])
        self.assertEqual(base64.b64encode(raw).decode("ascii"), prepared.visuals[0]["data"])

    def test_non_image_bytea_remains_text_output(self):
        actual = QueryRows([["\\x0102"]], ["payload"])
        prepared = self.tester.prepare_visual_rows(actual, ["bytea"])
        self.assertEqual([["\\x0102"]], prepared)
        self.assertEqual([], prepared.visuals)

    def test_scalar_row_value_labels_image_gallery_items(self):
        raw = b"\x89PNG\r\n\x1a\n" + b"\x00" * 8 + (4).to_bytes(4, "big") + (5).to_bytes(4, "big")
        actual = QueryRows([["weighted neighbors", "\\x" + raw.hex()]], ["title", "image"])
        prepared = self.tester.prepare_visual_rows(actual, ["text", "bytea"])
        self.assertEqual("weighted neighbors", prepared.visuals[0]["label"])

    def test_image_svg_embeds_buildtime_output_and_uses_alias_labels(self):
        svg = self.tester.visual_image_svg("raster-example", [{
            "source": "Output", "label": "grayscale", "media_type": "image/png",
            "width": 2, "height": 3, "data": "YWJj",
        }])
        self.assertIn('href="data:image/png;base64,YWJj"', svg)
        self.assertIn('preserveAspectRatio="none"', svg)
        self.assertIn('image-rendering:pixelated', svg)
        self.assertIn('data-postgis-geometry-id="raster-example-output-1"', svg)
        self.assertIn('>grayscale</text>', svg)

    def test_image_svg_uses_integer_upscale_for_pixel_art(self):
        svg = self.tester.visual_image_svg("pixels", [{
            "source": "Output", "label": "Output", "media_type": "image/png",
            "width": 150, "height": 90, "data": "YWJj",
        }])
        self.assertRegex(svg, r'<image[^>]+width="450"[^>]+height="270"')

    def test_image_gallery_uses_one_shared_pixel_scale(self):
        svg = self.tester.visual_image_svg("crop", [
            {"source": "Code", "label": "before", "media_type": "image/png",
             "width": 200, "height": 200, "data": "YQ=="},
            {"source": "Output", "label": "after", "media_type": "image/png",
             "width": 150, "height": 150, "data": "Yg=="},
        ])
        dimensions = re.findall(r'<image[^>]+width="([0-9.]+)"[^>]+height="([0-9.]+)"', svg)
        self.assertEqual([("200", "200"), ("150", "150")], dimensions)

    def test_multipolygon_ring_start_and_polygon_order_are_canonicalized(self):
        self.assertTrue(
            self.tester.values_equal(
                "MULTIPOLYGON(((0 0,1 0,1 1,0 0)),((10 10,11 10,11 11,10 10)))",
                "MULTIPOLYGON(((11 10,11 11,10 10,11 10)),((1 0,1 1,0 0,1 0)))",
            )
        )

    def test_geometrycollection_polygon_members_are_canonicalized(self):
        self.assertTrue(
            self.tester.values_equal(
                "GEOMETRYCOLLECTION(POLYGON((0 0,1 0,1 1,0 0)),POLYGON((10 10,11 10,11 11,10 10)))",
                "GEOMETRYCOLLECTION(POLYGON((11 10,11 11,10 10,11 10)),POLYGON((1 0,1 1,0 0,1 0)))",
            )
        )

    def test_polyhedral_surface_faces_and_ring_starts_are_canonicalized(self):
        self.assertTrue(
            self.tester.values_equal(
                "POLYHEDRALSURFACE(((0 0,1 0,0 1,0 0)),((10 10,11 10,10 11,10 10)))",
                "POLYHEDRALSURFACE(((11 10,10 11,10 10,11 10)),((1 0,0 1,0 0,1 0)))",
            )
        )

    def test_tin_faces_and_ring_starts_are_canonicalized(self):
        self.assertTrue(
            self.tester.values_equal(
                "TIN(((0 0,1 0,0 1,0 0)),((10 10,11 10,10 11,10 10)))",
                "TIN(((11 10,10 11,10 10,11 10)),((1 0,0 1,0 0,1 0)))",
            )
        )

    def test_multiline_expected_wkt_uses_canonical_comparison(self):
        actual = [["GEOMETRYCOLLECTION(POLYGON((1 0,1 1,0 0,1 0)),POLYGON((11 10,11 11,10 10,11 10)))"]]
        expected = [
            ["GEOMETRYCOLLECTION(POLYGON((0 0,1 0,1 1,0 0)),"],
            ["POLYGON((10 10,11 10,11 11,10 10)))"],
        ]
        self.assertTrue(self.tester.rows_equal(actual, expected))

    def test_version_comparison_checks_only_shape(self):
        self.assertTrue(self.tester.version_rows_equal([["3.7.0dev", "GEOS 3.13"]], [["3.6.0", "GEOS 3.12"]]))
        self.assertFalse(self.tester.version_rows_equal([["3.7.0dev"]], [["3.6.0", "GEOS 3.12"]]))

    def test_catalog_comparison_allows_empty_actual_catalog(self):
        expected = [["postgis", "3.7.0dev", "null"]]
        self.assertTrue(self.tester.catalog_rows_equal([], expected))
        self.assertTrue(self.tester.catalog_rows_equal([["postgis", "3.7.0dev", "null"]], expected))
        self.assertFalse(self.tester.catalog_rows_equal([["postgis", "3.7.0dev"]], expected))

    def test_typed_hexwkb_comparison_is_endian_independent(self):
        ndr = "0101000000000000000000F03F0000000000000040"
        xdr = "00000000013FF00000000000004000000000000000"
        queries = []
        self.tester.run_psql_scalar = lambda database, query: queries.append(query) or "t"
        actual = QueryRows([[xdr]], types=["geometry"])
        self.assertTrue(self.tester.typed_rows_equal("manual", actual, [[ndr]]))
        self.assertIn("ST_DumpPoints(actual)", queries[0])
        self.assertIn("greatest(1.0, abs(a.x), abs(e.x))", queries[0])
        self.assertIn("GeometryType(actual) = GeometryType(expected)", queries[0])
        self.assertIn("ST_CoordDim(actual) = ST_CoordDim(expected)", queries[0])
        self.assertIn("ST_Zmflag(actual) = ST_Zmflag(expected)", queries[0])
        self.assertIn("ST_SRID(actual) = ST_SRID(expected)", queries[0])

    def test_typed_geometry_comparison_canonicalizes_unordered_areal_rings(self):
        polygon = (
            "0103000000010000000400000000000000000000000000000000000000"
            "000000000000F03F00000000000000000000000000000000000000000000F03F"
            "00000000000000000000000000000000"
        )
        queries = []
        self.tester.run_psql_scalar = lambda database, query: queries.append(query) or "t"
        actual = QueryRows([[polygon]], types=["geometry"])
        self.assertTrue(self.tester.typed_rows_equal("manual", actual, [[polygon]]))
        self.assertIn("GeometryType(actual) IN ('POLYGON', 'MULTIPOLYGON', 'TRIANGLE')", queries[0])
        self.assertIn("ST_Normalize(actual)", queries[0])
        self.assertIn("generate_subscripts(points, 1)", queries[0])
        self.assertIn("jsonb_agg(face_key ORDER BY face_key)", queries[0])
        self.assertNotIn("array_agg(face_key ORDER BY face_key)", queries[0])

    def test_surface_comparison_canonicalizes_face_order_without_refentry_rules(self):
        query = self.tester.geometry_comparison_query(
            "00", expected_hex="00", actual_type="POLYHEDRALSURFACE"
        )
        self.assertEqual(2, query.count("jsonb_agg(face_key ORDER BY face_key)"))
        self.assertEqual(2, query.count("ST_Dump("))
        self.assertIn("ST_Normalize", query)
        self.assertEqual(2, query.count("ST_AsEWKT(point.geom, 12)"))
        self.assertNotIn("ST_GeomFromEWKT(ST_AsEWKT", query)
        self.assertNotIn("refentry", query.lower())

    def test_geometrycollection_comparison_normalizes_member_order(self):
        query = self.tester.geometry_comparison_query(
            "00", expected_wkt="GEOMETRYCOLLECTION EMPTY", actual_type="GEOMETRYCOLLECTION",
        )
        self.assertIn("ST_AsEWKT(ST_Normalize(actual))", query)
        self.assertIn("ST_AsEWKT(ST_Normalize(expected))", query)
        self.assertNotIn("actual_components AS", query)

    def test_areal_comparison_uses_documented_wkt_coordinate_keys(self):
        query = self.tester.geometry_comparison_query(
            "00", expected_wkt="POLYGON Z EMPTY", actual_type="POLYGON",
        )
        self.assertEqual(2, query.count("generate_subscripts(points, 1)"))
        self.assertEqual(2, query.count("ST_AsEWKT(point.geom, 0)"))
        self.assertIn("points[rotation:array_length(points, 1)]", query)
        self.assertIn("IS NOT DISTINCT FROM", query)
        self.assertNotIn("regexp_replace", query)
        self.assertNotIn("CG_3DDifference", query)

    def test_typed_hexwkb_comparison_accepts_readable_ewkt_expectation(self):
        ndr = "0101000020E6100000000000000000F03F0000000000000040"
        queries = []
        self.tester.run_psql_scalar = lambda database, query: queries.append(query) or "t"
        actual = QueryRows([[ndr]], types=["geometry"])
        self.assertTrue(
            self.tester.typed_rows_equal(
                "manual", actual, [["SRID=4326;POINT(1 2)"]]
            )
        )
        self.assertIn("ST_GeomFromEWKB", queries[0])
        self.assertIn("ST_GeomFromEWKT('SRID=4326;POINT(1 2)')", queries[0])
        self.assertNotIn("ST_AsText", queries[0])
        self.assertIn("ST_DumpPoints(actual)", queries[0])
        self.assertNotIn("regexp_replace", queries[0])

    def test_typed_hexwkb_comparison_accepts_readable_wkt_without_srid(self):
        ndr = "0102000020E610000002000000713D0AD7A38051C08FC2F5285C6F4540D7A3703D0A8751C0713D0AD7A3704540"
        queries = []
        def scalar(database, query):
            queries.append(query)
            if "ST_AsEWKT" in query:
                return "LINESTRING(-70.01 42.87,-70.11 42.88)"
            return "t"
        self.tester.run_psql_scalar = scalar
        actual = QueryRows([[ndr]], types=["geometry"])
        self.assertTrue(
            self.tester.typed_rows_equal(
                "manual", actual, [["LINESTRING(-70.01 42.87,-70.11 42.88)"]]
            )
        )
        self.assertIn("ST_AsEWKT", queries[0])

    def test_typed_hexwkb_comparison_preserves_explicit_srid_zero(self):
        ndr = "0101000000000000000000F03F0000000000000040"
        queries = []
        self.tester.run_psql_scalar = lambda database, query: queries.append(query) or "t"
        actual = QueryRows([[ndr]], types=["geometry"])
        self.assertTrue(
            self.tester.typed_rows_equal(
                "manual", actual, [["SRID=0;POINT(1 2)"]]
            )
        )
        self.assertIn("FALSE AND ST_SRID(expected) = 0", queries[0])

    def test_typed_hexwkb_comparison_accepts_legacy_3d_wkt_without_z_token(self):
        ndr = "010200008002000000010000000000F0BFFFFFFFFFFFFFFFBF0000000000000840020000000000F0BF00000000000010C00000000000000840"
        queries = []
        self.tester.run_psql_scalar = lambda database, query: queries.append(query) or "t"
        actual = QueryRows([[ndr]], types=["geometry"])
        self.assertTrue(
            self.tester.typed_rows_equal(
                "manual", actual, [["LINESTRING(-1 -2 3,-1 -4 3)"]]
            )
        )
        self.assertIn("ST_GeomFromEWKT('LINESTRING(-1 -2 3,-1 -4 3)')", queries[0])

    def test_typed_hexwkb_comparison_accepts_explicit_z_token(self):
        ndr = "010200008002000000010000000000F0BFFFFFFFFFFFFFFFBF0000000000000840020000000000F0BF00000000000010C00000000000000840"
        queries = []
        self.tester.run_psql_scalar = lambda database, query: queries.append(query) or "t"
        actual = QueryRows([[ndr]], types=["geometry"])
        self.assertTrue(
            self.tester.typed_rows_equal(
                "manual", actual, [["LINESTRING Z (-1 -2 3,-1 -4 3)"]]
            )
        )
        self.assertIn("ST_GeomFromEWKT('LINESTRING Z (-1 -2 3,-1 -4 3)')", queries[0])

    def test_rounded_typed_hexwkb_comparison_keeps_explicit_z_spelling(self):
        ndr = "010200008002000000010000000000F0BFFFFFFFFFFFFFFFBF0000000000000840020000000000F0BF00000000000010C00000000000000840"
        queries = []

        def scalar(database, query):
            queries.append(query)
            if "ST_AsText" in query:
                return "LINESTRING Z (-1.000 -2.000 3.000,-1.000 -4.000 3.000)"
            return "f"

        self.tester.run_psql_scalar = scalar
        actual = QueryRows([[ndr]], types=["geometry"])
        self.assertTrue(
            self.tester.typed_rows_equal(
                "manual",
                actual,
                [["LINESTRING Z (-1.000 -2.000 3.000,-1.000 -4.000 3.000)"]],
            )
        )
        self.assertIn("ST_AsText", queries[0])
        self.assertNotIn("ST_AsEWKT", queries[0])

    def test_typed_hexwkb_comparison_accepts_compact_ewkt_dimension_marker(self):
        point_m = "0101000040000000000000144000000000000028400000000000002640"
        queries = []
        self.tester.run_psql_scalar = lambda database, query: queries.append(query) or "t"
        actual = QueryRows([[point_m]], types=["geometry"])
        self.assertTrue(
            self.tester.typed_rows_equal(
                "manual", actual, [["POINTM(5 12 11)"]]
            )
        )
        self.assertIn("ST_GeomFromEWKT('POINTM(5 12 11)')", queries[0])
        self.assertNotIn("regexp_replace", queries[0])

    def test_typed_geometry_comparison_rejects_coordinate_mismatch(self):
        ndr = "010200000002000000000000000000F03F000000000000004000000000000008400000000000001040"
        queries = []
        self.tester.run_psql_scalar = lambda database, query: queries.append(query) or "f"
        actual = QueryRows([[ndr]], types=["geometry"])
        self.assertFalse(
            self.tester.typed_rows_equal(
                "manual", actual, [["LINESTRING(1 2,3 4)"]]
            )
        )
        self.assertEqual(1, len(queries))
        self.assertIn("ST_DumpPoints(actual)", queries[0])
        self.assertIn("greatest(1.0, abs(a.x), abs(e.x))", queries[0])

    def test_typed_nurbs_comparison_stays_native(self):
        nurbs = "011500000002000000030000000000000000000000000000000000000000000000000014400000000000002440000000000000000000000000"
        queries = []
        self.tester.run_psql_scalar = lambda database, query: queries.append(query) or "f"
        actual = QueryRows([[nurbs]], types=["geometry"])
        self.assertFalse(
            self.tester.typed_rows_equal(
                "manual", actual, [[nurbs]]
            )
        )
        self.assertEqual(1, len(queries))
        self.assertIn("ST_AsHEXEWKB(actual, 'XDR')", queries[0])
        self.assertIn("ST_AsHEXEWKB(expected, 'XDR')", queries[0])
        self.assertNotIn("ST_AsEWKT", queries[0])
        self.assertNotIn("ST_CurveToLine", queries[0])
        self.assertNotIn("ST_DumpPoints", queries[0])

    def test_typed_linestring_comparison_does_not_normalize_direction(self):
        ndr = "010200000002000000000000000000F03F000000000000004000000000000008400000000000001040"
        queries = []
        self.tester.run_psql_scalar = lambda database, query: queries.append(query) or "f"
        actual = QueryRows([[ndr]], types=["geometry"])
        self.assertFalse(
            self.tester.typed_rows_equal(
                "manual", actual, [["LINESTRING(3 4,1 2)"]]
            )
        )
        self.assertEqual(1, len(queries))
        self.assertIn(
            "CASE WHEN GeometryType(actual) IN "
            "('POLYGON', 'MULTIPOLYGON', 'TRIANGLE') THEN ST_Normalize(actual)",
            queries[0],
        )
        self.assertNotIn("ST_Reverse", queries[0])

    def test_typed_geometry_comparison_rejects_semantic_mismatch(self):
        ndr = "0101000020E6100000000000000000F03F0000000000000040"
        self.tester.run_psql_scalar = lambda database, query: "f"
        actual = QueryRows([[ndr]], types=["geometry"])
        self.assertFalse(
            self.tester.typed_rows_equal(
                "manual", actual, [["SRID=3857;POINT(1 2)"]]
            )
        )

    def test_typed_hexwkb_comparison_accepts_wrapped_readable_ewkt(self):
        ndr = "0101000020E6100000000000000000F03F0000000000000040"
        queries = []
        self.tester.run_psql_scalar = lambda database, query: queries.append(query) or "t"
        actual = QueryRows([[ndr]], types=["geometry"])

        self.assertTrue(
            self.tester.typed_rows_equal(
                "manual",
                actual,
                [["SRID=4326;POINT(1"], [" 2)"]],
            )
        )
        self.assertIn("ST_GeomFromEWKB", queries[0])
        self.assertIn("ST_DumpPoints(actual)", queries[0])

    def test_readable_geometry_expectation_is_detected(self):
        self.assertTrue(self.tester.expected_has_geometry_text([["POINT(1 2)"]]))
        self.assertTrue(self.tester.expected_has_geometry_text([["POLYGON EMPTY"]]))
        self.assertFalse(self.tester.expected_has_geometry_text([["not geometry"]]))

    def test_documented_visual_types_detect_readable_wkt_columns(self):
        self.assertEqual(
            ["geometry", ""],
            self.tester.documented_visual_types([
                ["POINT(1 2)", "label"],
                ["SRID=4326;POLYGON EMPTY", "still text"],
            ]),
        )


class ExampleTestRunnerTest(unittest.TestCase):
    def make_runner(self, actual_rows):
        tester = ExampleTester.__new__(ExampleTester)
        examples = []
        for label, query, expected in (
            ("first:1", "SELECT 1", [["1"]]),
            ("second:2", "SELECT 2", [["2"]]),
            ("third:3", "SELECT 3", [["3"]]),
        ):
            examples.append(
                {
                    "label": label,
                    "query": query,
                    "expected": expected,
                    "valid": True,
                    "catalog": False,
                    "version": False,
                    "visual_id": None,
                }
            )
        tester.examples = lambda: examples
        tester.check_environment = lambda database: None
        tester.run_psql_query = lambda database, query: actual_rows[query]
        return tester

    def test_run_examples_keep_going_reports_all_failures(self):
        tester = self.make_runner(
            {
                "SELECT 1": [["bad"]],
                "SELECT 2": [["2"]],
                "SELECT 3": [["also bad"]],
            }
        )

        stderr = io.StringIO()
        with self.assertRaisesRegex(RuntimeError, r"FAILED 2 example\(s\): first:1, third:3"):
            with contextlib.redirect_stderr(stderr):
                tester.run_examples("exampletest", keep_going=True)

        output = stderr.getvalue()
        self.assertIn("Example test failed: first:1", output)
        self.assertIn("Example test failed: third:3", output)
        self.assertIn("Actual:\nbad", output)
        self.assertIn("Actual:\nalso bad", output)

    def test_run_examples_without_keep_going_still_fails_fast(self):
        tester = self.make_runner(
            {
                "SELECT 1": [["bad"]],
                "SELECT 2": [["wrong but not reached"]],
                "SELECT 3": [["also not reached"]],
            }
        )

        with self.assertRaisesRegex(RuntimeError, "Example test failed: first:1") as raised:
            tester.run_examples("exampletest")

        self.assertNotIn("second:2", str(raised.exception))

    def test_visual_only_skips_unmarked_examples_and_global_environment_checks(self):
        tester = self.make_runner({"SELECT 1": [["1"]], "SELECT 2": [["2"]], "SELECT 3": [["3"]]})
        examples = tester.examples()
        examples[1]["visual_id"] = "second-visual"
        calls = []
        tester.check_environment = lambda database: self.fail("visual-only checked unrelated environment")
        tester.run_one_example = lambda database, example: calls.append(example["label"]) or example["expected"]
        tester.render_visual_example = lambda database, example, actual: {
            "id": example["visual_id"], "source": example["label"], "svg": "<svg/>\n"
        }
        tester.publish_visual_examples = lambda render_dir, visuals: calls.append(
            (render_dir, [visual["id"] for visual in visuals])
        )

        tester.run_examples("exampletest", render_dir="visuals", visual_only=True)

        self.assertEqual(["second:2", ("visuals", ["second-visual"])], calls)

    def test_visual_only_renders_typed_documented_output_at_any_job_count(self):
        for expected in (
            [["010300000000000000"]],
            [["POLYGON EMPTY"]],
        ):
            for jobs in (1, 2):
                with self.subTest(expected=expected, jobs=jobs):
                    tester = self.make_runner({"SELECT 1": expected})
                    example = tester.examples()[0]
                    example.update({
                        "expected": expected,
                        "documented_only": True,
                        "expected_headers": ["documented_geom"],
                        "visual_id": "documented-visual",
                        "visual_kind": "explicit",
                        "visual_preferred": True,
                    })
                    tester.examples = lambda: [example]
                    calls = []
                    tester.run_one_example = lambda database, selected: self.fail(
                        "documented output was executed"
                    )
                    tester.render_visual_example = lambda database, selected, actual: calls.append(
                        (selected["visual_id"], actual)
                    ) or {"id": selected["visual_id"], "source": selected["label"], "svg": "<svg/>\n"}
                    tester.publish_visual_examples = lambda render_dir, visuals: calls.append(
                        (render_dir, [visual["id"] for visual in visuals])
                    )

                    tester.run_examples("exampletest", render_dir="visuals", visual_only=True, jobs=jobs)

                    self.assertEqual("documented-visual", calls[0][0])
                    self.assertEqual(example["expected"], calls[0][1])
                    self.assertEqual(["documented_geom"], calls[0][1].headers)
                    self.assertEqual(["geometry"], calls[0][1].types)
                    self.assertEqual(("visuals", ["documented-visual"]), calls[1])

    def test_visual_only_preserves_documented_scalar_output_types(self):
        tester = self.make_runner({"SELECT 1": [["label", "POINT(1 2)"]]})
        example = tester.examples()[0]
        example.update({
            "expected": [["label", "POINT(1 2)"]],
            "documented_only": True,
            "expected_headers": ["caption", "documented_geom"],
            "visual_id": "documented-visual",
            "visual_kind": "explicit",
            "visual_preferred": True,
        })
        tester.examples = lambda: [example]
        calls = []
        tester.run_one_example = lambda database, selected: self.fail(
            "documented output was executed"
        )
        tester.render_visual_example = lambda database, selected, actual: calls.append(actual)
        tester.publish_visual_examples = lambda render_dir, visuals: None

        tester.run_examples("exampletest", render_dir="visuals", visual_only=True)

        self.assertEqual(["", "geometry"], calls[0].types)

    def test_documented_visual_render_failure_names_the_example(self):
        tester = self.make_runner({"SELECT 1": [["POLYGON EMPTY"]]})
        example = tester.examples()[0]
        example.update({
            "expected": [["POLYGON EMPTY"]],
            "documented_only": True,
            "expected_headers": ["documented_geom"],
            "visual_id": "documented-visual",
            "visual_kind": "explicit",
            "visual_preferred": True,
        })
        tester.examples = lambda: [example]

        def fail_render(database, selected, actual):
            raise RuntimeError("GEOS overlay failed")

        tester.render_visual_example = fail_render

        with self.assertRaisesRegex(
            RuntimeError,
            r"(?s)Documented visual failed to render: first:1 \(documented-visual\).*GEOS overlay failed",
        ):
            tester.run_examples("exampletest", render_dir="visuals", visual_only=True)

    def test_native_geometry_output_is_described_when_input_uses_serializer(self):
        hexwkb = "010200000002000000000000000000F03F0000000000000040" \
                 "00000000000008400000000000001040"
        tester = ExampleTester.__new__(ExampleTester)
        example = {
            "label": "ST_LineFromWKB:1",
            "query": "SELECT ST_LineFromWKB(ST_AsBinary(" \
                     "ST_GeomFromText('LINESTRING(1 2,3 4)'))) AS aline",
            "expected": [["LINESTRING(1 2,3 4)"]],
            "valid": True,
            "catalog": False,
            "version": False,
            "visual_id": "st-linefromwkb-1",
            "visual_kind": "geometry-output",
        }
        described = []
        tester.run_psql_query = lambda database, query: QueryRows(
            [[hexwkb]], ["aline"]
        )
        tester.query_output_types = lambda database, query: described.append(query) or ["geometry"]
        tester.run_psql_scalar = lambda database, query: "t"

        actual = tester.run_one_example("manual", example)

        self.assertEqual([example["query"]], described)
        self.assertEqual(["geometry"], actual.types)


class VisualExampleTest(unittest.TestCase):
    def setUp(self):
        self.tester = ExampleTester.__new__(ExampleTester)

    def point_radius(self, svg):
        match = re.search(r'class="point" d="M [^"]+ A ([0-9.]+) ([0-9.]+) 0 1 0', svg)
        self.assertIsNotNone(match)
        return float(match.group(1))

    def vertex_radius(self, svg):
        match = re.search(r'class="vertex"[^>]* r="([0-9.]+)"', svg)
        self.assertIsNotNone(match)
        return float(match.group(1))

    def test_svg_uses_postgis_fragments_and_stable_source_colors(self):
        svg = self.tester.visual_svg(
            "buffer-example",
            {
                "bounds": [0, 0, 10, 10],
                "parts": [
                    {"ord": 1, "source": "Code", "label": "Code 1", "type": "POINT", "x": 2, "y": -3},
                    {"ord": 2, "source": "Output", "label": "Output 1", "type": "POLYGON", "svg": "M 0 0 L 10 0 10 -10 Z"},
                ],
            },
        )
        self.assertIn('viewBox="0 0 600 380"', svg)
        self.assertIn('class="plot-grid"', svg)
        self.assertIn('pointer-events:stroke', svg)
        self.assertIn('pointer-events:all', svg)
        self.assertIn('data-postgis-geometry-id="buffer-example-code-1"', svg)
        self.assertRegex(
            svg,
            r'<path class="point" d="M [^\"]+ A [^\"]+ Z" fill="#2878b8" '
            r'stroke="white" stroke-width="[0-9.]+"/>',
        )
        self.assertIn(
            'd="M 0 0 L 10 0 10 -10 Z" stroke="#a62c2b" fill="#a62c2b"',
            svg,
        )
        self.assertIn('fill-opacity="0.24"', svg)
        self.assertIn('stroke-opacity="1"', svg)
        self.assertRegex(svg, r'fill-rule="evenodd" stroke-width="[0-9.]+"')
        self.assertIn('.geometry-layer.active{filter:brightness(.72)}', svg)
        self.assertIn('svg:has(.geometry-layer.active) .geometry-layer:not(.active){opacity:.18}', svg)
        self.assertIn('stroke="#dce2e7" stroke-width="1"', svg)
        self.assertIn('<text class="grid-label grid-label-x"', svg)
        self.assertIn('<text class="grid-label grid-label-y"', svg)
        self.assertIn('>0</text>', svg)
        self.assertIn('>10</text>', svg)
        self.assertIn('<rect width="11" height="11" rx="2" fill="#2878b8"/>', svg)
        self.assertIn('<rect width="11" height="11" rx="2" fill="#a62c2b"/>', svg)
        self.assertLess(
            svg.index('data-postgis-geometry-id="buffer-example-output-1"'),
            svg.index('data-postgis-geometry-id="buffer-example-code-1"'),
        )

    def test_svg_matrix_3x3_keeps_context_separate_from_output_grid(self):
        frames = [
            {"id": frame_id, "bounds": [0, 0, 1, 1]}
            for frame_id in ["a", "b"] + [f"cell-{index}" for index in range(9)]
        ]
        parts = [
            {"ord": 1, "source": "Code", "label": "a", "frame": "a", "type": "POINT", "x": 0, "y": 0},
            {"ord": 2, "source": "Code", "label": "b", "frame": "b", "type": "POINT", "x": 1, "y": 1},
        ]
        parts.extend(
            {
                "ord": index + 3,
                "source": "Output",
                "label": f"cell {index}",
                "frame": f"cell-{index}",
                "type": "POINT",
                "x": 0.5,
                "y": 0.5,
            }
            for index in range(9)
        )

        svg = self.tester.visual_svg(
            "matrix-cells",
            {"bounds": [0, 0, 1, 1], "frames": frames, "parts": parts},
            matrix_3x3=True,
        )

        backgrounds = re.findall(
            r'<rect class="plot-background" x="([^"]+)" y="([^"]+)"',
            svg,
        )
        self.assertEqual(11, len(backgrounds))
        row_counts = Counter(y for _x, y in backgrounds)
        self.assertEqual([2, 3, 3, 3], sorted(row_counts.values()))
        self.assertIn('viewBox="0 0 900 ', svg)

    def test_visual_candidate_includes_spatial_context_but_skips_trivial_points(self):
        self.assertIsNone(self.tester.visual_candidate("SELECT 'POINT(1 2)'", [["1"]]))
        self.assertIsNone(self.tester.visual_candidate(
            "SELECT ST_IsCollection('GEOMETRYCOLLECTION(POINT(0 0))')", [["t"]]
        ))
        self.assertIsNone(self.tester.visual_candidate(
            "SELECT ST_IsCollection("
            "'GEOMETRYCOLLECTION(MULTIPOINT((0 0),(1 1)),GEOMETRYCOLLECTION(POINT(2 2)))')",
            [["t"]],
        ))
        context = self.tester.visual_candidate(
            "SELECT ST_Intersects('LINESTRING(0 0,1 1)', 'POLYGON((0 0,2 0,2 2,0 0))')",
            [["t"]],
        )
        self.assertEqual({"kind": "input-context", "preferred": False}, context)
        collection_context = self.tester.visual_candidate(
            "SELECT ST_IsCollection("
            "'GEOMETRYCOLLECTION(POINT(0 0),LINESTRING(0 0,1 1))')",
            [["t"]],
        )
        self.assertEqual({"kind": "input-context", "preferred": False}, collection_context)
        output = self.tester.visual_candidate(
            "SELECT ST_Buffer('LINESTRING(0 0,1 1)', 1)",
            [["POLYGON((0 0,2 0,2 2,0 0))"]],
        )
        self.assertEqual({"kind": "geometry-output", "preferred": True}, output)
        point_output = self.tester.visual_candidate(
            "SELECT ST_Buffer('POINT(0 0)', 1)",
            [["POLYGON((0 0,1 0,1 1,0 0))"]],
        )
        self.assertEqual({"kind": "geometry-output", "preferred": True}, point_output)
        self.assertIsNone(self.tester.visual_candidate(
            "SELECT ST_AsText('LINESTRING(0 0,1 1)')", [["LINESTRING(0 0,1 1)"]]
        ))

    def test_visual_candidate_recognizes_natural_geometry_hex_output(self):
        polygon_hex = "010300000000000000"
        self.assertEqual(
            {"kind": "geometry-output", "preferred": True},
            self.tester.visual_candidate(
                "SELECT ST_PolygonFromText('POLYGON EMPTY')",
                [[polygon_hex]],
            ),
        )
        self.assertIsNone(self.tester.visual_candidate(
            "SELECT ST_AsHEXEWKB('POLYGON EMPTY'::geometry)",
            [[polygon_hex]],
        ))
        self.assertIsNone(self.tester.visual_candidate(
            "SELECT ST_PolygonFromText('POINT(1 2)')",
            [["0101000000000000000000F03F0000000000000040"]],
        ))

    def test_visual_candidate_does_not_draw_geography_inputs_as_planar(self):
        geography_context = self.tester.visual_candidate(
            "SELECT ST_Distance('LINESTRING(-122 47,0 51)'::geography, "
            "'POINT(-21 64)'::geography)",
            [["122235"]],
        )
        self.assertIsNone(geography_context)

        geometry_context = self.tester.visual_candidate(
            "SELECT ST_Distance('LINESTRING(-122 47,0 51)'::geometry, "
            "'POINT(-21 64)'::geometry)",
            [["13.3"]],
        )
        self.assertEqual({"kind": "input-context", "preferred": False}, geometry_context)

        geography_output = self.tester.visual_candidate(
            "SELECT ST_AsText(ST_Segmentize('LINESTRING(0 0,60 60)'::geography, 2000000))",
            [["LINESTRING(0 0,10 12,60 60)"]],
        )
        self.assertEqual({"kind": "geometry-output", "preferred": True}, geography_output)

    def test_code_geometry_candidates_include_constructed_context_in_sql_order(self):
        candidates = self.tester.code_geometry_candidates(
            "SELECT f('LINESTRING(0 0,10 10)', ST_Point(1, 2), "
            "ST_MakeEnvelope(2, 3, 8, 9, 4326))"
        )
        self.assertEqual([None, "Point", "Envelope"], [candidate.get("label") for candidate in candidates])
        self.assertEqual("LINESTRING(0 0,10 10)", candidates[0]["wkt"])
        self.assertEqual("POINT(1 2)", candidates[1]["wkt"])
        self.assertEqual(
            "SRID=4326;POLYGON((2 3,8 3,8 9,2 9,2 3))",
            candidates[2]["wkt"],
        )

    def test_code_geometry_candidates_use_cast_expression_aliases(self):
        candidates = self.tester.code_geometry_candidates(
            "SELECT 'LINESTRING(0 0,1 1)'::geometry AS line_a, "
            "'POINT(2 2)'::geography point_b, "
            "'POLYGON((0 0,1 0,0 0))'::geometry AS \"lake \"\"inner\"\" shore\""
        )
        self.assertEqual(
            ["line_a", "point_b", 'lake "inner" shore'],
            [candidate.get("label") for candidate in candidates],
        )

    def test_query_output_headers_use_only_outer_select_aliases(self):
        headers = self.tester.query_output_headers(
            "WITH data AS (SELECT 'POINT(0 0)'::geometry AS inner_geom), "
            "projected AS (SELECT inner_geom AS other_inner FROM data) "
            "SELECT ST_AsText(other_inner) AS input_geometry, "
            "ST_AsText(ST_Buffer(other_inner, 1)) AS result_geometry FROM projected"
        )
        self.assertEqual(["input_geometry", "result_geometry"], headers)

    def test_documented_visual_uses_query_aliases_for_layer_roles(self):
        captured = []
        self.tester.visual_payload = lambda database, layers: captured.extend(layers) or {
            "bounds": [0, 0, 1, 1],
            "parts": [{
                "ord": layer["ord"], "source": layer["source"], "label": layer["label"],
                "type": "POINT", "x": 0, "y": 0,
            } for layer in layers],
        }
        actual = QueryRows([["POINT(0 0)", "POINT(1 1)"]])
        self.tester.render_visual_example("manual", {
            "query": (
                "WITH data AS (SELECT 1) "
                "SELECT 'POINT(0 0)' AS input_center, 'POINT(1 1)' AS result FROM data"
            ),
            "visual_id": "documented-aliases", "label": "aliases:documented",
            "visual_refentry": "aliases", "visual_screen": 1,
            "visual_preferred": True, "visual_kind": "explicit",
        }, actual)
        self.assertEqual(["Code", "Output"], [layer["source"] for layer in captured])
        self.assertEqual(["center", "result"], [layer["label"] for layer in captured])

    def test_separate_output_visual_requests_one_frame_per_output_label(self):
        captured = []
        self.tester.visual_payload = lambda database, layers: captured.extend(layers) or {
            "bounds": [0, 0, 1, 1],
            "frames": [
                {"id": "Code", "bounds": [0, 0, 1, 1]},
                {"id": "cleaned", "bounds": [0, 0, 1, 1]},
                {"id": "invalid edges", "bounds": [0, 0, 1, 1]},
            ],
            "parts": [{
                "ord": layer["ord"], "source": layer["source"], "label": layer["label"],
                "frame": layer["requested_frame"], "type": "LINESTRING",
                "svg": "M 0 0 L 1 -1",
            } for layer in layers],
        }
        actual = QueryRows(
            [["LINESTRING(0 0,1 1)", "LINESTRING(0 1,1 0)"]],
            headers=["cleaned", "invalid_edges"],
        )
        self.tester.render_visual_example("manual", {
            "query": "SELECT 'POLYGON((0 0,1 0,1 1,0 0))'::geometry",
            "visual_id": "separate-output", "label": "comparison:documented",
            "visual_refentry": "comparison", "visual_screen": 1,
            "visual_preferred": True, "visual_kind": "explicit",
            "visual_separate_output": True,
        }, actual)
        self.assertEqual(
            ["Code", "cleaned", "invalid edges"],
            [layer["requested_frame"] for layer in captured],
        )

    def test_overlay_visual_requests_one_shared_frame(self):
        captured = []
        self.tester.visual_payload = lambda database, layers: captured.extend(layers) or {
            "bounds": [0, 0, 1, 1],
            "parts": [{
                "ord": layer["ord"], "source": layer["source"], "label": layer["label"],
                "type": "POINT", "x": 0, "y": 0,
            } for layer in layers],
        }
        actual = QueryRows([["POINT(1 1)"]], headers=["center"])
        self.tester.render_visual_example("manual", {
            "query": "SELECT 'POLYGON((0 0,1 0,1 1,0 0))'::geometry",
            "visual_id": "overlay", "label": "overlay:documented",
            "visual_refentry": "overlay", "visual_screen": 1,
            "visual_preferred": True, "visual_kind": "explicit",
            "visual_overlay": True,
        }, actual)
        self.assertEqual(["Overlay", "Overlay"], [layer["requested_frame"] for layer in captured])

    def test_visual_payload_separates_same_dimension_or_cross_srid_sources(self):
        queries = []
        self.tester.run_psql_scalar = lambda database, query: queries.append(query) or '{"frames":[],"parts":[]}'
        self.tester.visual_payload("manual", [])
        self.assertIn("source_dimensions", queries[0])
        self.assertIn("count(DISTINCT ST_SRID(geom))", queries[0])
        self.assertIn("AS shared_coordinate_space", queries[0])
        self.assertIn("bool_and(shared_coordinate_space)", queries[0])
        self.assertIn("CASE WHEN separate_sources THEN source ELSE 'Overlay' END", queries[0])
        self.assertIn("COALESCE(requested_frame", queries[0])
        self.assertIn("ST_Translate(geom, -min_x, -min_y)", queries[0])
        self.assertNotIn("frame_bounds_geoms AS", queries[0])
        self.assertNotIn("ST_CurveToLine(ST_Force2D(geom))", queries[0])
        self.assertIn("ST_Extent(ST_Force2D(all_geometries.geom))", queries[0])
        self.assertIn("ELSE ST_Extent(ST_Force2D(geom))", queries[0])
        self.assertIn("ELSE ST_AsSVG(ST_Force2D(render_geom), 0, 12)", queries[0])
        self.assertIn("WHEN ST_Dimension(render_geom) = 0 THEN NULL", queries[0])
        self.assertIn("'points'", queries[0])
        self.assertIn("ST_DumpPoints(render_geom) AS point", queries[0])
        self.assertIn("ST_Zmflag(raw_geom) IN (2, 3)", queries[0])
        self.assertIn("'has_z'", queries[0])
        self.assertIn("'points_xyz'", queries[0])
        self.assertIn("WITH RECURSIVE raw AS", queries[0])
        self.assertIn("FROM ST_Dump(parent.geom) AS dumped", queries[0])
        self.assertIn(
            "parent.part_type IN ('ST_PolyhedralSurface', 'ST_Tin')",
            queries[0],
        )
        self.assertIn("ST_GeometryN(parent.geom, part_num.n)", queries[0])
        self.assertIn("FROM geometry_tree", queries[0])
        self.assertIn("parent.part_path || child.path AS part_path", queries[0])
        self.assertIn("ORDER BY ord, part_path", queries[0])
        self.assertIn("'ST_MultiCurve', 'ST_MultiSurface'", queries[0])
        self.assertIn("'ST_Polygon', 'ST_CurvePolygon'", queries[0])
        self.assertIn("THEN ST_CurveToLine(raw_geom)", queries[0])
        self.assertNotIn("ST_Dump(framed.geom)", queries[0])
        self.assertNotIn("WHERE GeometryType(framed.geom) = 'CURVEPOLYGON'", queries[0])
        self.assertNotIn("WHERE GeometryType(framed.geom) <> 'CURVEPOLYGON'", queries[0])
        self.assertIn("geom AS render_geom", queries[0])
        self.assertNotIn("ST_HasArc(geom) OR GeometryType(geom) IN ('CURVEPOLYGON', 'MULTISURFACE')", queries[0])
        self.assertIn("ST_GeomFromEWKB(decode(hexwkb, 'hex'))", queries[0])
        self.assertIn("ELSE ST_GeomFromEWKT(wkt)", queries[0])
        self.assertNotIn("ST_AsEWKB(result.geom) = ST_AsEWKB(candidate.geom)", queries[0])
        self.assertIn("COALESCE(ST_Z(vertex3d.geom), 0)", queries[0])
        self.assertIn("ST_Z(vertex3d.geom)", queries[0])
        self.assertIn("COALESCE(total_points, ST_NPoints(framed.geom)) AS total_points", queries[0])
        self.assertIn(
            "upper(substring(ST_GeometryType(framed.geom) from 4)) AS root_type",
            queries[0],
        )
        self.assertIn(
            "'type', upper(substring(ST_GeometryType(geom) from 4))",
            queries[0],
        )
        self.assertIn("ST_GeometryType(geom) = 'ST_Point'", queries[0])
        self.assertIn(
            "ST_GeometryType(render_geom) IN ('ST_Polygon', 'ST_Triangle')",
            queries[0],
        )
        self.assertNotIn("GeometryType(framed.geom) AS root_type", queries[0])
        self.assertIn("geom AS raw_geom", queries[0])
        self.assertIn(") AS vertex3d", queries[0])
        self.assertIn("'root_type', root_type", queries[0])
        self.assertIn("'source_point_count', source_point_count", queries[0])
        self.assertIn("'marker_scale', marker_scale", queries[0])
        self.assertIn("root_type IN ('LINESTRING', 'POLYGON', 'TRIANGLE')", queries[0])
        self.assertIn("AND total_points <= 64", queries[0])
        self.assertNotIn("AND ST_NPoints(geom) <= 64", queries[0])
        self.assertNotIn("ST_Scale", queries[0])

    def test_visual_payload_sends_unicode_labels_without_json_escapes(self):
        queries = []
        self.tester.run_psql_scalar = lambda database, query: queries.append(query) or '{"frames":[],"parts":[]}'
        self.tester.visual_payload("manual", [{
            "ord": 1,
            "source": "Output",
            "label": "2 I(a) ∩ I(b)",
            "row_num": 1,
            "column_num": 1,
            "requested_frame": None,
            "wkt": "POINT(0 0)",
            "hexwkb": None,
            "source_point_count": 1,
            "total_points": 1,
            "marker_scale": None,
        }])
        match = re.search(r"decode\('([^']+)', 'base64'\)", queries[0])
        self.assertIsNotNone(match)
        payload = base64.b64decode(match.group(1)).decode("utf-8")
        self.assertIn("∩", payload)
        self.assertNotIn(r"\u2229", payload)

    def test_3d_edge_splitter_uses_postgis_overlay(self):
        queries = []
        self.tester.run_psql_scalar = lambda database, query: queries.append(query) or "[]"
        self.tester.split_3d_edge_segments("manual", [{
            "is_3d_face": True,
            "fill_3d_face": True,
            "face_3d": {
                "rings": [[[0, 0, 0], [2, 0, 0], [2, 2, 0], [0, 2, 0], [0, 0, 0]]],
                "depth_plane": (0, 0, 0),
                "edge_keys": [(("0", "0", "0"), ("2", "0", "0"))],
            },
            "edges_3d": [{
                "key": (("0", "0", "0"), ("2", "0", "0")),
                "start": [0, 0],
                "end": [2, 0],
                "start_depth": -1,
                "end_depth": -1,
                "visible_boundary": True,
            }],
        }], {
            0: [{
                "wkt": "POLYGON((0 0,2 0,2 2,0 2,0 0))",
                "svg": "M 0 0 L 2 0 L 2 -2 L 0 -2 Z",
                "draw": False,
                "edge_visibility": True,
            }, {
                "wkt": "POLYGON((0 0,1 0,1 1,0 1,0 0))",
                "svg": "M 0 0 L 1 0 L 1 -1 L 0 -1 Z",
                "draw": True,
                "edge_visibility": False,
            }],
        })
        encoded = re.search(r"decode\('([^']+)', 'base64'\)", queries[0]).group(1)
        edge_payload = json.loads(base64.b64decode(encoded))
        self.assertEqual(1, len(edge_payload["visible_patches"]))
        self.assertEqual(
            "POLYGON((0 0,2 0,2 2,0 2,0 0))",
            edge_payload["visible_patches"][0]["wkt"],
        )
        self.assertIn("ST_Split", queries[0])
        self.assertIn("ST_Intersection(e.geom, ST_Boundary(f.geom))", queries[0])
        self.assertIn("ST_CollectionExtract(geom, 1)", queries[0])
        self.assertIn("ST_Boundary(ST_CollectionExtract(geom, 2))", queries[0])
        self.assertIn("depth_crossings AS", queries[0])
        self.assertIn("delta.start_delta / (delta.start_delta - delta.end_delta)", queries[0])
        self.assertIn("ST_LineInterpolatePoint(e.geom, values.crossing_fraction)", queries[0])
        self.assertIn("visible_face_patches AS", queries[0])
        self.assertIn("data->'visible_patches'", queries[0])
        self.assertIn("visible_patch_edges AS", queries[0])
        self.assertIn("face_projection_union AS", queries[0])
        self.assertIn("face_projection_boundary AS MATERIALIZED", queries[0])
        self.assertIn("drawable_edges AS MATERIALIZED", queries[0])
        self.assertIn("ST_Intersection(e.geom, b.geom, p.grid_size)", queries[0])
        self.assertIn("WHERE (visible OR NOT silhouette_only)", queries[0])
        self.assertIn(
            "NOT silhouette_only OR ST_Length(geom) > 5 * p.cleanup_grid_size",
            queries[0],
        )
        self.assertIn("ST_LineLocatePoint(e.geom, ST_StartPoint((dumped).geom))", queries[0])
        self.assertIn("raw_faces AS", queries[0])
        self.assertIn("faces AS MATERIALIZED", queries[0])
        self.assertIn("visible_face_patches AS MATERIALIZED", queries[0])
        self.assertIn("ST_CollectionExtract(ST_MakeValid(f.geom), 3)", queries[0])
        self.assertIn("projected_extents AS", queries[0])
        self.assertIn("overlay_precision AS MATERIALIZED", queries[0])
        self.assertIn("scene_span * 1e-4 AS grid_size", queries[0])
        self.assertIn("scene_span * 4e-4 AS cleanup_grid_size", queries[0])
        self.assertIn("ST_Intersection(\n          s.geom,\n          f.geom,\n          p.grid_size", queries[0])
        self.assertIn("ON s.frame_id = f.frame_id", queries[0])
        self.assertIn("s.edge_key = f.edge_key", queries[0])
        self.assertNotIn("AND ST_Intersects(s.geom, f.geom)", queries[0])
        self.assertIn("ST_UnaryUnion(ST_Collect(v.geom), p.cleanup_grid_size)", queries[0])
        self.assertIn("JOIN face_projection_union AS f USING (frame_id, group_id)", queries[0])
        self.assertIn("ST_Difference(s.geom, v.geom, p.grid_size)", queries[0])
        self.assertIn("ST_Difference(h.geom, v.geom, p.cleanup_grid_size)", queries[0])
        self.assertNotIn("2e-2", queries[0])
        self.assertNotIn("8e-2", queries[0])
        self.assertNotIn("ST_SnapToGrid", queries[0])
        self.assertNotIn("ST_Buffer", queries[0])
        self.assertIn("straightened AS", queries[0])
        self.assertIn("ST_LineSubstring(", queries[0])
        self.assertIn(
            "start_fraction * ST_Length(original_geom) <= p.cleanup_grid_size",
            queries[0],
        )
        self.assertIn(
            "(1 - end_fraction) * ST_Length(original_geom) <= p.cleanup_grid_size",
            queries[0],
        )
        self.assertIn("ST_LineMerge(ST_UnaryUnion(ST_Collect(geom)))", queries[0])
        self.assertIn("FROM classified_segments", queries[0])
        self.assertIn("ST_LineLocatePoint(original_geom", queries[0])
        self.assertIn("plane_a * ST_X", queries[0])
        self.assertIn("'visible', visible", queries[0])
        self.assertNotIn("'visible', visible_boundary AND", queries[0])

    def test_ordinary_2d_svg_remains_byte_stable(self):
        svg = self.tester.visual_svg(
            "buffer-example",
            {
                "bounds": [0, 0, 10, 10],
                "parts": [
                    {"ord": 1, "source": "Code", "label": "Code 1",
                     "type": "POINT", "x": 2, "y": -3},
                    {"ord": 2, "source": "Output", "label": "Output 1",
                     "type": "POLYGON", "svg": "M 0 0 L 10 0 10 -10 Z",
                     "has_z": False,
                     "points_xyz": xyz_vertices(
                         [[0, 0, 0], [10, 0, 0], [10, 10, 0], [0, 0, 0]], ring=True
                     )},
                ],
            },
        )
        self.assertEqual(
            "ad9f03873ee58dbb970b5fab66598bf146db51f76331217814abeddbff05ae0b",
            hashlib.sha256(svg.encode()).hexdigest(),
        )

    def test_z_surface_faces_are_projected_back_to_front_without_culling_fill(self):
        payload = {
            "bounds": [0, 0, 2, 1],
            "parts": [
                {"ord": 1, "source": "Output", "label": "solid",
                 "root_type": "POLYHEDRALSURFACE", "type": "POLYGON",
                 "svg": "unused",
                 "has_z": True,
                 "points_xyz": xyz_vertices(
                     [[2, 0, 0], [2, 1, 0], [2, 1, 1], [2, 0, 1], [2, 0, 0]], ring=True
                 )},
                {"ord": 1, "source": "Output", "label": "solid",
                 "root_type": "POLYHEDRALSURFACE", "type": "POLYGON",
                 "svg": "unused",
                 "has_z": True,
                 "points_xyz": xyz_vertices(
                     [[0, 0, 0], [1, 0, 0], [1, 1, 0], [0, 1, 0], [0, 0, 0]], ring=True
                 )},
            ],
        }
        svg = self.tester.visual_svg("solid-3d", payload)
        self.assertEqual(svg, self.tester.visual_svg("solid-3d", payload))
        self.assertIn(">3D view</text>", svg)
        projected = self.tester.project_3d_payload(payload)
        self.assertEqual(
            [True, True],
            [part["fill_3d_face"] for part in projected["parts"]],
        )
        depths = [float(value) for value in re.findall(r'data-postgis-depth="([^"]+)"', svg)]
        self.assertEqual(sorted(depths), depths)
        fills = re.findall(
            r'data-postgis-face="[12]"[^>]+fill="(#[0-9a-f]{6})"', svg
        )
        self.assertEqual(2, len(fills))
        self.assertEqual(2, len(set(fills)))
        face_colors = re.findall(
            r'data-postgis-face="[12]"[^>]+stroke="(#[0-9a-f]{6})" fill="(#[0-9a-f]{6})"',
            svg,
        )
        self.assertEqual(2, len(face_colors))
        self.assertTrue(all(stroke != fill for stroke, fill in face_colors))
        self.assertEqual(2, svg.count('fill-opacity="1"'))
        self.assertNotIn('fill-opacity="0"', svg)
        self.assertNotIn('class="hidden-edge"', svg)
        self.assertEqual(2, svg.count('class="visible-edge"'))
        self.assertNotIn('fill-opacity="0.82"', svg)
        self.assertNotIn('stroke="#dce2e7" stroke-width="1"', svg)

    def test_wide_3d_scene_keeps_source_x_axis_readable(self):
        points = [
            (x, y, z)
            for x in (0, 280)
            for y in (0, 70)
            for z in (0, 10)
        ]
        view = self.tester.choose_3d_view(points)
        self.assertEqual((0.2, -1.5, 1.15), view)
        origin = self.tester.project_3d_point((0, 0, 0), view)
        x_axis = self.tester.project_3d_point((100, 0, 0), view)
        dx = x_axis[0] - origin[0]
        dy = x_axis[1] - origin[1]
        self.assertGreater(dx, 0)
        self.assertLess(abs(dy / dx), 0.1)

    def test_z_surface_draws_split_hidden_edges_between_fills_and_visible_edges(self):
        payload = {
            "bounds": [0, 0, 2, 1],
            "parts": [
                {"ord": 1, "source": "Output", "label": "solid",
                 "root_type": "POLYHEDRALSURFACE", "type": "POLYGON",
                 "svg": "unused", "has_z": True,
                 "points_xyz": xyz_vertices(points, ring=True)}
                for points in (
                    [[2, 0, 0], [2, 1, 0], [2, 1, 1], [2, 0, 1], [2, 0, 0]],
                    [[0, 0, 0], [1, 0, 0], [1, 1, 0], [0, 1, 0], [0, 0, 0]],
                )
            ],
        }
        self.tester.split_3d_edge_segments = lambda database, parts, face_patches_by_part=None: {
            0: [{"visible": False, "svg": "M 0 0 L 1 0"}],
            1: [{"visible": True, "svg": "M 0 0 L 0 1"}],
        }
        self.tester.split_3d_face_patches = lambda database, parts: {
            0: [{"depth": 0, "svg": "M 0 0 L 1 0 L 1 -1 Z"}],
            1: [{"depth": 0, "svg": "M 1 0 L 2 0 L 2 -1 Z"}],
        }
        svg = self.tester.visual_svg("hidden-edges", payload, database="manual")
        self.assertEqual(1, svg.count('class="hidden-edge"'))
        self.assertNotIn('class="hidden-edge-halo"', svg)
        self.assertEqual(1, svg.count('class="visible-edge"'))
        self.assertIn('stroke-dasharray=', svg)
        self.assertIn('stroke-dashoffset=', svg)
        self.assertLess(svg.index('fill-opacity="1"'), svg.index('class="hidden-edge"'))
        self.assertLess(svg.index('class="hidden-edge"'), svg.index('class="visible-edge"'))

    def test_separate_3d_layers_share_scene_depth_and_edge_phases(self):
        svg = self.tester.visual_svg("layered-solid", {
            "bounds": [0, 0, 2, 2],
            "parts": [
                {
                    "ord": 1, "source": "Output", "label": "front",
                    "root_type": "POLYGON", "type": "POLYGON",
                    "svg": "M 0 0 L 2 0 L 0 -2 Z",
                    "is_3d_face": True, "depth": 10, "shade": 1,
                    "fill_3d_face": True, "face_patches_3d_ready": True,
                    "face_patches_3d": [{
                        "depth": 10, "svg": "M 0 0 L 2 0 L 0 -2 Z",
                    }],
                    "edge_segments_3d_ready": True,
                    "edge_segments_3d": [{
                        "visible": False, "svg": "M 0 0 L 2 0",
                    }],
                },
                {
                    "ord": 2, "source": "Output", "label": "back",
                    "root_type": "POLYGON", "type": "POLYGON",
                    "svg": "M 0 -1 L 2 -1 L 0 -2 Z",
                    "is_3d_face": True, "depth": -10, "shade": 1,
                    "fill_3d_face": True, "face_patches_3d_ready": True,
                    "face_patches_3d": [{
                        "depth": -10, "svg": "M 0 -1 L 2 -1 L 0 -2 Z",
                    }],
                    "edge_segments_3d_ready": True,
                    "edge_segments_3d": [{
                        "visible": True, "svg": "M 0 -1 L 2 -1",
                    }],
                },
            ],
        })

        depths = [float(value) for value in re.findall(
            r'data-postgis-depth="([^"]+)"', svg
        )]
        self.assertEqual([-10, 10], depths)
        self.assertLess(svg.index('data-postgis-depth="10"'), svg.index('class="hidden-edge"'))
        self.assertLess(svg.index('class="hidden-edge"'), svg.index('class="visible-edge"'))
        self.assertEqual(2, svg.count('class="legend-row"'))

    def test_3d_output_scene_paints_after_code_scene(self):
        parts = []
        for ordinal, source, depth in (
            (1, "Code", 10),
            (2, "Output", -10),
        ):
            parts.append({
                "ord": ordinal, "source": source, "label": source,
                "root_type": "POLYGON", "type": "POLYGON",
                "svg": "M 0 0 L 1 0 L 0 -1 Z",
                "is_3d_face": True, "depth": depth, "shade": 1,
                "fill_3d_face": True, "face_patches_3d_ready": True,
                "face_patches_3d": [{
                    "depth": depth, "svg": "M 0 0 L 1 0 L 0 -1 Z",
                }],
                "edge_segments_3d_ready": True,
                "edge_segments_3d": ([{
                    "visible": True, "svg": "M 0 0 L 1 0",
                }] if source == "Code" else []),
            })

        svg = self.tester.visual_svg(
            "source-overlay", {"bounds": [0, 0, 1, 1], "parts": parts}
        )

        drawing = svg.split('<g class="legend-row"', 1)[0]
        code_groups = [match.start() for match in re.finditer(
            'data-postgis-geometry-id="source-overlay-code-1"', drawing
        )]
        output_groups = [match.start() for match in re.finditer(
            'data-postgis-geometry-id="source-overlay-output-1"', drawing
        )]
        self.assertEqual(2, len(code_groups))
        self.assertEqual(1, len(output_groups))
        self.assertLess(max(code_groups), min(output_groups))

    def test_z_surface_draws_visible_face_patches_instead_of_whole_faces(self):
        payload = {
            "bounds": [0, 0, 1, 1],
            "parts": [
                {"ord": 1, "source": "Output", "label": "solid",
                 "root_type": "POLYHEDRALSURFACE", "type": "POLYGON",
                 "svg": "unused", "has_z": True,
                 "points_xyz": xyz_vertices(
                     [[0, 0, 0], [1, 0, 0], [1, 1, 0], [0, 0, 0]], ring=True
                 )},
            ],
        }
        self.tester.split_3d_face_patches = lambda database, parts: {
            0: [
                {"depth": 0, "svg": "M 0 0 L 0.5 0 L 0 -0.5 Z"},
                {"depth": 0, "svg": "M 0.5 0 L 1 0 L 1 -0.5 Z"},
            ],
        }
        self.tester.split_3d_edge_segments = lambda database, parts, face_patches_by_part=None: {}

        svg = self.tester.visual_svg("face-patches", payload, database="manual")

        self.assertEqual(2, svg.count('data-postgis-face='))
        self.assertIn("M 0 0 L 0.5 0 L 0 -0.5 Z", svg)
        self.assertIn("M 0.5 0 L 1 0 L 1 -0.5 Z", svg)
        self.assertNotIn('d="unused"', svg)
        self.assertNotIn('class="visible-edge"', svg)
        self.assertNotIn('class="face-underpaint"', svg)
        face_colors = re.findall(
            r'data-postgis-face="1"[^>]+stroke="(#[0-9a-f]{6})" fill="(#[0-9a-f]{6})"[^>]+stroke-opacity="([^"]+)"',
            svg,
        )
        self.assertEqual(2, len(face_colors))
        self.assertTrue(all(stroke != fill for stroke, fill, _opacity in face_colors))
        self.assertTrue(all(opacity == "0" for _stroke, _fill, opacity in face_colors))

    def test_z_surface_does_not_fall_back_to_fully_occluded_face(self):
        payload = {
            "bounds": [0, 0, 1, 1],
            "parts": [
                {"ord": 1, "source": "Output", "label": "solid",
                 "root_type": "POLYHEDRALSURFACE", "type": "POLYGON",
                 "svg": "unused-front", "has_z": True,
                 "points_xyz": xyz_vertices(
                     [[0, 0, 1], [1, 0, 1], [1, 1, 1], [0, 0, 1]], ring=True
                 )},
                {"ord": 1, "source": "Output", "label": "solid",
                 "root_type": "POLYHEDRALSURFACE", "type": "POLYGON",
                 "svg": "unused-back", "has_z": True,
                 "points_xyz": xyz_vertices(
                     [[0, 0, 0], [1, 0, 0], [1, 1, 0], [0, 0, 0]], ring=True
                 )},
            ],
        }
        self.tester.split_3d_face_patches = lambda database, parts: {
            0: [{"depth": 1, "svg": "M 0 0 L 1 0 L 1 -1 Z"}],
        }
        self.tester.split_3d_edge_segments = (
            lambda database, parts, face_patches_by_part=None: {}
        )

        svg = self.tester.visual_svg("occluded-face", payload, database="manual")

        self.assertEqual(1, svg.count('class="area"'))
        self.assertIn('d="M 0 0 L 1 0 L 1 -1 Z"', svg)
        self.assertNotIn("unused-front", svg)
        self.assertNotIn("unused-back", svg)

    def test_z_surface_face_patches_split_on_depth_order_changes(self):
        queries = []
        self.tester.run_psql_scalar = lambda database, query: queries.append(query) or "[]"

        self.tester.split_3d_face_patches("manual", [
            {
                "is_3d_face": True,
                "fill_3d_face": False,
                "occlude_3d_face": True,
                "frame": "Overlay",
                "face_3d": {
                    "rings": [[[0, 0], [1, 0], [0, 1], [0, 0]]],
                    "depth_plane": (1, 0, 0),
                },
            },
            {
                "is_3d_face": True,
                "fill_3d_face": True,
                "frame": "Overlay",
                "face_3d": {
                    "rings": [[[0, 0], [1, 0], [0, 1], [0, 0]]],
                    "depth_plane": (-1, 0, 1),
                },
            },
        ])

        encoded = re.search(r"decode\('([^']+)', 'base64'\)", queries[0]).group(1)
        face_payload = json.loads(base64.b64decode(encoded))["faces"]
        self.assertEqual(2, len(face_payload))
        self.assertEqual(2, len({face["fill_group"] for face in face_payload}))
        self.assertEqual([False, True], [face["draw_fill"] for face in face_payload])
        self.assertIn("depth_split_lines AS", queries[0])
        self.assertIn("raw_faces AS", queries[0])
        self.assertIn("faces AS MATERIALIZED", queries[0])
        self.assertIn("ST_CollectionExtract(ST_MakeValid(f.geom), 3)", queries[0])
        self.assertIn("ST_Intersection(f.geom, o.geom, 1e-9)", queries[0])
        self.assertIn("ST_Polygonize", queries[0])
        self.assertIn("ST_UnaryUnion(ST_Collect(geom), 1e-9)", queries[0])
        self.assertIn("ST_Difference(p.geom, h.geom, 1e-9)", queries[0])
        self.assertIn("patch_precision AS MATERIALIZED", queries[0])
        self.assertIn("scene_span * scene_span * 1e-12 AS min_patch_area", queries[0])
        self.assertIn("ST_Area((dumped).geom) > p.min_patch_area", queries[0])
        self.assertIn("visible_face_unions AS", queries[0])
        self.assertIn("visible_fill_unions AS", queries[0])
        self.assertIn("rendered_and_edge_patches AS", queries[0])
        self.assertIn(
            "GROUP BY frame_id, group_id, fill_group",
            queries[0],
        )
        self.assertIn("'edge_visibility', edge_visibility", queries[0])
        self.assertNotIn("ST_ReducePrecision", queries[0])

    def test_coplanar_fill_groups_share_3d_render_style(self):
        queries = []
        self.tester.run_psql_scalar = lambda database, query: queries.append(query) or "[]"
        parts = []
        for root_type, visibility_group, ordinal in (
            ("POLYHEDRALSURFACE", "polyhedral", 1),
            ("POLYHEDRALSURFACE", "polyhedral", 1),
            ("TIN", "tin", 2),
            ("TIN", "tin", 2),
            ("POLYGON", "dump", 3),
            ("POLYGON", "dump", 4),
        ):
            parts.append({
                "is_3d_face": True,
                "fill_3d_face": True,
                "root_type": root_type,
                "frame": "Overlay",
                "visibility_group": visibility_group,
                "ord": ordinal,
                "face_3d": {
                    "rings": [[[0, 0], [1, 0], [0, 1], [0, 0]]],
                    "depth_plane": (0, 0, 1),
                },
            })

        self.tester.split_3d_face_patches("manual", parts)

        encoded = re.search(r"decode\('([^']+)', 'base64'\)", queries[0]).group(1)
        faces = json.loads(base64.b64decode(encoded))["faces"]
        self.assertEqual(faces[0]["fill_group"], faces[1]["fill_group"])
        self.assertEqual(faces[2]["fill_group"], faces[3]["fill_group"])
        self.assertNotEqual(faces[4]["fill_group"], faces[5]["fill_group"])

    def test_3d_output_layers_share_visibility_but_keep_paint_identity(self):
        payload = {
            "bounds": [0, 0, 1, 1],
            "parts": [
                {
                    "ord": ordinal, "source": source, "label": label,
                    "root_type": "POLYGON", "type": "POLYGON", "svg": "unused",
                    "has_z": True,
                    "points_xyz": xyz_vertices(points, ring=True),
                }
                for ordinal, source, label, points in (
                    (1, "Output", "face 1",
                     [[0, 0, 0], [1, 0, 0], [1, 1, 0], [0, 0, 0]]),
                    (2, "Output", "face 2",
                     [[0, 0, 0], [1, 1, 0], [0, 1, 0], [0, 0, 0]]),
                    (3, "Code", "input",
                     [[0, 0, 1], [1, 0, 1], [0, 1, 1], [0, 0, 1]]),
                )
            ],
        }

        output_1, output_2, code = self.tester.project_3d_payload(payload)["parts"]

        self.assertEqual(output_1["visibility_group"], output_2["visibility_group"])
        self.assertNotEqual(output_1["paint_group"], output_2["paint_group"])
        self.assertNotEqual(output_1["visibility_group"], code["visibility_group"])
        shared_edge = self.tester.edge_key_3d([0, 0, 0], [1, 1, 0])
        self.assertEqual(2, sum(
            edge["key"] == shared_edge
            for part in (output_1, output_2)
            for edge in part["edges_3d"]
        ))

    def test_2d_surface_does_not_draw_hidden_edges(self):
        svg = self.tester.visual_svg(
            "flat-area",
            {"bounds": [0, 0, 1, 1], "parts": [{
                "ord": 1, "source": "Output", "label": "area",
                "root_type": "POLYGON", "type": "POLYGON",
                "svg": "M 0 0 L 1 0 L 0 -1 Z", "has_z": False,
                "vertices": [[0, 0], [1, 0], [0, -1], [0, 0]],
            }]},
        )
        self.assertNotIn('class="hidden-edge"', svg)

    def test_tin_faces_keep_fills_when_their_projections_overlap(self):
        payload = {
            "bounds": [0, 0, 1, 1],
            "parts": [
                {"ord": 1, "source": "Code", "label": "TIN",
                 "root_type": "TIN", "type": "TRIANGLE", "svg": "unused",
                 "has_z": True, "points_xyz": xyz_vertices(points, ring=True)}
                for points in (
                    [[0, 0, 0], [0, 0, 1], [0, 1, 0], [0, 0, 0]],
                    [[0, 0, 0], [0, 1, 0], [1, 0, 0], [0, 0, 0]],
                    [[0, 0, 0], [1, 0, 0], [0, 0, 1], [0, 0, 0]],
                    [[1, 0, 0], [0, 1, 0], [0, 0, 1], [1, 0, 0]],
                )
            ],
        }
        svg = self.tester.visual_svg("tetrahedron", payload)
        self.assertEqual(4, svg.count('data-postgis-face='))
        self.assertEqual(4, svg.count('fill-opacity="1"'))
        self.assertNotIn('fill-opacity="0"', svg)
        self.assertNotIn('class="hidden-edge"', svg)
        self.assertEqual(4, svg.count('class="visible-edge"'))
        self.assertNotIn('fill-opacity="0.82"', svg)

    def test_output_tin_over_code_solid_renders_as_wire_overlay(self):
        payload = {
            "bounds": [0, 0, 1, 1],
            "parts": [
                {"ord": 1, "source": "Code", "label": "cube",
                 "root_type": "POLYHEDRALSURFACE", "type": "POLYGON",
                 "svg": "unused", "has_z": True,
                 "points_xyz": xyz_vertices(
                     [[0, 0, 0], [1, 0, 0], [0, 1, 0], [0, 0, 0]], ring=True
                 )}
            ] + [
                {"ord": 2, "source": "Output", "label": "tessellation",
                 "root_type": "TIN", "type": "TRIANGLE",
                 "svg": "unused", "has_z": True,
                 "points_xyz": xyz_vertices(points, ring=True)}
                for points in (
                    [[0, 0, 0], [1, 0, 0], [1, 1, 0], [0, 0, 0]],
                    [[0, 0, 0], [1, 1, 0], [0, 1, 0], [0, 0, 0]],
                )
            ],
        }
        ordinary = self.tester.project_3d_payload(payload)
        self.assertEqual([2, 2], [len(part["edges_3d"]) for part in ordinary["parts"][1:]])
        self.assertTrue(
            all(not part.get("tessellation_overlay") for part in ordinary["parts"][1:])
        )

        projected = self.tester.project_3d_payload(payload, show_tessellation=True)
        code_part, *output_parts = projected["parts"]

        self.assertTrue(code_part["fill_3d_face"])
        self.assertTrue(all(not part["fill_3d_face"] for part in output_parts))
        self.assertTrue(all(part["occlude_3d_face"] for part in output_parts))
        self.assertTrue(all(part["tessellation_overlay"] for part in output_parts))
        self.assertTrue(all(not part["draw_hidden_3d_edges"] for part in output_parts))
        seam_key = self.tester.edge_key_3d([0, 0, 0], [1, 1, 0])
        self.assertEqual(
            [[seam_key], [seam_key]],
            [[edge["key"] for edge in part["edges_3d"]] for part in output_parts],
        )

    def test_3d_wire_overlay_does_not_draw_hidden_segments_or_empty_faces(self):
        svg = self.tester.visual_svg("wire-overlay", {
            "bounds": [0, 0, 1, 1],
            "parts": [{
                "ord": 1, "source": "Output", "label": "tessellation",
                "root_type": "TIN", "type": "TRIANGLE",
                "svg": "M 0 0 L 1 0 L 1 -1 Z", "has_z": False,
                "is_3d_face": True, "depth": 0, "shade": 1,
                "fill_3d_face": False, "draw_hidden_3d_edges": False,
                "face_patches_3d_ready": True,
                "face_patches_3d": [{"svg": "M 0 0 L 1 0 L 1 -1 Z"}],
                "edge_segments_3d_ready": True,
                "edge_segments_3d": [
                    {"visible": True, "svg": "M 0 0 L 1 -1"},
                    {"visible": False, "svg": "M 0 0 L 1 0"},
                ],
            }],
        })

        self.assertEqual(1, svg.count('class="visible-edge"'))
        self.assertNotIn('class="hidden-edge"', svg)
        self.assertNotIn('class="area"', svg)

    def test_separate_frame_3d_output_keeps_face_fill(self):
        payload = {
            "frames": [
                {"id": "Code", "bounds": [0, 0, 5, 5]},
                {"id": "Output", "bounds": [0, 0, 5, 5]},
            ],
            "parts": [
                {"ord": 1, "source": "Code", "label": "Code", "frame": "Code",
                 "root_type": "POLYGON", "type": "POLYGON", "svg": "unused",
                 "has_z": True,
                 "points_xyz": xyz_vertices(
                     [[0, 0, 2], [5, 0, 2], [0, 5, 2], [0, 0, 2]],
                     ring=True,
                 )},
                {"ord": 2, "source": "Output", "label": "Output", "frame": "Output",
                 "root_type": "POLYGON", "type": "POLYGON", "svg": "unused",
                 "has_z": True,
                 "points_xyz": xyz_vertices(
                     [[0, 0, 2], [0, 5, 2], [5, 0, 2], [0, 0, 2]],
                     ring=True,
                 )},
            ],
        }

        projected = self.tester.project_3d_payload(payload)
        self.assertEqual(
            [True, True],
            [part["fill_3d_face"] for part in projected["parts"]],
        )

        svg = self.tester.visual_svg("separate-fill", payload)
        self.assertEqual(2, svg.count('fill-opacity="1"'))
        self.assertNotIn('fill-opacity="0"', svg)
        self.assertIn('data-postgis-frame="Code"', svg)
        self.assertIn('data-postgis-frame="Output"', svg)

    def test_separate_frame_2d_input_and_3d_output_keeps_output_fill(self):
        payload = {
            "frames": [
                {"id": "Code", "bounds": [0, 0, 5, 5]},
                {"id": "Output", "bounds": [0, 0, 5, 5]},
            ],
            "parts": [
                {"ord": 1, "source": "Code", "label": "Code", "frame": "Code",
                 "root_type": "POLYGON", "type": "POLYGON",
                 "svg": "M 0 0 L 0 -5 L 5 0 Z M 1 -1 L 3 -1 L 1 -3 Z",
                 "has_z": False,
                 "vertices": [[0, 0], [0, -5], [5, 0], [0, 0]]},
                {"ord": 2, "source": "Output", "label": "Output", "frame": "Output",
                 "root_type": "POLYGON", "type": "POLYGON", "svg": "unused",
                 "has_z": True,
                 "points_xyz": xyz_vertices(
                     [[0, 0, 0], [0, 5, 0], [5, 0, 0], [0, 0, 0]],
                     ring=True,
                 )},
            ],
        }

        projected = self.tester.project_3d_payload(payload)
        self.assertNotIn("fill_3d_face", projected["parts"][0])
        self.assertTrue(projected["parts"][1]["fill_3d_face"])

        svg = self.tester.visual_svg("separate-2d-3d-fill", payload)
        self.assertEqual(1, svg.count('fill-opacity="0.24"'))
        self.assertEqual(1, svg.count('fill-opacity="1"'))
        self.assertNotIn('fill-opacity="0"', svg)

    def test_mixed_frame_projects_2d_footprint_with_3d_roof(self):
        payload = {
            "frames": [{"id": "Overlay", "bounds": [0, 0, 4, 3]}],
            "parts": [
                {"ord": 1, "source": "Output", "label": "footprint",
                 "root_type": "POLYGON", "type": "POLYGON",
                 "svg": "M 0 0 L 4 0 4 -3 0 -3 Z", "has_z": False,
                 "points_xyz": xyz_vertices(
                     [[0, 0, 0], [4, 0, 0], [4, 3, 0], [0, 3, 0], [0, 0, 0]], ring=True
                 )},
                {"ord": 2, "source": "Output", "label": "roof",
                 "root_type": "POLYHEDRALSURFACE", "type": "POLYGON",
                 "svg": "unused", "has_z": True,
                 "points_xyz": xyz_vertices(
                     [[0, 0, 0], [4, 0, 0], [2, 1.5, 2], [0, 0, 0]], ring=True
                 )},
            ],
        }
        projected = self.tester.project_3d_payload(payload)
        self.assertTrue(projected["frames"][0]["three_dimensional"])
        self.assertTrue(all(part.get("is_3d_face") for part in projected["parts"]))
        self.assertNotEqual(payload["parts"][0]["svg"], projected["parts"][0]["svg"])
        svg = self.tester.visual_svg("mixed-roof", payload)
        self.assertEqual(2, svg.count("data-postgis-face="))
        self.assertNotIn('stroke="#dce2e7" stroke-width="1"', svg)

    def test_z_surface_draws_coplanar_fill_union_not_edge_visibility_patches(self):
        payload = {
            "bounds": [0, 0, 1, 1],
            "parts": [{
                "ord": 1, "source": "Output", "label": "solid",
                "root_type": "POLYHEDRALSURFACE", "type": "POLYGON",
                "svg": "unused", "has_z": True,
                "points_xyz": xyz_vertices(
                    [[0, 0, 0], [1, 0, 0], [1, 1, 0], [0, 0, 0]], ring=True
                ),
            }],
        }
        self.tester.split_3d_face_patches = lambda database, parts: {
            0: [
                {
                    "depth": 0,
                    "svg": "M 0 0 L 0.5 0 L 0 -0.5 Z",
                    "wkt": "POLYGON((0 0,0.5 0,0 0.5,0 0))",
                    "draw": False,
                    "edge_visibility": True,
                },
                {
                    "depth": 0,
                    "svg": "M 0 0 L 1 0 L 1 -1 Z",
                    "wkt": "POLYGON((0 0,1 0,1 1,0 0))",
                    "draw": True,
                    "edge_visibility": False,
                },
            ],
        }
        self.tester.split_3d_edge_segments = (
            lambda database, parts, face_patches_by_part=None: {}
        )

        svg = self.tester.visual_svg("coplanar-fill", payload, database="manual")

        self.assertEqual(1, svg.count('class="area"'))
        self.assertIn('d="M 0 0 L 1 0 L 1 -1 Z"', svg)
        self.assertNotIn('d="M 0 0 L 0.5 0 L 0 -0.5 Z"', svg)

    def test_coplanar_3d_triangle_seams_are_not_drawn_as_edges(self):
        payload = {
            "bounds": [0, 0, 1, 1],
            "parts": [
                {"ord": 1, "source": "Output", "label": "surface",
                 "root_type": "TIN", "type": "TRIANGLE", "svg": "unused",
                 "has_z": True, "points_xyz": xyz_vertices(points, ring=True)}
                for points in (
                    [[0, 0, 0], [1, 0, 0], [1, 1, 0], [0, 0, 0]],
                    [[0, 0, 0], [1, 1, 0], [0, 1, 0], [0, 0, 0]],
                )
            ],
        }
        projected = self.tester.project_3d_payload(payload)
        edge_keys = [
            edge["key"]
            for part in projected["parts"]
            for edge in part.get("edges_3d") or []
        ]
        self.assertEqual(4, len(edge_keys))
        self.assertNotIn(self.tester.edge_key_3d([0, 0, 0], [1, 1, 0]), edge_keys)

    def test_coplanar_3d_triangle_faces_share_the_geometry_palette(self):
        payload = {
            "bounds": [0, 0, 1, 1],
            "parts": [
                {"ord": 1, "source": "Output", "label": "surface",
                 "root_type": "TIN", "type": "TRIANGLE", "svg": "unused",
                 "has_z": True, "points_xyz": xyz_vertices(points, ring=True)}
                for points in (
                    [[0, 0, 0], [1, 0, 0], [1, 1, 0], [0, 0, 0]],
                    [[0, 0, 0], [1, 1, 0], [0, 1, 0], [0, 0, 0]],
                )
            ],
        }

        svg = self.tester.visual_svg("coplanar-tin", payload)

        fills = re.findall(r'data-postgis-face="[12]"[^>]+fill="(#[0-9a-f]{6})"', svg)
        self.assertEqual(2, len(fills))
        self.assertEqual(1, len(set(fills)))

    def test_non_coplanar_3d_triangle_shared_edges_are_kept(self):
        payload = {
            "bounds": [0, 0, 1, 1],
            "parts": [
                {"ord": 1, "source": "Output", "label": "surface",
                 "root_type": "TIN", "type": "TRIANGLE", "svg": "unused",
                 "has_z": True, "points_xyz": xyz_vertices(points, ring=True)}
                for points in (
                    [[0, 0, 0], [1, 0, 0], [1, 1, 0], [0, 0, 0]],
                    [[0, 0, 0], [1, 1, 0], [0, 1, 1], [0, 0, 0]],
                )
            ],
        }
        projected = self.tester.project_3d_payload(payload)
        edge_keys = [
            edge["key"]
            for part in projected["parts"]
            for edge in part.get("edges_3d") or []
        ]
        self.assertIn(self.tester.edge_key_3d([0, 0, 0], [1, 1, 0]), edge_keys)

    def test_dense_capped_surface_marks_side_facets_as_silhouette_only(self):
        segment_count = 32
        ring = [
            [
                round(10 * math.cos(2 * math.pi * index / segment_count), 12),
                round(10 * math.sin(2 * math.pi * index / segment_count), 12),
            ]
            for index in range(segment_count)
        ]
        faces = []
        for index in range(segment_count):
            start = ring[index]
            end = ring[(index + 1) % segment_count]
            faces.append([
                [start[0], start[1], 0],
                [end[0], end[1], 0],
                [end[0], end[1], 100],
                [start[0], start[1], 100],
                [start[0], start[1], 0],
            ])
        sharp_edge = [[20, 0, 0], [20, 0, 0.1]]
        faces.extend([
            [
                [20, 0, 0], [21, 0, 0], [21, 0, 0.1],
                [20, 0, 0.1], [20, 0, 0],
            ],
            [
                [20, 0, 0], [20, 0, 0.1], [20, 1, 0.1],
                [20, 1, 0], [20, 0, 0],
            ],
        ])
        faces.extend([
            [[x, y, 0] for x, y in ring] + [[ring[0][0], ring[0][1], 0]],
            [[x, y, 100] for x, y in reversed(ring)]
            + [[ring[-1][0], ring[-1][1], 100]],
        ])
        payload = {
            "bounds": [-10, -10, 10, 10],
            "parts": [
                {
                    "ord": 1,
                    "source": "Output",
                    "label": "tube",
                    "root_type": "POLYHEDRALSURFACE",
                    "type": "POLYGON",
                    "svg": "unused",
                    "has_z": True,
                    "points_xyz": xyz_vertices(points, ring=True),
                }
                for points in faces
            ],
        }

        projected = self.tester.project_3d_payload(payload)
        self.assertEqual(
            [True] * (segment_count + 2) + [False, False],
            [part["smooth_3d_side"] for part in projected["parts"]],
        )
        self.assertEqual(
            [False] * (segment_count + 2) + [True, True],
            [part["smooth_3d_cap"] for part in projected["parts"]],
        )
        flags_by_key = {}
        for part in projected["parts"]:
            self.assertTrue(part["smooth_3d_surface"])
            for edge in part["edges_3d"]:
                flags_by_key.setdefault(edge["key"], set()).add(edge["silhouette_only"])

        side_keys = {
            self.tester.edge_key_3d([x, y, 0], [x, y, 100])
            for x, y in ring
        }
        cap_keys = {
            self.tester.edge_key_3d(
                [ring[index][0], ring[index][1], z],
                [
                    ring[(index + 1) % segment_count][0],
                    ring[(index + 1) % segment_count][1],
                    z,
                ],
            )
            for index in range(segment_count)
            for z in (0, 100)
        }
        self.assertEqual({True}, set().union(*(flags_by_key[key] for key in side_keys)))
        self.assertEqual({False}, set().union(*(flags_by_key[key] for key in cap_keys)))
        self.assertEqual(
            {False}, flags_by_key[self.tester.edge_key_3d(*sharp_edge)]
        )

    def test_smooth_face_fill_overlaps_without_white_seams(self):
        payload_part = {
            "ord": 1,
            "source": "Output",
            "label": "shell",
            "root_type": "POLYHEDRALSURFACE",
            "type": "POLYGON",
            "svg": "unused",
            "has_z": True,
            "is_3d_face": True,
            "depth": 0,
            "shade": 0.8,
            "fill_3d_face": True,
            "smooth_3d_surface": True,
            "smooth_3d_side": True,
            "face_patches_3d_ready": True,
            "face_patches_3d": [{
                "svg": "M 0 0 L 1 0 L 0 -1 Z",
                "depth": 0,
                "draw": True,
            }],
            "edge_segments_3d_ready": True,
            "edge_segments_3d": [],
        }
        svg = self.tester.visual_svg("smooth-shell", {
            "bounds": [0, 0, 1, 1],
            "parts": [payload_part],
        })

        face = re.search(r'<path class="area"[^>]+>', svg).group(0)
        fill = re.search(r'fill="(#[0-9a-f]{6})"', face).group(1)
        self.assertEqual(
            self.tester.shade_color("#a62c2b", self.tester.SMOOTH_3D_SIDE_SHADE),
            fill,
        )
        self.assertIn(f'stroke="{fill}"', face)
        self.assertIn('stroke-opacity="1"', face)
        self.assertNotIn('stroke="white"', face)

        payload = {
            "bounds": [0, 0, 1, 1],
            "parts": [{
                **payload_part,
                "shade": 0.2,
                "smooth_3d_side": False,
                "smooth_3d_cap": True,
            }],
        }
        cap_svg = self.tester.visual_svg("smooth-cap", payload)
        cap_face = re.search(r'<path class="area"[^>]+>', cap_svg).group(0)
        cap_fill = re.search(r'fill="(#[0-9a-f]{6})"', cap_face).group(1)
        self.assertEqual(
            self.tester.shade_color(
                "#a62c2b", self.tester.SMOOTH_3D_CAP_MIN_SHADE
            ),
            cap_fill,
        )

    def test_dense_triangular_surface_keeps_only_silhouette_edges(self):
        segment_count = 32
        ring = [
            [
                round(10 * math.cos(2 * math.pi * index / segment_count), 12),
                round(10 * math.sin(2 * math.pi * index / segment_count), 12),
            ]
            for index in range(segment_count)
        ]
        faces = []
        for index in range(segment_count):
            start = ring[index]
            end = ring[(index + 1) % segment_count]
            start_bottom = [start[0], start[1], 0]
            end_bottom = [end[0], end[1], 0]
            start_top = [start[0], start[1], 100]
            end_top = [end[0], end[1], 100]
            faces.extend([
                [start_bottom, end_bottom, end_top, start_bottom],
                [start_bottom, end_top, start_top, start_bottom],
            ])
        payload = {
            "bounds": [-10, -10, 10, 10],
            "parts": [
                {
                    "ord": 1,
                    "source": "Output",
                    "label": "smooth mesh",
                    "root_type": "POLYHEDRALSURFACE",
                    "type": "POLYGON",
                    "svg": "unused",
                    "has_z": True,
                    "points_xyz": xyz_vertices(points, ring=True),
                }
                for points in faces
            ],
        }

        projected = self.tester.project_3d_payload(payload)
        self.assertTrue(all(part["smooth_3d_surface"] for part in projected["parts"]))
        self.assertTrue(all(not part["smooth_3d_side"] for part in projected["parts"]))
        flags_by_key = {}
        for part in projected["parts"]:
            for edge in part["edges_3d"]:
                flags_by_key.setdefault(edge["key"], set()).add(edge["silhouette_only"])

        side_keys = {
            self.tester.edge_key_3d([x, y, 0], [x, y, 100])
            for x, y in ring
        }
        boundary_keys = {
            self.tester.edge_key_3d(
                [ring[index][0], ring[index][1], z],
                [
                    ring[(index + 1) % segment_count][0],
                    ring[(index + 1) % segment_count][1],
                    z,
                ],
            )
            for index in range(segment_count)
            for z in (0, 100)
        }
        self.assertEqual({True}, set().union(*(flags_by_key[key] for key in side_keys)))
        self.assertEqual(
            {False}, set().union(*(flags_by_key[key] for key in boundary_keys))
        )

    def test_faceted_shell_hides_short_bevels_and_shallow_mesh_seams(self):
        vertices = (
            (0.0, 0.0, 0.0),
            (4.0, 0.0, 0.0),
            (0.0, 4.0, 0.0),
            (0.0, 0.0, 4.0),
        )
        fraction = 0.04

        def truncated(start, end):
            return [
                (1 - fraction) * vertices[start][axis]
                + fraction * vertices[end][axis]
                for axis in range(3)
            ]

        faces = []
        for vertex in range(4):
            cap = [truncated(vertex, other) for other in range(4) if other != vertex]
            faces.append(cap + [cap[0]])

        shallow_edge = None
        for first, second, third in ((0, 1, 2), (0, 1, 3), (0, 2, 3), (1, 2, 3)):
            boundary = [
                truncated(first, second),
                truncated(second, first),
                truncated(second, third),
                truncated(third, second),
                truncated(third, first),
                truncated(first, third),
            ]
            first_vector = [
                vertices[second][axis] - vertices[first][axis]
                for axis in range(3)
            ]
            second_vector = [
                vertices[third][axis] - vertices[first][axis]
                for axis in range(3)
            ]
            normal = [
                first_vector[1] * second_vector[2] - first_vector[2] * second_vector[1],
                first_vector[2] * second_vector[0] - first_vector[0] * second_vector[2],
                first_vector[0] * second_vector[1] - first_vector[1] * second_vector[0],
            ]
            normal_length = math.sqrt(sum(value ** 2 for value in normal))
            center = [
                sum(point[axis] for point in boundary) / len(boundary)
                + 0.02 * normal[axis] / normal_length
                for axis in range(3)
            ]
            if shallow_edge is None:
                shallow_edge = self.tester.edge_key_3d(boundary[0], center)
            for index, start in enumerate(boundary):
                end = boundary[(index + 1) % len(boundary)]
                faces.append([start, end, center, start])

        payload = {
            "bounds": [0, 0, 4, 4],
            "parts": [
                {
                    "ord": 1,
                    "source": "Output",
                    "label": "faceted shell",
                    "root_type": "POLYHEDRALSURFACE",
                    "type": "POLYGON",
                    "svg": "unused",
                    "has_z": True,
                    "points_xyz": xyz_vertices(points, ring=True),
                }
                for points in faces
            ],
        }

        projected = self.tester.project_3d_payload(payload)
        flags_by_key = {}
        for part in projected["parts"]:
            self.assertTrue(part["smooth_3d_surface"])
            for edge in part["edges_3d"]:
                flags_by_key.setdefault(edge["key"], set()).add(
                    edge["silhouette_only"]
                )

        short_bevel = self.tester.edge_key_3d(truncated(0, 1), truncated(0, 2))
        long_crease = self.tester.edge_key_3d(truncated(0, 1), truncated(1, 0))
        self.assertEqual({True}, flags_by_key[short_bevel])
        self.assertEqual({True}, flags_by_key[shallow_edge])
        self.assertEqual({False}, flags_by_key[long_crease])

        tetrahedron = {
            "bounds": [0, 0, 4, 4],
            "parts": [
                {
                    "ord": 1,
                    "source": "Code",
                    "label": "tetrahedron",
                    "root_type": "POLYHEDRALSURFACE",
                    "type": "POLYGON",
                    "svg": "unused",
                    "has_z": True,
                    "points_xyz": xyz_vertices(
                        [vertices[index] for index in face] + [vertices[face[0]]],
                        ring=True,
                    ),
                }
                for face in ((0, 1, 2), (0, 1, 3), (0, 2, 3), (1, 2, 3))
            ],
        }
        tetrahedron_parts = self.tester.project_3d_payload(tetrahedron)["parts"]
        self.assertTrue(all(
            not part["smooth_3d_surface"] for part in tetrahedron_parts
        ))
        self.assertTrue(all(
            not edge["silhouette_only"]
            for part in tetrahedron_parts
            for edge in part["edges_3d"]
        ))

    def test_faceted_shell_keeps_short_screen_space_crease_connector(self):
        vertices = [
            (-0.193132, 0.126622, 4.000118),
            (-0.034518, -0.228316, 4.003695),
            (0.087369, 0.087369, 4.195106),
            (0.083342, 0.081596, 3.800677),
            (0.047923, 3.840255, 0.159745),
            (3.866667, 0.133333, 0.133333),
            (0.091495, 3.926447, -0.198877),
            (3.919597, 0.080403, -0.201008),
            (0.037318, -0.006838, -0.227802),
            (3.877427, -0.195387, 0.011527),
            (4.22786, -0.006276, 0.037067),
            (0.012489, -0.196238, 0.121111),
            (-0.202457, 0.079493, 0.077623),
            (0.101221, 4.198699, 0.060053),
            (-0.228261, 4.028106, -0.020978),
        ]
        faces = [
            (0, 1, 2), (3, 2, 1), (2, 3, 0), (4, 5, 6),
            (7, 8, 6), (5, 7, 6), (5, 9, 10), (5, 10, 7),
            (9, 7, 10), (11, 1, 12), (12, 8, 11), (5, 3, 9),
            (6, 13, 4), (13, 14, 4), (6, 14, 13), (5, 4, 3),
            (9, 8, 7), (11, 8, 9), (3, 4, 0), (14, 0, 4),
            (14, 12, 0), (12, 1, 0), (9, 3, 1), (11, 9, 1),
            (6, 8, 14), (14, 8, 12),
        ]
        payload = {
            "frames": [
                {"id": "faceted", "bounds": [0, 0, 4, 4]},
                {"id": "labels", "bounds": [0, 0, 4, 4]},
            ],
            "bounds": [0, 0, 4, 4],
            "parts": [
                {
                    "ord": 1,
                    "frame": "faceted",
                    "source": "Output",
                    "label": "faceted shell",
                    "root_type": "POLYHEDRALSURFACE",
                    "type": "POLYGON",
                    "svg": "unused",
                    "has_z": True,
                    "points_xyz": xyz_vertices(
                        [vertices[index] for index in face]
                        + [vertices[face[0]]],
                        ring=True,
                    ),
                }
                for face in faces
            ] + [{
                "ord": 2,
                "frame": "labels",
                "source": "Output",
                "label": "2D control",
                "root_type": "POINT",
                "type": "POINT",
                "svg": "unused",
                "has_z": False,
                "x": 0,
                "y": 0,
                "points": [],
                "vertices": [],
            }],
        }

        projected = self.tester.project_3d_payload(payload)
        flags_by_key = {}
        for part in projected["parts"]:
            if not part.get("is_3d_face"):
                continue
            self.assertTrue(part["smooth_3d_surface"])
            for edge in part["edges_3d"]:
                flags_by_key.setdefault(edge["key"], set()).add(
                    edge["silhouette_only"]
                )

        connector = self.tester.edge_key_3d(vertices[0], vertices[3])
        long_crease = self.tester.edge_key_3d(vertices[3], vertices[5])
        rear_left_crease = self.tester.edge_key_3d(vertices[8], vertices[7])
        rear_right_crease = self.tester.edge_key_3d(vertices[12], vertices[14])
        short_bevel = self.tester.edge_key_3d(vertices[1], vertices[2])
        self.assertEqual({False}, flags_by_key[connector])
        self.assertEqual({False}, flags_by_key[long_crease])
        self.assertEqual({False}, flags_by_key[rear_left_crease])
        self.assertEqual({False}, flags_by_key[rear_right_crease])
        self.assertEqual({True}, flags_by_key[short_bevel])

    def test_faceted_skinny_strip_uses_broader_smoothing_angle(self):
        skinny_strip_faces = [
            {"edge_to_other_length_ratio": 0.91},
            {"edge_to_other_length_ratio": 0.88},
        ]
        supported_crease_faces = [
            {"edge_to_other_length_ratio": 0.91},
            {"edge_to_other_length_ratio": 0.45},
        ]

        self.assertEqual(
            self.tester.SMOOTH_3D_EDGE_MAX_ANGLE,
            self.tester.faceted_smooth_edge_max_angle(skinny_strip_faces),
        )
        self.assertEqual(
            self.tester.FACETED_3D_SMOOTH_EDGE_MAX_ANGLE,
            self.tester.faceted_smooth_edge_max_angle(supported_crease_faces),
        )

    def test_dense_smooth_surface_keeps_an_individual_sharp_crease(self):
        heights = [round(0.01 * index * index, 12) for index in range(32)]
        faces = [
            [
                [index, 0, heights[index]],
                [index + 1, 0, heights[index + 1]],
                [index + 1, 1, heights[index + 1]],
                [index, 1, heights[index]],
                [index, 0, heights[index]],
            ]
            for index in range(31)
        ]
        faces.append([
            [31, 0, heights[31]],
            [32, 0, 100],
            [32, 1, 100],
            [31, 1, heights[31]],
            [31, 0, heights[31]],
        ])
        payload = {
            "bounds": [0, 0, 32, 1],
            "parts": [
                {
                    "ord": 1,
                    "source": "Output",
                    "label": "mostly smooth mesh",
                    "root_type": "POLYHEDRALSURFACE",
                    "type": "POLYGON",
                    "svg": "unused",
                    "has_z": True,
                    "points_xyz": xyz_vertices(points, ring=True),
                }
                for points in faces
            ],
        }

        projected = self.tester.project_3d_payload(payload)
        flags_by_key = {}
        for part in projected["parts"]:
            self.assertTrue(part["smooth_3d_surface"])
            for edge in part["edges_3d"]:
                flags_by_key.setdefault(edge["key"], set()).add(edge["silhouette_only"])

        shallow_key = self.tester.edge_key_3d(
            [15, 0, heights[15]], [15, 1, heights[15]]
        )
        sharp_key = self.tester.edge_key_3d(
            [31, 0, heights[31]], [31, 1, heights[31]]
        )
        self.assertEqual({True}, flags_by_key[shallow_key])
        self.assertEqual({False}, flags_by_key[sharp_key])

    def test_mixed_frame_projects_3d_line_with_3d_polygon(self):
        line_points = [[0, 0, 0], [4, 0, 1], [4, 3, 2]]
        payload = {
            "frames": [{"id": "Overlay", "bounds": [0, 0, 4, 3]}],
            "parts": [
                {"ord": 1, "source": "Code", "label": "ring",
                 "root_type": "LINESTRING", "type": "LINESTRING",
                 "svg": "M 0 0 L 4 0 4 -3", "closed": False, "has_z": True,
                 "points_xyz": xyz_vertices(line_points)},
                {"ord": 2, "source": "Output", "label": "polygon",
                 "root_type": "POLYGON", "type": "POLYGON",
                 "svg": "unused", "has_z": True,
                 "points_xyz": xyz_vertices(line_points + [line_points[0]], ring=True)},
            ],
        }
        projected = self.tester.project_3d_payload(payload)
        line = projected["parts"][0]
        self.assertNotEqual(payload["parts"][0]["svg"], line["svg"])
        self.assertFalse(line.get("is_3d_face", False))
        min_x, min_y, max_x, max_y = projected["frames"][0]["bounds"]
        for x, svg_y in line["vertices"]:
            self.assertGreaterEqual(x, min_x)
            self.assertLessEqual(x, max_x)
            self.assertGreaterEqual(-svg_y, min_y)
            self.assertLessEqual(-svg_y, max_y)
        svg = self.tester.visual_svg("makepolygon", payload)
        self.assertIn('data-postgis-geometry-id="makepolygon-code-1"', svg)
        self.assertIn('data-postgis-geometry-id="makepolygon-output-1"', svg)

    def test_separate_frames_share_3d_projection_policy(self):
        ring_points = [[0, 0, 1], [4, 0, 2], [4, 3, 2], [0, 3, 1], [0, 0, 1]]
        payload = {
            "frames": [
                {"id": "Code", "bounds": [0, 0, 4, 3]},
                {"id": "Output", "bounds": [0, 0, 4, 3]},
            ],
            "parts": [
                {"ord": 1, "source": "Code", "label": "ring", "frame": "Code",
                 "root_type": "LINESTRING", "type": "LINESTRING",
                 "svg": "M 0 0 L 4 0 4 -3 0 -3 0 0", "closed": True,
                 "has_z": True, "points_xyz": xyz_vertices(ring_points)},
                {"ord": 2, "source": "Output", "label": "polygon", "frame": "Output",
                 "root_type": "POLYGON", "type": "POLYGON", "svg": "unused",
                 "has_z": True, "points_xyz": xyz_vertices(ring_points, ring=True)},
            ],
        }
        projected = self.tester.project_3d_payload(payload)
        self.assertTrue(all(frame["three_dimensional"] for frame in projected["frames"]))
        self.assertNotEqual(payload["parts"][0]["svg"], projected["parts"][0]["svg"])
        self.assertEqual(
            projected["parts"][0]["svg"] + " Z",
            projected["parts"][1]["svg"],
        )
        svg = self.tester.visual_svg("separate-3d", payload)
        self.assertNotIn('stroke="#dce2e7" stroke-width="1"', svg)
        self.assertIn(">Code (3D view)</text>", svg)
        self.assertIn(">Output (3D view)</text>", svg)

    def test_separate_3d_curvepolygon_keeps_surface_fill_and_hole(self):
        outer = [[0, 0, 2], [5, 0, 2], [0, 5, 2], [0, 0, 2]]
        inner = [[1, 1, 2], [1, 3, 2], [3, 1, 2], [1, 1, 2]]
        points_xyz = [
            {"path": [ring_number, point_number], "point": point}
            for ring_number, ring in enumerate((outer, inner), 1)
            for point_number, point in enumerate(ring, 1)
        ]
        payload = {
            "frames": [
                {"id": "Code", "bounds": [0, 0, 5, 5]},
                {"id": "Output", "bounds": [0, 0, 5, 5]},
            ],
            "parts": [
                {"ord": 1, "source": "Code", "label": "Code", "frame": "Code",
                 "root_type": "POLYGON", "type": "POLYGON", "svg": "unused",
                 "has_z": True, "points_xyz": points_xyz},
                {"ord": 2, "source": "Output", "label": "Output", "frame": "Output",
                 "root_type": "CURVEPOLYGON", "type": "CURVEPOLYGON", "svg": "unused",
                 "has_z": True, "points_xyz": points_xyz},
            ],
        }

        projected = self.tester.project_3d_payload(payload)
        self.assertTrue(all(part["is_3d_face"] for part in projected["parts"]))
        self.assertEqual(2, projected["parts"][1]["svg"].count("M "))

        svg = self.tester.visual_svg("forcecurve", payload)
        self.assertEqual(2, svg.count('class="area"'))
        self.assertEqual(2, svg.count('fill-opacity="1"'))
        self.assertNotIn('class="start-vertex"', svg)

    def test_separate_2d_output_is_not_projected_with_3d_input(self):
        ring_3d = [[0, 0, 2], [0, 5, 2], [5, 0, 2], [0, 0, 2]]
        ring_2d = [[0, 0, 0], [0, 5, 0], [5, 0, 0], [0, 0, 0]]
        payload = {
            "frames": [
                {"id": "Code", "bounds": [0, 0, 5, 5]},
                {"id": "Output", "bounds": [0, 0, 5, 5]},
            ],
            "parts": [
                {"ord": 1, "source": "Code", "label": "input", "frame": "Code",
                 "root_type": "POLYGON", "type": "POLYGON", "svg": "input-2d-placeholder",
                 "has_z": True, "points_xyz": xyz_vertices(ring_3d, ring=True)},
                {"ord": 2, "source": "Output", "label": "output", "frame": "Output",
                 "root_type": "POLYGON", "type": "POLYGON", "svg": "M 0 0 L 0 -5 5 0 Z",
                 "has_z": False, "points_xyz": xyz_vertices(ring_2d, ring=True)},
            ],
        }
        projected = self.tester.project_3d_payload(payload)
        self.assertTrue(projected["frames"][0]["three_dimensional"])
        self.assertNotIn("three_dimensional", projected["frames"][1])
        self.assertNotEqual(payload["parts"][0]["svg"], projected["parts"][0]["svg"])
        self.assertEqual(payload["parts"][1]["svg"], projected["parts"][1]["svg"])
        svg = self.tester.visual_svg("force2d", payload)
        self.assertIn(">Code (3D view)</text>", svg)
        self.assertIn(">Output</text>", svg)
        self.assertNotIn(">Output (3D view)</text>", svg)
        self.assertIn('stroke="#dce2e7" stroke-width="1"', svg)

    def test_3d_frame_projects_zero_z_point_and_line(self):
        payload = {
            "frames": [{"id": "Overlay", "bounds": [0, 0, 3, 2]}],
            "parts": [
                {"ord": 1, "source": "Output", "label": "surface",
                 "root_type": "POLYGON", "type": "POLYGON", "svg": "unused",
                 "has_z": True, "points_xyz": xyz_vertices(
                     [[0, 0, 0], [3, 0, 0], [0, 2, 2], [0, 0, 0]], ring=True
                 )},
                {"ord": 2, "source": "Code", "label": "origin",
                 "root_type": "POINT", "type": "POINT", "x": 2, "y": -1,
                 "has_z": False, "points_xyz": xyz_vertices([[2, 1, 0]])},
                {"ord": 3, "source": "Code", "label": "edge",
                 "root_type": "LINESTRING", "type": "LINESTRING",
                 "svg": "M 0 0 L 3 -2", "closed": False, "has_z": False,
                 "points_xyz": xyz_vertices([[0, 0, 0], [3, 2, 0]])},
            ],
        }
        projected = self.tester.project_3d_payload(payload)
        frame_points = [
            vertex["point"]
            for part in payload["parts"]
            for vertex in (part.get("points_xyz") or [])
        ]
        view = self.tester.choose_3d_view(frame_points)
        expected_x, expected_y, _ = self.tester.project_3d_point([2, 1, 0], view)
        self.assertEqual(expected_x, projected["parts"][1]["x"])
        self.assertEqual(-expected_y, projected["parts"][1]["y"])
        self.assertEqual([[expected_x, -expected_y]], projected["parts"][1]["points"])
        self.assertNotEqual(payload["parts"][2]["svg"], projected["parts"][2]["svg"])
        svg = self.tester.visual_svg("mixed-zero-z", payload)
        self.assertIn('data-postgis-geometry-id="mixed-zero-z-code-1"', svg)
        self.assertIn('data-postgis-geometry-id="mixed-zero-z-code-2"', svg)

    def test_3d_projection_chooses_non_edge_on_view_for_skinny_footprints(self):
        payload = {
            "frames": [{"id": "Overlay", "bounds": [1, 10, 175, 155]}],
            "parts": [
                {"ord": 1, "source": "Code", "label": "poly",
                 "root_type": "POLYGON", "type": "POLYGON", "svg": "unused",
                 "has_z": True, "points_xyz": xyz_vertices(
                     [[175, 150, 5], [20, 40, 5], [35, 45, 5],
                      [50, 60, 5], [100, 100, 5], [175, 150, 5]],
                     ring=True,
                 )},
                {"ord": 2, "source": "Code", "label": "mline",
                 "root_type": "MULTILINESTRING", "type": "LINESTRING",
                 "svg": "unused", "closed": False, "has_z": True,
                 "points_xyz": xyz_vertices(
                     [[175, 155, 2], [20, 40, 20], [50, 60, -2],
                      [125, 100, 1], [175, 155, 1]]
                 )},
                {"ord": 3, "source": "Output", "label": "lol3d",
                 "root_type": "LINESTRING", "type": "LINESTRING",
                 "svg": "unused", "closed": False, "has_z": True,
                 "points_xyz": xyz_vertices([[175, 150, 5], [1, 10, 2]])},
                {"ord": 4, "source": "Output", "label": "lol2d",
                 "root_type": "LINESTRING", "type": "LINESTRING",
                 "svg": "unused", "closed": False, "has_z": False,
                 "points_xyz": xyz_vertices([[175, 150, 0], [1, 10, 0]])},
            ],
        }
        projected = self.tester.project_3d_payload(payload)
        min_x, min_y, max_x, max_y = projected["frames"][0]["bounds"]
        self.assertGreater(max_x - min_x, 0)
        self.assertGreater(max_y - min_y, 0)
        self.assertGreaterEqual(
            min(max_x - min_x, max_y - min_y) / max(max_x - min_x, max_y - min_y),
            0.55,
        )

    def test_3d_projection_preserves_vertical_separation(self):
        ground = self.tester.project_3d_point([4, 7, 0])
        raised = self.tester.project_3d_point([4, 7, 3])
        self.assertEqual(ground[0], raised[0])
        self.assertGreater(raised[1], ground[1])
        self.assertGreater(raised[2], ground[2])

    def test_3d_projection_does_not_collapse_opposite_cube_corners(self):
        origin = self.tester.project_3d_point([0, 0, 0])
        opposite = self.tester.project_3d_point([1, 1, 1])
        self.assertNotEqual(origin[:2], opposite[:2])

    def test_svg_renders_separate_source_frames(self):
        svg = self.tester.visual_svg(
            "comparison",
            {
                "frames": [
                    {"id": "Code", "label": "Code (SRID 2249)", "bounds": [0, 0, 10, 10]},
                    {"id": "Output", "label": "Output (SRID 4326)", "bounds": [0, 0, 10, 10]},
                ],
                "parts": [
                    {"ord": 1, "source": "Code", "label": "Code", "frame": "Code",
                     "type": "LINESTRING", "svg": "M 0 0 L 10 -10", "closed": False},
                    {"ord": 2, "source": "Output", "label": "Output", "frame": "Output",
                     "type": "LINESTRING", "svg": "M 0 0 L 10 -10", "closed": False},
                ],
            },
        )
        self.assertIn('class="frame-label"', svg)
        self.assertIn('>Code (SRID 2249)</text>', svg)
        self.assertIn('>Output (SRID 4326)</text>', svg)
        self.assertIn('data-postgis-frame="Code"', svg)
        self.assertIn('data-postgis-frame="Output"', svg)
        self.assertEqual(2, svg.count('class="plot-background"'))
        self.assertNotIn("font:", svg)

    def test_svg_renders_four_frames_in_two_column_grid(self):
        frames = [
            {"id": label, "label": label, "bounds": [0, 0, 10, 10]}
            for label in ("(a) simple", "(b) not simple", "(c) simple", "(d) not simple")
        ]
        svg = self.tester.visual_svg(
            "four-frames",
            {
                "frames": frames,
                "parts": [
                    {"ord": index, "source": "Code", "label": frame["label"],
                     "frame": frame["id"], "type": "LINESTRING",
                     "svg": "M 0 0 L 10 -10", "closed": False}
                    for index, frame in enumerate(frames, 1)
                ],
            },
        )
        self.assertEqual(4, svg.count('class="plot-background"'))
        self.assertEqual(2, svg.count('y="30" width="270" height="220"'))
        self.assertEqual(2, svg.count('y="270" width="270" height="220"'))
        self.assertIn('x="155" y="20"', svg)
        self.assertIn('x="445" y="260"', svg)
        self.assertIn('height="548" viewBox="0 0 600 548"', svg)

    def test_svg_renders_three_frames_in_one_wide_row(self):
        frames = [
            {"id": label, "label": label, "bounds": [0, 0, 10, 10]}
            for label in ("Code", "linework", "structure")
        ]
        svg = self.tester.visual_svg(
            "three-frames",
            {
                "frames": frames,
                "parts": [
                    {"ord": index, "source": "Output" if index > 1 else "Code",
                     "label": frame["label"], "frame": frame["id"],
                     "type": "LINESTRING", "svg": "M 0 0 L 10 -10", "closed": False}
                    for index, frame in enumerate(frames, 1)
                ],
            },
        )
        self.assertIn('width="900"', svg)
        self.assertIn('viewBox="0 0 900 380"', svg)
        self.assertEqual(3, svg.count('class="plot-background"'))
        self.assertEqual(3, svg.count('y="30" width="273.333333333" height="292"'))
        self.assertIn('>Code</text>', svg)
        self.assertIn('>linework</text>', svg)
        self.assertIn('>structure</text>', svg)

    def test_separate_explicit_inputs_use_their_labels_as_frames(self):
        captured = []
        self.tester.visual_payload = lambda database, layers: captured.extend(layers) or {
            "frames": [
                {"id": layer["requested_frame"], "label": layer["requested_frame"],
                 "bounds": [0, 0, 1, 1]}
                for layer in layers
            ],
            "parts": [{
                "ord": layer["ord"], "source": layer["source"], "label": layer["label"],
                "frame": layer["requested_frame"], "type": "POINT", "x": 0, "y": 0,
            } for layer in layers],
        }
        actual = QueryRows(
            [["POINT(0 0)", "POINT(1 1)"]],
            ["input_(a) simple", "input_(b) not simple"],
        )
        self.tester.render_visual_example("manual", {
            "query": "SELECT 'POINT(9 9)'",
            "visual_id": "labeled-inputs", "label": "labels:labeled-inputs",
            "visual_refentry": "labels", "visual_screen": 3,
            "visual_preferred": True, "visual_kind": "geometry-output",
            "visual_separate_output": True,
        }, actual)
        self.assertEqual(["(a) simple", "(b) not simple"], [
            layer["requested_frame"] for layer in captured
        ])

    def test_svg_uses_one_color_for_polyhedral_surface_faces(self):
        svg = self.tester.visual_svg(
            "solid",
            {"bounds": [0, 0, 2, 1], "parts": [
                {"ord": 1, "source": "Output", "label": "solid",
                 "root_type": "POLYHEDRALSURFACE", "type": "POLYGON",
                 "svg": "M 0 0 L 1 0 0 -1 Z"},
                {"ord": 1, "source": "Output", "label": "solid",
                 "root_type": "POLYHEDRALSURFACE", "type": "POLYGON",
                 "svg": "M 1 0 L 2 0 2 -1 Z"},
            ]},
        )
        self.assertEqual(2, svg.count('stroke="#a62c2b" fill="#a62c2b"'))
        self.assertNotIn('stroke="#d95f3d"', svg)

    def test_svg_draws_bounded_vertices_for_linear_geometries(self):
        svg = self.tester.visual_svg(
            "vertices",
            {"bounds": [0, 0, 1, 1], "parts": [{
                "ord": 1, "source": "Output", "label": "Output", "type": "POLYGON",
                "svg": "M 0 0 L 1 0 0 -1 Z",
                "vertices": [[0, 0], [1, 0], [0, -1], [0, 0]],
            }]},
        )
        self.assertEqual(3, svg.count('class="vertex"'))

    def test_svg_can_hide_output_area_vertices_for_visual_clarity(self):
        payload = {"bounds": [0, 0, 1, 1], "parts": [{
            "ord": 1, "source": "Output", "label": "circle", "type": "POLYGON",
            "svg": "M 0 0 L 1 0 0 -1 Z",
            "vertices": [[0, 0], [1, 0], [0, -1], [0, 0]],
        }]}
        default_svg = self.tester.visual_svg("vertices-default", payload)
        hidden_svg = self.tester.visual_svg(
            "vertices-hidden", payload, hide_output_area_vertices=True
        )
        self.assertEqual(3, default_svg.count('class="vertex"'))
        self.assertNotIn('class="vertex"', hidden_svg)

    def test_svg_suppresses_vertices_for_dense_output_areas(self):
        vertices = [[index, index % 2] for index in range(17)]
        vertices.append(vertices[0])
        svg = self.tester.visual_svg(
            "dense-output-area",
            {"bounds": [0, 0, 16, 1], "parts": [{
                "ord": 1, "source": "Output", "label": "Output", "type": "POLYGON",
                "svg": "M 0 0 L 1 -1 L 2 0 L 3 -1 L 16 0 Z",
                "vertices": vertices,
            }]},
        )
        self.assertNotIn('class="vertex"', svg)

    def test_svg_keeps_vertices_for_dense_input_areas(self):
        vertices = [[index, index % 2] for index in range(17)]
        vertices.append(vertices[0])
        svg = self.tester.visual_svg(
            "dense-input-area",
            {"bounds": [0, 0, 16, 1], "parts": [{
                "ord": 1, "source": "Code", "label": "Code", "type": "POLYGON",
                "svg": "M 0 0 L 1 -1 L 2 0 L 3 -1 L 16 0 Z",
                "vertices": vertices,
            }]},
        )
        self.assertEqual(17, svg.count('class="vertex"'))

    def test_svg_uses_input_underlay_and_direction_arrows_for_open_lines(self):
        svg = self.tester.visual_svg(
            "directions",
            {"bounds": [0, 0, 2, 1], "parts": [
                {"ord": 1, "source": "Code", "label": "Code", "type": "LINESTRING",
                 "svg": "M 0 0 L 2 -1", "closed": False, "vertices": [[0, 0], [2, -1]]},
                {"ord": 2, "source": "Output", "label": "Output", "type": "LINESTRING",
                 "svg": "M 0 0 L 2 -1", "closed": False, "vertices": [[0, 0], [2, -1]]},
            ]},
        )
        self.assertEqual(2, svg.count('class="direction-arrow"'))
        self.assertIn('fill="none" fill-opacity="0"', svg)
        self.assertRegex(svg, r'class="direction-arrow"[^>]*fill="#2878b8"')
        self.assertRegex(svg, r'class="direction-arrow"[^>]*fill="#a62c2b"')
        widths = re.findall(r'class="line"[^>]*stroke-width="([0-9.]+)"', svg)
        self.assertGreater(float(widths[0]), float(widths[1]))
        self.assertNotIn('class="start-vertex"', svg)

    def test_svg_uses_direction_arrows_for_area_rings_when_requested(self):
        payload = {"bounds": [0, 0, 4, 4], "parts": [{
            "ord": 1, "source": "Code", "label": "Code", "type": "POLYGON",
            "svg": "M 0 0 L 4 0 L 4 -4 L 0 -4 Z M 1 -1 L 1 -3 L 3 -3 L 3 -1 Z",
            "vertices": [
                [0, 0], [4, 0], [4, -4], [0, -4], [0, 0],
                [1, -1], [1, -3], [3, -3], [3, -1], [1, -1],
            ],
            "direction_paths": [
                [[0, 0], [4, 0], [4, -4], [0, -4], [0, 0]],
                [[1, -1], [1, -3], [3, -3], [3, -1], [1, -1]],
            ],
        }]}
        default_svg = self.tester.visual_svg("area-directions", payload)
        directed_svg = self.tester.visual_svg(
            "area-directions", payload, show_direction=True
        )
        self.assertNotIn("ring-direction-arrow", default_svg)
        self.assertEqual(4, directed_svg.count("ring-direction-arrow"))
        self.assertRegex(
            directed_svg,
            r'class="direction-arrow ring-direction-arrow"[^>]*fill="#2878b8"',
        )

    def test_svg_marks_start_vertex_for_closed_lines(self):
        svg = self.tester.visual_svg(
            "scroll",
            {"bounds": [0, 0, 2, 1], "parts": [
                {"ord": 1, "source": "Code", "label": "input", "type": "LINESTRING",
                 "svg": "M 0 0 L 2 0 1 -1 0 0", "closed": True,
                 "vertices": [[0, 0], [2, 0], [1, -1], [0, 0]]},
                {"ord": 2, "source": "Output", "label": "scrolled", "type": "LINESTRING",
                 "svg": "M 1 -1 L 0 0 2 0 1 -1", "closed": True,
                 "vertices": [[1, -1], [0, 0], [2, 0], [1, -1]]},
            ]},
        )
        self.assertEqual(2, svg.count('class="start-vertex"'))
        self.assertIn("<title>Start vertex for input</title>", svg)
        self.assertIn("<title>Start vertex for scrolled</title>", svg)
        self.assertIn('data-postgis-start-vertex="true"', svg)
        self.assertNotIn('class="direction-arrow"', svg)

    def test_svg_points_have_contrasting_halo(self):
        svg = self.tester.visual_svg(
            "point-halo",
            {"bounds": [0, 0, 1, 1], "parts": [{
                "ord": 1, "source": "Output", "label": "Output", "type": "POINT",
                "x": 0.5, "y": -0.5,
            }]},
        )
        self.assertRegex(svg, r'class="point"[^>]*stroke="white" stroke-width="[0-9.]+"')

    def test_svg_renders_multipoint_coordinates_as_point_markers(self):
        svg = self.tester.visual_svg(
            "multipoint-markers",
            {"bounds": [0, 0, 10, 10], "parts": [{
                "ord": 1, "source": "Code", "label": "Code", "type": "MULTIPOINT",
                "dimension": 0, "points": [[1, -1], [2, -2], [3, -3]],
                "svg": 'cx="1" cy="-1",cx="2" cy="-2"',
            }]},
        )
        self.assertEqual(3, svg.count('class="point"'))
        self.assertNotIn('d="cx=', svg)

    def test_svg_scales_point_markers_down_for_dense_payloads(self):
        sparse = self.tester.visual_svg(
            "sparse-points",
            {"bounds": [0, 0, 10, 10], "parts": [{
                "ord": 1, "source": "Output", "label": "Output", "type": "POINT",
                "x": 5, "y": -5, "total_points": 8, "source_point_count": 8,
            }]},
        )
        dense = self.tester.visual_svg(
            "dense-points",
            {"bounds": [0, 0, 10, 10], "parts": [{
                "ord": 1, "source": "Output", "label": "Output", "type": "POINT",
                "x": 5, "y": -5, "total_points": 512, "source_point_count": 512,
            }]},
        )
        self.assertLess(self.point_radius(dense), self.point_radius(sparse))

    def test_svg_keeps_medium_density_code_points_below_four_pixel_radius(self):
        svg = self.tester.visual_svg(
            "delaunay-points",
            {"bounds": [0, 0, 200, 200], "parts": [{
                "ord": 1, "source": "Code", "label": "Code", "type": "POINT",
                "x": 5, "y": -5, "total_points": 42, "source_point_count": 42,
            }]},
        )
        self.assertLess(self.point_radius(svg), 4.0)

    def test_svg_scales_vertices_down_for_dense_payloads(self):
        sparse = self.tester.visual_svg(
            "sparse-vertices",
            {"bounds": [0, 0, 4, 1], "parts": [{
                "ord": 1, "source": "Output", "label": "Output", "type": "LINESTRING",
                "svg": "M 0 0 L 4 -1", "closed": False,
                "vertices": [[0, 0], [4, -1]],
                "total_points": 2, "source_point_count": 2,
            }]},
        )
        dense = self.tester.visual_svg(
            "dense-vertices",
            {"bounds": [0, 0, 4, 1], "parts": [{
                "ord": 1, "source": "Output", "label": "Output", "type": "LINESTRING",
                "svg": "M 0 0 L 1 -1 L 2 0 L 3 -1 L 4 0", "closed": False,
                "vertices": [[0, 0], [1, -1], [2, 0], [3, -1], [4, 0]],
                "total_points": 512, "source_point_count": 512,
            }]},
        )
        self.assertLess(self.vertex_radius(dense), self.vertex_radius(sparse))

    def test_svg_scales_point_markers_with_marker_scale_but_caps_them(self):
        smaller = self.tester.visual_svg(
            "marker-scale-small",
            {"bounds": [0, 0, 10, 10], "parts": [{
                "ord": 1, "source": "Output", "label": "Output", "type": "POINT",
                "x": 5, "y": -5, "marker_scale": 0.5, "total_points": 1, "source_point_count": 1,
            }]},
        )
        larger = self.tester.visual_svg(
            "marker-scale-large",
            {"bounds": [0, 0, 10, 10], "parts": [{
                "ord": 1, "source": "Output", "label": "Output", "type": "POINT",
                "x": 5, "y": -5, "marker_scale": 999, "total_points": 1, "source_point_count": 1,
            }]},
        )
        self.assertLess(self.point_radius(smaller), self.point_radius(larger))
        self.assertLess(self.point_radius(larger), self.point_radius(smaller) * 4)

    def test_svg_does_not_arrow_closed_curves(self):
        svg = self.tester.visual_svg(
            "closed-curve",
            {"bounds": [0, 0, 1, 1], "parts": [{
                "ord": 1, "source": "Output", "label": "Output", "type": "CIRCULARSTRING",
                "svg": "M 0 0 C 0 1 1 1 0 0", "closed": True,
            }]},
        )
        self.assertNotIn('class="direction-arrow"', svg)

    def test_render_visual_example_uses_geometry_column_names(self):
        self.tester.visual_payload = lambda database, layers: {
            "bounds": [0, 0, 1, 1],
            "parts": [{
                "ord": layer["ord"], "source": layer["source"], "label": layer["label"],
                "type": "POINT", "x": 0, "y": 0,
            } for layer in layers],
        }
        actual = QueryRows([["POINT(0 0)", "POINT(1 1)"]], ["center", "nearest"])
        visual = self.tester.render_visual_example("manual", {
            "query": "SELECT 1", "visual_id": "labels", "label": "labels:1",
            "visual_refentry": "labels", "visual_screen": 1,
            "visual_preferred": True, "visual_kind": "geometry-output",
        }, actual)
        self.assertIn(">center</text>", visual["svg"])
        self.assertIn(">nearest</text>", visual["svg"])

    def test_render_visual_example_groups_multirow_outputs_by_cluster_column(self):
        captured = []
        self.tester.visual_payload = lambda database, layers: captured.extend(layers) or {
            "bounds": [0, 0, 2, 2],
            "parts": [{
                "ord": layer["ord"], "source": layer["source"], "label": layer["label"],
                "type": "POINT", "x": layer["row_num"] or 0, "y": -(layer["row_num"] or 0),
            } for layer in layers],
        }
        actual = QueryRows(
            [["POINT(0 0)", "A"], ["POINT(1 1)", "A"], ["POINT(2 2)", "B"]],
            ["geom", "cluster"],
        )
        visual = self.tester.render_visual_example("manual", {
            "query": "SELECT 1", "visual_id": "clustered", "label": "clustered:1",
            "visual_refentry": "clustered", "visual_screen": 1,
            "visual_preferred": True, "visual_kind": "geometry-output",
        }, actual)
        self.assertEqual(["cluster A", "cluster A", "cluster B"], [layer["label"] for layer in captured])
        self.assertEqual([1, 1, 2], [layer["ord"] for layer in captured])
        self.assertEqual(2, len(visual["layers"]))
        self.assertEqual(2, visual["svg"].count('class="legend-row"'))

    def test_render_visual_example_uses_text_companion_for_cluster_labels(self):
        captured = []
        self.tester.visual_payload = lambda database, layers: captured.extend(layers) or {
            "bounds": [0, 0, 2, 2],
            "parts": [{
                "ord": layer["ord"], "source": layer["source"], "label": layer["label"],
                "type": "POINT", "x": layer["row_num"] or 0, "y": -(layer["row_num"] or 0),
            } for layer in layers],
        }
        actual = QueryRows(
            [
                ["0", "Melbourne, Sydney", "MULTIPOINT((144.96 -37.81),(151.21 -33.87))"],
                ["1", "Boston, New York", "MULTIPOINT((-71.06 42.36),(-74.01 40.71))"],
            ],
            ["cid", "places", "cluster_geom"],
        )
        visual = self.tester.render_visual_example("manual", {
            "query": "SELECT 1", "visual_id": "clustered-places", "label": "clustered:2",
            "visual_refentry": "clustered", "visual_screen": 2,
            "visual_preferred": True, "visual_kind": "geometry-output",
        }, actual)
        self.assertEqual(
            ["cid 0: Melbourne, Sydney", "cid 1: Boston, New York"],
            [layer["label"] for layer in captured],
        )
        self.assertIn(">cid 0: Melbourne, Sydney</text>", visual["svg"])

    def test_render_visual_example_reads_marker_scale_companion_column(self):
        captured = []
        self.tester.visual_payload = lambda database, layers: captured.extend(layers) or {
            "bounds": [0, 0, 1, 1],
            "parts": [{
                "ord": layer["ord"], "source": layer["source"], "label": layer["label"],
                "type": "POINT", "x": 0, "y": 0, "marker_scale": layer["marker_scale"],
            } for layer in layers],
        }
        actual = QueryRows([["POINT(0 0)", "3.5"]], ["geom", "marker_scale"])
        self.tester.render_visual_example("manual", {
            "query": "SELECT 1", "visual_id": "scaled", "label": "scaled:1",
            "visual_refentry": "scaled", "visual_screen": 1,
            "visual_preferred": True, "visual_kind": "geometry-output",
        }, actual)
        self.assertEqual(3.5, captured[0]["marker_scale"])

    def test_render_visual_example_defers_geometry_point_counts_to_postgis(self):
        captured = []
        self.tester.visual_payload = lambda database, layers: captured.extend(layers) or {
            "bounds": [0, 0, 1, 1],
            "parts": [{
                "ord": layer["ord"], "source": layer["source"], "label": layer["label"],
                "type": "MULTIPOINT", "svg": "M 0 0",
            } for layer in layers],
        }
        actual = QueryRows([["MULTIPOINT((0 0),(1 1))"]], ["points"])
        self.tester.render_visual_example("manual", {
            "query": "SELECT 'MULTIPOINT((0 0),(1 1))'::geometry",
            "visual_id": "point-counts", "label": "point-counts:1",
            "visual_refentry": "point-counts", "visual_screen": 1,
            "visual_preferred": True, "visual_kind": "geometry-output",
        }, actual)
        self.assertTrue(captured)
        self.assertTrue(all(layer["total_points"] is None for layer in captured))
        self.assertTrue(all(layer["source_point_count"] == 1 for layer in captured))

    def test_render_visual_example_uses_typed_hexwkb_without_text_roundtrip(self):
        captured = []
        self.tester.visual_payload = lambda database, layers: captured.extend(layers) or {
            "bounds": [0, 0, 1, 1],
            "parts": [{
                "ord": layer["ord"], "source": layer["source"], "label": layer["label"],
                "type": "POLYGON", "svg": "M 0 0 L 1 0 L 0 -1 Z",
            } for layer in layers],
        }
        actual = QueryRows(
            [["010300000000000000"]],
            ["st_polygonfromtext"],
            ["geometry"],
        )
        visual = self.tester.render_visual_example("manual", {
            "query": "SELECT ST_PolygonFromText('POLYGON EMPTY')",
            "visual_id": "hexwkb", "label": "ST_PolygonFromText:1",
            "visual_refentry": "ST_PolygonFromText", "visual_screen": 1,
            "visual_preferred": True, "visual_kind": "geometry-output",
        }, actual)
        self.assertEqual("010300000000000000", captured[-1]["hexwkb"])
        self.assertIsNone(captured[-1]["wkt"])
        self.assertEqual("Output", captured[-1]["label"])
        self.assertIn(">Output</text>", visual["svg"])

    def test_native_geometry_visual_exposes_raw_hex_alongside_readable_ewkt(self):
        self.tester.visual_payload = lambda database, layers: {
            "bounds": [0, 0, 1, 1],
            "parts": [{
                "ord": layer["ord"], "source": layer["source"], "label": layer["label"],
                "type": "POINT", "x": 1, "y": 2,
            } for layer in layers],
        }
        hexwkb = "0101000000000000000000F03F0000000000000040"
        actual = QueryRows([[hexwkb]], ["point"], ["geometry"])
        visual = self.tester.render_visual_example("manual", {
            "query": "SELECT 'POINT(1 2)'::geometry",
            "expected": [["POINT(1 2)"]],
            "visual_id": "native-output", "label": "native-output:1",
            "visual_refentry": "native-output", "visual_screen": 1,
            "visual_preferred": True, "visual_kind": "geometry-output",
        }, actual)
        self.assertEqual(hexwkb, visual["native_output"])

    def test_render_visual_example_skips_provisional_hex_text(self):
        hexwkb = "010300000000000000"
        actual = QueryRows([[hexwkb]], ["payload"], ["text"])
        visual = self.tester.render_visual_example("manual", {
            "query": f"SELECT '{hexwkb}'::text",
            "expected": [[hexwkb]],
            "visual_id": "hex-text", "label": "hex-text:1",
            "visual_refentry": "hex-text", "visual_screen": 1,
            "visual_preferred": True, "visual_kind": "geometry-output",
        }, actual)
        self.assertIsNone(visual)

    def test_render_visual_example_skips_empty_only_geometry(self):
        self.tester.visual_payload = lambda database, layers: {
            "bounds": None, "frames": None, "parts": None,
        }
        hexwkb = "010300000000000000"
        actual = QueryRows([[hexwkb]], ["empty_polygon"], ["geometry"])
        visual = self.tester.render_visual_example("manual", {
            "query": "SELECT 'POLYGON EMPTY'::geometry",
            "expected": [[hexwkb]],
            "visual_id": "empty", "label": "empty:1",
            "visual_refentry": "empty", "visual_screen": 1,
            "visual_preferred": True, "visual_kind": "geometry-output",
        }, actual)
        self.assertIsNone(visual)

    def test_render_visual_example_exposes_text_when_empty_output_leaves_input_context(self):
        def payload(database, layers):
            code = next(layer for layer in layers if layer["source"] == "Code")
            return {
                "bounds": [0, 0, 1, 1],
                "parts": [{
                    "ord": code["ord"], "source": "Code", "label": code["label"],
                    "type": "LINESTRING", "svg": "M 0 0 L 1 -1",
                }],
            }

        self.tester.visual_payload = payload
        hexwkb = "010300000000000000"
        actual = QueryRows([[hexwkb]], ["empty_result"], ["geometry"])
        visual = self.tester.render_visual_example("manual", {
            "query": "SELECT f('LINESTRING(0 0,1 1)'::geometry)",
            "expected": [[hexwkb]],
            "visual_id": "empty-result", "label": "empty-result:1",
            "visual_refentry": "empty-result", "visual_screen": 1,
            "visual_preferred": True, "visual_kind": "geometry-output",
        }, actual)
        self.assertTrue(visual["output_omitted"])
        self.assertFalse(visual["preferred"])
        self.assertEqual("input-context-fallback", visual["kind"])
        self.assertEqual(["Code"], [layer["source"] for layer in visual["layers"]])

    def test_explicit_input_layers_replace_inferred_code_context(self):
        captured = []
        self.tester.visual_payload = lambda database, layers: captured.extend(layers) or {
            "bounds": [0, 0, 1, 1],
            "parts": [{
                "ord": layer["ord"], "source": layer["source"], "label": layer["label"],
                "type": "POINT", "x": 0, "y": 0,
            } for layer in layers],
        }
        actual = QueryRows([[
            "POINT(0 0)", "LINESTRING(0 0,1 1)", "POLYGON((0 0,1 0,0 1,0 0))",
        ]], ["input_point", "input_line", "result"])

        self.tester.render_visual_example("manual", {
            "query": "SELECT 'POINT(0 0)', 'LINESTRING(0 0,1 1)'",
            "visual_id": "complete-layers", "label": "labels:complete",
            "visual_refentry": "labels", "visual_screen": 2,
            "visual_preferred": True, "visual_kind": "geometry-output",
        }, actual)

        self.assertEqual(["Code", "Code", "Output"], [layer["source"] for layer in captured])
        self.assertEqual(["point", "line", "result"], [layer["label"] for layer in captured])

    def test_named_outputs_retain_inferred_code_context(self):
        captured = []
        self.tester.visual_payload = lambda database, layers: captured.extend(layers) or {
            "bounds": [0, 0, 1, 1],
            "parts": [{
                "ord": layer["ord"], "source": layer["source"], "label": layer["label"],
                "type": "POINT", "x": 0, "y": 0,
            } for layer in layers],
        }
        actual = QueryRows(
            [["POINT(160 40)", "POINT(125.75 115.34)"]],
            ["cp_pt_line", "cp_line_pt"],
        )

        self.tester.render_visual_example("manual", {
            "query": (
                "SELECT ST_ClosestPoint(pt, line), ST_ClosestPoint(line, pt) "
                "FROM (SELECT 'POINT (160 40)'::geometry AS pt, "
                "'LINESTRING (10 30,50 50,30 110,70 90,180 140,130 190)'::geometry AS line) t"
            ),
            "visual_id": "closest-point", "label": "ST_ClosestPoint:1",
            "visual_refentry": "ST_ClosestPoint", "visual_screen": 1,
            "visual_preferred": True, "visual_kind": "geometry-output",
        }, actual)

        self.assertEqual(["Code", "Code", "Output", "Output"], [
            layer["source"] for layer in captured
        ])
        self.assertEqual(["pt", "line", "cp pt line", "cp line pt"], [
            layer["label"] for layer in captured
        ])

    def test_render_visual_example_hides_single_output_column_header(self):
        self.tester.visual_payload = lambda database, layers: {
            "bounds": [0, 0, 1, 1],
            "parts": [{
                "ord": layer["ord"], "source": layer["source"], "label": layer["label"],
                "type": "LINESTRING", "svg": "M 0 0 L 1 -1", "closed": False,
            } for layer in layers],
        }
        actual = QueryRows([["LINESTRING(0 0,1 1)"]], ["polysnapped"])
        visual = self.tester.render_visual_example("manual", {
            "query": "SELECT 1", "visual_id": "default-label", "label": "labels:2",
            "visual_refentry": "labels", "visual_screen": 2,
            "visual_preferred": True, "visual_kind": "geometry-output",
        }, actual)
        self.assertIn(">Output</text>", visual["svg"])
        self.assertNotIn(">polysnapped</text>", visual["svg"])

    def test_render_visual_example_keeps_single_output_alias_in_comparison(self):
        self.tester.visual_payload = lambda database, layers: {
            "bounds": [0, 0, 1, 1],
            "parts": [{
                "ord": layer["ord"], "source": layer["source"], "label": layer["label"],
                "type": "LINESTRING", "svg": "M 0 0 L 1 -1", "closed": False,
            } for layer in layers],
        }
        actual = QueryRows([["LINESTRING(0 0,1 1)"]], ["distance_segment"])
        visual = self.tester.render_visual_example("manual", {
            "query": "SELECT 'LINESTRING(0 0,1 1)'::geometry AS input",
            "visual_id": "comparison-label", "label": "labels:3",
            "visual_refentry": "labels", "visual_screen": 3,
            "visual_preferred": True, "visual_kind": "geometry-output",
        }, actual)
        self.assertIn(">input</text>", visual["svg"])
        self.assertIn(">distance segment</text>", visual["svg"])

    def test_svg_distinguishes_parts_of_multipart_areas(self):
        svg = self.tester.visual_svg(
            "triangulation",
            {"bounds": [0, 0, 2, 1], "parts": [
                {"ord": 1, "source": "Output", "label": "Output", "type": "POLYGON",
                 "svg": "M 0 0 L 1 0 0 -1 Z"},
                {"ord": 1, "source": "Output", "label": "Output", "type": "POLYGON",
                 "svg": "M 1 0 L 2 0 2 -1 Z"},
            ]},
        )
        self.assertIn('stroke="#a62c2b" fill="#a62c2b"', svg)
        self.assertIn('stroke="#d95f3d" fill="#d95f3d"', svg)

    def test_svg_draws_lines_above_area_fills(self):
        svg = self.tester.visual_svg(
            "line-over-area",
            {"bounds": [0, 0, 2, 1], "parts": [
                {"ord": 1, "source": "Code", "label": "observer", "type": "LINESTRING",
                 "svg": "M 0 0 L 2 0", "closed": False},
                {"ord": 2, "source": "Output", "label": "visibility", "type": "POLYGON",
                 "svg": "M 0 0 L 2 0 1 -1 Z"},
            ]},
        )
        self.assertLess(
            svg.index('data-postgis-geometry-id="line-over-area-output-1"'),
            svg.index('data-postgis-geometry-id="line-over-area-code-1"'),
        )

    def test_svg_rejects_empty_non_point_paths(self):
        with self.assertRaisesRegex(RuntimeError, "empty SVG path"):
            self.tester.visual_svg(
                "broken-curve",
                {"bounds": [0, 0, 1, 1], "parts": [
                    {"ord": 1, "source": "Output", "label": "Output", "type": "MULTICURVE", "svg": ""}
                ]},
            )

    def test_svg_rejects_negative_arc_radii(self):
        with self.assertRaisesRegex(RuntimeError, "invalid SVG path"):
            self.tester.visual_svg(
                "bad-arc",
                {"bounds": [0, 0, 1, 1], "parts": [{
                    "ord": 1, "source": "Output", "label": "Output",
                    "type": "CIRCULARSTRING", "svg": "M 0 0 A -1 -1 0 0 0 1 -1",
                    "closed": False,
                }]},
            )

    def test_grid_is_bounded_for_large_coordinate_offsets(self):
        svg = self.tester.visual_svg(
            "large-offset",
            {"bounds": [1e16, 1e16, 1e16 + 2, 1e16 + 2], "parts": [
                {"ord": 1, "source": "Code", "label": "Code", "type": "POINT", "x": 1e16, "y": -1e16}
            ]},
        )
        self.assertLessEqual(svg.count("<line "), 64)

    def test_multi_frame_3d_payload_uses_shared_projected_bounds(self):
        def triangle_part(ord_, label, frame, size):
            points = [
                [0, 0, 0],
                [size, 0, 0],
                [0, size, 0],
                [0, 0, size],
            ]
            return {
                "ord": ord_,
                "source": "Output",
                "label": label,
                "frame": frame,
                "type": "POLYHEDRALSURFACE",
                "has_z": True,
                "points_xyz": [
                    {"path": [face, index], "point": point}
                    for face, ring in enumerate((
                        [points[0], points[1], points[2], points[0]],
                        [points[0], points[1], points[3], points[0]],
                        [points[0], points[2], points[3], points[0]],
                        [points[1], points[2], points[3], points[1]],
                    ), 1)
                    for index, point in enumerate(ring, 1)
                ],
            }

        projected = self.tester.project_3d_payload({
            "frames": [
                {"id": "outer", "bounds": [0, 0, 4, 4]},
                {"id": "inner", "bounds": [0, 0, 2, 2]},
            ],
            "parts": [
                triangle_part(1, "outer", "outer", 4),
                triangle_part(2, "inner", "inner", 2),
            ],
        })
        self.assertEqual(
            projected["frames"][0]["bounds"],
            projected["frames"][1]["bounds"],
        )

    def test_unparseable_documented_output_falls_back_to_input_context(self):
        calls = []
        def payload(_database, layers):
            calls.append(layers)
            if any(layer["source"] == "Output" for layer in layers):
                raise RuntimeError("geometry contains non-closed rings; parse error")
            return {"bounds": [0, 0, 1, 1], "parts": [
                {"ord": 1, "source": "Code", "label": "Code", "type": "LINESTRING", "svg": "M 0 0 L 1 -1"}
            ]}
        self.tester.visual_payload = payload
        visual = self.tester.render_visual_example("db", {
            "query": "SELECT ST_AsText(f('LINESTRING(0 0,1 1)'))",
            "label": "example:1", "visual_id": "example", "visual_refentry": "example",
            "visual_screen": 1, "visual_preferred": True, "visual_kind": "geometry-output",
        }, [["POLYGON((0 0,1 0,1 1))"]])
        self.assertEqual(2, len(calls))
        self.assertTrue(visual["output_omitted"])
        self.assertFalse(visual["preferred"])
        self.assertEqual("input-context-fallback", visual["kind"])

    def test_publish_visual_examples_replaces_stale_assets(self):
        with tempfile.TemporaryDirectory() as directory:
            target = Path(directory) / "visual-examples"
            target.mkdir()
            (target / "stale.svg").write_text("stale", encoding="utf-8")
            self.tester.publish_visual_examples(
                target,
                [{
                    "id": "one", "source": "ST_Buffer:1", "svg": "<svg/>\n",
                    "kind": "geometry-output", "preferred": True,
                    "refentry": "ST_Buffer", "screen": 1, "output_omitted": False,
                    "native_output": "01010000000000000000000000000000000000000000",
                    "layers": [{
                        "id": "one-code-1", "source": "Code", "label": "center",
                        "row": None, "column": None, "dimension": 0,
                        "closed": None, "srid": 0, "frame": "Overlay",
                    }],
                }],
            )
            self.assertFalse((target / "stale.svg").exists())
            self.assertEqual("<svg/>\n", (target / "one.svg").read_text(encoding="utf-8"))
            self.assertEqual(
                [{
                    "id": "one", "kind": "geometry-output", "preferred": True,
                    "refentry": "ST_Buffer", "screen": 1, "source": "ST_Buffer:1",
                    "output_omitted": False,
                    "width": None, "height": None,
                    "native_output": "01010000000000000000000000000000000000000000",
                    "layers": [{
                        "id": "one-code-1", "source": "Code", "label": "center",
                        "row": None, "column": None, "dimension": 0,
                        "closed": None, "srid": 0, "frame": "Overlay",
                    }],
                }],
                json.loads((target / "manifest.json").read_text(encoding="utf-8")),
            )
            manifest_xml = (target / "manifest.xml").read_text(encoding="utf-8")
            self.assertIn('id="one"', manifest_xml)
            self.assertIn('preferred="true"', manifest_xml)
            self.assertIn('<native-output format="hexewkb">01010000', manifest_xml)
            self.assertIn('<layer id="one-code-1" source="Code" label="center"', manifest_xml)
            self.assertIn('dimension="0" srid="0" frame="Overlay"', manifest_xml)

    def test_publish_visual_examples_records_svg_dimensions(self):
        with tempfile.TemporaryDirectory() as directory:
            target = Path(directory) / "visual-examples"
            self.tester.publish_visual_examples(
                target,
                [{
                    "id": "one", "source": "ST_Buffer:1",
                    "svg": '<svg xmlns="http://www.w3.org/2000/svg" width="600" height="402"></svg>\n',
                    "kind": "geometry-output", "preferred": True,
                    "refentry": "ST_Buffer", "screen": 1, "output_omitted": False,
                    "layers": [],
                }],
            )
            manifest_json = json.loads((target / "manifest.json").read_text(encoding="utf-8"))
            self.assertEqual("600", manifest_json[0]["width"])
            self.assertEqual("402", manifest_json[0]["height"])
            manifest_xml = (target / "manifest.xml").read_text(encoding="utf-8")
            self.assertIn('width="600"', manifest_xml)
            self.assertIn('height="402"', manifest_xml)


if __name__ == "__main__":
    unittest.main()
