2020/09/20: We released
torchsparsev1.1, which is significantly faster than our
torchsparsev1.0 and is also achieves 1.9x speedup over MinkowskiEngine v0.5 alpha when running MinkUNet18C!
2020/08/30: We released
torchsparsev1.0.
We release
torchsparse, a high-performance computing library for efficient 3D sparse convolution. This library aims at accelerating sparse computation in 3D, in particular the Sparse Convolution operation.
The major advantage of this library is that we support all computation on the GPU, especially the kernel map construction (which is done on the CPU in latest MinkowskiEngine V0.4.3).
You may run the following command to install torchsparse.
pip install --upgrade git+https://github.com/mit-han-lab/torchsparse.git
Note that this library depends on Google's sparse hash map project. In order to install this library, you may run
sudo apt-get install libsparsehash-dev
on Ubuntu servers. If you are not sudo, please clone Google's codebase, compile it and install locally. Finally, add the path to this library to your
CPLUS_INCLUDE_PATHenvironmental variable.
For GPU server users, we currently support PyTorch 1.6.0 + CUDA 10.2 + CUDNN 7.6.2. For CPU users, we support PyTorch 1.6.0 (CPU version), MKLDNN backend is optional.
Our SPVNAS project (ECCV2020) is built with torchsparse. You may navigate to this project and follow the instructions in that codebase to play around.
Here, we also provide a walk-through on some important concepts in torchsparse.
In torchsparse, we have two data structures for point cloud storage, namely
torchsparse.SparseTensorand
torchsparse.PointTensor. Both structures has two data fields
C(coordinates) and
F(features). In
SparseTensor, we assume that all coordinates are integer and do not duplicate. However, in
PointTensor, all coordinates are floating-point and can duplicate.
The way to convert a point cloud to
SparseTensorso that it can be consumed by networks built with Sparse Convolution or Sparse Point-Voxel Convolution is to use the function
torchsparse.utils.sparse_quantize. An example is given here:
inds, labels, inverse_map = sparse_quantize(pc, feat, labels, return_index=True, return_invs=True)
where
pc,
feat,
labelscorresponds to point cloud (coordinates, should be integer), feature and ground-truth. The
indsdenotes unique indices in the point cloud coordinates, and
inverse_mapdenotes the unique index each point is corresponding to. The
inverse mapis used to restore full point cloud prediction from downsampled prediction.
To combine a list of
SparseTensors to a batch, you may want to use the
torchsparse.utils.sparse_collate_fnfunction.
Detailed results are given in SemanticKITTI dataset preprocessing code in our SPVNAS project.
The computation interface in torchsparse is straightforward and very similar to original PyTorch. An example here defines a basic convolution block:
class BasicConvolutionBlock(nn.Module): def __init__(self, inc, outc, ks=3, stride=1, dilation=1): super().__init__() self.net = nn.Sequential( spnn.Conv3d(inc, outc, kernel_size=ks, dilation=dilation, stride=stride), spnn.BatchNorm(outc), spnn.ReLU(True) )def forward(self, x): out = self.net(x) return out
where
spnndenotes
torchsparse.nn, and
spnn.Conv3dmeans 3D sparse convolution operation,
spnn.BatchNormand
spnn.ReLUdenotes 3D sparse tensor batchnorm and activations, respectively. We also support direct convolution kernel call via
torchsparse.nn.functional, for example:
outputs = torchsparse.nn.functional.conv3d(inputs, kernel, stride=1, dilation=1, transpose=False)
where we need to define
inputs(SparseTensor),
kernel(of shape k^3 x OC x IC when k > 1, or OC x IC when k = 1, where k denotes the kernel size and IC, OC means input / output channels). The
outputsis still a SparseTensor.
Detailed examples are given in here, where we use the
torchsparse.nn.functionalinterfaces to implement weight-shared 3D-NAS modules.
Sparse hash map query is important in 3D sparse computation. It is mainly used to infer a point's memory location (i.e. index) given its coordinates. For example, we use this operation in kernel map construction part of 3D sparse convolution, and also sparse voxelization / devoxelization in Sparse Point-Voxel Convolution. Here, we provide the following example for hash map API:
source_hash = torchsparse.nn.functional.sphash(torch.floor(source_coords).int()) target_hash = torchsparse.nn.functional.sphash(torch.floor(target_coords).int()) idx_query = torchsparse.nn.functional.sphashquery(source_hash, target_hash)
In this example,
sphashis the function converting integer coordinates to hashing. The
sphashquery(source_hash, target_hash)performs the hash table lookup. Here, the hash map has key
target_hashand value corresponding to point indices in the target point cloud tensor. For each point in the
source_coords, we find the point index in
target_coordswhich has the same coordinate as it.
We here provides an entire training example with dummy input here. In this example, we cover
You are also welcomed to check out our SPVNAS project to implement training / inference with real data.
We benchmark the performance of our torchsparse and latest MinkowskiEngine V0.4.3 here, latency is measured on NVIDIA GTX 1080Ti GPU:
| Network | Latency (ME V0.4.3) | Latency (torchsparse V1.0.0) | | :----------------------: | :-----------------: | :--------------------------: | | MinkUNet18C (MACs / 10) | 224.7 | 124.3 | | MinkUNet18C (MACs / 4) | 244.3 | 160.9 | | MinkUNet18C (MACs / 2.5) | 269.6 | 214.3 | | MinkUNet18C | 323.5 | 294.0 |
If you find this code useful, please consider citing:
@inproceedings{ tang2020searching, title = {Searching Efficient 3D Architectures with Sparse Point-Voxel Convolution}, author = {Tang, Haotian* and Liu, Zhijian* and Zhao, Shengyu and Lin, Yujun and Lin, Ji and Wang, Hanrui and Han, Song}, booktitle = {European Conference on Computer Vision}, year = {2020} }
This library is inspired by MinkowskiEngine, SECOND and SparseConvNet.