# darts_differentiable_architecture_search__794ac17a.pdf

Published as a conference paper at ICLR 2019

DARTS: DIFFERENTIABLE ARCHITECTURE SEARCH

Hanxiao Liu CMU hanxiaol@cs.cmu.com

Karen Simonyan Deep Mind simonyan@google.com

Yiming Yang CMU yiming@cs.cmu.edu

This paper addresses the scalability challenge of architecture search by formulating the task in a differentiable manner. Unlike conventional approaches of applying evolution or reinforcement learning over a discrete and non-differentiable search space, our method is based on the continuous relaxation of the architecture representation, allowing efficient search of the architecture using gradient descent. Extensive experiments on CIFAR-10, Image Net, Penn Treebank and Wiki Text-2 show that our algorithm excels in discovering high-performance convolutional architectures for image classification and recurrent architectures for language modeling, while being orders of magnitude faster than state-of-the-art non-differentiable techniques. Our implementation has been made publicly available to facilitate further research on efficient architecture search algorithms.

1 INTRODUCTION

Discovering state-of-the-art neural network architectures requires substantial effort of human experts. Recently, there has been a growing interest in developing algorithmic solutions to automate the manual process of architecture design. The automatically searched architectures have achieved highly competitive performance in tasks such as image classification (Zoph & Le, 2017; Zoph et al., 2018; Liu et al., 2018b;a; Real et al., 2018) and object detection (Zoph et al., 2018).

The best existing architecture search algorithms are computationally demanding despite their remarkable performance. For example, obtaining a state-of-the-art architecture for CIFAR-10 and Image Net required 2000 GPU days of reinforcement learning (RL) (Zoph et al., 2018) or 3150 GPU days of evolution (Real et al., 2018). Several approaches for speeding up have been proposed, such as imposing a particular structure of the search space (Liu et al., 2018b;a), weights or performance prediction for each individual architecture (Brock et al., 2018; Baker et al., 2018) and weight sharing/inheritance across multiple architectures (Elsken et al., 2017; Pham et al., 2018b; Cai et al., 2018; Bender et al., 2018), but the fundamental challenge of scalability remains. An inherent cause of inefficiency for the dominant approaches, e.g. based on RL, evolution, MCTS (Negrinho & Gordon, 2017), SMBO (Liu et al., 2018a) or Bayesian optimization (Kandasamy et al., 2018), is the fact that architecture search is treated as a black-box optimization problem over a discrete domain, which leads to a large number of architecture evaluations required.

In this work, we approach the problem from a different angle, and propose a method for efficient architecture search called DARTS (Differentiable ARchi Tecture Search). Instead of searching over a discrete set of candidate architectures, we relax the search space to be continuous, so that the architecture can be optimized with respect to its validation set performance by gradient descent. The data efficiency of gradient-based optimization, as opposed to inefficient black-box search, allows DARTS to achieve competitive performance with the state of the art using orders of magnitude less computation resources. It also outperforms another recent efficient architecture search method, ENAS (Pham et al., 2018b). Notably, DARTS is simpler than many existing approaches as it does not involve controllers (Zoph & Le, 2017; Baker et al., 2017; Zoph et al., 2018; Pham et al., 2018b; Zhong et al., 2018), hypernetworks (Brock et al., 2018) or performance predictors (Liu et al., 2018a), yet it is generic enough handle both convolutional and recurrent architectures.

The idea of searching architectures within a continuous domain is not new (Saxena & Verbeek, 2016; Ahmed & Torresani, 2017; Veniat & Denoyer, 2017; Shin et al., 2018), but there are several major

Current affiliation: Google Brain.

Published as a conference paper at ICLR 2019

distinctions. While prior works seek to fine-tune a specific aspect of an architecture, such as filter shapes or branching patterns in a convolutional network, DARTS is able to learn high-performance architecture building blocks with complex graph topologies within a rich search space. Moreover, DARTS is not restricted to any specific architecture family, and is applicable to both convolutional and recurrent networks.

In our experiments (Sect. 3) we show that DARTS is able to design a convolutional cell that achieves 2.76 0.09% test error on CIFAR-10 for image classification using 3.3M parameters, which is competitive with the state-of-the-art result by regularized evolution (Real et al., 2018) obtained using three orders of magnitude more computation resources. The same convolutional cell also achieves 26.7% top-1 error when transferred to Image Net (mobile setting), which is comparable to the best RL method (Zoph et al., 2018). On the language modeling task, DARTS efficiently discovers a recurrent cell that achieves 55.7 test perplexity on Penn Treebank (PTB), outperforming both extensively tuned LSTM (Melis et al., 2018) and all the existing automatically searched cells based on NAS (Zoph & Le, 2017) and ENAS (Pham et al., 2018b).

Our contributions can be summarized as follows:

We introduce a novel algorithm for differentiable network architecture search based on bilevel optimization, which is applicable to both convolutional and recurrent architectures.

Through extensive experiments on image classification and language modeling tasks we show that gradient-based architecture search achieves highly competitive results on CIFAR-10 and outperforms the state of the art on PTB. This is a very interesting result, considering that so far the best architecture search methods used non-differentiable search techniques, e.g. based on RL (Zoph et al., 2018) or evolution (Real et al., 2018; Liu et al., 2018b).

We achieve remarkable efficiency improvement (reducing the cost of architecture discovery to a few GPU days), which we attribute to the use of gradient-based optimization as opposed to non-differentiable search techniques.

We show that the architectures learned by DARTS on CIFAR-10 and PTB are transferable to Image Net and Wiki Text-2, respectively.

The implementation of DARTS is available at https://github.com/quark0/darts

2 DIFFERENTIABLE ARCHITECTURE SEARCH

We describe our search space in general form in Sect. 2.1, where the computation procedure for an architecture (or a cell in it) is represented as a directed acyclic graph. We then introduce a simple continuous relaxation scheme for our search space which leads to a differentiable learning objective for the joint optimization of the architecture and its weights (Sect. 2.2). Finally, we propose an approximation technique to make the algorithm computationally feasible and efficient (Sect. 2.3).

2.1 SEARCH SPACE

Following Zoph et al. (2018); Real et al. (2018); Liu et al. (2018a;b), we search for a computation cell as the building block of the final architecture. The learned cell could either be stacked to form a convolutional network or recursively connected to form a recurrent network.

A cell is a directed acyclic graph consisting of an ordered sequence of N nodes. Each node x(i) is a latent representation (e.g. a feature map in convolutional networks) and each directed edge (i, j) is associated with some operation o(i,j) that transforms x(i). We assume the cell to have two input nodes and a single output node. For convolutional cells, the input nodes are defined as the cell outputs in the previous two layers (Zoph et al., 2018). For recurrent cells, these are defined as the input at the current step and the state carried from the previous step. The output of the cell is obtained by applying a reduction operation (e.g. concatenation) to all the intermediate nodes.

Each intermediate node is computed based on all of its predecessors:

i<j o(i,j)(x(i)) (1)

Published as a conference paper at ICLR 2019

Figure 1: An overview of DARTS: (a) Operations on the edges are initially unknown. (b) Continuous relaxation of the search space by placing a mixture of candidate operations on each edge. (c) Joint optimization of the mixing probabilities and the network weights by solving a bilevel optimization problem. (d) Inducing the final architecture from the learned mixing probabilities.

A special zero operation is also included to indicate a lack of connection between two nodes. The task of learning the cell therefore reduces to learning the operations on its edges.

2.2 CONTINUOUS RELAXATION AND OPTIMIZATION

Let O be a set of candidate operations (e.g., convolution, max pooling, zero) where each operation represents some function o( ) to be applied to x(i). To make the search space continuous, we relax the categorical choice of a particular operation to a softmax over all possible operations:

o(i,j)(x) = X

exp(α(i,j) o ) P

o O exp(α(i,j) o ) o(x) (2)

where the operation mixing weights for a pair of nodes (i, j) are parameterized by a vector α(i,j) of dimension |O|. The task of architecture search then reduces to learning a set of continuous variables α = α(i,j) , as illustrated in Fig. 1. At the end of search, a discrete architecture can be obtained by

replacing each mixed operation o(i,j) with the most likely operation, i.e., o(i,j) = argmaxo O α(i,j) o . In the following, we refer to α as the (encoding of the) architecture.

After relaxation, our goal is to jointly learn the architecture α and the weights w within all the mixed operations (e.g. weights of the convolution filters). Analogous to architecture search using RL (Zoph & Le, 2017; Zoph et al., 2018; Pham et al., 2018b) or evolution (Liu et al., 2018b; Real et al., 2018) where the validation set performance is treated as the reward or fitness, DARTS aims to optimize the validation loss, but using gradient descent.

Denote by Ltrain and Lval the training and the validation loss, respectively. Both losses are determined not only by the architecture α, but also the weights w in the network. The goal for architecture search is to find α that minimizes the validation loss Lval(w , α ), where the weights w associated with the architecture are obtained by minimizing the training loss w = argminw Ltrain(w, α ).

This implies a bilevel optimization problem (Anandalingam & Friesz, 1992; Colson et al., 2007) with α as the upper-level variable and w as the lower-level variable:

min α Lval(w (α), α) (3)

s.t. w (α) = argminw Ltrain(w, α) (4)

The nested formulation also arises in gradient-based hyperparameter optimization (Maclaurin et al., 2015; Pedregosa, 2016; Franceschi et al., 2018), which is related in a sense that the architecture α could be viewed as a special type of hyperparameter, although its dimension is substantially higher than scalar-valued hyperparameters such as the learning rate, and it is harder to optimize.

Published as a conference paper at ICLR 2019

Algorithm 1: DARTS Differentiable Architecture Search

Create a mixed operation o(i,j) parametrized by α(i,j) for each edge (i, j) while not converged do

1. Update architecture α by descending αLval(w ξ w Ltrain(w, α), α) (ξ = 0 if using first-order approximation) 2. Update weights w by descending w Ltrain(w, α) Derive the final architecture based on the learned α.

2.3 APPROXIMATE ARCHITECTURE GRADIENT

Evaluating the architecture gradient exactly can be prohibitive due to the expensive inner optimization. We therefore propose a simple approximation scheme as follows:

αLval(w (α), α) (5) αLval(w ξ w Ltrain(w, α), α) (6)

where w denotes the current weights maintained by the algorithm, and ξ is the learning rate for a step of inner optimization. The idea is to approximate w (α) by adapting w using only a single training step, without solving the inner optimization (equation 4) completely by training until convergence. Related techniques have been used in meta-learning for model transfer (Finn et al., 2017), gradientbased hyperparameter tuning (Luketina et al., 2016) and unrolled generative adversarial networks (Metz et al., 2017). Note equation 6 will reduce to αLval(w, α) if w is already a local optimum for the inner optimization and thus w Ltrain(w, α) = 0.

The iterative procedure is outlined in Alg. 1. While we are not currently aware of the convergence guarantees for our optimization algorithm, in practice it is able to reach a fixed point with a suitable choice of ξ1. We also note that when momentum is enabled for weight optimisation, the one-step unrolled learning objective in equation 6 is modified accordingly and all of our analysis still applies.

Applying chain rule to the approximate architecture gradient (equation 6) yields

αLval(w , α) ξ 2 α,w Ltrain(w, α) w Lval(w , α) (7)

where w = w ξ w Ltrain(w, α) denotes the weights for a one-step forward model. The expression above contains an expensive matrix-vector product in its second term. Fortunately, the complexity can be substantially reduced using the finite difference approximation. Let ϵ be a small scalar2 and w = w ϵ w Lval(w , α). Then:

2 α,w Ltrain(w, α) w Lval(w , α) αLtrain(w+, α) αLtrain(w , α)

Evaluating the finite difference requires only two forward passes for the weights and two backward passes for α, and the complexity is reduced from O(|α||w|) to O(|α| + |w|).

First-order Approximation When ξ = 0, the second-order derivative in equation 7 will disappear. In this case, the architecture gradient is given by αLval(w, α), corresponding to the simple heuristic of optimizing the validation loss by assuming the current w is the same as w (α). This leads to some speed-up but empirically worse performance, according to our experimental results in Table 1 and Table 2. In the following, we refer to the case of ξ = 0 as the first-order approximation, and refer to the gradient formulation with ξ > 0 as the second-order approximation.

2.4 DERIVING DISCRETE ARCHITECTURES

To form each node in the discrete architecture, we retain the top-k strongest operations (from distinct nodes) among all non-zero candidate operations collected from all the previous nodes. The strength of an operation is defined as exp(α(i,j) o ) P

o O exp(α(i,j) o ). To make our derived architecture comparable with

1A simple working strategy is to set ξ equal to the learning rate for w s optimizer. 2We found ϵ = 0.01/ w Lval(w , α) 2 to be sufficiently accurate in all of our experiments.

Published as a conference paper at ICLR 2019

0.5 1.0 1.5 2.0 2.5 3.0 3.5 Architecture ( )

Weights (w)

= 0.7 = 0.5

Figure 2: Learning dynamics of our iterative algorithm when Lval(w, α) = αw 2α + 1 and Ltrain(w, α) = w2 2αw+α2, starting from (α(0), w(0)) = (2, 2). The analytical solution for the corresponding bilevel optimization problem is (α , w ) = (1, 1), which is highlighted in the red circle. The dashed red line indicates the feasible set where constraint equation 4 is satisfied exactly (namely, weights in w are optimal for the given architecture α). The example shows that a suitable choice of ξ helps to converge to a better local optimum.

those in the existing works, we use k = 2 for convolutional cells (Zoph et al., 2018; Liu et al., 2018a; Real et al., 2018) and k = 1 for recurrent cells (Pham et al., 2018b).

The zero operations are excluded in the above for two reasons. First, we need exactly k non-zero incoming edges per node for fair comparison with the existing models. Second, the strength of the zero operations is underdetermined, as increasing the logits of zero operations only affects the scale of the resulting node representations, and does not affect the final classification outcome due to the presence of batch normalization (Ioffe & Szegedy, 2015).

3 EXPERIMENTS AND RESULTS

Our experiments on CIFAR-10 and PTB consist of two stages, architecture search (Sect. 3.1) and architecture evaluation (Sect. 3.2). In the first stage, we search for the cell architectures using DARTS, and determine the best cells based on their validation performance. In the second stage, we use these cells to construct larger architectures, which we train from scratch and report their performance on the test set. We also investigate the transferability of the best cells learned on CIFAR-10 and PTB by evaluating them on Image Net and Wiki Text-2 (WT2) respectively.

3.1 ARCHITECTURE SEARCH

3.1.1 SEARCHING FOR CONVOLUTIONAL CELLS ON CIFAR-10

We include the following operations in O: 3 3 and 5 5 separable convolutions, 3 3 and 5 5 dilated separable convolutions, 3 3 max pooling, 3 3 average pooling, identity, and zero. All operations are of stride one (if applicable) and the convolved feature maps are padded to preserve their spatial resolution. We use the Re LU-Conv-BN order for convolutional operations, and each separable convolution is always applied twice (Zoph et al., 2018; Real et al., 2018; Liu et al., 2018a).

Our convolutional cell consists of N = 7 nodes, among which the output node is defined as the depthwise concatenation of all the intermediate nodes (input nodes excluded). The rest of the setup follows Zoph et al. (2018); Liu et al. (2018a); Real et al. (2018), where a network is then formed by stacking multiple cells together. The first and second nodes of cell k are set equal to the outputs of cell k 2 and cell k 1, respectively, and 1 1 convolutions are inserted as necessary. Cells located at the 1/3 and 2/3 of the total depth of the network are reduction cells, in which all the operations adjacent to the input nodes are of stride two. The architecture encoding therefore is (αnormal, αreduce), where αnormal is shared by all the normal cells and αreduce is shared by all the reduction cells.

Detailed experimental setup for this section can be found in Sect. A.1.1.

3.1.2 SEARCHING FOR RECURRENT CELLS ON PENN TREEBANK

Our set of available operations includes linear transformations followed by one of tanh, relu, sigmoid activations, as well as the identity mapping and the zero operation. The choice of these candidate operations follows Zoph & Le (2017); Pham et al. (2018b).

Our recurrent cell consists of N = 12 nodes. The very first intermediate node is obtained by linearly transforming the two input nodes, adding up the results and then passing through a tanh activation

Published as a conference paper at ICLR 2019

0 10 20 GPU Hours

Best Valid Error So Far (%)

DARTS (run 1) DARTS (run 2) DARTS (run 3) DARTS (run 4) Amoeba Net-A NASNet-A

0 5 10 15 20 25 GPU Hours

Valid Error (%)

DARTS Amoeba Net-A NASNet-A

0 1 2 3 4 GPU Hours

Best Valid Perplexity So Far

Penn Treebank

DARTS (run 1) DARTS (run 2) DARTS (run 3) DARTS (run 4) ENAS

0 1 2 3 4 GPU Hours

Valid Perplexity

Penn Treebank

Figure 3: Search progress of DARTS for convolutional cells on CIFAR-10 and recurrent cells on Penn Treebank. We keep track of the most recent architectures over time. Each architecture snapshot is re-trained from scratch using the training set (for 100 epochs on CIFAR-10 and for 300 epochs on PTB) and then evaluated on the validation set. For each task, we repeat the experiments for 4 times with different random seeds, and report the median and the best (per run) validation performance of the architectures over time. As references, we also report the results (under the same evaluation setup; with comparable number of parameters) of the best existing cells discovered using RL or evolution, including NASNet-A (Zoph et al., 2018) (2000 GPU days), Amoeba Net-A (3150 GPU days) (Real et al., 2018) and ENAS (0.5 GPU day) (Pham et al., 2018b).

function, as done in the ENAS cell (Pham et al., 2018b). The rest of the cell is learned. Other settings are similar to ENAS, where each operation is enhanced with a highway bypass (Zilly et al., 2016) and the cell output is defined as the average of all the intermediate nodes. As in ENAS, we enable batch normalization in each node to prevent gradient explosion during architecture search, and disable it during architecture evaluation. Our recurrent network consists of only a single cell, i.e. we do not assume any repetitive patterns within the recurrent architecture.

Detailed experimental setup for this section can be found in Sect. A.1.2.

c_{k-2} 0 sep_conv_3x3

1 sep_conv_3x3

2 skip_connect

3 skip_connect

sep_conv_3x3

sep_conv_3x3

sep_conv_3x3

dil_conv_3x3

Figure 4: Normal cell learned on CIFAR-10.

max_pool_3x3

2 max_pool_3x3

max_pool_3x3

max_pool_3x3

max_pool_3x3

skip_connect

skip_connect

skip_connect c_{k}

Figure 5: Reduction cell learned on CIFAR-10.

Figure 6: Recurrent cell learned on PTB.

3.2 ARCHITECTURE EVALUATION

To determine the architecture for final evaluation, we run DARTS four times with different random seeds and pick the best cell based on its validation performance obtained by training from scratch for

Published as a conference paper at ICLR 2019

a short period (100 epochs on CIFAR-10 and 300 epochs on PTB). This is particularly important for recurrent cells, as the optimization outcomes can be initialization-sensitive (Fig. 3).

To evaluate the selected architecture, we randomly initialize its weights (weights learned during the search process are discarded), train it from scratch, and report its performance on the test set. We note the test set is never used for architecture search or architecture selection.

Detailed experimental setup for architecture evaluation on CIFAR-10 and PTB can be found in Sect. A.2.1 and Sect. A.2.2, respectively. Besides CIFAR-10 and PTB, we further investigated the transferability of our best convolutional cell (searched on CIFAR-10) and recurrent cell (searched on PTB) by evaluating them on Image Net (mobile setting) and Wiki Text-2, respectively. More details of the transfer learning experiments can be found in Sect. A.2.3 and Sect. A.2.4.

Table 1: Comparison with state-of-the-art image classifiers on CIFAR-10 (lower error rate is better). Note the search cost for DARTS does not include the selection cost (1 GPU day) or the final evaluation cost by training the selected architecture from scratch (1.5 GPU days).

Architecture Test Error Params Search Cost #ops Search (%) (M) (GPU days) Method

Dense Net-BC (Huang et al., 2017) 3.46 25.6 manual

NASNet-A + cutout (Zoph et al., 2018) 2.65 3.3 2000 13 RL NASNet-A + cutout (Zoph et al., 2018) 2.83 3.1 2000 13 RL Block QNN (Zhong et al., 2018) 3.54 39.8 96 8 RL Amoeba Net-A (Real et al., 2018) 3.34 0.06 3.2 3150 19 evolution Amoeba Net-A + cutout (Real et al., 2018) 3.12 3.1 3150 19 evolution Amoeba Net-B + cutout (Real et al., 2018) 2.55 0.05 2.8 3150 19 evolution Hierarchical evolution (Liu et al., 2018b) 3.75 0.12 15.7 300 6 evolution PNAS (Liu et al., 2018a) 3.41 0.09 3.2 225 8 SMBO ENAS + cutout (Pham et al., 2018b) 2.89 4.6 0.5 6 RL ENAS + cutout (Pham et al., 2018b)* 2.91 4.2 4 6 RL

Random search baseline + cutout 3.29 0.15 3.2 4 7 random DARTS (first order) + cutout 3.00 0.14 3.3 1.5 7 gradient-based DARTS (second order) + cutout 2.76 0.09 3.3 4 7 gradient-based

* Obtained by repeating ENAS for 8 times using the code publicly released by the authors. The cell for final evaluation is chosen according to the same selection protocol as for DARTS. Obtained by training the corresponding architectures using our setup. Best architecture among 24 samples according to the validation error after 100 training epochs.

Table 2: Comparison with state-of-the-art language models on PTB (lower perplexity is better). Note the search cost for DARTS does not include the selection cost (1 GPU day) or the final evaluation cost by training the selected architecture from scratch (3 GPU days).

Architecture Perplexity Params Search Cost #ops Search valid test (M) (GPU days) Method

Variational RHN (Zilly et al., 2016) 67.9 65.4 23 manual LSTM (Merity et al., 2018) 60.7 58.8 24 manual LSTM + skip connections (Melis et al., 2018) 60.9 58.3 24 manual LSTM + 15 softmax experts (Yang et al., 2018) 58.1 56.0 22 manual

NAS (Zoph & Le, 2017) 64.0 25 1e4 CPU days 4 RL ENAS (Pham et al., 2018b)* 68.3 63.1 24 0.5 4 RL ENAS (Pham et al., 2018b) 60.8 58.6 24 0.5 4 RL

Random search baseline 61.8 59.4 23 2 4 random DARTS (first order) 60.2 57.6 23 0.5 4 gradient-based DARTS (second order) 58.1 55.7 23 1 4 gradient-based

* Obtained using the code (Pham et al., 2018a) publicly released by the authors. Obtained by training the corresponding architecture using our setup. Best architecture among 8 samples according to the validation perplexity after 300 training epochs.

Published as a conference paper at ICLR 2019

Table 3: Comparison with state-of-the-art image classifiers on Image Net in the mobile setting.

Architecture Test Error (%) Params + Search Cost Search top-1 top-5 (M) (M) (GPU days) Method

Inception-v1 (Szegedy et al., 2015) 30.2 10.1 6.6 1448 manual Mobile Net (Howard et al., 2017) 29.4 10.5 4.2 569 manual Shuffle Net 2 (g = 3) (Zhang et al., 2017) 26.3 5 524 manual

NASNet-A (Zoph et al., 2018) 26.0 8.4 5.3 564 2000 RL NASNet-B (Zoph et al., 2018) 27.2 8.7 5.3 488 2000 RL NASNet-C (Zoph et al., 2018) 27.5 9.0 4.9 558 2000 RL Amoeba Net-A (Real et al., 2018) 25.5 8.0 5.1 555 3150 evolution Amoeba Net-B (Real et al., 2018) 26.0 8.5 5.3 555 3150 evolution Amoeba Net-C (Real et al., 2018) 24.3 7.6 6.4 570 3150 evolution PNAS (Liu et al., 2018a) 25.8 8.1 5.1 588 225 SMBO

DARTS (searched on CIFAR-10) 26.7 8.7 4.7 574 4 gradient-based

3.3 RESULTS ANALYSIS

The CIFAR-10 results for convolutional architectures are presented in Table 1. Notably, DARTS achieved comparable results with the state of the art (Zoph et al., 2018; Real et al., 2018) while using three orders of magnitude less computation resources (i.e. 1.5 or 4 GPU days vs 2000 GPU days for NASNet and 3150 GPU days for Amoeba Net). Moreover, with slightly longer search time, DARTS outperformed ENAS (Pham et al., 2018b) by discovering cells with comparable error rates but less parameters. The longer search time is due to the fact that we have repeated the search process four times for cell selection. This practice is less important for convolutional cells however, because the performance of discovered architectures does not strongly depend on initialization (Fig. 3).

Alternative Optimization Strategies To better understand the necessity of bilevel optimization, we investigated a simplistic search strategy, where α and w are jointly optimized over the union of the training and validation sets using coordinate descent. The resulting best convolutional cell (out of 4 runs) yielded 4.16 0.16% test error using 3.1M parameters, which is worse than random search. In the second experiment, we optimized α simultaneously with w (without alteration) using SGD, again over all the data available (training + validation). The resulting best cell yielded 3.56 0.10% test error using 3.0M parameters. We hypothesize that these heuristics would cause α (analogous to hyperparameters) to overfit the training data, leading to poor generalization. Note that α is not directly optimized on the training set in DARTS.

Table 2 presents the results for recurrent architectures on PTB, where a cell discovered by DARTS achieved the test perplexity of 55.7. This is on par with the state-of-the-art model enhanced by a mixture of softmaxes (Yang et al., 2018), and better than all the rest of the architectures that are either manually or automatically discovered. Note that our automatically searched cell outperforms the extensively tuned LSTM (Melis et al., 2018), demonstrating the importance of architecture search in addition to hyperparameter search. In terms of efficiency, the overall cost (4 runs in total) is within 1 GPU day, which is comparable to ENAS and significantly faster than NAS (Zoph & Le, 2017).

It is also interesting to note that random search is competitive for both convolutional and recurrent models, which reflects the importance of the search space design. Nevertheless, with comparable or less search cost, DARTS is able to significantly improve upon random search in both cases (2.76 0.09 vs 3.29 0.15 on CIFAR-10; 55.7 vs 59.4 on PTB).

Results in Table 3 show that the cell learned on CIFAR-10 is indeed transferable to Image Net. It is worth noticing that DARTS achieves competitive performance with the state-of-the-art RL method (Zoph et al., 2018) while using three orders of magnitude less computation resources.

Table 4 shows that the cell identified by DARTS transfers to WT2 better than ENAS, although the overall results are less strong than those presented in Table 2 for PTB. The weaker transferability between PTB and WT2 (as compared to that between CIFAR-10 and Image Net) could be explained by the relatively small size of the source dataset (PTB) for architecture search. The issue of transferability could potentially be circumvented by directly optimizing the architecture on the task of interest.

Published as a conference paper at ICLR 2019

Table 4: Comparison with state-of-the-art language models on WT2.

Architecture Perplexity Params Search Cost Search valid test (M) (GPU days) Method

LSTM + augmented loss (Inan et al., 2017) 91.5 87.0 28 manual LSTM + continuous cache pointer (Grave et al., 2016) 68.9 manual LSTM (Merity et al., 2018) 69.1 66.0 33 manual LSTM + skip connections (Melis et al., 2018) 69.1 65.9 24 manual LSTM + 15 softmax experts (Yang et al., 2018) 66.0 63.3 33 manual

ENAS (Pham et al., 2018b) (searched on PTB) 72.4 70.4 33 0.5 RL

DARTS (searched on PTB) 71.2 69.6 33 1 gradient-based

Obtained by training the corresponding architecture using our setup.

4 CONCLUSION

We presented DARTS, a simple yet efficient architecture search algorithm for both convolutional and recurrent networks. By searching in a continuous space, DARTS is able to match or outperform the state-of-the-art non-differentiable architecture search methods on image classification and language modeling tasks with remarkable efficiency improvement by several orders of magnitude.

There are many interesting directions to improve DARTS further. For example, the current method may suffer from discrepancies between the continuous architecture encoding and the derived discrete architecture. This could be alleviated, e.g., by annealing the softmax temperature (with a suitable schedule) to enforce one-hot selection. It would also be interesting to explore performance-aware architecture derivation schemes based on the one-shot model learned during the search process.

ACKNOWLEDGEMENTS

The authors would like to thank Zihang Dai, Hieu Pham and Zico Kolter for useful discussions.

Karim Ahmed and Lorenzo Torresani. Connectivity learning in multi-branch networks. ar Xiv preprint ar Xiv:1709.09582, 2017.

G Anandalingam and TL Friesz. Hierarchical optimization: An introduction. Annals of Operations Research, 34(1):1 11, 1992.

Bowen Baker, Otkrist Gupta, Nikhil Naik, and Ramesh Raskar. Designing neural network architectures using reinforcement learning. ICLR, 2017.

Bowen Baker, Otkrist Gupta, Ramesh Raskar, and Nikhil Naik. Accelerating neural architecture search using performance prediction. ICLR Workshop, 2018.

Gabriel Bender, Pieter-Jan Kindermans, Barret Zoph, Vijay Vasudevan, and Quoc Le. Understanding and simplifying one-shot architecture search. In International Conference on Machine Learning, pp. 549 558, 2018.

Andrew Brock, Theodore Lim, James M Ritchie, and Nick Weston. Smash: one-shot model architecture search through hypernetworks. ICLR, 2018.

Han Cai, Tianyao Chen, Weinan Zhang, Yong Yu, and Jun Wang. Efficient architecture search by network transformation. AAAI, 2018.

Benoît Colson, Patrice Marcotte, and Gilles Savard. An overview of bilevel optimization. Annals of operations research, 153(1):235 256, 2007.

Terrance De Vries and Graham W Taylor. Improved regularization of convolutional neural networks with cutout. ar Xiv preprint ar Xiv:1708.04552, 2017.

Published as a conference paper at ICLR 2019

Thomas Elsken, Jan-Hendrik Metzen, and Frank Hutter. Simple and efficient architecture search for convolutional neural networks. ar Xiv preprint ar Xiv:1711.04528, 2017.

Chelsea Finn, Pieter Abbeel, and Sergey Levine. Model-agnostic meta-learning for fast adaptation of deep networks. In ICML, pp. 1126 1135, 2017.

Luca Franceschi, Paolo Frasconi, Saverio Salzo, and Massimilano Pontil. Bilevel programming for hyperparameter optimization and meta-learning. ICML, 2018.

Yarin Gal and Zoubin Ghahramani. A theoretically grounded application of dropout in recurrent neural networks. In NIPS, pp. 1019 1027, 2016.

Edouard Grave, Armand Joulin, and Nicolas Usunier. Improving neural language models with a continuous cache. ar Xiv preprint ar Xiv:1612.04426, 2016.

Andrew G Howard, Menglong Zhu, Bo Chen, Dmitry Kalenichenko, Weijun Wang, Tobias Weyand, Marco Andreetto, and Hartwig Adam. Mobilenets: Efficient convolutional neural networks for mobile vision applications. ar Xiv preprint ar Xiv:1704.04861, 2017.

Gao Huang, Zhuang Liu, Kilian Q Weinberger, and Laurens van der Maaten. Densely connected convolutional networks. In CVPR, 2017.

Hakan Inan, Khashayar Khosravi, and Richard Socher. Tying word vectors and word classifiers: A loss framework for language modeling. ICLR, 2017.

Sergey Ioffe and Christian Szegedy. Batch normalization: Accelerating deep network training by reducing internal covariate shift. ar Xiv preprint ar Xiv:1502.03167, 2015.

Kirthevasan Kandasamy, Willie Neiswanger, Jeff Schneider, Barnabas Poczos, and Eric Xing. Neural architecture search with bayesian optimisation and optimal transport. NIPS, 2018.

Diederik P Kingma and Jimmy Ba. Adam: A method for stochastic optimization. ar Xiv preprint ar Xiv:1412.6980, 2014.

Ben Krause, Emmanuel Kahembwe, Iain Murray, and Steve Renals. Dynamic evaluation of neural sequence models. ar Xiv preprint ar Xiv:1709.07432, 2017.

Chenxi Liu, Barret Zoph, Jonathon Shlens, Wei Hua, Li-Jia Li, Li Fei-Fei, Alan Yuille, Jonathan Huang, and Kevin Murphy. Progressive neural architecture search. ECCV, 2018a.

Hanxiao Liu, Karen Simonyan, Oriol Vinyals, Chrisantha Fernando, and Koray Kavukcuoglu. Hierarchical representations for efficient architecture search. ICLR, 2018b.

Ilya Loshchilov and Frank Hutter. Sgdr: Stochastic gradient descent with warm restarts. ar Xiv preprint ar Xiv:1608.03983, 2016.

Jelena Luketina, Mathias Berglund, Klaus Greff, and Tapani Raiko. Scalable gradient-based tuning of continuous regularization hyperparameters. In ICML, pp. 2952 2960, 2016.

Dougal Maclaurin, David Duvenaud, and Ryan Adams. Gradient-based hyperparameter optimization through reversible learning. In ICML, pp. 2113 2122, 2015.

Gábor Melis, Chris Dyer, and Phil Blunsom. On the state of the art of evaluation in neural language models. ICLR, 2018.

Stephen Merity, Nitish Shirish Keskar, and Richard Socher. Regularizing and optimizing lstm language models. ICLR, 2018.

Luke Metz, Ben Poole, David Pfau, and Jascha Sohl-Dickstein. Unrolled generative adversarial networks. ICLR, 2017.

Renato Negrinho and Geoff Gordon. Deeparchitect: Automatically designing and training deep architectures. ar Xiv preprint ar Xiv:1704.08792, 2017.

Published as a conference paper at ICLR 2019

Adam Paszke, Sam Gross, Soumith Chintala, Gregory Chanan, Edward Yang, Zachary De Vito, Zeming Lin, Alban Desmaison, Luca Antiga, and Adam Lerer. Automatic differentiation in pytorch. In NIPS-W, 2017.

Fabian Pedregosa. Hyperparameter optimization with approximate gradient. In ICML, 2016.

Hieu Pham, Melody Y Guan, Barret Zoph, Quoc V Le, and Jeff Dean. Authors implementation of Efficient Neural Architecture Search via Parameter Sharing . https://github.com/ melodyguan/enas/tree/2734eb2657847f090e1bc5c51c2b9cbf0be51887, 2018a. Accessed: 2018-04-05.

Hieu Pham, Melody Y Guan, Barret Zoph, Quoc V Le, and Jeff Dean. Efficient neural architecture search via parameter sharing. ICML, 2018b.

Boris T Polyak and Anatoli B Juditsky. Acceleration of stochastic approximation by averaging. SIAM Journal on Control and Optimization, 30(4):838 855, 1992.

Esteban Real, Alok Aggarwal, Yanping Huang, and Quoc V Le. Regularized evolution for image classifier architecture search. ar Xiv preprint ar Xiv:1802.01548, 2018.

Shreyas Saxena and Jakob Verbeek. Convolutional neural fabrics. In NIPS, pp. 4053 4061, 2016.

Richard Shin, Charles Packer, and Dawn Song. Differentiable neural network architecture search. In ICLR Workshop, 2018.

Christian Szegedy, Wei Liu, Yangqing Jia, Pierre Sermanet, Scott Reed, Dragomir Anguelov, Dumitru Erhan, Vincent Vanhoucke, Andrew Rabinovich, Jen-Hao Rick Chang, et al. Going deeper with convolutions. In CVPR, 2015.

Tom Veniat and Ludovic Denoyer. Learning time/memory-efficient deep architectures with budgeted super networks. ar Xiv preprint ar Xiv:1706.00046, 2017.

Zhilin Yang, Zihang Dai, Ruslan Salakhutdinov, and William W Cohen. Breaking the softmax bottleneck: a high-rank rnn language model. ICLR, 2018.

Xiangyu Zhang, Xinyu Zhou, Mengxiao Lin, and Jian Sun. Shufflenet: An extremely efficient convolutional neural network for mobile devices. ar Xiv preprint ar Xiv:1707.01083, 2017.

Zhao Zhong, Junjie Yan, Wei Wu, Jing Shao, and Cheng-Lin Liu. Practical block-wise neural network architecture generation. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, pp. 2423 2432, 2018.

Julian Georg Zilly, Rupesh Kumar Srivastava, Jan Koutník, and Jürgen Schmidhuber. Recurrent highway networks. ar Xiv preprint ar Xiv:1607.03474, 2016.

Barret Zoph and Quoc V Le. Neural architecture search with reinforcement learning. ICLR, 2017.

Barret Zoph, Vijay Vasudevan, Jonathon Shlens, and Quoc V Le. Learning transferable architectures for scalable image recognition. CVPR, 2018.

Published as a conference paper at ICLR 2019

A EXPERIMENTAL DETAILS

A.1 ARCHITECTURE SEARCH

A.1.1 CIFAR-10

Since the architecture will be varying throughout the search process, we always use batch-specific statistics for batch normalization rather than the global moving average. Learnable affine parameters in all batch normalizations are disabled during the search process to avoid rescaling the outputs of the candidate operations.

To carry out architecture search, we hold out half of the CIFAR-10 training data as the validation set. A small network of 8 cells is trained using DARTS for 50 epochs, with batch size 64 (for both the training and validation sets) and the initial number of channels 16. The numbers were chosen to ensure the network can fit into a single GPU. We use momentum SGD to optimize the weights w, with initial learning rate ηw = 0.025 (annealed down to zero following a cosine schedule without restart (Loshchilov & Hutter, 2016)), momentum 0.9, and weight decay 3 10 4. We use zero initialization for architecture variables (the α s in both the normal and reduction cells), which implies equal amount of attention (after taking the softmax) over all possible ops. At the early stage this ensures weights in every candidate op to receive sufficient learning signal (more exploration). We use Adam (Kingma & Ba, 2014) as the optimizer for α, with initial learning rate ηα = 3 10 4, momentum β = (0.5, 0.999) and weight decay 10 3. The search takes one day on a single GPU3.

For architecture search, both the embedding and the hidden sizes are set to 300. The linear transformation parameters across all incoming operations connected to the same node are shared (their shapes are all 300 300), as the algorithm always has the option to focus on one of the predecessors and mask away the others. Tying the weights leads to memory savings and faster computation, allowing us to train the continuous architecture using a single GPU. Learnable affine parameters in batch normalizations are disabled, as we did for convolutional cells. The network is then trained for 50 epochs using SGD without momentum, with learning rate ηw = 20, batch size 256, BPTT length 35, and weight decay 5 10 7. We apply variational dropout (Gal & Ghahramani, 2016) of 0.2 to word embeddings, 0.75 to the cell input, and 0.25 to all the hidden nodes. A dropout of 0.75 is also applied to the output layer. Other training settings are identical to those in Merity et al. (2018); Yang et al. (2018). Similarly to the convolutional architectures, we use Adam for the optimization of α (initialized as zeros), with initial learning rate ηα = 3 10 3, momentum β = (0.9, 0.999) and weight decay 10 3. The search takes 6 hours on a single GPU.

A.2 ARCHITECTURE EVALUATION

A.2.1 CIFAR-10

A large network of 20 cells is trained for 600 epochs with batch size 96. The initial number of channels is increased from 16 to 36 to ensure our model size is comparable with other baselines in the literature (around 3M). Other hyperparameters remain the same as the ones used for architecture search. Following existing works (Pham et al., 2018b; Zoph et al., 2018; Liu et al., 2018a; Real et al., 2018), additional enhancements include cutout (De Vries & Taylor, 2017), path dropout of probability 0.2 and auxiliary towers with weight 0.4. The training takes 1.5 days on a single GPU with our implementation in Py Torch (Paszke et al., 2017). Since the CIFAR results are subject to high variance even with exactly the same setup (Liu et al., 2018b), we report the mean and standard deviation of 10 independent runs for our full model.

To avoid any discrepancy between different implementations or training settings (e.g. the batch sizes), we incorporated the NASNet-A cell (Zoph et al., 2018) and the Amoeba Net-A cell (Real et al., 2018) into our training framework and reported their results under the same settings as our cells.

3All of our experiments were performed using NVIDIA GTX 1080Ti GPUs.

Published as a conference paper at ICLR 2019

A single-layer recurrent network with the discovered cell is trained until convergence with batch size 64 using averaged SGD (Polyak & Juditsky, 1992) (ASGD), with learning rate ηw = 20 and weight decay 8 10 7. To speedup, we start with SGD and trigger ASGD using the same protocol as in Yang et al. (2018); Merity et al. (2018). Both the embedding and the hidden sizes are set to 850 to ensure our model size is comparable with other baselines. The token-wise dropout on the embedding layer is set to 0.1. Other hyperparameters remain exactly the same as those for architecture search. For fair comparison, we do not finetune our model at the end of the optimization, nor do we use any additional enhancements such as dynamic evaluation (Krause et al., 2017) or continuous cache (Grave et al., 2016). The training takes 3 days on a single 1080Ti GPU with our Py Torch implementation. To account for implementation discrepancies, we also incorporated the ENAS cell (Pham et al., 2018b) into our codebase and trained their network under the same setup as our discovered cells.

A.2.3 IMAGENET

We consider the mobile setting where the input image size is 224 224 and the number of multiply-add operations in the model is restricted to be less than 600M.

A network of 14 cells is trained for 250 epochs with batch size 128, weight decay 3 10 5 and initial SGD learning rate 0.1 (decayed by a factor of 0.97 after each epoch). Other hyperparameters follow Zoph et al. (2018); Real et al. (2018); Liu et al. (2018a)4. The training takes 12 days on a single GPU.

A.2.4 WIKITEXT-2

We use embedding and hidden sizes 700, weight decay 5 10 7, and hidden-node variational dropout 0.15. Other hyperparameters remain the same as in our PTB experiments.

B SEARCH WITH INCREASED DEPTH

To better understand the effect of depth for architecture search, we conducted architecture search on CIFAR-10 by increasing the number of cells in the stack from 8 to 20. The initial number of channels is reduced from 16 to 6 due to memory budget of a single GPU. All the other hyperparameters remain the same. The search cost doubles and the resulting cell achieves 2.88 0.09% test error, which is slightly worse than 2.76 0.09% obtained using a shallower network. This particular setup may have suffered from the enlarged discrepancy of the number of channels between architecture search and final evaluation. Moreover, searching with a deeper model might require different hyperparameters due to the increased number of layers to back-prop through.

C COMPLEXITY ANALYSIS

In this section, we analyze the complexity of our search space for convolutional cells.

Each of our discretized cell allows Q4 k=1 (k+1)k

2 (72) 109 possible DAGs without considering graph isomorphism (recall we have 7 non-zero ops, 2 input nodes, 4 intermediate nodes with 2 predecessors each). Since we are jointly learning both normal and reduction cells, the total number of architectures is approximately (109)2 = 1018. This is greater than the 5.6 1014 of PNAS (Liu et al., 2018a) which learns only a single type of cell.

Also note that we retained the top-2 predecessors per node only at the very end, and our continuous search space before this final discretization step is even larger. Specifically, each relaxed cell (a fully connected graph) contains 2 + 3 + 4 + 5 = 14 learnable edges, allowing (7 + 1)14 4 1012 possible configurations (+1 to include the zero op indicating a lack of connection). Again, since we are learning both normal and reduction cells, the total number of architectures covered by the continuous space before discretization is (4 1012)2 1025.

4We did not conduct extensive hyperparameter tuning.