Set SQLite constant_rows_max, and use it in join_data_coordinates.

With QG generation now relying heavily on join_data_coordinates, we
were hitting this limit sometimes.

python/lsst/daf/butler/direct_query_driver/_driver.py

@@ -156,7 +156,7 @@ def __init__(
         self._exit_stack: ExitStack | None = None
         self._raw_page_size = raw_page_size
         self._postprocessing_filter_factor = postprocessing_filter_factor
-        self._constant_rows_limit = constant_rows_limit
+        self._constant_rows_limit = min(constant_rows_limit, db.get_constant_rows_max())
         self._cursors: set[_Cursor] = set()
 
     def __enter__(self) -> None:

python/lsst/daf/butler/registry/databases/sqlite.py

@@ -398,6 +398,12 @@ def constant_rows(
         ]
         return sqlalchemy.sql.union_all(*selects).alias(name)
 
+    def get_constant_rows_max(self):
+        # Docstring inherited.
+        # This is the default SQLITE_MAX_COMPOUND_SELECT (see
+        # https://www.sqlite.org/limits.html):
+        return 500
+
     @property
     def has_distinct_on(self) -> bool:
         # Docstring inherited.