2. TensorLy’s backend system
In short, you can write your code using TensorLy and you can transparently combine it and execute with any of the backends. Currently we support NumPy PyTorch, MXNet, JAX, TensorFlow and CuPy as backends.
To represent tensors and for numerical computation, TensorLy supports several backends transparently: the ubiquitous NumPy (the default), MXNet, and PyTorch. For the end user, the interface is exactly the same, but under the hood, a different library is used to represent multi-dimensional arrays and perform computations on these.
In other words, you write your code using TensorLy and can then decide whether the computation is done using NumPy, PyTorch or MXNet.
2.2. Why backends?
The goal of TensorLy is to make tensor methods accessible. While NumPy needs no introduction, other backends such as MXNet and PyTorch backends are especially useful as they allows to perform transparently computation on CPU or GPU. Last but not least, using MXNet or PyTorch as a backend, we are able to combine tensor methods and deep learning easily!
2.3. How do I change the backend?
To change the backend, e.g. to NumPy, you can change the value of
default_backend in tensorly/__init__.
Alternatively during the execution, assuming you have imported TensorLy as
import tensorly as tl, you can change the backend in your code by calling
NumPy is installed by default with TensorLy if you haven’t already installed it. However, to keep dependencies as minimal as possible, and to not complexify installation, neither MXNet nor PyTorch are installed. If you want to use them as backend, you will have to install them first. It is easy however, simply refer to their respective installation instructions:
Once you change the backend, all the computation is done using that backend.
2.4. Context of a tensor
Different backends have different parameters associated with the tensors. For instance, in NumPy we traditionally set the dtype when creating an ndarray, while in mxnet we also have to change the context (GPU or CPU), with the ctx argument. Similarly, in PyTorch, we might want to create a FloatTensor for CPU and a cuda.FloatTensor for GPU.
To handle this difference, we implemented a context function, that, given a tensor, returns a dictionary of values characterising that tensor. A function getting a tensor as input and creating a new tensor should use that context to create the new tensor.
import tensorly as tl def trivial_fun(tensor): """ Trivial function that takes a tensor and create a new one with value tensor + 2... """ # context is a dict of values context = tl.context(tensor) # when creating a new tensor we use these as parameters new_tensor = tl.tensor(tensor + 2, **context) return new_tensor
2.5. Basic functions
We have isolated the basic functions required for tensor methods in the backend, and provide a uniform API using wrappers when necessary. In practice, this means that function like min, max, reshape, etc, are accessible from the backend:
import tensorly as tl import numpy as np tl.set_backend('pytorch') # or any other backend tensor = tl.tensor(np.random.random((10, 10, 10))) # This will call the correct function depending on the backend min_value = tl.min(tensor) unfolding = tl.unfold(tensor, mode=0) U, S, V = tl.partial_svd(unfolding, n_eigenvecs=5)
This will allow your code to work transparently with any of the backend.
2.6. Case study: TensorLy and PyTorch
Let’s go through the creation and decomposition of a tensor, using PyTorch.
2.6.1. On CPU
First, we import tensorly and set the backend:
import tensorly as tl tl.set_backend('pytorch')
Now, let’s create a random tensor using the
from tensorly import random tensor = random.random_tensor((10, 10, 10)) # tensor is a PyTorch Tensor!
We can decompose it easily, here using a Tucker decomposition: First, we reate a decomposition instance, which keeps the number of parameters the same and with a random initialization. We then fit it to our tensor.
from tensorly.decomposition import Tucker decomp = Tucker(rank='same', init='random') cp_tensor = decomp.fit_transform(tensor)
You can reconstruct the full tensor and measure the reconstruction error:
rec = cp_tensor.to_tensor() error = tl.norm(tensor - rec)/tl.norm(tensor)
2.6.2. On GPU
Now, imaging you want everything to run on GPU: this is very easy using TensorLy and the PyTorch backend: you simply send the tensor to the GPU!
There are to main ways to do this: either you specify the context during the creation of the tensor or you use pytorch tensors’ properties to send them to the desired device post-creation.
# Specify context during creation tensor = random.random_tensor(shape=(10, 10, 10), device='cuda', dtype=tl.float32) # Posthoc tensor = random.random_tensor(shape=(10, 10, 10)) tensor = tensor.to('cuda')
The rest is exactly the same, nothing more to do!
decomp = Tucker(rank='same', init='random') cp_tensor = decomp.fit_transform(tensor) # Runs on GPU!