from d2l import mxnet as d2l
from mxnet import autograd, gluon, init, np, npx
from mxnet.gluon import nn
npx.set_np()The previous section did data-parallel training the hard way — manual all_reduce, manual replica management. In practice, every framework wraps it in a one-liner:
nn.DataParallel(net) (multi-GPU on one host) or nn.parallel.DistributedDataParallel (multi-host).gluon.Trainer(..., kvstore='device').tf.distribute.MirroredStrategy().Same numerical result; orders of magnitude less boilerplate; NCCL all-reduce under the hood.
We use a small ResNet for these experiments — the speedup from data parallelism only matters once the per-GPU compute is non-trivial:
def resnet18(num_classes):
"""A slightly modified ResNet-18 model."""
def resnet_block(num_channels, num_residuals, first_block=False):
blk = nn.Sequential()
for i in range(num_residuals):
if i == 0 and not first_block:
blk.add(d2l.Residual(
num_channels, use_1x1conv=True, strides=2))
else:
blk.add(d2l.Residual(num_channels))
return blk
net = nn.Sequential()
# This model uses a smaller convolution kernel, stride, and padding and
# removes the max-pooling layer
net.add(nn.Conv2D(64, kernel_size=3, strides=1, padding=1),
nn.BatchNorm(), nn.Activation('relu'))
net.add(resnet_block(64, 2, first_block=True),
resnet_block(128, 2),
resnet_block(256, 2),
resnet_block(512, 2))
net.add(nn.GlobalAvgPool2D(), nn.Dense(num_classes))
return netWrap the model in the framework’s data-parallel container. Parameters are replicated to each GPU automatically:
The wrapper also handles inference — splits the input minibatch across replicas, gathers outputs:
def evaluate_accuracy_gpus(net, data_iter, split_f=d2l.split_batch):
"""Compute the accuracy for a model on a dataset using multiple GPUs."""
# Query the list of devices
devices = list(net.collect_params().values())[0].list_ctx()
# No. of correct predictions, no. of predictions
metric = d2l.Accumulator(2)
for features, labels in data_iter:
X_shards, y_shards = split_f(features, labels, devices)
# Run in parallel
pred_shards = [net(X_shard) for X_shard in X_shards]
metric.add(sum(float(d2l.accuracy(pred_shard, y_shard)) for
pred_shard, y_shard in zip(
pred_shards, y_shards)), labels.size)
return metric[0] / metric[1]The loop looks like ordinary single-GPU training because the wrapper owns the distributed work:
The important lesson is the interface: after wrapping the model, most training code should not need to know how many GPUs are present.
Use this as the throughput baseline before the data-parallel wrapper adds replication and gradient averaging.
The training loop is unchanged; the wrapper splits the minibatch and synchronizes gradients under the hood.
DataParallel, MirroredStrategy) reduce data-parallel SGD to one line of setup.DistributedDataParallel / MultiWorkerMirroredStrategy — same idea, NCCL/Gloo across the network.