PointFlow: 3D Point Cloud Generation with Continuous Normalizing Flows
Abstract
As 3D point clouds become the representation of choice for multiple vision and graphics applications, the ability to synthesize or reconstruct highresolution, highfidelity point clouds becomes crucial. Despite the recent success of deep learning models in discriminative tasks of point clouds, generating point clouds remains challenging. This paper proposes a principled probabilistic framework to generate 3D point clouds by modeling them as a distribution of distributions. Specifically, we learn a twolevel hierarchy of distributions where the first level is the distribution of shapes and the second level is the distribution of points given a shape. This formulation allows us to both sample shapes and sample an arbitrary number of points from a shape. Our generative model, named PointFlow, learns each level of the distribution with a continuous normalizing flow. The invertibility of normalizing flows enables the computation of the likelihood during training and allows us to train our model in the variational inference framework. Empirically, we demonstrate that PointFlow achieves stateoftheart performance in point cloud generation. We additionally show that our model can faithfully reconstruct point clouds and learn useful representations in an unsupervised manner. The code is available at https://github.com/stevenygd/PointFlow.
mygraygray0.6
1 Introduction
Point clouds are becoming popular as a 3D representation because they can capture a much higher resolution than voxel grids and are a stepping stone to more sophisticated representations such as meshes. Learning a generative model of point clouds could benefit a wide range of point cloud synthesis tasks such as reconstruction and superresolution, by providing a better prior of point clouds. However, a major roadblock in generating point clouds is the complexity of the space of point clouds. A cloud of points corresponding to a chair is best thought of as samples from a distribution that corresponds to the surface of the chair, and the chair itself is best thought of as a sample from a distribution of chair shapes. As a result, in order to generate a chair according to this formulation, we need to characterize a distribution of distributions, which is underexplored by existing generative models.
In this paper, we propose PointFlow, a principled generative model for 3D point clouds that learns a distribution of distributions: the former being the distribution of shapes and the latter being the distribution of points given a shape. Our key insight is that instead of directly parametrizing the distribution of points in a shape, we model this distribution as an invertible parameterized transformation of 3D points from a prior distribution (e.g., a 3D Gaussian). Intuitively, under this model, generating points for a given shape involves sampling points from a generic Gaussian prior, and then moving them according to this parameterized transformation to their new location in the target shape, as illustrated in Figure LABEL:fig:teaser. In this formulation, a given shape is then simply the variable that parametrizes such transformation, and a category is simply a distribution of this variable. Interestingly, we find that representing this distribution too as a transformation of a prior distribution leads to a more expressive model of shapes. In particular, we use the recently proposed continuous normalizing flow framework to model both kinds of transformations [NormalizingFlow, neuralODE, ffjord].
This parameterization confers several advantages. The invertibility of these transformations allows us to not just sample but also estimate probability densities. The ability to estimate probability densities in turn allows us to train these models in a principled manner using the variational inference framework [VAE], where we maximize a variational lower bound on the loglikelihood of a training set of point clouds. This probabilistic framework for training further lets us avoid the complexities of training GANs or handcrafting good distance metrics for measuring the difference between two sets of points. Experiments show that PointFlow outperforms previous stateoftheart generative models of point clouds, and achieves compelling results in point cloud reconstruction and unsupervised feature learning.
2 Related work
Deep learning for point clouds. Deep learning has been introduced to improve performance in various point cloud discriminative tasks including classification [pointnet, pointnet++, foldingnet, deepsets], segmentation [pointnet, shoef2019pointwise], and critical points sampling [learning2sample]. Recently, substantial progress has been made in point cloud synthesis tasks such as autoencoding [achlioptas_L3DP, foldingnet, atlasnet], singleview 3D reconstruction [pointsetgen, GALGA, DeformNetFD, EfficientDP, MultiresolutionTN], stereo reconstruction [usenko2015reconstructing], and point cloud completion [PUnet, ECNet]. Many point cloud synthesis works convert a point distribution to a matrix by sampling ( is predefined) points from the distribution so that existing generative models are readily applicable. For example, Gadelha \etal [MultiresolutionTN] apply variational autoencoders (VAEs) [VAE] and Zamorski \etal [3DAAE] apply adversarial autoencoders (AAEs) [AAE] to point cloud generation. Achlioptas \etal [achlioptas_L3DP] explore generative adversarial networks (GANs) [GAN, WGAN, iWGAN] for point clouds in both raw data space and latent space of a pretrained autoencoder. In the above methods, the autoencoders are trained with heuristic loss functions that measure the distance between two point sets, such as Chamfer distance (CD) or earth mover’s distance (EMD). Sun \etal [pointgrow] apply autoregressive models [pixelCNN] with a discrete point distribution to generate one point at a time, also using a fixed number of points per shape.
However, treating a point cloud as a fixeddimensional matrix has several drawbacks. First, the model is restricted to generate a fixed number of points. Getting more points for a particular shape requires separate upsampling models such as [PUnet, ECNet, patchbasedupsampling]. Second, it ignores the permutation invariance property of point sets, which might lead to suboptimal parameter efficiency. Heuristic set distances are also far from ideal objectives from a generative modeling perspective since they make the original probabilistic interpretation of VAE/AAE no longer applicable when used as the reconstruction objective. In addition, exact EMD is slow to compute while approximations could lead to biased or noisy gradients. CD has been shown to incorrectly favor point clouds that are overly concentrated in the mode of the marginal point distribution [achlioptas_L3DP].
Some recent works introduce sophisticated decoders consisting of a cascade [foldingnet] or a mixture [atlasnet] of smaller decoders to map one (or a mixture of) 2D uniform distribution(s) to the target point distribution, overcoming the shortcomings of using a fixed number of points. However, they still rely on heuristic set distances that lack a probabilistic guarantee. Also, their methods only learn the distribution of points for each shape, but not the distribution of shapes. Li \etal [PCGAN] propose a “sandwiching” reconstruction objective that combines a variant of WGAN [WGAN] loss with EMD. They also train another GAN in the latent space to learn shape distribution, similar to Achlioptas \etal [achlioptas_L3DP]. In contrast, our method is simply trained endtoend by maximizing a variational lower bound on the loglikelihood, does not require multistage training, and does not have any instability issues common for GANbased methods.
Generative models. There are several popular frameworks of deep generative models, including generative adversarial networks [GAN, WGAN, StyleGAN], variational autoencoders [VAE, stochastic_bp], autoregressive models [PixelRNN, pixelCNN], and flowbased models [Nice, NormalizingFlow, RealNVP, GLOW]. In particular, flowbased models and autoregressive models can both perform exact likelihood evaluation, while flowbased models are much more efficient to sample from. Flowbased models have been successfully applied to a variety of generation tasks such as image generation [GLOW, RealNVP, Nice], video generation [VideoFlow], and voice synthesis [WaveGlow]. Also, there has been recent work that combines flows with other generative models, such as GANs [FlowGAN, danihelka2017comparison], autoregressive models [NAF, MAF, IAF], and VAEs [IAF, SylvesterNF, VLAE, NormalizingFlow, SylvesterNF, neuralODE, ffjord].
Most existing deep generative models aim at learning the distribution of fixeddimensional variables. Learning the distribution of distributions, where the data consists of a set of sets, is still underexplored. Edwards and Storkey [nerual_stats] propose a hierarchical VAE named Neural Statistician that consumes a set of sets. They are mostly interested in the fewshot case where each set only has a few samples. Also, they are focused on classifying sets or generating new samples from a given set. While our method is also applicable to these tasks, our focus is on learning the distribution of sets and generating new sets (point clouds in our case). In addition, our model employs a tighter lower bound on the loglikelihood, thanks to the use of normalizing flow in modeling both the reconstruction likelihood and the prior.
3 Overview
Consider a set of shapes from a particular class of object, where each shape is represented as a set of 3D points . As discussed in Section 1, each point is best thought of as being sampled from a point distribution , usually a uniform distribution over the surface of an object . Each shape is itself a sample from a distribution over shapes that captures what shapes in this category look like.
Our goal is to learn the distribution of shapes, each shape itself being a distribution of points. In other words, our generative model should be able to both sample shapes and sample an arbitrary number of points from a shape.
We propose to use continuous normalizing flows to model the distribution of points given a shape. A continuous normalizing flow can be thought of as a vector field in the D Euclidean space, which induces a distribution of points through transforming a generic prior distribution (\eg, a standard Gaussian). To sample points from the induced distribution, we simply sample points from the prior and move them according to the vector field. Moreover, the continuous normalizing flow is invertible, which means we can move data points back to the prior distribution to compute the exact likelihood. This model is highly intuitive and interpretable, allowing a close inspection of the generative process as shown in Figure LABEL:fig:teaser.
We parametrize each continuous normalizing flow with a latent variable that represents the shape. As a result, modeling the distribution of shapes can be reduced to modeling the distribution of the latent variable. Interestingly, we find continuous normalizing flow also effective in modeling the latent distribution. Our full generative model thus consists of two levels of continuous normalizing flows, one modeling the shape distribution by modeling the distribution of the latent variable, and the other modeling the point distribution given a shape.
In order to optimize the generative model, we construct a variational lower bound on the loglikelihood by introducing an inference network that infers a latent variable distribution from a point cloud. Here, we benefit from the fact that the invertibility of the continuous normalizing flow enables likelihood computation. This allows us to train our model endtoend in a stable manner, unlike previous work based on GANs that requires twostage training [achlioptas_L3DP, PCGAN]. As a side benefit, we find the inference network learns a useful representation of point clouds in an unsupervised manner.
4 Background
4.1 Continuous normalizing flow
A normalizing flow [NormalizingFlow] is a series of invertible mappings that transform an initial known distribution to a more complicated one. Formally, let denote a series of invertible transformations we want to apply to a latent variable with a distribution . is the output variable. Then the probability density of the output variable is given by the change of variables formula:
(1) 
where can be computed from using the inverse flow: . In practice, are usually instantiated as neural networks with an architecture that makes the determinant of the Jacobian easy to compute. The normalizing flow has been generalized from a discrete sequence to a continuous transformation [ffjord, neuralODE] by defining the transformation using a continuoustime dynamic , where is a neural network that has an unrestricted architecture. The continuous normalizing flow (CNF) model for with a prior distribution at the start time can be written as:
(2) 
and can be computed using the inverse flow . A blackbox ordinary differential equation (ODE) solver can been applied to estimate the outputs and the input gradients of a continuous normalizing flow [ffjord, neuralODE].
4.2 Variational autoencoder
Suppose we have a random variable that we are building generative models for. The variational autoencoder (VAE) is a framework that allows one to learn from a dataset of observations of [VAE, stochastic_bp]. The VAE models the data distribution via a latent variable with a prior distribution , and a decoder which captures the (hopefully simpler) distribution of given . During training, it additionally learns an inference model (or encoder) . The encoder and decoder are jointly trained to maximize a lower bound on the loglikelihood of the observations
(3) 
which is also called the evidence lower bound (ELBO). One can interpret ELBO as the sum of the negative reconstruction error (the first term) and a latent space regularizer (the second term). In practice, is usually modeled as a diagonal Gaussian whose mean and standard deviation are predicted by a neural network with parameters . To efficiently optimize the ELBO, sampling from is done by reparametrizing as , where .
5 Model
We now have the paraphernalia needed to define our generative model of point clouds. Using the terminology of the VAE, we need three modules: the encoder that encodes a point cloud into a shape representation , a prior over shape representations, and a decoder that models the distribution of points given the shape representation. We use a simple permutationinvariant encoder to predict , following the architecture in Achlioptas \etal [achlioptas_L3DP]. We use continuous normalizing flows for both the prior and the generator , which are described below.
5.1 Flowbased point generation from shape representations
We first decompose the reconstruction loglikelihood of a point set into the sum of loglikelihood of each point
(4) 
We propose to model using a conditional extension of CNF. Specifically, a point in the point set is the result of transforming some point in the prior distribution using a CNF conditioned on :
where defines the continuoustime dynamics of the flow conditioned on . Note that the inverse of is given by with . The reconstruction likelihood of follows equation (2):
(5) 
Note that can be computed in closed form with the Gaussian prior.
5.2 Flowbased prior over shapes
Although it is possible to use a simple Gaussian prior over shape representations, it has been shown that a restricted prior tends to limit the performance of VAEs [VLAE]. To alleviate this problem, we use another CNF to parametrize a learnable prior. Formally, we rewrite the KL divergence term in Equation 3 as
(6) 
where is the entropy and is the prior distribution with learnable parameters , obtained by transforming a simple Gaussian with a CNF:
where defines the continuoustime dynamics of the flow . Similarly as described above, the inverse of is given by with . The log probability of the prior distribution can be computed by:
(7) 
5.3 Final training objective
Plugging Equation 4, 5, 6, 7 into Equation 3, the ELBO of a point set can be finally written as
(8) 
Our model is trained endtoend by maximizing the ELBO of all point sets in the dataset
(9) 
We can interpret this objective as the sum of three parts:

Prior: encourages the encoded shape representation to have a high probability under the prior, which is modeled by a CNF as described in Section 5.2. We use the reparameterization trick [VAE] to enable a differentiable Monte Carlo estimate of the expectation:
where and are mean and standard deviation of the isotropic Gaussian posterior and is simply set to . is sampled from the standard Gaussian distribution .

Reconstruction likelihood: is the reconstruction loglikelihood of the input point set, computed as described in Section 5.1. The expectation is also estimated using Monte Carlo sampling.

Posterior Entropy: is the entropy of the approximated posterior:
All the training details (\eg, hyperparameters, model architectures) are included in Section B of the appendix.
5.4 Sampling
To sample a shape representation, we first draw then pass it through to get . To generate a point given a shape representation , we first sample a point from , then pass through conditioned on to produce a point on the shape : . To sample a point cloud with size , we simply repeat it for times. Combining these two steps allows us to sample a point cloud with points from our model:
6 Experiments
In this section, we first introduce existing metrics for evaluating point cloud generation, discuss their limitations, and introduce a new metric that overcomes these limitations. We then compare the proposed method with previous stateoftheart generative models of point clouds, using both previous metrics and the proposed one. We additionally evaluate the reconstruction and representation learning ability of the autoencoder part of our model.
6.1 Evaluation metrics
Following prior work, we use Chamfer distance (CD) and earth mover’s distance (EMD) to measure the similarity between point clouds. Formally, they are defined as follows:
where and are two point clouds with the same number of points and is a bijection between them. Note that most previous methods use either CD or EMD in their training objectives, which tend to be favored if evaluated under the same metric. Our method, however, does not use CD or EMD during training.
# Parameters (M)  JSD ()  MMD ()  COV (%, )  1NNA (%, )  

Category  Model  Full  Gen  CD  EMD  CD  EMD  CD  EMD  
Airplane  rGAN  7.22  6.91  7.44  0.261  5.47  42.72  18.02  93.58  99.51 
lGAN (CD)  1.97  1.71  4.62  0.239  4.27  43.21  21.23  86.30  97.28  
lGAN (EMD)  1.97  1.71  3.61  0.269  3.29  47.90  50.62  87.65  85.68  
PCGAN  9.14  1.52  4.63  0.287  3.57  36.46  40.94  94.35  92.32  
PointFlow (ours)  1.61  1.06  4.92  0.217  3.24  46.91  48.40  75.68  75.06  
Training set      \textcolormygray6.61  \textcolormygray0.226  3.08  \textcolormygray42.72  \textcolormygray49.14  70.62  67.53  
Chair  rGAN  7.22  6.91  11.5  2.57  12.8  33.99  9.97  71.75  99.47 
lGAN (CD)  1.97  1.71  4.59  2.46  8.91  41.39  25.68  64.43  85.27  
lGAN (EMD)  1.97  1.71  2.27  2.61  7.85  40.79  41.69  64.73  65.56  
PCGAN  9.14  1.52  3.90  2.75  8.20  36.50  38.98  76.03  78.37  
PointFlow (ours)  1.61  1.06  1.74  2.42  7.87  46.83  46.98  60.88  59.89  
Training set      1.50  1.92  7.38  57.25  55.44  59.67  58.46  
Car  rGAN  7.22  6.91  12.8  1.27  8.74  15.06  9.38  97.87  99.86 
lGAN (CD)  1.97  1.71  4.43  1.55  6.25  38.64  18.47  63.07  88.07  
lGAN (EMD)  1.97  1.71  2.21  1.48  5.43  39.20  39.77  69.74  68.32  
PCGAN  9.14  1.52  5.85  1.12  5.83  23.56  30.29  92.19  90.87  
PointFlow (ours)  1.61  1.06  0.87  0.91  5.22  44.03  46.59  60.65  62.36  
Training set      0.86  \textcolormygray1.03  \textcolormygray5.33  48.30  51.42  57.39  53.27 
Let be the set of generated point clouds and be the set of reference point clouds with . To evaluate generative models, we first consider the three metrics introduced by Achlioptas \etal [achlioptas_L3DP]:

JensenShannon Divergence (JSD) are computed between the marginal point distributions:
where . and are marginal distributions of points in the reference and generated sets, approximated by discretizing the space into voxels and assigning each point to one of them. However, it only considers the marginal point distributions but not the distribution of individual shapes. A model that always outputs the “average shape” can obtain a perfect JSD score without learning any meaningful shape distributions.

Coverage (COV) measures the fraction of point clouds in the reference set that are matched to at least one point cloud in the generated set. For each point cloud in the generated set, its nearest neighbor in the reference set is marked as a match:
where can be either CD or EMD. While coverage is able to detect mode collapse, it does not evaluate the quality of generated point clouds. In fact, it is possible to achieve a perfect coverage score even if the distances between generated and reference point clouds are arbitrarily large.

Minimum matching distance (MMD) is proposed to complement coverage as a metric that measures quality. For each point cloud in the reference set, the distance to its nearest neighbor in the generated set is computed and averaged:
where can be either CD or EMD. However, MMD is actually very insensitive to lowquality point clouds in , since they are unlikely to be matched to real point clouds in . In the extreme case, one can imagine that consists of mostly very lowquality point clouds with one additional point cloud in each mode of , yet has a reasonably good MMD score.
As discussed above, all existing metrics have their limitations. As will be shown later, we also empirically find all these metrics sometimes give generated point clouds even better scores than real point clouds, further casting doubt on whether they can ensure a fair model comparison. We therefore introduce another metric that we believe is better suited for evaluating generative models of point clouds:

1nearest neighbor accuracy (1NNA) is proposed by LopezPaz and Oquab [tst] for twosample tests, assessing whether two distributions are identical. It has also been explored as a metric for evaluating GANs [1nnacc]. Let and be the nearest neighbor of in . NNA is the leaveoneout accuracy of the 1NN classifier:
where is the indicator function. For each sample, the 1NN classifier classifies it as coming from or according to the label of its nearest sample. If and are sampled from the same distribution, the accuracy of such a classifier should converge to given a sufficient number of samples. The closer the accuracy is to , the more similar and are, and therefore the better the model is at learning the target distribution. In our setting, the nearest neighbor can be computed using either CD or EMD. Unlike JSD, 1NNA considers the similarity between shape distributions rather than between marginal point distributions. Unlike COV and MMD, 1NNA directly measures distributional similarity and takes both diversity and quality into account.
6.2 Generation
We compare our method with three existing generative models for point clouds: rawGAN [achlioptas_L3DP], latentGAN [achlioptas_L3DP], and PCGAN [PCGAN], using their official implementations that are either publicly available or obtained by contacting the authors. We train each model using point clouds from one of the three categories in the ShapeNet [shapenet] dataset: airplane, chair, and car. The point clouds are obtained by sampling points uniformly from the mesh surface. All points in each category are normalized to have zeromean per axis and unitvariance globally. Following prior convention [achlioptas_L3DP], we use points for each shape during both training and testing, although our model is able to sample an arbitrary number of points. We additionally report the performance of point clouds sampled from the training set, which is considered as an upper bound since they are from the target distribution.
In Table 1, we report the performance of different models, as well as their number of parameters in total (full) or in the generative pathways (gen). We first note that all the previous metrics (JSD, MMD, and COV) sometimes assign a better score to point clouds generated by models than those from the training set (marked in \textcolormygraygray). The 1NNA metric does not seem to have this problem and always gives a better score to shapes from the training set. Our model outperforms all baselines across all three categories according to 1NNA and also obtains the best score in most cases as evaluated by other metrics. Besides, our model has the fewest parameters among compared models. In Section C of the appendix, we perform additional ablation studies to show the effectiveness of different components of our model. Figure 2 shows some examples of novel point clouds generated by our model. Figure 3 shows examples of point clouds reconstructed from given inputs.
Method  MN40 (%)  MN10 (%) 
SPH [SPH]  68.2  79.8 
LFD [LFD]  75.5  79.9 
TL Network [TLNet]  74.4   
VConvDAE [VConvDAEDV]  75.5  80.5 
3DGAN [3dgan]  83.3  91.0 
lGAN (EMD) [achlioptas_L3DP]  84.0  95.4 
lGAN (CD) [achlioptas_L3DP]  84.5  95.4 
PointGrow [pointgrow]  85.7   
MRTNetVAE [MultiresolutionTN]  86.4   
FoldingNet [foldingnet]  88.4  94.4 
lGAN (CD) [achlioptas_L3DP] ^{†}  87.0  92.8 
lGAN (EMD) [achlioptas_L3DP] ^{†}  86.7  92.2 
PointFlow (ours)  86.8  93.7 

We run the official code of lGAN on our preprocessed dataset using the same encoder architecture as our model.
6.3 Autoencoding
We further quantitatively compare the reconstruction ability of our flowbased autoencoder with lGAN [achlioptas_L3DP] and AtlasNet [atlasnet]. Following the setting of AtlasNet, the stateoftheart in this task, we train our autoencoder on all shapes in the ShapeNet dataset. The autoencoder is trained with the reconstruction likelihood objective only. At test time, we sample points per shape and split them into an input set and a reference set, each consisting of points. We then compute the distance (CD or EMD) between the reconstructed input set and the reference set ^{1}^{1}1We use a separate reference set because we expect the autoencoder to learn the point distribution. Exactly reproducing the input points is acceptable behavior, but should not be given a higher score than randomly sampling points from the underlying point distribution.. Although our model is not directly trained with EMD, it obtains the best EMD score, even higher than lGAN trained with EMD and AtlasNet which has more than times more parameters.
Model  # Parameters (M)  CD  EMD 

lGAN (CD) [achlioptas_L3DP]  1.77  7.12  7.95 
lGAN (EMD) [achlioptas_L3DP]  1.77  8.85  5.26 
AtlasNet [atlasnet]  44.9  5.13  5.97 
PointFlow (ours)  1.30  7.54  5.18 
6.4 Unsupervised representation learning
We finally evaluate the representation learning ability of our autoencoders. Specifically, we extract the latent representations of our autoencoder trained in the full ShapeNet dataset and train a linear SVM classifier on top of it on ModelNet10 or ModelNet40 [modelnet]. Only for this task, we normalize each individual point cloud to have zeromean per axis and unitvariance globally, following prior works [deepsets, achlioptas_L3DP]. We also apply random rotations along the gravity axis when training the autoencoder.
A problem with this task is that different authors have been using different encoder architectures with a different number of parameters, making it hard to perform an applestoapples comparison. In addition, different authors may use different preprocessing protocols (as also noted by Yang \etal [foldingnet]), which could also affect the numbers.
In Table 2, we still show the numbers reported by previous papers, but also include a comparison with lGAN [achlioptas_L3DP] trained using the same encoder architecture and the exact same data as our model. On ModelNet10, the accuracy of our model is and higher than lGAN (EMD) and lGAN (CD), respectively. On ModelNet40, the performance of the three models is very close.
7 Conclusion and future works
In this paper, we propose PointFlow, a generative model for point clouds consisting of two levels of continuous normalizing flows trained with variational inference. Future work includes applications to other tasks such as point cloud reconstruction from a single image.
8 Acknowledgment
This work was supported in part by a research gift from Magic Leap. Xun Huang was supported by NVIDIA Graduate Fellowship.
References
Appendix A Overview
In the appendix, we first describe the detailed hyperparameters and model architectures for our experiments in Section B. We then compare our model with additional baselines to understand the effect of different model components in Section C. Limitations and typical failure cases are discussed in Section D. Finally, additional visualizations of latent space tSNE, interpolations and flow transformations are presented in Section E, Section F, and Section G respectively.
Appendix B Training details
In this section, we provide details about our network architectures and training hyperparameters. We will release the code to reproduce our experiments. Please refer to algorithm 1 for the detailed training procedure.
Encoder. The architecture of our encoder follows that of Achlioptas \etal [achlioptas_L3DP]. Specifically, we first use 1D Convolution with filter size , , , and to process each point independently and then use max pooling to create a dimension feature as done in PointNet [pointnet]. Such a feature is invariant to the permutation of points due to the maxpooling. Finally, we apply a threelayer MLP with and hidden dimensions to convert the permutation invariant feature to a dimension one. For the unsupervised representation learning experiment, we set following convention. For all other experiments, is set to .
CNF prior. The CNF prior models the distribution . We follow FFJORD [ffjord]’s released code to use three concatsquash layers to model the dynamics . A concatsquash layer is defined as:
(10) 
where , , , , , and are all trainable parameters and is the sigmoid function. uses three concatsquash layers with a hidden dimension . Tanh is used as the nonlinearity between layers.
We use a Moving Batch Normalization layer to learn the scale of each dimension before and after the CNF, following FFJORD’s released code [ffjord]. Specifically, Moving Batch Normalization is defined as
(11) 
where and are trainable parameters, Different from batch normalization proposed by Ioffe and Szegedy [bn], and are running averages of the batch mean and standard deviation. MovingBatchNorm is invertible : . Its log determinant is given as:
(12) 
CNF decoder. The CNF decoder models the reconstruction likelihood . We extend the concatsquash layer to condition on the latent vector :
(13) 
where are all learnable parameters. The CNF decoder uses four conditional concatsquash layers with a hidden dimension to model the dynamic . The nonlinearity between layers is Tanh. Similar to the CNF prior model, we also add a Moving Batch Normalization layer before and after the CNF. In this case, all 3D points (from different shapes) from a batch are used to compute the batch statistics.
Hyperparameters. We use an Adam optimizer with an initial learning rate , , and . The learning rate decays linearly to starting at the 2000 epoch and ends at the epoch. We do not use any weight decay. We also learn the integration time during training by backpropogation [neuralODE].
Appendix C Additional comparisons
# Parameters (M)  JSD ()  MMD ()  COV (%, )  1NNA (%, )  

Category  Model  Full  Gen  CD  EMD  CD  EMD  CD  EMD  
Airplane  NS [nerual_stats]  2.29  1.00  1.74  0.655  4.51  7.81  4.51  99.61  99.61 
VAECNF  1.47  0.92  6.30  0.261  3.35  41.98  46.17  88.64  82.72  
WGANCNF  1.75  1.06  4.29  0.254  3.23  42.47  48.40  75.80  75.68  
PointFlow (ours)  1.61  1.06  4.92  0.217  3.24  46.91  48.40  75.68  75.06  
Training set      6.61  0.226  3.08  42.72  49.14  70.62  67.53 
In this section, we compare our model to more baselines to show the effectiveness of the model design. The first baseline is Neural Statistician (NS) [nerual_stats], a stateoftheart generative model for sets. We modify its official code for generating 2D spatial coordinates of MNIST digits to make it work with 3D point cloud coordinates. We use the same encoder architecture as our model, and use the VAE decoder provided by authors with the input dimension changed from to . It differs from our model mainly in 1) using VAEs instead of CNFs to model the reconstruction likelihood, and 2) using a simple Gaussian prior instead of a flowbased one. The second baseline is VAECNF, where we use the CNF to model the reconstruction likelihood but not prior. Specifically, the VAECNF optimizes ELBO in the following form:
(14) 
where is a standard Gaussian and is the KLdivergence. As another baseline, we follow lGAN [achlioptas_L3DP] to train a WGAN [iWGAN] in the latent space of our pretrained autoencoder. Both the discriminator and the generator are MLP with batch normalization between layers. The generator has three layers with hidden dimensions 256. The discriminator has three layers with hidden dimensions 512.
The results are presented in Table 4. Neural Statistician [nerual_stats] is able to learn the marginal point distribution but fails to learn the correct shape distribution, as it obtains the best marginal JSD but very poor scores according to metrics that measure similarities between shape distributions. Also, using a flexible prior parameterized by a CNF (PointFlow) is better than using a simple Gaussian prior (VAECNF) or a prior learned with a latent GAN (WGANCNF) that requires twostage training.
Appendix D Limitation and failure cases
In this section, we discuss the limitation of our model and present visualizations of difficult cases where our model fails. As mentioned in FFJORD [ffjord], each integration requires evaluating the neural networks modeling the dynamics multiple times. The number of function evaluations tends to increase as the training proceeds since the dynamic becomes more complex and more function evaluations are needed to achieve the same numerical precision. This issue limits our model size and makes the convergence slow. Grathwohl \etalindicate that using regularization such as weight decay could alleviate such an issue, but we empirically find that using regularization tends to hurt performance. Future advances in invertible models like CNF might help improve this issue. Typical failure case appears when reconstructing or generating the rare shape or shapes with many thin structures as presented in Figure 4.
Appendix E Latent space visualizations
We provide visualization of the sampled latent vectors in Figure 5. We sample latent vectors and run tSNE [maaten2008visualizing] to visualize these latent vectors in 2D. Shapes with similar styles are close in the latent space.
Appendix F Interpolation
In this section, we present interpolation between two different shapes using our model. For two shapes and , we first compute the mean of the posterior distribution using . Let and be the means of the posterior distribution for and respectively. We use and as the latent representation for these two shapes. We then use the inverse prior flow to transform and back to the prior space. Let and be the corresponding vectors for and in the prior space. We use spherical interpolation between and to retrieve a series of vectors . For each , we use the CNF prior and the CNF decoder to generate the corresponding shape . Figure 6 contains examples of the interpolation.
Appendix G More flow transformation
Figure 7 presents more examples of flow transformations from the Gaussian prior to different shapes.