diff --git a/README.md b/README.md
index e259682..f61f7af 100644
--- a/README.md
+++ b/README.md
@@ -6,8 +6,15 @@ This is a repository for code implementing the **smoothed subpixel projection (S
 * G. Romano, R. Arrieta, and S. G. Johnson, [“Differentiating through binarized topology changes: Second-order subpixel-smoothed projection,”](http://arxiv.org/abs/2601.10737) arXiv.org e-Print archive, 2601.10737, January 2026.
 * R. Arrieta, G. Romano, and S. G. Johnson, [“Hyperparameter-free minimum-lengthscale constraints for topology optimization,”](http://arxiv.org/abs/2507.16108) arXiv.org e-Print archive, 2507.16108, July 2025.
 
+## Documentation
+
+* [Supported features](docs/features.md) — a running list of what the Julia and Python
+  packages each implement, along with known limitations and gaps.
+
 ## Installation
 
+### Python
+
 Install the PyPI distribution:
 
 ```bash
@@ -25,3 +32,15 @@ For local development:
 ```bash
 python -m pip install -e ".[dev]"
 ```
+
+### Julia
+
+The Julia package is not registered yet, so install it from this repository:
+
+```julia
+using Pkg
+Pkg.develop(path="src/julia/SSP")
+```
+
+See [`src/julia/SSP/README.md`](src/julia/SSP/README.md) for usage of both the high-level
+and low-level Julia APIs.
diff --git a/docs/features.md b/docs/features.md
new file mode 100644
index 0000000..a8c2266
--- /dev/null
+++ b/docs/features.md
@@ -0,0 +1,43 @@
+# Supported Features
+
+The Julia package lives in [`src/julia/SSP`](../src/julia/SSP) and the Python package in
+[`src/python/ssp_topopt`](../src/python/ssp_topopt). The table below tracks what each one
+currently implements; **please keep it up to date when adding or removing functionality.**
+
+| Feature | Julia (`SSP`) | Python (`ssp_topopt`) |
+| --- | --- | --- |
+| Conic ("hat") filter | ✅ `conic_filter` | ✅ `conic_filter` |
+| Filter radius from an eroded threshold point | ❌ | ✅ `get_conic_radius_from_eta_e` |
+| Plain tanh projection | ❌ (internal only) | ✅ `tanh_projection` |
+| First-order subpixel smoothing (SSP1), linear interpolation | ✅ `ssp1_linear` | ✅ `ssp1_bilinear` |
+| First-order subpixel smoothing (SSP1), cubic interpolation | ✅ `ssp1` | ❌ |
+| Second-order subpixel smoothing (SSP2), differentiable through topology changes | ✅ `ssp2` | ✅ `ssp2` |
+| Finite and infinite projection strength (0 ≤ β ≤ ∞) | ✅ | ✅ |
+| Dilation/erosion of the projected contour | ✅ `dilation_distance` argument | ❌ |
+| Minimum-lengthscale constraints for solid and void | ✅ `constraint_solid`, `constraint_void` | ❌ |
+| Lengthscale constraints compatible with any SSP order | ✅ (constraints act on `rho_filtered`/`rho_projected`) | ❌ |
+| Reverse-mode automatic differentiation | ✅ hand-written adjoints, exposed to Zygote.jl and friends through a ChainRulesCore.jl extension | ✅ through JAX (`grad`, `jit`, `vmap`) |
+| Dimensionality | N-dimensional code paths (only 2D is currently tested) | 2D only |
+| Periodic filter axes | ❌ | ✅ `periodic_axes` argument of `conic_filter` |
+| Low-level `init`/`solve!`/`adjoint_solve!` API with reduced allocations | ✅ | ❌ |
+| Explicit control over padding/boundary conditions, kernels, interpolation, and projection target points | ✅ (low-level API) | ❌ |
+
+## Known limitations
+
+* The subpixel fill factor is the analytic expression for a *circular* smoothing kernel, so
+  both implementations assume an isotropic grid (`dx == dy`). Julia asserts that all grid
+  steps are equal; Python takes a single scalar `resolution` in the projection routines
+  (`conic_filter` does accept an anisotropic `resolution`).
+* Only 2D usage is covered by the tests and examples in this repository, even though the
+  Julia routines are written generically over the number of dimensions.
+* The high-level Julia `conic_filter` always pads by replicating the boundary values.
+  Other padding styles (`FillPadding`, `Inner`) are only reachable through the low-level API.
+
+## Not yet supported
+
+Contributions welcome — these are known gaps rather than fundamental limitations:
+
+* Python: cubic-interpolation SSP1, dilation/erosion, and minimum-lengthscale constraints.
+* Python: a low-level API with reusable workspaces.
+* Julia: `get_conic_radius_from_eta_e`-style helpers and periodic filter axes.
+* Both: validated 3D usage and anisotropic grid spacings in the projection.
diff --git a/src/julia/SSP/README.md b/src/julia/SSP/README.md
index 2c40f38..e747cf9 100644
--- a/src/julia/SSP/README.md
+++ b/src/julia/SSP/README.md
@@ -3,6 +3,9 @@
 A Smoothed Subpixel Projection (SSP) package for topology optimization in Julia.
 Supports N-dimensional data and reverse-mode automatic differentiation with minimal allocations.
 
+For a feature-by-feature comparison of this package with its Python cousin, see
+[`docs/features.md`](../../../docs/features.md).
+
 ## Usage
 
 This package provides a high-level API nearly identical to its python cousin.