Kernels in Jax.
Project description
This project has now been incorporated into GPJax.
Kernels in Jax.
Introduction
JaxKern is Python library for working with kernel functions in JAX. We currently support the following kernels:
- Stationary
- Radial basis function (Squared exponential)
- Matérn
- Powered exponential
- Rational quadratic
- White noise
- Periodic
- Non-stationary
- Linear
- Polynomial
- Non-Euclidean
- Graph kernels
In addition to this, we implement kernel approximations using the Random Fourier feature approach.
Example
The following code snippet demonstrates how the first order Matérn kernel can be computed and, subsequently, approximated using random Fourier features.
import jaxkern as jk
import jax.numpy as jnp
import jax.random as jr
key = jr.PRNGKey(123)
# Define the points on which we'll evaluate the kernel
X = jr.uniform(key, shape = (10, 1), minval=-3., maxval=3.)
Y = jr.uniform(key, shape = (20, 1), minval=-3., maxval=3.)
# Instantiate the kernel and its parameters
kernel = jk.Matern32()
params = kernel.init_params(key)
# Compute the 10x10 Gram matrix
Kxx = kernel.gram(params, X)
# Compute the 10x20 cross-covariance matrix
Kxy = kernel.cross_covariance(params, X, Y)
# Build a RFF approximation
approx = RFF(kernel, num_basis_fns = 5)
rff_params = approx.init_params(key)
# Build an approximation to the Gram matrix
Qff = approx.gram(rff_params, X)
Code Structure
All kernels are supplied with a gram
and cross_covariance
method. When computing a Gram matrix, there is often some structure in the data (e.g., Markov) that can be exploited to yield a sparse matrix. To instruct JAX how to operate on this, the return type of gram
is a Linear Operator from JaxLinOp.
Within GPJax, all kernel computations are handled using JaxKern.
Documentation
A full set of documentation is a work in progress. However, many of the details in JaxKern can be found in the GPJax documentation.
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
Hashes for jaxkern-nightly-0.0.5.dev20230630.tar.gz
Algorithm | Hash digest | |
---|---|---|
SHA256 | 35796fc16a2e4997ec02b9dbcc3204ebed283ef8814fe09d6bf8509c02577a3c |
|
MD5 | c10297c97f736a6144bf5f98d42141d2 |
|
BLAKE2b-256 | 66b302a466757f6ef11cbc0d431ed5ecfa3e1c527c8122f22ad98940c873131f |
Hashes for jaxkern_nightly-0.0.5.dev20230630-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | b14d704caffc6785142333970953369df79efd72b1e74786f18925838d63ff89 |
|
MD5 | 4d88bcc75cac981605b753acf1823e3e |
|
BLAKE2b-256 | fe66daa932200292ef2d83b9650fb7a9afb2745bdb2e63ab6c39c95f59de75a8 |