geometric_kernels 0.4

A Python Package offering geometric kernels in NumPy, TensorFlow, PyTorch, and Jax.

GeometricKernels

GeometricKernels is a library that implements kernels — most importantly, heat and Matérn kernels — on non-Euclidean spaces such as Riemannian manifolds, graphs and meshes. This enables kernel methods — in particular Gaussian process models — to be deployed on such spaces.

Installation

0. [Optionally] create and activate a new virtual environment.

   You can use uv

   uv venv
   

   or conda

   conda create -n [env_name] python=3.[version]
   conda activate [env_name]
   

   or virtualenv

   virtualenv [env_name]
   source [env_name]/bin/activate
   

1. Install the library in the active environment by running

   pip install geometric_kernels
   

   NOTE: If you use uv, swap pip with uv pip everywhere. If you initialized a project with uv init, use uv add instead.

   If you want to install specific GitHub branch called [branch], run

   pip install "git+https://github.com/geometric-kernels/GeometricKernels@[branch]"
   

2. Install a backend of your choice

   We use LAB to support multiple backends (e.g., TensorFlow, Jax, PyTorch). However, you are not required to install all of them on your system to use the GeometricKernels package. Simply install the backend (and (optionally) a GP package) of your choice. For example,

   • Tensorflow

     pip install tensorflow tensorflow-probability
     

     Optionally, you can install the Tensorflow-based Gaussian processes library GPflow, for which we provide a frontend.

     pip install gpflow
     

   • PyTorch

     pip install torch
     

     Optionally, you can install the PyTorch-based Gaussian processes library GPyTorch, for which we provide a frontend.

     pip install gpytorch
     

   • JAX (the cpu version)—the gpu and tpu versions can be installed similarly.

     pip install "jax[cpu]"
     

     Optionally, you can install the JAX-based Gaussian processes library GPJax, for which we provide a frontend.

     pip install gpjax
     

A basic example

This example shows how to compute a 3x3 kernel matrix for the Matern52 kernel on the standard two-dimensional sphere. It relies on the numpy-based backend. Look up the information on how to use other backends in the documentation.

# Import a backend.
import numpy as np
# Import the geometric_kernels backend.
import geometric_kernels
# Import a space and an appropriate kernel.
from geometric_kernels.spaces import Hypersphere
from geometric_kernels.kernels import MaternGeometricKernel

# Create a manifold (2-dim sphere).
hypersphere = Hypersphere(dim=2)

# Define 3 points on the sphere.
xs = np.array([[0., 0., 1.], [0., 1., 0.], [1., 0., 0.]])

# Initialize kernel.
kernel = MaternGeometricKernel(hypersphere)
params = kernel.init_params()
params["nu"] = np.array([5/2])
params["lengthscale"] = np.array([1.])

# Compute and print out the 3x3 kernel matrix.
print(np.around(kernel.K(params, xs), 2))

This should output

[[1.   0.36 0.36]
 [0.36 1.   0.36]
 [0.36 0.36 1.  ]]

Notebooks

You can find numerous example notebooks for GeometricKernels here in the docs or here in the tree. These cover every space (e.g., sphere, graph), backend (e.g., Torch, JAX), frontend (e.g., GPyTorch, GPflow), and more.

Application Examples

Looking for practical use cases of GeometricKernels? Check out

• PeMS Regression: A benchmark suite for graph node regression with uncertainty. This project employs GeometricKernels among other tools, and offers processed data, baseline models, and an example notebook for experiments on graph-structured data. Notably, in this benchmark, geometric Gaussian processes built with GeometricKernels have been shown to outperform various alternative methods, including ensembles of graph neural networks and Bayesian graph neural networks.

• Bayesian optimization demonstration: A minimal notebook illustrating the use of GeometricKernels with the botorch library for Bayesian optimization. This is a simple, self-contained example designed to demonstrate core concepts rather than to reflect a real-world scenario.

Documentation

The documentation for GeometricKernels is available on a separate website.

For development and running the tests

If you want to contribute to the library, thank you! Follow these instructions to set up the environment and run the tests and the code formatting.

Initialize the virtual environment.

make venv

You can change the python version and the venv directory like this:

make venv UV_PYTHON=3.11 VENV_DIR=.venv

NOTE: We use uv for development. It is not strictly necessary, and if you don't want to use it, you can still run make lint and make test if you set UV_RUN= to be empty. You will need to set up the environment yourself.

Install all backends and the dev requirements (Pytest, black, etc.). This will install all the backends.

make install

NOTE: If not using uv, you can still install the dev requirements via pip install -e .[dev].

Run the style checks

make lint

Run the tests

make test

If you want to run Jupyter with your uv development environment, check out this page.

Example: If you want to learn how to implement your own space or kernel component, checkout the CustomSpacesAndKernels.ipynb notebook.

If you have a question

Post it in issues using the "How do I do ..." and other issues template and the "question" label.

This link chooses the right template and label for you.

Citation

If you are using GeometricKernels, please cite the following paper:

@article{mostowsky2024,
      title = {The GeometricKernels Package: Heat and Matérn Kernels for Geometric Learning on Manifolds, Meshes, and Graphs},
      author = {Peter Mostowsky and Vincent Dutordoir and Iskander Azangulov and Noémie Jaquier and Michael John Hutchinson and Aditya Ravuri and Leonel Rozo and Alexander Terenin and Viacheslav Borovitskiy},
      year = {2024},
      journal = {arXiv:2407.08086},
}

Furthermore, please consider citing the theoretical papers the package is based on. You can find the relevant references for any space in

• the docstring of the respective space class,
• at the end of the respective tutorial notebook.