Mesh-TensorFlow: Deep Learning for Supercomputers

11/05/2018 ∙ by Noam Shazeer, et al. ∙ 8

Batch-splitting (data-parallelism) is the dominant distributed Deep Neural Network (DNN) training strategy, due to its universal applicability and its amenability to Single-Program-Multiple-Data (SPMD) programming. However, batch-splitting suffers from problems including the inability to train very large models (due to memory constraints), high latency, and inefficiency at small batch sizes. All of these can be solved by more general distribution strategies (model-parallelism). Unfortunately, efficient model-parallel algorithms tend to be complicated to discover, describe, and to implement, particularly on large clusters. We introduce Mesh-TensorFlow, a language for specifying a general class of distributed tensor computations. Where data-parallelism can be viewed as splitting tensors and operations along the "batch" dimension, in Mesh-TensorFlow, the user can specify any tensor-dimensions to be split across any dimensions of a multi-dimensional mesh of processors. A Mesh-TensorFlow graph compiles into a SPMD program consisting of parallel operations coupled with collective communication primitives such as Allreduce. We use Mesh-TensorFlow to implement an efficient data-parallel, model-parallel version of the Transformer sequence-to-sequence model. Using TPU meshes of up to 512 cores, we train Transformer models with up to 5 billion parameters, surpassing state of the art results on WMT'14 English-to-French translation task and the one-billion-word language modeling benchmark. Mesh-Tensorflow is available at .



There are no comments yet.


This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

Batch-splitting (data-parallelism) is the dominant distributed Deep Neural Network (DNN) training strategy, due to its universal applicability and its amenability to Single-Program-Multiple-Data (SPMD) programming. However, batch-splitting suffers from several major problems when training very large models. The memory required to store parameters and/or activations and the time necessary to synchronize parameters can make purely-data-parallel algorithms impossible or inefficient. Different distribution strategies (model-parallelism Dean:2012:LSD:2999134.2999271) can solve these issues, but specifying these strategies can be complicated, and the current MIMD implementations generate very large programs which can be difficult to compile and to optimize.

We solve this problem by introducing Mesh-TensorFlow, a language for specifying a general class of distributed tensor computations. Where data-parallelism can be viewed as splitting tensors and operations along the "batch" dimension, in Mesh-TensorFlow, the user can specify any tensor-dimensions to be split across any dimensions of a multi-dimensional mesh of processors. A Mesh-TensorFlow graph compiles into a SPMD program consisting of parallel operations coupled with collective communication primitives such as Allreduce. We use Mesh-TensorFlow to implement an efficient data-parallel, model-parallel version of the Transformer Vaswani17 sequence-to-sequence model. Using TPU meshes of up to 512 cores, we train Transformer models with up to 5 billion parameters, surpassing state-of-the-art results on WMT’14 English-to-French translation task and the one-billion-word Language modeling benchmark.

2 Hardware Assumptions

While much work deals with heterogeneous and/or unreliable hardware, we focus on clusters of identical, reliable processors, each with a local memory. We define a mesh as an n-dimensional array of such processors. The mesh is only a naming abstraction and does not imply a physical network topology. As such, different meshes can be defined over the same set of physical processors. For example, a 512-core TPU cluster with a 16x16x2 toroidal network interconnect could be represented by a 3-dimensional mesh with shape [16, 16, 2], a two-dimensional mesh with shape [32, 16], a one-dimensional mesh with shape [512], etc. The physical network topology does affect performance; particularly important is the performance of MPI Allreduce, grouped by splitting the mesh by a subset of the dimensions, which can be very efficient Patarasuk2009 Jain10y.:optimal if each such group is physically connected.

3 Inspiration: Single-Program-Multiple-Data (SPMD) Batch-Splitting

We first review a commonly-used variant of synchronous data-parallelism where each processor keeps an identical copy of all parameters (Algorithm 1). For each step, the batch of training examples is split into sub-batches, one for each processor. Each processor computes the forward and backward passes on its sub-batch, resulting in gradients on the model parameters. These gradients are then summed across all processors and the results broadcast to all processors (MPI-allreduce). Finally, each processor updates its own copy of the parameters.

1:Compute partial parameter gradients Local computation
3: Local computation
Algorithm 1 Synchronous data-parallelism with replicated parameters. Each processor maintains a complete copy of all weights . The batch of training examples for timestep is partitioned among the set of processors: . Below is the computation performed on one processor .

This algorithm is typically implemented using Single-Program-Multiple-Data (SPMD) programming, with every processor running the same program of local operations and MPI-allreduce primitives.

One way to see this algorithm is that every tensor and every operation in the computation is either split across all processors (if it has a "batch" dimension), or fully replicated across all processors (if it does not have a "batch" dimension). Operations which reduce out the "batch" dimension require an additional MPI-allreduce to produce the correct result. We can describe this as splitting the computation across the "batch" dimension. Mesh-TensorFlow generalizes this idea to splitting computations across arbitrary dimensions.

4 Mesh-TensorFlow: Beyond Batch Splitting

Mesh-Tensorflow generalizes from the batch-splitting algorithm described in section 3 to allow for splitting across different Tensor dimensions. The similarities are as follows:

  • Each tensor in the computation is represented by one (not-necessarily-distinct) slice of the tensor on each processor.

  • Each operation in the computation is implemented as one operation on each processor. Most operations require no communication, with each processor producing its slice of the output from its slices of the inputs. Some operations additionally require collective communication primitives such as MPI-allreduce.

  • Based on the above, the computation can be implemented as a SPMD program.

The new elements in Mesh-TensorFlow are as follows:

  • Tensors have named dimensions. This allows for the idea of a logical dimension (like "batch") which will be split in the same way for different tensors and operations. It is illegal for a tensor to have two identically-named dimensions.

  • Rather than an unstructured set of processors, Mesh-Tensorflow allows for an n-dimensional mesh of processors (section 2). The mesh also has named dimensions.

  • A global "computation layout" is a partial map from tensor-dimension to mesh-dimension specifying which tensor-dimensions are split across which dimensions of the processor-mesh. For example, batch-splitting (data-parallelism) would be expressed by using a one-dimensional mesh with dimension "all_processors" and using the computation layout [("batch", "all_processors")]. This means that all tensors with a "batch" dimension are split along that dimension across all processors, while all other tensors are fully replicated.

5 Tensor Representations

A tensor is represented as one slice of the tensor per processor. The layout of a tensor is an injective partial map from the tensor’s dimensions to dimensions of the mesh, and is computed as the restriction of the global computation layout to that tensor’s dimensions. It is illegal for two dimensions of the same tensor to map to the same mesh dimension. If a tensor’s layout is empty, it is fully replicated on each processor. For every (tensor-dimension, mesh-dimension) pair in the tensor’s layout, the slice on a processor is restricted along that tensor-dimension to a stripe corresponding to that processor’s coordinate along that mesh-dimension. The current implementation of Mesh-TensorFlow requires the size of the tensor-dimension to be evenly divisible by the size of the mesh-dimension.

6 Operation Implementation

Each operation is implemented by parallel computation on every processor, and sometimes collective communication. We describe the implementations of some important operations here:

Component-wise Operations

Mesh-TensorFlow supports component-wise operations where the shapes (and hence the layouts) of the input and output tensors are identical. These are trivially implemented by parallel operations on each processor to compute that processor’s slice of the output from that processor’s slice(s) of the input(s).

Reduction (reduce_sum(), reduce_max(), etc.)

Mesh-TensorFlow supports reductions where the output dimensions are a subset of the input dimensions. These can be implemented by local reductions of each slice, followed by MPI-allreduce across any mesh dimensions corresponding to reduced-out Tensor dimensions. The allreduce operation is necessary because the local reduction only sums across a subset of the split tensor-dimension. Bandwidth-efficient implementations of allreduce exist when the processors for each group are connected in any type of tree. Patarasuk2009 Jain10y.:optimal

Einstein Summation (matrix multiplication, etc.)

Einstein-summation (einsum) notation (as defined in numpy, TensorFlow, etc.) is a way of expressing a class of operations including (batch) matrix multiplication, reductions and broadcasts, where the operation is defined by the names of the dimensions of the input and output tensors. Mesh-TensorFlow’s use of named dimensions makes using einsum particularly convenient. Einsum can be defined as broadcasting all inputs to a shape consisting the union of all their dimensions, multiplying them component-wise, then reducing out all dimensions not in the specified output shape. Einsum is implemented by parallel einsum operations on each processor of that processor’s input slices, followed by MPI-allreduce across any mesh dimensions corresponding to reduced-out Tensor dimensions.

6.1 Reshape

While reshape is simple in the non-distributed case, Mesh-TensorFlow reshape can require network communication, since the layout of the output tensor may differ from that of the input tensor. Even keeping the same dimension sizes, changing the dimension names (and hence the layout) can result in several different communication patterns: If a dimension is split in the input but not in the output, the implementation involves MPI-allgather communication across the corresponding mesh-dimension. If a dimension is split in the output but not in the input, the implementation involves no communication, just slicing on each processor. MPI-alltoall is used in the case where different dimensions in the input and the output are split across the same mesh dimension, as might be the case when switching between data-parallelism and model-parallelism for different layers of the same model, as in Shazeer17.

7 Mesh-TensorFlow syntax

The Mesh-TensorFlow language is nearly identical to TensorFlow tensorflow2015-whitepaper, with the familiar notions of graphs, tensors, operations, variables, devices (called meshes), and automatic gradient computation. The principal difference is that in Mesh-TensorFlow, tensor-dimensions have a name as well as a size. The shape of each tensor is a statically-known tuple of such dimensions. Shapes are inferred automatically when possible, as they are in TensorFlow. Binary component-wise operations like addition employ implicit broadcasting in the case where the shape of one operand is a subset of the shape of the other.

The initial implementation of Mesh-TensorFlow is a Python library. The user builds a Mesh-TensorFlow graph in python, which the library "lowers" to generate part of a TensorFlow graph. As of the writing of this paper, implementations exist for generating SPMD TensorFlow code for TPUs, or MIMD code (using device placement) for multi-CPU/GPU configurations.

8 Example: Two Fully-Connected Layers

We consider a simple example of two fully-connected layers in the middle of a neural network. The input layer and the output layer each have units, and the hidden layer has units. The hidden layer also has a bias and activation.


This Mesh-TensorFlow code fragment runs these layers on a batch of inputs.

... batch = mtf.Dimension("batch", b) io = mtf.Dimension("io", d_io) hidden = mtf.Dimension("hidden", d_h) # x.shape == [batch, io] w = mtf.get_variable("w", shape=[io, hidden]) bias = mtf.get_variable("bias", shape=[hidden]) v = mtf.get_variable("v", shape=[hidden, io]) h = mtf.relu(mtf.einsum(x, w, output_shape=[batch, hidden]) + bias) y = mtf.einsum(h, v, output_shape=[batch, io]) ...

The code above defines only the mathematical model. We now discuss several different computation layouts. Each will produce identical results, but will have different performance characteristics. We also provide illustrations of the layouts in Appendix A.

8.1 Data-Parallel Layout

To train the above model in data-parallel mode on a mesh of processors, we would define:

  mesh_shape = [("all", n)]
  computation_layout = [("batch", "all")]

When the Mesh-TensorFlow graph is compiled with this layout, the parameter tensors , , and are replicated on all processors, but the activation matrices , , , etc. are split across the batch dimension. For example, each processor keeps a slice of with shape .

There is no inter-processor communication in the forward pass. However, the gradient computations for the parameters are mtf.einsum operations which reduce out the batch dimension, and hence produce operations when they are compiled. The number of values allreduced per processor is equal to the number of parameters, approximately .

8.2 Model-Parallel Layout

Rather than splitting the batch, we can split the units in the hidden layer:

  mesh_shape = [("all", n)]
  computation_layout = [("hidden", "all")]

When the Mesh-TensorFlow graph is compiled with this layout, the input and output layers , and are replicated on all processors, but the hidden activations and the parameter tensors , and are all split across the hidden dimension. For example, each processor keeps a slice of with shape and a slice of with shape .

When computing , the split hidden dimension is reduced out. Consequently, the results of that computation get allreduced across all processors. A similar allreduce happens in computing the gradients on . In all, the number of values allreduced per processor is .

8.3 Data-Parallel, Model-Parallel Layouts

On a two-dimensional mesh of processors, we can employ both data-parallelism and model-parallelism:

  mesh_shape = [("rows", r), ("cols", c)]
  computation_layout = [("batch", "rows"), ("hidden", "cols")]

In this layout, each row of processors handles a fraction of the batch, while each column of processors handles a fraction of the hidden units. Each processor keeps a slice of x with shape , with processors in the same row having identical slices. The hidden activation tensor is tiled in two dimensions, with each processor keeping a slice with shape .

This layout causes partitioned-allreduce operations in several places. For example, in computing , we reduce out the hidden dimension, which is split over the cols dimension of the mesh, so the results of the operation need to be summed up by processor-column, as opposed to over the entire mesh. In all, the number of values allreduced per processor is

If we have a three-dimensional mesh of processors, we can even split the computation in three dimensions:

  mesh_shape = [("rows", r), ("cols", c), ("planes", p)]
  computation_layout = [
    ("batch", "rows"), ("hidden", "cols"), ("io", "planes"])

In this case, every matrix in the computation is tiled across two mesh dimensions and replicated in the third, and every einsum requires an allreduce across one mesh dimension.

8.4 Inefficient Layouts

For a computation layout to be efficient, all expensive operations need to be split (as opposed to replicated) across all mesh dimensions. For example, the empty layout below produces correct results, but since it replicates all computation on every processor, it saves no time or memory. A general rule is that any expensive einsum operation should have one input dimension that is split across each batch dimension.

  mesh_shape = [("all", n)]
  computation_layout = []

8.5 Illegal Layouts

The computation layout below is illegal, because it causes the tensor to have two dimensions which are split across the same dimension of the mesh.

  mesh_shape = [("all", n)]
  computation_layout = [("batch", "all"), ("hidden", "all")]

8.6 Performance Comparison

Layout Comp. Time Comm. Time Memory/Processor
[("batch", "all")]
[("hidden", "all")]
[("batch", "rows"),
("hidden", "cols")]
[("batch", "rows"),
("hidden", "cols"),
("io", "planes")]

Table 1: Computation, communication and memory costs for different layouts of the computation in Algorithm 1. Constant factors and lower-order terms are dropped.

Table 1 shows the computational costs associated with our example computation layouts. The computation time is dominated by that of einsum operations. The communication time comes from the Allreduce operations, which are necessary whenever the inner dimension of einsum is split. Assuming that the mesh has physical links between all pairs of logically adjacent processors, each Allreduce operations can be done in time proportional to the size of one slice divided by the per-link network bandwidth Jain10y.:optimal.

The network-boundedness of the computation is proportional to the value shown in the table column marked , with the constant of proportionality depending on the ratio of communication and computation speeds on the given hardware. In the data-parallel layout, the value is , the inverse of the per-processor batch size. Performance suffers if the per-processor batch is too small. In the model-parallel layout, the value is , the inverse of the number of hidden units per processor. Performance suffers if the hidden layer is sliced too finely. For good performance, batch size is irrelevant, but we need the hidden layer to get larger as we increase the number of processors. In the first data-parallel, model-parallel layout, the value is . In this layout, we can quadratically increase the number of processors while only linearly increasing the batch size and hidden layer sizes necessary to maintain good efficiency. The final layout lets us cubically increase the number of processors in a 3-dimensional mesh, while only linearly increasing the batch size and the layer sizes.

9 Model-Parallel "Transformer"

We implemented a model-parallel layout of the Transformer attention-based sequence-to-sequence model described in Vaswani17. The complete implementation is available in the tensor2tensor library on github. The layout is given by:

  mesh_shape = [("all", n)]
  computation_layout = [
    ("vocab", "all"), ("d_ff", "all"), ("heads", "all")]

That is, the dimensions representing the vocabulary size, the size of the feed-forward hidden layer, and the number of attention heads are each split across all processors. This layout works because every expensive operation in the model has exactly one of these dimensions, and no tensor in the model has more than one. Similarly to the model-parallel layout for our example network (Section 8.2), network-boundedness and memory usage per processor remain constant if we scale all of these dimensions proportionally to the number of processors. We did just this, training transformer models with ever larger hidden layers and numbers of attention heads on ever larger TPU clusters (we did not increase the vocabulary size). As expected, we saw very similar performance characteristics between the models. This scaling turns out to be highly beneficial to model quality (Section 9.1).

To use even more processors, we combined this model-parallelism with data parallelism, splitting the batch across one dimension of a 2-dimensional TPU mesh and the dimensions described above across the other dimension of the mesh:

  mesh_shape = [("rows", r), ("cols", c")]
  computation_layout = [("batch", "rows"), ("vocab", "cols"),
                        ("d_ff", "cols"), ("heads", "cols")]

This layout maintains constant performance if the batch size is scaled proportionally to r and the mentioned model dimensions are scaled proportionally to c. Using this layout, we trained Transformer models with feed-forward hidden dimensions up to 262144 and up to 256 attention heads on 2-dimensional TPUv2 meshes of up to 16x32=512 cores, maintaining computational efficiency of over 50% (6 PFLOP/s out of a maximum 11.5 PFLOP/s) on the largest models.

9.1 Experiments and Results

To examine the benefit of scaling the Transformer model in the manner suggested by the previous section, we trained such models on machine translation and language modeling tasks. Results are given in Tables 2 and 3.

For the billion-word language modeling benchmark, we trained the models for 10 epochs. The largest model (4.9B parameters) took 13 hours to train on a 512-core TPUv2 cluster. Batch size for all models was 256 sequences of 256 tokens each (each sequence was the concatenation of multiple training sentences). The batch was split along the mesh dimension of size 16 and the model dimensions were split along the mesh dimension of size 32. Per-word dev-perplexity for the largest model was 24.0, but dropped to 23.5 when the model was evaluated with the logits multiplied by 0.9 (likely due to overfitting). This represents the best published result on this dataset. As expected, perplexity was lower for larger models. We have included random samples from these models in Appendix

C. On the languagemodel_wiki_noref_v128k_l1k dataset from the Tensor2Tensor library111No published results exist for this dataset., consisting of over 5 billion tokens of text from Wikipedia, perplexity continued to improve significantly with a model size of 5 billion parameters.

On the WMT14 En-Fr translation tasks (3), we trained the models for 3 epochs. The largest model (2.9B parameters) was trained for 22 hours on a 128-core TPUv2 cluster. Quality improved with model size, with the largest model achieved BLEU score 43.9 (evaluated using sacrebleu), the best published result to date. For the WMT14 En-De dataset, gains from model size were smaller, presumably due to the small size of the training data.

Additional details about the configurations for these experiments are available as part of the tensor2tensor library on github.

Parameters Billion-Word Benchmark Wikipedia
(Billions) Word-Perplexity Subword-Perplexity
4096 4 0.14 35.0 8.74
8192 8 0.22 31.7 8.03
16384 16 0.37 28.9 7.44
32768 32 0.67 26.8 6.99
65516 64 1.28 25.1 6.55
131072 128 2.48 24.1 6.24
262144 256 4.90 24.0(23.5) 6.01
Prev Best DNN Shazeer17  6.5 28.0
Best DNN Ensemble Jozefowicz16 26.1
Best Ensemble (different methods)Jozefowicz16 23.7
Table 2: Transformer-Decoder Language Models: ,
Parameters WMT14 EN-DE WMT14 EN-FR
(Billions) BLEU BLEU
2048 4 128 0.15 25.5 41.8
4096 8 128 0.24 26.5 42.5
8192 16 128 0.42 27.1 43.3
16384 32 128 0.77 27.5 43.5
32768 64 128 1.48 27.5 43.8
65536 128 128 2.89 26.7 43.9
4096 16 64 0.21 28.4 41.8 Vaswani17
Table 3: Transformer Machine-Translation Results. ,

10 Related Work

A large part of deep learning computations is a series of matrix multiplications and tensor contractions (Einsums). Distributed matrix multiplication is a well-studied problem in high performance computing. Efficient algorithms partition the computational space, instead of partitioning work by the output matrix/tensor (owners compute), to minimize communication. This technique is sometimes called iteration space tiling Wolfe:1989:MIS:76263.76337, replication SD_EUROPAR_2011, or task parallelism Calvin:2015:STA:2833179.2833186. Mesh-TensorFlow can express a wide range of uniform partitionings of the iteration space and therefore can adopt many best known mappings, e.g., 3D Aggarwal:1990:CCP:77831.77836; berntsen1989communication and 2.5D SD_EUROPAR_2011 algorithms for square matrices, CARMA demmel2013communication for rectangular matrices, 1.5D SpDM3 algorithm for matrices with different sparsities, best tile sizes for direct convolutions demmel2018communication, etc., although sometimes with higher memory requirements. Furthermore, in most existing work, when multiple multiplications are composed together, the user has to specify the data layout for each matrix separately pmlr-v84-koanantakool18a. Mesh-TensorFlow lets the user name the dimension to split, simplifying the process and allowing for much easier mapping explorations. Feature-wise, Mesh-TensorFlow shares many similarities with the Cyclops Tensor Framework solomonik2014massively, a distributed tensor contraction library originally developed for quantum chemistry applications, which also supports replication and arbitrary mappings.

In the context of deep learning, partitioning the iteration space, e.g., interpolating between data and model parallelism, is relatively new. Gholami et al. 

gholami2017integrated analytically showed that using both data and model parallelism at the same time can be more beneficial than using just one of them. Building on top of 1.5D matrix multiplication algorithms, their algorithm can support replication and arbitrary processor grid shapes. However, they only explored the parallelization of AlexNet krizhevsky2012imagenet and they have not implemented the algorithm. Jia et al. jia2018exploring; jia2018beyond implemented a framework that uses cost modeling to pick the best parallelization strategy, including how to partition work for each operation. Their parallelizable dimensions are defined as the set of all divisible dimensions in the output tensor (owners compute), and therefore their mapping can be suboptimal in terms of communication. We expand on this in Appendix B.

11 Future Work

The Mesh-TensorFlow library is available at and is under active development. Some potential areas for development are:

  • Automated search for optimal computation layout.

  • Implementations of different models and operations. For example, convolutions on spatially-partitioned tensors will require the communication of "halo" regions, as described in Jin18.

  • Implementation of SPMD programming on CPU/GPU clusters.

12 Conclusion

In this paper, we introduce the Mesh-TensorFlow language, facilitating a broad class of SPMD distributed tensor computations. Applying Mesh-TensorFlow to the Transformer model, we are able to train models with 5 billion parameters on up to 512-core clusters, establishing new state-of-the-art results for WMT14 En-Fr translation task and the One Billion Word language modeling benchmark.

Appendix A Illustrations for the Two Fully-Connected Layers Example

This section provides the illustrations of the four layouts mentioned in the Two Fully-Connected Layers example in Section 8. The overall computation is shown in Figure 1. We draw a matrix multiplication by putting to the left of and putting above . For each matrix, we put its name inside, its number of rows on the left or right side, and its number of columns above or below it. We omit the numbers of rows or columns that can be implied from adjacent matrices, i.e., knowing that the multiplication dimensions must match.

Figure 1: The overall computation of the Two Fully-Connected Layers example. First, is multiplied with . The temporary result is summed with , applied component-wise , and then stored in . Lastly, we multiply with to get . The component-wise operations and the intermediate matrix are not shown in the figure.

Figure 2 presents the purely data-parallel layout with processors. We number the processors 0 and 1, respectively. The ranks of the processors that store each matrix part are written in blue on the matrix part. The matrix names are moved to the bottom-left corners. The whole and are labeled with both 0 and 1 because both of them are fully replicated between the two processors. The purely model-parallel layout are drawn similarly in Figure 3.

Figures 4 and 5 show the mixed data-and-model-parallel layout with a 2-by-2 and a 2-by-2-by-2 processor meshes, respectively. We give each processor a serialized rank as shown in the figure, and use the serialized rank to label matrix slices.

Figure 2: The data-parallel layout for the Two Fully-Connected Layers example, with processors . Blue numbers on matrices indicate the ranks of the processors the matrix slices reside on. The dimension is split among all processors. and are fully replicated.
Figure 3: The model-parallel layout for the Two Fully-Connected Layers example, with processors . Blue numbers on matrices indicate the ranks of the processors the matrix slices reside on. The dimension is split among all processors. and are fully replicated.
Figure 4: The mixed data-and-model-parallel layout for the Two Fully-Connected Layers example. There are 4 processors, arranged into a 2-by-2 mesh. Each processor is assigned a serialized rank which is used to label matrix slices that it owns.
Figure 5: The mixed data-and-model-parallel layout for the Two Fully-Connected Layers example. There are 8 processors, arranged into a 2-by-2-by-2 mesh. Each processor is assigned a serialized rank which is used to label matrix slices that it owns.

Appendix B Per-operation Parallelizations and Communication Costs

Communication is much more expensive than computation and is usually the bottleneck in a parallel program, especially in distributed setting. The section shows how the more common owner-compute parallelization strategies can be communication-suboptimal. We start with a simplified overview of the parallelization schemes used in distributed matrix multiplication and tensor contractions (Einsums), focusing on their communication bandwidth costs. (See irony2004communication; ballard2011minimizing; hbl for rigorous communication lower bounds analyses.) We only discuss distributed matrix multiplication here since the concept is trivially generalizable to its tensor counterpart.

Iteration Space.

The iteration space is the set of all index tuples required to compute a problem. For example, the matrix multiplication problem computes . Its iteration space consists of all possible tuples , where is -by-, is -by-, and is -by-, as shown in Figure 6.

Figure 6: The iteration space of the matrix multiplication problem . Each voxel represents the computation .


Let be the number of processors. Parallelization corresponds to partitioning the set of voxels into (not necessarily) equal subsets for each processor to compute. For matrix multiplication, the most widely-used partitionings are grouped into three categories ballard2013communication: 1D, 2D, and 3D, based on the number of dimensions of the iteration space that are split. Figure 7 shows an example for each category. The left image is 1D partitioning () because only the axis is split. The middle image splits axes and so it is 2D partitioning (). The right image is 3D partitioning () because all three axes are split.

Figure 7: Example partitionings of the iteration space for the matrix multiplication problem. The names 1D, 2D, and 3D comes from the number of dimensions that are split.

Owner computes.

Owner-compute strategies split a matrix (or matrices) equally among processors and each processor is responsible for all computations related to the matrix chunk it owns. The 1D and 2D partitionings in Figure 7 are owner-compute strategies. In 1D case, each processor owns a slice of matrices and each, and computes the whole slab requires for its slice of . In 2D case, each processor owns a patch of and computes a pencil corresponding to it. The 3D partitioning goes beyond owner-compute rule, since no processor is responsible for all computations associated with the data it has. Owner-compute schemes are more common because they are the most intuitive as we often view the output data as the unit of work. We will show why its communication costs are usually suboptimal in the next paragraph.


Here, we focus on the number of elements that have to be transferred by a processor. Let be the voxel subset assigned to a processor, and , and be the projections of onto the , , and planes, respectively. The total number of elements a processor has to access to complete its computation is simply , where denotes set cardinality. Since a processor can only hold a limited amount of data in memory, the rest of the elements must come through communication. The volume of the subset designates the computational workload. As mentioned in the paper, we would like to maximize the computation-to-communication ratio, therefore we want to have as low surface-to-volume ratio as possible. Assuming only takes cuboid shapes, then the best shape is a cube.

Owner-compute methods fall short when it cannot partition the space into cubes. To illustrate, we compare 2D and 3D partitionings for processors in Figure 7. When , each pencil in the 2D partitioning has a computation-to-communication ratio,

Each cube in the 3D partitioning has a higher computation-to-communication ratio,


Mesh-TensorFlow can express more fine-grained parallelism than owner-compute, even though all we have to do is just specifying the data layout for each tensor. That is because our layout allows the tensor to be replicated. This, combines with multiple layouts from multiple tensors involved in an operation, can split the iteration space in as many dimensions as necessary (up to the rank of the iteration space).

Appendix C Random Samples from Transformer Language Models Trained on Billion-Word Language Modeling Benchmark

The following samples were randomly generated from the Transformer language models described in the paper. All sentences were seeded with the initial words "According to Ray Kurzweil", and continued randomly by the model. While all the models produce mostly grammatical sentences, the larger models exhibit more world knowledge.

c.1 Model with 0.14B Parameters, PPL=35.0

According to Ray Kurzweil …

  • … , professor of musical economics at the University of California , Berkeley , the military ’s " biggest challenge " might have been growing the region ’s growing urban sprawl .

  • … , the chief executive of ITC Emerging Markets , the metals sector should contribute about $ 210 billion by year-end 2008 – a figure that has now risen to $ 3 billion .

  • … , the head of PMI , he has now written off all his employees fairly .

  • … , senior vice president of sales and marketing at Aspen , there is a percentage of the online response for the auto company escalating to 18 customers and increasing to 30 days a year .

  • … , director of the Economic Policy Institute , India ranks 21st out of 221 nations in the current budget deficit , which ranks 13th globally .

  • … , Victor Hesketh ’s food agency that produces , sends and distributes the raw beef and pork through the state ’s Food and Drug Administration , and they ’re also allowed to take the lamb by their own account alongside the chicken and veal .

  • … , the author of " Smoking : The Letters of Hysteria and Reclining the State of South Carolina " ( 2007 ) , 30 percent of liquor ’s coconut is sold on the first batch .

  • … , an MIT student who is not involved in anything more than a stock-market move , the latest system of debt and bankruptcy parallels the completely unregulated collection of debt that emerged in the early 1990s , when a financial industry boom brought a financial crisis that spilled almost everywhere from the United States .

c.2 Model with 0.37B Parameters, PPL=28.9

According to Ray Kurzweil …

  • … , the owner and web guru at Stanford University , even a single person might fall in love with the internet simultaneously .

  • … , PhD , a professor of computer science at Princeton who has staked his reputation on the ability of code technicians to simultaneously amplify a cursor ’s legal weight , machines will go digital by using meta-dimensions instead of a brick , and design schemes closer to the core of GPUs .

  • … , chief executive of the company , very few people work through chronology , and most people can ’t use tables on their machines .

  • … , we are saving the most jobs in the world .

  • … , creator of the Star Wars creator and creator of the popular game , " Indiana Jones and the Kingdom of the Crystal Skull , " here comes Martin Schafer ( " Spider-Man 3 " ) to describe the businessman as a grand master .

  • … , a technology expert at the Massachusetts Institute of Technology and a host of from academia , Ardipithecus ramidus was frequently incubated with solar-powered cells to frighten away and frighten out the predators , and housed in them so they could damage Lego ’s control panels .

  • … , the famously headstrong and egocentric stack-on-stack , a keyboard is more than a part of the laptop ’s muscle when it is standing upright rather than crouching , and its outlet at the right end of the screen was probably the key to how the phone turned into its belly .

  • … , founder of the Stanford-funded company , similar " recycled book " concepts are also more of an obsession , as this lets authors track their emotions in print product forms and refuse to look at one all- , -certain identity .

c.3 Model with 1.28B Parameters, PPL=25.1

According to Ray Kurzweil …

  • … , a transplant surgeon and Harvard astrophysicist , the results are illuminated in Honolulu , Vietnam .

  • … of UNK , a software company in California , computer games are alive and well , with millions of gamers , but most games do not make money .

  • … , the James Watson and James Watson Professor of Physics at MIT , if we all assume that the project will continue to go ahead , we will eventually be in an era when human beings are capable of figuring out just what our genetic make-up is focused on .

  • … , the physicist who has kept many of these principles alive , the world has vanished into " more and more of the static world " – and in many ways we ’re losing the earth ’s ability to appreciate water .

  • … , creator and the only original idea from the series , the tablet is expected to be a device that combines a USB 2.0 card with an iPad , a laptop and an iPod Touch – with UNK sharp-wired connections to the Internet .

  • … and a panel of experts testifying in Los Angeles , six years in Congress attempts to improve " brain " of Americans by hitting them with a $ 50 annual fee , then having them likely pay the company $ 3,000 for every additional year in life they are in a job .

  • … , creator of the Review of Medical UNK , the organisation could be the " holy grail " of degenerative brain disease .

  • … , music analyst and co-founder of zero-carbon quantum computing firm Redpoint , if you listen carefully , you ’ll hear a far more universal " idea " of how this new car is supposed to be or what it should not be , or how much going to cost .

c.4 Model with 4.9B Parameters, PPL=24.0

According to Ray Kurzweil …

  • … , chief technology officer for the US Department of Energy , aviation has " potential to be the largest and fastest growing source of consumer and commercial emissions . "

  • … , the futurist son of the futurist who wrote The Singularity is Near , the " early days " of Google are not the right time to push forward with the next great leap .

  • … , creator of the first modern computer , the sexy cyborg was the brainchild of an MIT professor , Thomas Harris , and a former banking entrepreneur , Henry Lee , who was looking for an UNK , a light that could be recovered and used to light up the Internet .

  • … , the inventor of the modern personal computer , the shrinking human brain could eventually replace the Internet as a tool of human intelligence and imagination .

  • … , the expert and co-author of " The Singularity is Near : Comprehending the Technological Future of Engineering , " people are looking for ways to protect and make their lives better .

  • … , creator of the Gaia hypothesis , earlier computer systems should become not just more efficient , but more-efficient , increasing their efficiency by reducing human errors ( the unexpected , but often the regrettable ) and even the number of errors .

  • … , who will make an appearance at this year ’s Consumer Electronics Show in Las Vegas next week , these mobile gadgets will be able to " talk " to each other .

  • … , the futurist turned futurist , the onset of Alzheimer ’s coincided precisely with the rate of unemployment in America .