arXiv:1702.07825v2 [cs.CL] 7 Mar 2017

Deep Voice: Real-time Neural Text-to-Speech

Sercan Ö. Arık†
SERCANARIK @ BAIDU . COM
Mike Chrzanowski†
MIKECHRZANOWSKI @ BAIDU . COM
Adam Coates†
ADAMCOATES @ BAIDU . COM
Gregory Diamos†
GREGDIAMOS @ BAIDU . COM
Andrew Gibiansky†
GIBIANSKYANDREW @ BAIDU . COM
Yongguo Kang†
KANGYONGGUO @ BAIDU . COM
Xian Li†
LIXIAN 05@ BAIDU . COM
John Miller†
MILLERJOHN @ BAIDU . COM
Andrew Ng†
ANDREWNG @ BAIDU . COM
Jonathan Raiman†
JONATHANRAIMAN @ BAIDU . COM
Shubho Sengupta†
SSENGUPTA @ BAIDU . COM
MOHAMMAD @ BAIDU . COM
Mohammad Shoeybi†
Baidu Silicon Valley Artificial Intelligence Lab, 1195 Bordeaux Dr. Sunnyvale, CA 94089

Abstract
We present Deep Voice, a production-quality
text-to-speech system constructed entirely from
deep neural networks. Deep Voice lays the
groundwork for truly end-to-end neural speech
synthesis.
The system comprises five major building blocks: a segmentation model for
locating phoneme boundaries, a grapheme-tophoneme conversion model, a phoneme duration
prediction model, a fundamental frequency prediction model, and an audio synthesis model.
For the segmentation model, we propose a novel
way of performing phoneme boundary detection
with deep neural networks using connectionist
temporal classification (CTC) loss. For the audio synthesis model, we implement a variant
of WaveNet that requires fewer parameters and
trains faster than the original. By using a neural network for each component, our system is
simpler and more flexible than traditional text-tospeech systems, where each component requires
laborious feature engineering and extensive domain expertise. Finally, we show that inference
with our system can be performed faster than real
time and describe optimized WaveNet inference
kernels on both CPU and GPU that achieve up to
400x speedups over existing implementations.
†

Authors are listed alphabetically by last name.

Submitted February 24, 2017 for review for the International Conference on Machine Learning (ICML) 2017. Copyright 2017 by
the author(s).

1. Introduction
Synthesizing artificial human speech from text, commonly
known as text-to-speech (TTS), is an essential component
in many applications such as speech-enabled devices, navigation systems, and accessibility for the visually-impaired.
Fundamentally, it allows human-technology interaction
without requiring visual interfaces. Modern TTS systems
are based on complex, multi-stage processing pipelines,
each of which may rely on hand-engineered features and
heuristics. Due to this complexity, developing new TTS
systems can be very labor intensive and difficult.
Deep Voice is inspired by traditional text-to-speech
pipelines and adopts the same structure, while replacing all
components with neural networks and using simpler features: first we convert text to phoneme and then use an
audio synthesis model to convert linguistic features into
speech (Taylor, 2009). Unlike prior work (which uses
hand-engineered features such as spectral envelope, spectral parameters, aperiodic parameters, etc.), our only features are phonemes with stress annotations, phoneme durations, and fundamental frequency (F0). This choice of
features makes our system more readily applicable to new
datasets, voices, and domains without any manual data annotation or additional feature engineering. We demonstrate
this claim by retraining our entire pipeline without any hyperparameter changes on an entirely new dataset that contains solely audio and unaligned textual transcriptions and
generating relatively high quality speech. In a conventional
TTS system this adaptation requires days to weeks of tuning, whereas Deep Voice allows you to do it in only a few
hours of manual effort and the time it takes models to train.

Deep Voice: Real-time Neural TTS

Real-time inference is a requirement for a productionquality TTS system; without it, the system is unusable for
most applications of TTS. Prior work has demonstrated that
a WaveNet (van den Oord et al., 2016) can generate close to
human-level speech. However, WaveNet inference poses a
daunting computational problem due to the high-frequency,
autoregressive nature of the model, and it has been hitherto
unknown whether such models can be used in a production system. We answer this question in the affirmative and
demonstrate efficient, faster-than-real-time WaveNet inference kernels that produce high-quality 16 kHz audio and
realize a 400X speedup over previous WaveNet inference
implementations (Paine et al., 2016).

2. Related Work
Previous work uses neural networks as substitutes for
several TTS system components, including grapheme-tophoneme conversion models (Rao et al., 2015; Yao &
Zweig, 2015), phoneme duration prediction models (Zen
& Sak, 2015), fundamental frequency prediction models
(Pascual & Bonafonte, 2016; Ronanki et al., 2016), and
audio synthesis models (van den Oord et al., 2016; Mehri
et al., 2016). Unlike Deep Voice, however, none of these
systems solve the entire problem of TTS and many of them
use specialized hand-engineered features developed specifically for their domain.
Most recently, there has been a lot of work in parametric audio synthesis, notably WaveNet, SampleRNN, and
Char2Wav (van den Oord et al., 2016; Mehri et al., 2016;
Sotelo et al., 2017). While WaveNet can be used for
both conditional and unconditional audio generation, SampleRNN is only used for unconditional audio generation.
Char2Wav extends SampleRNN with an attention-based
phoneme duration model and the equivalent of an F0 prediction model, effectively providing local conditioning information to a SampleRNN-based vocoder.
Deep Voice differs from these systems in several key aspects that notably increase the scope of the problem. First,
Deep Voice is completely standalone; training a new Deep
Voice system does not require a pre-existing TTS system,
and can be done from scratch using a dataset of short audio clips and corresponding textual transcripts. In contrast,
reproducing either of the aforementioned systems requires
access and understanding of a pre-existing TTS system, because they use features from another TTS system either at
training or inference time.
Second, Deep Voice minimizes the use of hand-engineered
features; it uses one-hot encoded characters for grapheme
to phoneme conversion, one-hot encoded phonemes and
stresses, phoneme durations in milliseconds, and normalized log fundamental frequency that can be computed from

waveforms using any F0 estimation algorithm. All of these
can easily be obtained from audio and transcripts with minimal effort. In contrast, prior works use a much more complex feature representation, that effectively makes reproducing the system impossible without a pre-existing TTS
system. WaveNet uses several features from a TTS system
(Zen et al., 2013), that include values such as the number
of syllables in a word, position of syllables in the phrase,
position of the current frame in the phoneme, and dynamic
features of the speech spectrum like spectral and excitation
parameters, as well as their time derivatives. Char2Wav
relies on vocoder features from the WORLD TTS system
(Morise et al., 2016) for pre-training their alignment module which include F0, spectral envelope, and aperiodic parameters.
Finally, we focus on creating a production-ready system,
which requires that our models run in real-time for inference. Deep Voice can synthesize audio in fractions of a
second, and offers a tunable trade-off between synthesis
speed and audio quality. In contrast, previous results with
WaveNet require several minutes of runtime to synthesize
one second of audio. We are unaware of similar benchmarks for SampleRNN, but the 3-tier architecture as described in the original publication requires approximately
4-5X as much compute during inference as our largest
WaveNet models, so running the model in real-time may
prove challenging.

3. TTS System Components
As shown in Fig. 1, the TTS system consists of five major
building blocks:
• The grapheme-to-phoneme model converts from
written text (English characters) to phonemes (encoded using a phonemic alphabet such as ARPABET).
• The segmentation model locates phoneme boundaries in the voice dataset. Given an audio file and a
phoneme-by-phoneme transcription of the audio, the
segmentation model identifies where in the audio each
phoneme begins and ends.
• The phoneme duration model predicts the temporal
duration of every phoneme in a phoneme sequence (an
utterance).
• The fundamental frequency model predicts whether
a phoneme is voiced. If it is, the model predicts the fundamental frequency (F0) throughout the
phoneme’s duration.
• The audio synthesis model combines the outputs
of the grapheme-to-phoneme, phoneme duration, and
fundamental frequency prediction models and synthesizes audio at a high sampling rate, corresponding to

Deep Voice: Real-time Neural TTS
(a)

Training

Audio

Text

Segmentation

Grapheme-to-Phoneme
Phoneme
Dictionary

Durations

Audio Synthesis
Duration
Prediction

Phonemes
Audio

Fundamental Frequency
(F0) Prediction

F0 Extraction

(b)

Inference
Duration Prediction

Du

ra

tio

ns

Text

Grapheme-to-Phoneme
Phoneme
Dictionary

Phonemes

Audio Synthesis

Fundamental Frequency
(F0) Prediction

Audio

ile

rof

P
F0

Figure 1. System diagram depicting (a) training procedure and (b) inference procedure, with inputs on the left and outputs on the right.
In our system, the duration prediction model and the F0 prediction model are performed by a single neural network trained with a joint
loss. The grapheme-to-phoneme model is used as a fallback for words that are not present in a phoneme dictionary, such as CMUDict.
Dotted lines denote non-learned components.

the desired text.

During inference, text is fed through the grapheme-tophoneme model or a phoneme dictionary to generate
phonemes. Next, the phonemes are provided as inputs to
the phoneme duration model and F0 prediction model to
assign durations to each phoneme and generate an F0 contour. Finally, the phonemes, phoneme durations, and F0
are used as local conditioning input features to the audio
synthesis model, which generates the final utterance.
Unlike the other models, the segmentation model is not
used during inference. Instead, it is used to annotate
the training voice data with phoneme boundaries. The
phoneme boundaries imply durations, which can be used
to train the phoneme duration model. The audio, annotated with phonemes and phoneme durations as well as
fundamental frequency, is used to train the audio synthesis model.
In the following sections, we describe all the building
blocks in detail.
3.1. Grapheme-to-Phoneme Model
Our grapheme-to-phoneme model is based on the encoderdecoder architecture developed by (Yao & Zweig, 2015).
However, we use a multi-layer bidirectional encoder with
a gated recurrent unit (GRU) nonlinearity and an equally
deep unidirectional GRU decoder (Chung et al., 2014). The
initial state of every decoder layer is initialized to the final

hidden state of the corresponding encoder forward layer.
The architecture is trained with teacher forcing and decoding is performed using beam search. We use 3 bidirectional
layers with 1024 units each in the encoder and 3 unidirectional layers of the same size in the decoder and a beam
search with a width of 5 candidates. During training, we
use dropout with probability 0.95 after each recurrent layer.
For training, we use the Adam optimization algorithm with
β1 = 0.9, β2 = 0.999, ε = 10−8 , a batch size of 64, a
learning rate of 10−3 , and an annealing rate of 0.85 applied
every 1000 iterations (Kingma & Ba, 2014).
3.2. Segmentation Model
Our segmentation model is trained to output the alignment between a given utterance and a sequence of target
phonemes. This task is similar to the problem of aligning
speech to written output in speech recognition. In that domain, the connectionist temporal classification (CTC) loss
function has been shown to focus on character alignments
to learn a mapping between sound and text (Graves et al.,
2006). We adapt the convolutional recurrent neural network architecture from a state-of-the-art speech recognition system (Amodei et al., 2015) for phoneme boundary
detection.
A network trained with CTC to generate sequences of
phonemes will produce brief peaks for every output
phoneme. Although this is sufficient to roughly align the
phonemes to the audio, it is insufficient to detect precise

Deep Voice: Real-time Neural TTS

phoneme boundaries. To overcome this, we train to predict
sequences of phoneme pairs rather than single phonemes.
The network will then tend to output phoneme pairs at
timesteps close to the boundary between two phonemes in
a pair.
To illustrate our label encoding, consider the string
“Hello!”. To convert this to a sequence of phoneme pair
labels, convert the utterance to phonemes (using a pronunciation dictionary such as CMUDict or a grapheme-tophoneme model) and pad the phoneme sequence on either
end with the silence phoneme to get “sil HH EH L OW sil”.
Finally, construct consecutive phoneme pairs and get “(sil,
HH), (HH, EH), (EH, L), (L, OW), (OW, sil)”.
Input audio is featurized by computing 20 Mel-frequency
cepstral coefficients (MFCCs) with a ten millisecond stride.
On top of the input layer, there are two convolution layers (2D convolutions in time and frequency), three bidirectional recurrent GRU layers, and finally a softmax output
layer. The convolution layers use kernels with unit stride,
height nine (in frequency bins), and width five (in time)
and the recurrent layers use 512 GRU cells (for each direction). Dropout with a probability of 0.95 is applied
after the last convolution and recurrent layers. To compute the phoneme-pair error rate (PPER), we decode using
beam search. To decode phoneme boundaries, we perform
a beam search with width 50 with the constraint that neighboring phoneme pairs overlap by at least one phoneme and
keep track of the positions in the utterance of each phoneme
pair.
For training, we use the Adam optimization algorithm with
β1 = 0.9, β2 = 0.999, ε = 10−8 , a batch size of 128, a
learning rate of 10−4 , and an annealing rate of 0.95 applied
every 500 iterations (Kingma & Ba, 2014).
3.3. Phoneme Duration and Fundamental Frequency
Model
We use a single architecture to jointly predict phoneme duration and time-dependent fundamental frequency. The input to the model is a sequence of phonemes with stresses,
with each phoneme and stress being encoded as a one-hot
vector. The architecture comprises two fully connected layers with 256 units each followed by two unidirectional recurrent layers with 128 GRU cells each and finally a fullyconnected output layer. Dropout with a probability of 0.8 is
applied after the initial fully-connected layers and the last
recurrent layer.
The final layer produces three estimations for every input
phoneme: the phoneme duration, the probability that the
phoneme is voiced (i.e. has a fundamental frequency), and
20 time-dependent F0 values, which are sampled uniformly
over the predicted duration.

The model is optimized by minimizing a joint loss that
combines phoneme duration error, fundamental frequency
error, the negative log likelihood of the probability that
the phoneme is voiced, and a penalty term proportional to
the absolute change of F0 with respect to time to impose
smoothness. The specific functional form of the loss function is described in Appendix B.
For training, we use the Adam optimization algorithm with
β1 = 0.9, β2 = 0.999, ε = 10−8 , a batch size of 128, a
learning rate of 3 × 10−4 , and an annealing rate of 0.9886
applied every 400 iterations (Kingma & Ba, 2014).
3.4. Audio Synthesis Model
Our audio synthesis model is a variant of WaveNet.
WaveNet consists of a conditioning network, which upsamples linguistic features to the desired frequency, and
an autoregressive network, which generates a probability distribution P(y) over discretized audio samples y ∈
{0, 1, . . . , 255}. We vary the number of layers `, the number of residual channels r (dimension of the hidden state of
every layer), and the number of skip channels s (the dimension to which layer outputs are projected prior to the output
layer).
WaveNet consists of an upsampling and conditioning network, followed by ` 2×1 convolution layers with r residual
output channels and gated tanh nonlinearities. We break
the convolution into two matrix multiplies per timestep
with Wprev and Wcur . These layers are connected with
residual connections. The hidden state of every layer is
concatenated to an `r vector and projected to s skip channels with Wskip , followed by two layers of 1 × 1 convolutions (with weights Wrelu and Wout ) with relu nonlinearities.
WaveNet uses transposed convolutions for upsampling and
conditioning. We find that our models perform better, train
faster, and require fewer parameters if we instead first encode the inputs with a stack of bidirectional quasi-RNN
(QRNN) layers (Bradbury et al., 2016) and then perform
upsampling by repetition to the desired frequency.
Our highest-quality final model uses ` = 40 layers, r = 64
residual channels, and s = 256 skip channels. For training, we use the Adam optimization algorithm with β1 =
0.9, β2 = 0.999, ε = 10−8 , a batch size of 8, a learning
rate of 10−3 , and an annealing rate of 0.9886 applied every
1,000 iterations (Kingma & Ba, 2014).
Please refer to Appendix A for full details of our WaveNet
architecture and the QRNN layers we use.

Deep Voice: Real-time Neural TTS

4. Results
We train our models on an internal English speech database
containing approximately 20 hours of speech data segmented into 13,079 utterances. In addition, we present
audio synthesis results for our models trained on a subset
of the Blizzard 2013 data (Prahallad et al., 2013). Both
datasets are spoken by a professional female speaker.
All of our models are implemented using the TensorFlow
framework (Abadi et al., 2015).
4.1. Segmentation Results
We train on 8 TitanX Maxwell GPUs, splitting each batch
equally among the GPUs and using a ring all-reduce to average gradients computed on different GPUs, with each
iteration taking approximately 1300 milliseconds. After
approximately 14,000 iterations, the model converges to a
phoneme pair error rate of 7%. We also find that phoneme
boundaries do not have to be precise, and randomly shifting phoneme boundaries by 10-30 milliseconds makes no
difference in the audio quality, and so suspect that audio
quality is insensitive to the phoneme pair error rate past a
certain point.
4.2. Grapheme-to-Phoneme Results
We train a grapheme-to-phoneme model on data obtained
from CMUDict (Weide, 2008). We strip out all words that
do not start with a letter, contain numbers, or have multiple
pronunciations, which leaves 124,978 out of the original
133,854 grapheme-phoneme sequence pairs.
We train on a single TitanX Maxwell GPU with each iteration taking approximately 150 milliseconds. After approximately 20,000 iterations, the model converges to a
phoneme error rate of 5.8% and a word error rate of 28.7%,
which are on par with previous reported results (Yao &
Zweig, 2015). Unlike prior work, we do not use a language
model during decoding and do not include words with multiple pronunciations in our data set.
4.3. Phoneme Duration and Fundamental Frequency
Results
We train on a single TitanX Maxwell GPU with each iteration taking approximately 120 milliseconds. After approximately 20,000 iterations, the model converges to a mean
absolute error of 38 milliseconds (for phoneme duration)
and 29.4 Hz (for fundamental frequency).
4.4. Audio Synthesis Results
We divide the utterances in our audio dataset into one
second chunks with a quarter second of context for each
chunk, padding each utterance with a quarter second of si-

lence at the beginning. We filter out chunks that are predominantly silence and end up with 74,348 total chunks.
We trained models with varying depth, including 10, 20,
30, and 40 layers in the residual layer stack. We find that
models below 20 layers result in poor quality audio. The
20, 30, and 40 layer models all produce high quality recognizable speech, but the 40 layer models have less noise
than the 20 layer models, which can be detected with highquality over-ear headphones.
Previous work has emphasized the importance of receptive
field size in determining model quality. Indeed, the 20 layer
models have half the receptive field as the 40 layer models. However, when run at 48 kHz, models with 40 layers
have only 83 milliseconds of receptive field, but still generate high quality audio. This suggests the receptive field
of the 20 layer models is sufficient, and we conjecture the
difference in audio quality is due to some other factor than
receptive field size.
We train on 8 TitanX Maxwell GPUs with one chunk per
GPU, using a ring allreduce to average gradients computed
on different GPUs. Each iteration takes approximately 450
milliseconds. Our model converges after approximately
300,000 iterations. We find that a single 1.25s chunk is sufficient to saturate the compute on the GPU and that batching does not increase training efficiency.
As is common with high-dimensional generative models
(Theis et al., 2015), model loss is somewhat uncorrelated
with perceptual quality of individual samples. While models with unusually high loss sound distinctly noisy, models
that optimize below a certain threshold do not have a loss
indicative of their quality. In addition, changes in model
architecture (such as depth and output frequency) can have
a significant impact on model loss while having a small effect on audio quality.
To estimate perceptual quality of the individual stages of
our TTS pipeline, we crowdsourced mean opinion score
(MOS) ratings (ratings between one and five, higher values
being better) from Mechanical Turk using the CrowdMOS
toolkit and methodology (Ribeiro et al., 2011). In order to
separate the effect of the audio preprocessing, the WaveNet
model quality, and the phoneme duration and fundamental
frequency model quality, we present MOS scores for a variety of utterance types, including synthesis results where
the WaveNet inputs (duration and F0) are extracted from
ground truth audio rather than synthesized by other models. The results are presented in Table 1. We purposefully
include ground truth samples in every batch of samples that
raters evaluate to highlight the delta from human speech
and allow raters to distinguish finer grained differences between models; the downside of this approach is that the resulting MOS scores will be significantly lower than if raters

Deep Voice: Real-time Neural TTS

are presented only with synthesized audio samples.
First of all, we find a significant drop in MOS when simply
downsampling the audio stream from 48 kHz to 16 kHz, especially in combination with µ-law companding and quantization, likely because a 48 kHz sample is presented to the
raters as a baseline for a 5 score, and a low quality noisy
synthesis result is presented as a 1. When used with ground
truth durations and F0, our models score highly, with the
95% confidence intervals of our models intersecting those
of the ground truth samples. However, using synthesized
frequency reduces the MOS, and further including synthesized durations reduces it significantly. We conclude that
the main barrier to progress towards natural TTS lies with
duration and fundamental frequency prediction, and our
systems have not meaningfully progressed past the state of
the art in that regard. Finally, our best models run slightly
slower than real-time (see Table 2), so we demonstrate that
synthesis quality can be traded for inference speed by adjusting model size by obtaining scores for models that run
1X and 2X faster than real-time.
We also tested WaveNet models trained on the full set of
features from the original WaveNet publication, but found
no perceptual difference between those models and models
trained on our reduced feature set.
4.5. Blizzard Results
To demonstrate the flexibility of our system, we retrained
all of our models with identical hyperparameters on the
Blizzard 2013 dataset (Prahallad et al., 2013). For our experiments, we used a 20.5 hour subset of the dataset segmented into 9,741 utterances. We evaluated the model using the procedure described in Section 4.4, which encourages raters to compare synthesized audio directly with the
ground truth. On the held out set, 16 kHz companded and
expanded audio receives a MOS score of 4.65±0.13, while
our synthesized audio received a MOS score of 2.67±0.37.

5. Optimizing Inference
Although WaveNet has shown promise in generating highquality synthesized speech, initial experiments reported
generation times of many minutes or hours for short utterances. WaveNet inference poses an incredibly challenging computational problem due to the high-frequency, autoregressive nature of the model, which requires orders of
magnitude more timesteps than traditional recurrent neural
networks. When generating audio, a single sample must
be generated in approximately 60 µs (for 16 kHz audio) or
20 µs (for 48 kHz audio). For our 40 layer models, this
means that a single layer (consisting of several matrix multiplies and nonlinearities) must complete in approximately
1.5 µs. For comparison, accessing a value that resides

in main memory on a CPU can take 0.1 µs. In order to
perform inference at real-time, we must take great care to
never recompute any results, store the entire model in the
processor cache (as opposed to main memory), and optimally utilize the available computational units. These same
techniques could be used to accelerate image synthesis with
PixelCNN (Oord et al., 2016) to fractions of a second per
image.
Synthesizing one second of audio with our 40 layer
WaveNet model takes approximately 55×109 floating point
operations (FLOPs). The activations in any given layer depend on the activations in the previous layer and the previous timestep, so inference must be done one timestep
and one layer at a time. A single layer requires only
42 × 103 FLOPs, which makes achieving meaningful parallelism difficult. In addition to the compute requirements,
the model has approximately 1.6 × 106 parameters, which
equate to about 6.4 MB if represented in single precision.
(See Appendix E for a complete performance model.)
On CPU, a single Haswell or Broadwell core has a peak
single-precision throughput of approximately 77 × 109
FLOPs and an L2-to-L1 cache bandwidth of approximately
140 GB/s 1 . The model must be loaded from cache once per
timestep, which requires a bandwidth of 100 GB/s. Even
if the model were to fit in L2 cache, the implementation
would need to utilize 70% of the maximum bandwidth and
70% of the peak FLOPS in order to do inference in realtime on a single core. Splitting the calculations across
multiple cores reduces the difficulty of the problem, but
nonetheless it remains challenging as inference must operate at a significant fraction of maximum memory bandwidth and peak FLOPs and while keeping threads synchronized.
A GPU has higher memory bandwidth and peak FLOPs
than a CPU but provides a more specialized and hence
restrictive computational model. A naive implementation
that launches a single kernel for every layer or timestep is
untenable, but an implementation based on the persistent
RNN technique (Diamos et al., 2016) may be able to take
advantage of the throughput offered by GPUs.
We implement high-speed optimized inference kernels for
both CPU and GPU and demonstrate that WaveNet inference at faster-than-real-time speeds is achievable. Table 2
lists the CPU and GPU inference speeds for different models. In both cases, the benchmarks include only the autoregressive, high-frequency audio generation and do not
include the generation of linguistic conditioning features
(which can be done in parallel for the entire utterance). Our
CPU kernels run at real-time or faster-than-real-time for a
1

Assuming two 8-wide AVX FMA instructions every cycle
and an L2-to-L1 bandwidth of 64 bytes per cycle.

Deep Voice: Real-time Neural TTS

Type
Ground Truth (48 kHz)
Ground Truth
Ground Truth (companded and expanded)
Synthesized
Synthesized (48 kHz)
Synthesized (Synthesized F0)
Synthesized (Synthesized Duration and F0)
Synthesized (2X real-time inference)
Synthesized (1X real-time inference)

Model Size
None
None
None
` = 40, r = 64, s = 256
` = 40, r = 64, s = 256
` = 40, r = 64, s = 256
` = 40, r = 64, s = 256
` = 20, r = 32, s = 128
` = 20, r = 64, s = 128

MOS±CI
4.75 ± 0.12
4.45 ± 0.16
4.34 ± 0.18
3.94 ± 0.26
3.84 ± 0.24
2.76 ± 0.31
2.00 ± 0.23
2.74 ± 0.32
3.35 ± 0.31

Table 1. Mean Opinion Scores (MOS) and 95% confidence intervals (CIs) for utterances. This MOS score is a relative MOS score
obtained by showing raters the same utterance across all the model types (which encourages comparative rating and allows the raters
to distinguish finer grained differences). Every batch of samples also includes the ground truth 48 kHz recording, which makes all our
ratings comparative to natural human voices. 474 ratings were collected for every sample. Unless otherwise mentioned, models used
phoneme durations and F0 extracted from the ground truth, rather than synthesized by the duration prediction and frequency prediction
models, as well as a 16384 Hz audio sampling rate.

Model
` = 20, r = 32, s = 128
` = 20, r = 32, s = 128
` = 20, r = 64, s = 128
` = 20, r = 64, s = 128
` = 20, r = 64, s = 128
` = 40, r = 64, s = 256
` = 40, r = 64, s = 256
` = 40, r = 64, s = 256
` = 20, r = 32, s = 128
` = 20, r = 64, s = 128
` = 40, r = 32, s = 128
` = 40, r = 64, s = 128

Platform
CPU
CPU
CPU
CPU
CPU
CPU
CPU
CPU
GPU
GPU
GPU
GPU

Data Type
float32
float32
int16
float32
float32
int16
float32
float32
float32
float32
float32
float32

Number of Threads
6
2
2
6
2
2
6
2
N/A
N/A
N/A
N/A

Speed-up Over Real-time
2.7
2.05
1.2
1.11
0.79
0.67
0.61
0.35
0.39
0.29
0.23
0.17

Table 2. CPU and GPU inference kernel benchmarks for different models in float32 and int16. At least one main and one auxiliary
thread were used for all CPU kernels. These kernels operate on a single utterance with no batching. CPU results are from a Intel Xeon
E5-2660 v3 Haswell processor clocked at 2.6 GHz and GPU results are from a GeForce GTX Titan X Maxwell GPU.

subset of models, while the GPU models do not yet match
this performance.

5.1. CPU Implementation
We achieve real-time CPU inference by avoiding any recomputation, doing cache-friendly memory accesses, parallelizing work via multithreading with efficient synchronization, minimizing nonlinearity FLOPs, avoiding cache
thrashing and thread contention via thread pinning, and using custom hardware-optimized routines for matrix multiplication and convolution.
For the CPU implementation, we split the computation into
the following steps:
1. Sample Embedding: Compute the WaveNet input
causal convolution by doing two sample embeddings,
one for the current timestep and one for the previous
timestep, and summing them with a bias. That is,
x(0) = Wemb,prev · yi−1 + Wemb,cur · yi + Bembed (1)

Deep Voice: Real-time Neural TTS
(1)

Main Threads

(2)

acur

acur

x(1)

x(1)

(`−1)

...

(1)

(`)

acur

acur

zs → za → p

x(`−1) x(`)

(2)

(`−1)

acur → acur → . . . → acur

x(0)

(`)

→ acur

x(1) → x(2) → . . . → x(`−1) → x(`)

Sync Points
timestep t
Auxiliary Threads

q

(1)

q

(2)

...

q

(`−1)

q

(`)

(1)
aprev

→

(2)
aprev

→ ... →

(`−1)
aprev

→

(`)
aprev

timestep t + 1

Figure 2. Two groups of threads run in parallel. Computation of the Wskip is offloaded to the auxiliary threads while the main threads
progress through the stack of WaveNet layers. While the main threads are computing the output layer, the auxiliary threads prepare the
left Wprev half of the WaveNet layer convolutions for the upcoming timestep. Arrows indicate where one thread group waits on results
from the other thread group, and are implemented as spinlocks.

2. Layer Inference: For every layer j from j = 1 to `
with dilation width d:
(a) Compute the left half of the width-two dilated
convolution via a matrix-vector multiply:
(j−1)

(j)
a(j)
prev = Wprev · xi−d

(2)

(b) Compute the right half of the dilated convolution:
(j−1)

(j)
a(j)
cur = Wcur · xi

(3)

(c) Compute the hidden state h(j) given the condi(j)
tioning vector Lh :
(j)

a

h(j)

=

a(j)
prev

a(j)
cur

(j)
Bh

(j)
Lh

+
+
+

 

(j)
(j)
= tanh a0:r · σ ar:2r ,

(4)
(5)

where v0:r denotes the first r elements of the vector v and vr:2r denotes the next r elements. Then,
compute the input to the next layer via a matrixvector multiply:
(j)
(j)
· h(j) + Bres
x(j) = Wres

(6)

(d) Compute the contribution to the skip-channel
matrix multiply from this layer, accumulating
over all layers, with q (0) = Bskip :
(j)

q (j) = q (j−1) + Wskip · h(j)

(7)

3. Output: Compute the two output 1 × 1 convolutions:
 
zs = relu q (`)
(8)
za = relu (Wrelu · zs + Brelu )
p = softmax (Wout · za + Bout )

(9)
(10)

Finally, sample yi+1 randomly from the distribution p.

We parallelize these across two groups of threads as depicted in Figure 2. A group of main threads computes
(j)
x(0) , acur , h(j) , and x(j) , za , and p. A group of auxiliary
(j)
(j)
threads computes aprev , q (j) , and zs , with the aprev being
computed for the next upcoming timestep while the main
threads compute za and p. Each of these groups can consist of a single thread or of multiple threads; if there are
multiple threads, each thread computes one block of each
matrix-vector multiply, binary operation, or unary operation, and thread barriers are inserted as needed. Splitting
the model across multiple threads both splits up the compute and can also be used to ensure that the model weights
fit into the processor L2 cache.
Pinning threads to physical cores (or disabling hyperthreading) is important for avoiding thread contention and
cache thrashing and increases performance by approximately 30%.
Depending on model size, the nonlinearities (tanh,
sigmoid, and softmax) can also take a significant fraction of inference time, so we replace all nonlinearities with
high-accuracy approximations, which are detailed in Appendix C. The maximum absolute error arising from these
approximations is 1.5 × 10−3 for tanh, 2.5 × 10−3 for
sigmoid, and 2.4 × 10−5 for ex . With approximate instead
of exact nonlinearities, performance increases by roughly
30%.
We also implement inference with weight matrices quantized to int16 and find no change in perceptual quality
when using quantization. For larger models, quantization
offers a significant speedup when using fewer threads, but
overhead of thread synchronization prevents it from being
useful with a larger number of threads.
Finally, we write custom AVX assembly kernels for matrixvector multiplication using PeachPy (Dukhan, 2015) specialized to our matrix sizes. Inference using our custom
assembly kernels is up to 1.5X faster than Intel MKL and
3.5X faster than OpenBLAS when using float32. Nei-

Deep Voice: Real-time Neural TTS

ther library provides the equivalent int16 operations.
5.2. GPU Implementation
Due to their computational intensity, many neural models
are ultimately deployed on GPUs, which can have a much
higher computational throughput than CPUs. Since our
model is memory bandwidth and FLOP bound, it may seem
like a natural choice to run inference on a GPU, but it turns
out that comes with a different set of challenges.
Usually, code is run on the GPU in a sequence of kernel
invocations, with every matrix multiply or vector operation
being its own kernel. However, the latency for a CUDA
kernel launch (which may be up to 50 µs) combined with
the time needed to load the entire model from GPU memory are prohibitively large for an approach like this. An
inference kernel in this style ends up being approximately
1000X slower than real-time.
To get close to real-time on a GPU, we instead build a kernel using the techniques of persistent RNNs (Diamos et al.,
2016) which generates all samples in the output audio in a
single kernel launch. The weights for the model are loaded
to registers once and then used without unloading them for
the entire duration of inference. Due to the mismatch between the CUDA programming model and such persistent
kernels, the resulting kernels are specialized to particular
model sizes and are incredibly labor-intensive to write. Although our GPU inference speeds are not quite real-time
(Table 2), we believe that with these techniques and a better implementation we can achieve real-time WaveNet inference on GPUs as well as CPUs. Implementation details
for the persistent GPU kernels are available in Appendix D.

6. Conclusion
In this work, we demonstrate that current Deep Learning
approaches are viable for all the components of a highquality text-to-speech engine by building a fully neural system. We optimize inference to faster-than-real-time speeds,
showing that these techniques can be applied to generate audio in real-time in a streaming fashion. Our system
is trainable without any human involvement, dramatically
simplifying the process of creating TTS systems.
Our work opens many new possible directions for exploration. Inference performance can be further improved
through careful optimization, model quantization on GPU,
and int8 quantization on CPU, as well as experimenting with other architectures such as the Xeon Phi. Another natural direction is removing the separation between
stages and merging the segmentation, duration prediction,
and fundamental frequency prediction models directly into
the audio synthesis model, thereby turning the problem into
a full sequence-to-sequence model, creating a single end-

to-end trainable TTS system, and allowing us to train the
entire system with no intermediate supervision. In lieu of
fusing the models, improving the duration and frequency
models via larger training datasets or generative modeling
techniques may have an impact on voice naturalness.

Deep Voice: Real-time Neural TTS

References
Abadi, Martı́n, Agarwal, Ashish, Barham, Paul, Brevdo,
Eugene, Chen, Zhifeng, Citro, Craig, Corrado, Greg S.,
Davis, Andy, Dean, Jeffrey, Devin, Matthieu, Ghemawat, Sanjay, Goodfellow, Ian, Harp, Andrew, Irving, Geoffrey, Isard, Michael, Jia, Yangqing, Jozefowicz,
Rafal, Kaiser, Lukasz, Kudlur, Manjunath, Levenberg,
Josh, Mané, Dan, Monga, Rajat, Moore, Sherry, Murray,
Derek, Olah, Chris, Schuster, Mike, Shlens, Jonathon,
Steiner, Benoit, Sutskever, Ilya, Talwar, Kunal, Tucker,
Paul, Vanhoucke, Vincent, Vasudevan, Vijay, Viégas,
Fernanda, Vinyals, Oriol, Warden, Pete, Wattenberg,
Martin, Wicke, Martin, Yu, Yuan, and Zheng, Xiaoqiang.
TensorFlow: Large-scale machine learning on heterogeneous systems, 2015. URL http://tensorflow.
org/. Software available from tensorflow.org.
Amodei, Dario, Anubhai, Rishita, Battenberg, Eric, Case,
Carl, Casper, Jared, Catanzaro, Bryan, Chen, Jingdong,
Chrzanowski, Mike, Coates, Adam, Diamos, Greg, et al.
Deep speech 2: End-to-end speech recognition in english
and mandarin. arXiv preprint arXiv:1512.02595, 2015.
Boersma, Paulus Petrus Gerardus et al. Praat, a system
for doing phonetics by computer. Glot international, 5,
2002.
Bradbury, James, Merity, Stephen, Xiong, Caiming, and
Socher, Richard. Quasi-recurrent neural networks. arXiv
preprint arXiv:1611.01576, 2016.
Chung, Junyoung, Gulcehre, Caglar, Cho, KyungHyun,
and Bengio, Yoshua. Empirical evaluation of gated recurrent neural networks on sequence modeling. arXiv
preprint arXiv:1412.3555, 2014.
Diamos, Greg, Sengupta, Shubho, Catanzaro, Bryan,
Chrzanowski, Mike, Coates, Adam, Elsen, Erich, Engel,
Jesse, Hannun, Awni, and Satheesh, Sanjeev. Persistent
rnns: Stashing recurrent weights on-chip. In Proceedings of The 33rd International Conference on Machine
Learning, pp. 2024–2033, 2016.
Dukhan, Marat. Peachpy meets opcodes: direct machine
code generation from python. In Proceedings of the 5th
Workshop on Python for High-Performance and Scientific Computing, pp. 3. ACM, 2015.

Mehri, Soroush, Kumar, Kundan, Gulrajani, Ishaan, Kumar, Rithesh, Jain, Shubham, Sotelo, Jose, Courville,
Aaron, and Bengio, Yoshua. Samplernn: An unconditional end-to-end neural audio generation model. arXiv
preprint arXiv:1612.07837, 2016.
Morise, Masanori, Yokomori, Fumiya, and Ozawa, Kenji.
World: a vocoder-based high-quality speech synthesis
system for real-time applications. IEICE TRANSACTIONS on Information and Systems, 99(7):1877–1884,
2016.
Oord, Aaron van den, Kalchbrenner, Nal, and
Kavukcuoglu, Koray. Pixel recurrent neural networks.
arXiv preprint arXiv:1601.06759, 2016.
Paine, Tom Le, Khorrami, Pooya, Chang, Shiyu, Zhang,
Yang, Ramachandran, Prajit, Hasegawa-Johnson,
Mark A, and Huang, Thomas S. Fast wavenet generation algorithm. arXiv preprint arXiv:1611.09482,
2016.
Pascual, Santiago and Bonafonte, Antonio. Multi-output
rnn-lstm for multiple speaker speech synthesis with αinterpolation model. way, 1000:2, 2016.
Prahallad, Kishore, Vadapalli, Anandaswarup, Elluru,
Naresh, et al. The blizzard challenge 2013indian language task. In In Blizzard Challenge Workshop 2013,
2013.
Rao, Kanishka, Peng, Fuchun, Sak, Haşim, and Beaufays, Françoise. Grapheme-to-phoneme conversion using long short-term memory recurrent neural networks.
In Acoustics, Speech and Signal Processing (ICASSP),
2015 IEEE International Conference on, pp. 4225–4229.
IEEE, 2015.
Ribeiro, Flávio, Florêncio, Dinei, Zhang, Cha, and Seltzer,
Michael. Crowdmos: An approach for crowdsourcing
mean opinion score studies. In Acoustics, Speech and
Signal Processing (ICASSP), 2011 IEEE International
Conference on, pp. 2416–2419. IEEE, 2011.
Ronanki, Srikanth, Henter, Gustav Eje, Wu, Zhizheng, and
King, Simon. A template-based approach for speech
synthesis intonation generation using lstms. Interspeech
2016, pp. 2463–2467, 2016.

Graves, Alex, Fernández, Santiago, Gomez, Faustino, and
Schmidhuber, Jürgen. Connectionist temporal classification: Labelling unsegmented sequence data with recurrent neural networks. In Proceedings of the 23rd International Conference on Machine Learning, ICML ’06,
pp. 369–376, New York, NY, USA, 2006. ACM.

Sotelo, Jose, Mehri, Soroush, Kumar, Kundan, Santos,
Joao Felipe, Kastner, Kyle, Courville, Aaron, and Bengio, Yoshua. Char2wav: End-to-end speech synthesis. In
ICLR 2017 workshop submission, 2017. URL https:
//openreview.net/forum?id=B1VWyySKx.

Kingma, D. and Ba, J. Adam: A method for stochastic
optimization. arXiv preprint arXiv:1412.6980, 2014.

Stephenson, Ian. Production Rendering, Design and Implementation. Springer, 2005.

Deep Voice: Real-time Neural TTS

Taylor, Paul. Text-to-Speech Synthesis. Cambridge University Press, New York, NY, USA, 1st edition, 2009. ISBN
0521899273, 9780521899277.
Theis, Lucas, Oord, Aäron van den, and Bethge, Matthias.
A note on the evaluation of generative models. arXiv
preprint arXiv:1511.01844, 2015.
van den Oord, Aäron, Dieleman, Sander, Zen, Heiga, Simonyan, Karen, Vinyals, Oriol, Graves, Alex, Kalchbrenner, Nal, Senior, Andrew, and Kavukcuoglu, Koray. Wavenet: A generative model for raw audio. CoRR
abs/1609.03499, 2016.
Weide, R.
The CMU pronunciation dictionary 0.7.
Carnegie Mellon University, 2008.
Yao, Kaisheng and Zweig, Geoffrey.
Sequence-tosequence neural net models for grapheme-to-phoneme
conversion. arXiv preprint arXiv:1506.00196, 2015.
Zen, Heiga and Sak, Haşim. Unidirectional long short-term
memory recurrent neural network with recurrent output
layer for low-latency speech synthesis. In Acoustics,
Speech and Signal Processing (ICASSP), 2015 IEEE International Conference on, pp. 4470–4474. IEEE, 2015.
Zen, Heiga, Senior, Andrew, and Schuster, Mike. Statistical
parametric speech synthesis using deep neural networks.
In Proceedings of the IEEE International Conference on
Acoustics, Speech, and Signal Processing (ICASSP), pp.
7962–7966, 2013.

Deep Voice: Real-time Neural TTS

Appendices
A. WaveNet Architecture and Details
The WaveNet consists of a conditioning network c = C(v), which converts low-frequency linguistic features v to the
native audio frequency, and an auto-regressive process P (yi |c, yi−1 , . . . , yi−R ) which predicts the next audio sample given
the conditioning for the current timestep c and a context of R audio samples. R is the receptive field size, and is a property
determined by the structure of the network. A sketch of the WaveNet architecture is shown in Figure 3. The network details
are described in the following subsections.

+

L layers

softmax

1x1 conv,
r outputs
1x1 conv,
s outputs
tanh

x

+

+

ReLU

1x1 conv,
a outputs

ReLU

1x1 conv,
a outputs

σ

+

Split
channels

nearest neighbor
upsampling to 16kHz

interleave
stack

2x1 dilated conv,
2r outputs
Forward QRNN L x r

stack

tanh

Forward QRNN q0 /2
2x1 conv,
r outputs

Backward QRNN L x r

16kHz audio input,
a one-hot features

Backward QRNN q0 /2

256Hz phoneme and F0 features,
227 dimensions

Figure 3. The modified WaveNet architecture. Components are colored according to function: teal inputs, green convolutions and
QRNNs, yellow unary operations and softmax, pink binary operations, and indigo reshapes, transposes, and slices.

A.1. Auto-regressive WaveNet
The structure of the auto-regressive network is parameterized by the number of layers `, the number of skip channels s,
and the number of residual channels r.
Audio is quantized to a = 256 values using µ-law companding, as described in Section 2.2 of WaveNet. The one-hot
encoded values go through an initial 2x1 convolution which generates the input x(0) ∈ Rr for the first layer in the residual
stack:
x(0) = Wembed ∗ y + Bembed ,

(11)

Deep Voice: Real-time Neural TTS

where ∗ is the one-dimensional convolution operator. Since the input audio y is a one-hot vector, this convolution can be
done via embeddings instead of matrix multiplies. Each subsequent layer computes a hidden state vector h(i) and then (due
to the residual connections between layers) adds to its input x(i−1) to generate its output x(i) :


 
(i)
(i)
(i)
(12)
h(i) = tanh Wh ∗ x(i−1) + Bh + Lh · σ Wg(i) ∗ x(i−1) + Bg(i) + L(i)
g
x(i) = x(i−1) + Wr(i) · h(i) + Br(i) ,

(13)

where L(i) is the output for that layer of the conditioning network. Since each layer adds its output to its input, the
dimensionality of the layers must remain fixed to the number of residual channels, r. Although here this is written as two
convolutions, one for Wh and one for Wg , it is actually done more efficiently with a single convolution with r input and
2r output channels. During inference, this convolution is replaced with two matrix-vector multiplies with matrices Wprev
(the left half of the convolution) and Wcur (the right half). Thus we can reformulate the computation of h(i) for a specific
timestep t as follows:
(i−1)

(i−1)

(i)
(i)
h0(i) = Wprev
· xt−d + Wcur
· xt

 

0(i)
0(i)
h(i) = tanh h0:r · σ hr:2r ,
(i)

+ B (i) + L(i)

(i)

(14)
(15)

(i)

(i)

where L(i) and is a concatenation of Lh and Lg and B (i) and is a concatenation of Bh and Bg .
The hidden state h(i) from each of the layers 1 through ` is concatenated and projected with a learned Wskip down to the
number of skip channels s:
 (1) 
h
h(2) 


h ∈ R`r
(16)
h= . 
 .. 
h(`)
zs = relu (Wskip · h + Bskip ) ,

zs ∈ Rs

(17)

where relu(x) = max(0, x).
zs is then fed through two fully connected relu layers to generate the output distribution p ∈ Ra :
za = relu (Wrelu · zs + Brelu ) ,
p = softmax (Wout · za + Bout )

za ∈ Ra

(18)
(19)

A.2. Conditioning Network
When trained without conditioning information, WaveNet models produce human-like “babbling sounds”, as they lack sufficient long-range information to reproduce words. In order to generate recognizable speech, every timestep is conditioned
by an associated set of linguistic features. This is done by biasing every layer with a per-timestep conditioning vector
generated from a lower-frequency input signal containing phoneme, stress, and fundamental frequency features.
The frequency of the audio is significantly higher than the frequency of the linguistic conditioning information, so an
upsampling procedure is used to convert from lower-frequency linguistic features to higher-frequency conditioning vectors
for each WaveNet layer.
The original WaveNet does upsampling by repetition or through a transposed convolution. Instead, we first pass our input
features through two bidirectional quasi-RNN layers (Bradbury et al., 2016) with fo-pooling and 2x1 convolutions. A
unidirectional QRNN layer with fo-pooling is defined by the following equations:
h̃ = tanh (Wh ∗ x + Bh )

(20)

o = σ (Wo ∗ x + Bo )

(21)

f = σ (Wf ∗ x + Bf )

(22)

ht = ft · ht−1 + (1 − ft ) · h̃t

(23)

zt = ot · ht

(24)

Deep Voice: Real-time Neural TTS

A bidirectional QRNN layer is computed by running two unidirectional QRNNs, one on the input sequence and one on
a reversed copy of the input sequence, and then stacking their output channels. After both QRNN layers, we interleave
the channels, so that the tanh and the sigmoid in the WaveNet both get channels generated by the forward QRNN and
backward QRNN.
Following the bidirectional QRNN layers, we upsample to the native audio frequency by repetition2 .
We find that the model is very sensitive to the upsampling procedure: although many variations of the conditioning network
converge, they regularly produce phoneme mispronunciations.
A.3. Input Featurization
Our WaveNet is trained with 8-bit µ-law companded audio which is downsampled to 16384 Hz from 16-bit dual-channel
PCM audio at 48000 Hz. It is conditioned on a 256 Hz phoneme signal. The conditioning feature vector has 227 dimensions. Of these two are for fundamental frequency. One of these indicates whether the current phoneme is voiced (and thus
has an F0) and the other is normalized log-frequency, computed by normalizing the log of F0 to minimum observed F0 to
be approximately between -1 and 1. The rest of the features describe the current phoneme, the two previous phonemes,
and the two next phonemes, with each phoneme being encoded via a 40-dimensional one-hot vector for phoneme identity
(with 39 phonemes for ARPABET phonemes and 1 for silence) and a 5-dimensional one-hot vector for phoneme stress (no
stress, primary stress, secondary stress, tertiary stress, and quaternary stress). Not all of the datasets we work with have
tertiary or quaternary stress, and those features are always zero for the datasets that do not have those stress levels.
In our experiments, we found that including the phoneme context (two previous and two next phonemes) is crucial for
upsampling via transposed convolution and less critical but still important for our QRNN-based upsampling. Although
sound quality without the phoneme context remains high, mispronunciation of a subset of the utterances becomes an issue.
We also found that including extra prosody features such as word and syllable breaks, pauses, phoneme and syllable counts,
frame position relative to phoneme, etc, were unhelpful and did not result in higher quality synthesized samples.
In order to convert from phonemes annotated with durations to a fixed-frequency phoneme signal, we sample the phonemes
at regular intervals, effectively repeating each phoneme (with context and F0) a number proportional to its duration. As a
result, phoneme duration is effectively quantized to 1/256 sec ≈ 4ms.
We use Praat (Boersma et al., 2002) in batch mode to compute F0 at the appropriate frequency, with a minimum F0
of 75 and a maximum F0 of 500. The Praat batch script used to generate F0 is available at https://github.com/baiduresearch/deep-voice/blob/master/scripts/f0-script.praat and can be run with praat --run f0-script.praat.
A.4. Sampling from Output Distribution
At every timestep, the synthesis model produces a distribution over samples, P (s), conditioned on the previous samples
and the linguistic features. In order to produce the samples, there are a variety of ways you could choose to use this
distribution:
• Direct Sampling: Sample randomly from P (y).
• Temperature Sampling: Sample randomly from a distribution adjusted by a temperature t
P̃t (y) =

1
1
P (y) /t ,
Z

(25)

where Z is a normalizing constant.
• Mean: Take the mean of the distribution EP [y].
• Mode: Take the most likely sample, argmax P (y).
2

Upsampling using bilinear interpolation slowed convergence and reduced generation quality by adding noise or causing mispronunciations, while bicubic upsampling led to muffled sounds. Upsampling by repetition is done by computing the ratio of the output
frequency to the input frequency and repeating every element in the input signal an appropriate number of times.

Deep Voice: Real-time Neural TTS

• Top k: Sample from an adjusted distribution that only permits the top k samples
(
0
if y < kth(P (y))
P̃k (y) =
P (y)/Z otherwise,

(26)

where Z is a normalizing constant.
We find that out of these different sampling methods, only direct sampling produces high quality outputs. Temperature sampling produces acceptable quality results, and indeed outperforms direct sampling early on in training, but for converged
models is significantly worse. This observation indicates that the generative audio model accurately learns a conditional
sample distribution and that modifying this distribution through the above heuristics is worse than just using the learned
distribution.
A.5. Training
We observed several tendencies of the models during training. As expected, the randomly initialized model produces white
noise. Throughout training, the model gradually increases the signal to noise ratio, and the volume of the white noise
dies down while the volume of the speech signal increases. The speech signal can be inaudible for tens of thousands of
iterations before it dominates the white noise.
In addition, because the model is autoregressive, rare mistakes can produce very audible disturbances. For example, a
common failure mode is to produce a small number of incorrect samples during sampling, which then results in a large
number incorrect samples due to compounding errors. This is audible as a brief period of loud noise before the model
stabilizes. The likelihood of this happening is higher early on in training, and does not happen in converged models.

B. Phoneme Model Loss
The loss for the nth phoneme is
Ln = |t̂n − tn | + λ1 CE(p̂n , pn ) + λ2

T
−1
X

c0n,t − F 0n,t | + λ3
|F

t=0

T
−2
X

c0n,t+1 − F
c0n,t |,
|F

(27)

t=0

where λi ’s are tradeoff constants, b
tn and tn are the estimated and ground-truth durations of the nth phoneme, p̂n and pn
are the estimated and ground-truth probabilities that the nth phoneme is voiced, CE is the cross-entropy function, Fbn,t and
Fn,t are the estimated and ground-truth values of the fundamental frequency of the nth phoneme at time t. T time samples
are equally spaced along the phoneme duration.

C. Nonlinearity Approximation Details
During inference, we replace exact implementations of the neural network nonlinearities with high-accuracy rational approximations. In this appendix, we detail the derivation of these approximations.
C.1. tanh and sigmoid approximation
Denoting ẽ(x) as an approximation to e|x| , we use the following approximations for tanh and σ:
ẽ(x) − 1/ẽ(x)
tanh(x) ≈ sign(x)
ẽ(x) + 1/ẽ(x)
( ẽ(x)
x≥0
σ(x) ≈ 1+ẽ(x)
1
x≤0
1+ẽ(x)

(28)
(29)

We choose a forth-order polynomial to represent ẽ(x). The following fit produces accurate values for both tanh(x) and
σ(x):
ẽ(x) = 1 + |x| + 0.5658x2 + 0.143x4
|x|

(30)

By itself, ẽ(x) is not a very good approximate function for e , but it yields good approximations when used to approximate
tanh and σ as described in Equations 28 and 29.

Deep Voice: Real-time Neural TTS

C.2. ex approximation
We follow the approach of (Stephenson, 2005) to calculate an approximate ex function. Instead of approximating ex
directly, we approximate 2x and use the identity ex = 2x/ln 2 .
Let bxc to be the floor of x ∈ R. Then,
2x = 2bxc · 2x−bxc


= 2bxc · 1 + (2x−bxc − 1)
where 0 ≤ 2x−bxc − 1 < 1 since 0 ≤ x − bxc < 1. If we use a 32-bit float to represent 2x , then bxc + 127 and 2x−bxc − 1
are represented by the exponent and fraction bits of 2x . Therefore, if we interpret the bytes pattern of 2x as a 32-bits integer
(represented by I2x ), we have
I2x = (bxc + 127) · 223 + (2x−bxc − 1) · 223 .
(31)
Rearranging the Equation 31 and using z = x − bxc results to
I2x = (x + 126 + {2z − z}) · 223

(32)

If we can accurately approximate g(z) = 2z − z over z ∈ [0, 1), then interpreting back the byte representation of I2x in
Equation 32 as a 32-bits float, we can accurately approximate 2x . We use a rational approximation as
g(z) ≈ −4.7259162 +

27.7280233
− 1.49012907z,
4.84252568 − z

(33)

which gives are maximum error 2.4 × 10−5 for x ∈ (−∞, 0].

D. Persistent GPU Kernels
A NVIDIA GPU has multiple Streaming Multiprocessors (SMs), each of which has a register file and a L1 cache. There
is also a coherent L2 cache that is shared by all SMs. The inference process needs to generate one sample every 61 µs.
Due to the high latency of a CUDA kernel launch and of reading small matrices from GPU memory, the entire audio
generation process must be done by a single kernel with the weights loaded into the register file across all SMs. This raises
two challenges—how to split the model across registers in a way to minimize communication between SMs and how to
communicate between SMs given the restrictions imposed by the CUDA programming model.
We split the model across the register file of 24 SMs, numbered SM1 · SM24, of a TitanX GPU. We do not use SM24.
SM1 to SM20 store two adjacent layers of the residual stack. This means SM1 stores layers 1 and 2, SM2 stores layers
3 and 4 and so on and so forth. Each layer has three matrices and three bias vectors—Wprev , Bprev , Wcur , Bcur , that are
for the dilated convolutions and Wr , Br . Thus SMi generates two hidden states h(2i) and h(2i+1) and an output x(2i) .
Each SM also stores the rows of the Wskip matrix that will interact with the generated hidden state vectors. Thus Wskip is
partitioned across 20 SMs. Only SM20 needs to store Bskip . SM21 stores Wrelu and Brelu . Finally, Wout is split across
two SMs—SM22 and SM23 because of register file limitations and SM23 stores Bout .
The next challenge is to coordinate the data transfer between SMs, since the CUDA programming model executes one
kernel across all SMs in parallel. However we want execution to go sequentially in a round robin fashion from SM1 to
SM23 and back again from SM1 as we generate one audio sample at a time. We launch our CUDA kernel with 23 thread
blocks and simulate such sequential execution by spinning on locks, one for each SM, that are stored in global memory and
cached in L2. First SM1 executes two layers of the WaveNet model to generate h(1) , h(2) and x(2) . It then unlocks the lock
that SM2 is spinning on and sets its own lock. It does this by bypassing the L1 cache to write to global memory so that all
SMs have a coherent view of the locks. Then SM2 does the same for SM3 and this sequential locking and unlocking chain
continues for each SM. Finally SM23 generates the output distribution p for timestep t and unlocks SM1 so that entire
process can repeat to generate p for timestep t + 1.
Just like locks, we pass data between SMs, by reading and writing to global memory by bypassing the L1 cache. Since
NVIDIA GPUs have a coherent L2 cache, a global memory write bypassing the L1, followed by a memory fence results
in a coherent view of memory across SMs.

Deep Voice: Real-time Neural TTS

This partitioning scheme however is quite inflexible and only works for specific values of l, r and s shown in Table 2. This
is because each SM has a fixed sized register file and combined with the relatively inflexible and expensive communication mechanism between SMs implies that splitting weight matrices between SMs is challenging. Any change in those
parameters means a new kernel has to be written, which is a very time consuming process.
There are two main reasons why the GPU kernels are slower than CPU kernels. Firstly, synchronization between SMs in a
GPU is expensive since it is done by busy waiting on locks in L2 cache. Secondly even though we divide the model in a
way that will fit in the register file of each SM, the CUDA compiler still spills to L1 cache. We hope that with handcrafted
assembly code, we will be able to match the performance of CPU kernels. However, the lack of parallelism in WaveNet
inference makes it difficult to hide the latencies inherent in reading and writing small matrices from GPU memory which
are exposed in the absence of a rich cache hierarchy in GPUs.

E. Performance model
We present a performance model for the autoregressive WaveNet architecture described in Appendix A.1. In our model a
dot product between two vectors of dimension r takes 2r FLOPS—r multiplications and r additions. This means that a
matrix-vector multiply between W , an r × r matrix and x, a r × 1 vector takes 2r · r = 2r2 FLOPs. Thus calculating h0(i)
uses


Cost h0(i) = (2r · 2r) + (2r · 2r) + 2r + 2r + 2r FLOPs
(34)
Let division and exponentiaton take fd and fe FLOPs respectively. This means tanh and σ takes (fd + 2fe + 1) FLOPs.
Thus calculating h(i) takes 2r · (fd + 2fe + 1) + r FLOPs. Finally calculating x(i) for each layer takes r + (2r · r) + r
FLOPs. This brings the total FLOPs for calculating one layer to
Cost (layer) = 10r2 + 11r + 2r(fd + fe ) FLOPs

(35)

Under the same model, calculating zs takes (` · 2r) · s + s + s FLOPs, where we assume that relu takes 1 FLOP. Similarly,
calculating za takes 2s · a + a + a FLOPs and Wout · za + Bout takes 2a · a + a FLOPs.
Calculating the numerically stable softmax takes one max, one subtract, one exponentiation, one sum and one division per
element of a vector. Hence calculating p takes 3a + a(fd + fe ) FLOPs.
Adding it all up, our final performance model to generate each audio sample is as follows:

Cost (sample) = ` 10r2 + 11r + 2r(fd + fe ) + s(2r · ` + 2) + a(2s + 2a + 3) + a(3 + fd + fe ) FLOPS

(36)

If we let ` = 40, r = 64, and s = a = 256, and assume that fd = 10 and fe = 10, with a sampling frequency of 16384Hz,
we have approximately 55 × 109 FLOPs for every second of synthesis.

