doc: update documentation links and improve clarity on API usage

fix: enforce 1-D input requirement for histogram functions and update tests

Author: ElouenGinat

README.md

@@ -136,7 +136,13 @@ uv run pytest
 
 ## Documentation
 
-See the [documentation index](docs/index.rst), [API comparison guide](docs/api_comparison.rst), and [demo notebook](docs/demo.ipynb) for the current API reference and usage examples.
+Full documentation is hosted at **[khiops.github.io/khisto-python](https://khiops.github.io/khisto-python/)**.
+
+- [API Reference](https://khiops.github.io/khisto-python/array/histogram/index.html) — NumPy-like histogram API
+- [Matplotlib Integration](https://khiops.github.io/khisto-python/matplotlib/index.html) — `hist` plotting function
+- [Core API](https://khiops.github.io/khisto-python/core/index.html) — full access to histogram granularity levels
+- [API Comparison](https://khiops.github.io/khisto-python/api_comparison.html) — side-by-side with NumPy and Matplotlib
+- [Demo Notebook](https://khiops.github.io/khisto-python/demo.html) — interactive walkthrough
 
 ## License
 

docs/api_comparison.rst

@@ -197,7 +197,7 @@ Key Differences
      - ``bins``, ``stacked``, and ``weights`` raise a ``TypeError``
    * - **Multiple datasets**
      - Supported
-     - Sequences are concatenated into one dataset. Sequences of different lengths are not supported.
+     - Not supported. Only 1-D arrays are accepted.
 
 Usage Comparison
 """"""""""""""""

pyproject.toml

@@ -7,6 +7,45 @@ license = "BSD-3-Clause-Clear"
 license-files = ["LICENSE"]
 requires-python = ">=3.10"
 dependencies = ["numpy>=2.0"]
+authors = [
+    {name = "Elouen Ginat", email = "elouen.ginat@orange.com"},
+]
+keywords = [
+    "histogram",
+    "binning",
+    "optimal",
+    "khiops",
+    "statistics",
+    "data-analysis",
+    "density-estimation",
+    "distribution",
+    "numpy",
+    "matplotlib",
+    "visualization",
+]
+classifiers = [
+    "Development Status :: 4 - Beta",
+    "Intended Audience :: Developers",
+    "Intended Audience :: Science/Research",
+    "License :: OSI Approved :: BSD License",
+    "Operating System :: OS Independent",
+    "Programming Language :: Python :: 3",
+    "Programming Language :: Python :: 3.10",
+    "Programming Language :: Python :: 3.11",
+    "Programming Language :: Python :: 3.12",
+    "Programming Language :: Python :: 3.13",
+    "Topic :: Scientific/Engineering",
+    "Topic :: Scientific/Engineering :: Mathematics",
+    "Topic :: Scientific/Engineering :: Visualization",
+    "Typing :: Typed",
+]
+
+[project.urls]
+Homepage = "https://github.com/khiops/khisto-python"
+Documentation = "https://khiops.github.io/khisto-python/"
+Repository = "https://github.com/khiops/khisto-python"
+Issues = "https://github.com/khiops/khisto-python/issues"
+Changelog = "https://github.com/khiops/khisto-python/blob/main/CHANGELOG.md"
 
 [tool.khiops]
 khiops-version = "11.0.1-a.3"

src/khisto/array/histogram/api.py

@@ -62,7 +62,7 @@ def histogram(
     Parameters
     ----------
     a : array_like
-        Input data. The histogram is computed over the flattened array.
+        Input data. Must be 1-dimensional.
     range : tuple of (float, float), optional
         The lower and upper range of the bins. Values outside the range are
         ignored. If not provided, the range is ``(a.min(), a.max())``.
@@ -105,7 +105,13 @@ def histogram(
        histograms for large-scale data sets. Computational Statistics & Data
        Analysis, 180:0-0, 2023.
     """
-    arr = np.asarray(a, dtype=np.float64).flatten()
+    arr = np.asarray(a, dtype=np.float64)
+
+    if arr.ndim != 1:
+        raise ValueError(
+            f"Expected 1-D array, got {arr.ndim}-D array instead. "
+            "Reshape your data or flatten it before calling histogram."
+        )
 
     if max_bins is not None and max_bins <= 0:
         raise ValueError("max_bins must be a positive integer or None.")

src/khisto/matplotlib/hist.py

@@ -72,10 +72,8 @@ def hist(
 
     Parameters
     ----------
-    x : array_like or sequence of array_like
-        Input data. Nested arrays are concatenated and histogrammed as a
-        single dataset. This behavior is different from Matplotlib's hist, which treats a sequence of
-        arrays as multiple datasets and plots them separately.
+    x : array_like
+        Input data. Must be 1-dimensional.
     range : tuple of (float, float), optional
         Lower and upper range of the bins. Values outside the range are
         ignored.

tests/array/test_histogram.py

@@ -172,10 +172,8 @@ def test_bin_edges_length(self, normal_data):
         hist, bin_edges = histogram(normal_data)
         assert len(bin_edges) == len(hist) + 1
 
-    def test_2d_array_flattening(self):
-        """Test that 2D arrays are flattened."""
+    def test_2d_array_raises(self):
+        """Test that 2D arrays raise ValueError."""
         data_2d = np.array([[1, 2, 3], [4, 5, 6]])
-        hist, bin_edges = histogram(data_2d)
-
-        # Should process all 6 values
-        assert np.sum(hist) == 6
+        with pytest.raises(ValueError, match="Expected 1-D array"):
+            histogram(data_2d)

tests/plot/test_matplotlib_histogram.py

@@ -151,15 +151,6 @@ def test_reverse_cumulative_frequency_histogram(self, normal_data):
         assert np.isclose(n[0], len(normal_data))
         plt.close(fig)
 
-    def test_sequence_of_arrays_is_combined(self, normal_data):
-        """Test that sequences of arrays are combined into one histogram."""
-        fig, ax = plt.subplots()
-        datasets = [normal_data[:500], normal_data[500:]]
-        n, bins, patches = hist(datasets, ax=ax)
-
-        assert np.sum(n) == len(normal_data)
-        plt.close(fig)
-
     def test_unsupported_bins_parameter(self, normal_data):
         """Test that bins raises a clear error message."""
         fig, ax = plt.subplots()