U NDERSTANDING AND R OBUSTIFYING D IFFERENTIABLE A RCHITECTURE S EARCH
Anonymous authors
Paper under double-blind review
A BSTRACT
Differentiable Architecture Search (DARTS) has attracted a lot of attention due to its simplicity and small search costs achieved by a continuous relaxation and an approximation of the resulting bi level optimization problem
However,
DARTS does not work robustly for new problems: we identify a wide range of search spaces for which DARTS yields degenerate architectures with very poor test performance.
We study this failure mode and show that, while DARTS successfully minimizes validation loss, the found solutions generalize poorly when they coincide with high validation loss curvature in the architecture space.
We show that by adding one of various types of regularization we can robustify DARTS to ﬁnd solutions with less curvature and better generalization properties.
Based on these observations, we propose several simple variations of DARTS that perform substantially more robustly in practice.
Our observations are robust across ﬁve search spaces on three image classiﬁcation tasks and also hold for the very different domains of disparity estimation (a dense regression task) and language modelling
I NTRODUCTION
Neural Architecture Search (NAS), the process of automatically designing neural network architectures, has recently attracted attention by achieving state-of-the-art performance on a variety of tasks (Zoph & Le, 2017; Real et al., 2019).
Differentiable architecture search (DARTS)
(Liu et al., 2019) signiﬁcantly improved the efﬁciency of NAS over prior work, reducing its costs to the same order of magnitude as training a single neural network.
This expanded the scope of NAS substantially, allowing it to also be applied on more expensive problems, such as semantic segmentation (Chenxi et al., 2019) or disparity estimation (Saikia et al., 2019).
However, several researchers have also reported DARTS to not work well, in some cases even no better than random search (Li & Talwalkar, 2019; Sciuto et al., 2019).
Why is this?
How can these seemingly contradicting results be explained?
The overall goal of this paper is to understand and overcome such failure modes of DARTS.
To this end, we make the following contributions:
We identify 12 NAS benchmarks based on four search spaces in which standard DARTS yields degenerate architectures with poor test performance across several datasets (Section 3).
By computing the eigenspectrum of the Hessian of the validation loss with respect to the architectural parameters, we show that there is a strong correlation between its dominant eigenvalue and the architecture’s generalization error.
Based on this ﬁnding, we propose a simple variation of DARTS with early stopping that performs substantially more robustly (Section 4).
We show that, related to previous work on sharp/ﬂat local minima, regularizing the inner objective of DARTS more strongly allows it to ﬁnd solutions with smaller Hessian spectrum and better generalization properties.
Based on these insights, we propose two practical robustiﬁcations of DARTS that overcome its failure modes in all our 12 NAS benchmarks (Section 5).
Our ﬁndings are robust across a wide range of NAS benchmarks based on image recognition and also hold for the very different domains of language modelling (PTB) and disparity estimation.
They consolidate the ﬁndings of the various results in the literature and lead to a substantially more robust version of DARTS.
We provide our implementation and scripts to facilitate reproducibility 1 .
B ACKGROUND AND R ELATED W ORK 2. R ELATION BETWEEN FLAT / SHARP MINIMA AND GENERALIZATION PERFORMANCE
Already Hochreiter & Schmidhuber (1997) observed that ﬂat minima of the training loss yield better generalization performance than sharp minima.
Recent work (Keskar et al., 2016; Yao et al., 2018) focuses more on the settings of large/small batch size training, where observations show that small batch training tends to get attracted to ﬂatter minima and generalizes better.
Similarly, Nguyen et al.
(2018) observed that this phenomenon manifests also in the hyperparameter space.
They showed that whenever the hyperparameters overﬁt the validation data, the minima lie in a sharper region of the space.
This motivated us to conduct a similar analysis in the context of differentiable architecture search later in Section 4.1, where we see the same effect in the space of neural network architectures.
2. B I - LEVEL O PTIMIZATION
We start by a short introduction of the bi-level optimization problem (Colson et al., 2007).
These are problems which contain two optimization tasks, nested within each other.
Deﬁnition 2.1.
Given the outer objective function F : R P × R N → R and the inner objective function f : R P × R N → R , the bi-level optimization problem is given by
[Equation]
where y ∈ R P and θ ∈ R N are the outer and inner variables, respectively.
One may also see the bi-level problem as a constrained optimization problem, with the inner problem as a constraint.
In general, even in the case when the inner objective (2) is strongly convex and has an unique minimizer θ ∗ ( y ) =
arg min θ ∈ R N f
( y, θ ) , it is not possible to directly optimize the outer objective (1).
A possible method around this issue is to use the implicit function theorem to retrieve the derivative of the solution map (or response map) θ ∗ (
y ) ∈ F ⊆ R N w.r.t.
y (Bengio, 2000; Pedregosa, 2016; Beirami et al., 2017).
Another strategy is to approximate the inner problem with a dynamical system (Domke, 2012; Maclaurin et al., 2015; Franceschi et al., 2017; 2018), where the optimization dynamics could, e.g., describe gradient descent.
In the case that the minimizer of the inner problem is unique, under some conditions the set of minimizers of this approximate problem will indeed converge to the minimizers of the bilevel problem (1) (see Franceschi et al.
2. N EURAL A RCHITECTURE S EARCH
Neural Architecture Search (NAS) denotes the process of automatically designing neural network architectures in order to overcome the cumbersome trial-and-error process when designing architectures manually.
We brieﬂy review NAS here and refer to the recent survey by Elsken et al.
(2019b)
for a more thorough overview.
Prior work mostly employs either reinforcement learning techniques (Baker et al., 2017a;
Zoph & Le, 2017; Zhong et al., 2018; Zoph et al., 2018) or evolutionary algorithms (Stanley & Miikkulainen, 2002; Liu et al., 2018b; Miikkulainen et al., 2017;
Real et al., 2017; 2019) to optimize the discrete architecture space.
As these methods are often very expensive, various works focus on reducing the search costs by, e.g., employing network morphisms (Cai et al., 2018a;b; Elsken et al., 2017; 2019a), weight sharing within search models (Saxena & Verbeek, 2016; Bender et al., 2018; Pham et al., 2018) or multi-ﬁdelity optimization (Baker et al., 2017b;
Falkner et al., 2018; Li et al., 2017; Zela et al., 2018), but their applicability still often remains restricted to rather simple tasks and small datasets.
2. D IFFERENTIABLE A RCHITECTURE S EARCH (DARTS)
A recent line of work focuses on relaxing the discrete neural architecture search problem to a continuous one that can be solved by gradient descent (Liu et al., 2019; Xie et al., 2019; Casale et al., 2019;
Cai et al., 2019).
In DARTS (Liu et al., 2019), this is achieved by simply using a weighted sum of possible candidate operations for each layer, whereas the real-valued weights then effectively parametrize the network’s architecture.
We will now review DARTS in more detail, as our work builds directly upon it.
Continuous relaxation of the search space.
In agreement with prior work (Zoph et al., 2018;
Real et al., 2019), DARTS optimizes only substructures called cells that are stacked to deﬁne the full network architecture.
Each cell contains N nodes organized in a directed acyclic graph.
The graph contains two inputs nodes (given by the outputs of the previous two cells), a set of intermediate nodes, and one output node (given by concatenating all intermediate nodes).
Each intermediate node x ( j ) represents a feature map.
See Figure 1 for an illustration of such a cell.
Instead of applying a single operation to a speciﬁc node during architecture search, Liu et al.
(2019) relax the decision which operation to choose by computing the intermediate node as a mixture of candidate operations, applied to predecessor nodes x ( i ) , i < j , x ( j )
= (cid:80) i<
j (cid:80) o ∈O
[Equation]
,jo
[Equation]
jo
[Equation]
x
[Equation] [Equation]
, where O denotes the set of all candidate operations (e.g., 3 × 3 convolution, skip connection, 3 × 3 max pooling, etc.)
and α = ( α i,jo ) i,j,o serves as a real-valued parameterization of the architecture.
Gradient-based optimization of the search space.
DARTS then optimizes both the weights of the search network (often called the weight-sharing or one-shot model, since the weights of all individual subgraphs/architectures are shared) and architectural parameters by alternating gradient descent.
The network weights and the architecture parameters are optimized on the training and validation set, respectively.
This can be interpreted as solving the bi-level optimization problem (1), (2), where F and f are the validation and training loss, L valid and L train , respectively, while y and θ denote the architectural parameters α and network weights w , respectively.
Note that DARTS only approximates the lower-level solution by a single gradient step (see Appendix A for more details).
At the end of the search phase, a discrete cell is obtained by choosing the k most important incoming operations for each intermediate node while all others are pruned.
Importance is measured by the operation weighting factor
[Equation]
W HEN DARTS FAILS
We now describe various search spaces and demonstrate that standard DARTS fails on them.
We start with four search spaces similar to the original CIFAR-10 search space but simpler, and evaluate across three different datasets (CIFAR-10, CIFAR-100 and SVHN).
They are quite standard in that they use the same macro architecture as the original DARTS paper (Liu et al., 2018a), consisting of normal and reduction cells; however, they only allow a subset of operators for the cell search space:
S1: This search space uses a different set of only two operators per edge, which we identiﬁed using an ofﬂine process that iteratively dropped the operations from the original DARTS search space with the least importance.
This pre-optimized space has the advantage of being quite small while still including many strong architectures.
We refer to Appendix
B
for details on its construction and an illustration (Figure 8).
S2:
In this space, the set of candidate operations per edge is { 3 × 3 SepConv , SkipConnect } .
We choose these operations since they are the most frequent ones in the discovered cells reported by Liu et al.
S3:
In this space, the set of candidate operations per edge is { 3 × 3 SepConv , SkipConnect , Zero } , where the Zero operation simply replaces every value in the input feature map by zeros.
S4:
In this space, the set of candidate operations per edge is { 3 × 3 SepConv , Noise } , where the Noise operation simply replaces every value from the input feature map by noise (cid:15) ∼ N (0 , 1) .
This is
the only space out of S1-S4 that is not a strict subspace of the original DARTS space; we intentionally added the Noise operation, which actively harms performance and should therefore not be selected by DARTS.
We ran DARTS on each of these spaces, using exactly the same setup as Liu et al.
Figure 1 shows the poor cells DARTS selected on these search spaces for CIFAR-10 (see Appendix G for c_{k-2}
skip_connect c_{k-1} skip_connect 1 skip_connect 2 skip_connect 3 skip_connect
skip_connect
skip_connect
skip_connect c_{k} (a) Spacec_{k-2}
skip_connect 2 skip_connect c_{k-1} skip_connect 1 skip_connect skip_connect
skip_connect 3 skip_connect c_{k} sep_conv_3x (b) Spacec_{k-2}
skip_connect 1 skip_connect 2 skip_connect 3 skip_connect c_{k-1} skip_connect skip_connect
skip_connect
skip_connect c_{k} (c) Spacec_{k-2} 0 sep_conv_3x3sep_conv_3x3 2 noise 3 noise c_{k-1} sep_conv_3x3 noise noise c_{k} noise (d) Space analogous results on the other datasets).
Already visually, one suspect that the found cells are suboptimal: the parameter-less skip connections dominate in almost all the edges for spaces S1-S3, and for S4 even the harmful Noise operation was selected for ﬁve out of eight operations.
Table 1 (ﬁrst column) conﬁrms the very poor performance standard DARTS yields on all of these search spaces and on different datasets.
We note that Liu et al.
and Xie et al.
(2019) argue that the
Zero operation can help to search for the architecture topology and choice of operators jointly, but in our experiments it did not help to reduce the importance weight of the skip connection (compare
Figure 1b vs. Figure 1c).
We emphasize that search spaces S1-S3 are very natural, and, as strict subspaces of the original space, should merely be easier to search than that.
They are in no way special or constructed in an adversarial manner.
Only S4 was constructed speciﬁcally to show-case the failure mode of DARTS selecting the obviously suboptimal Noise operator.
S5: Very small search space with known global optimum.
Knowing the global minimum has the advantage that one can benchmark the performance of algorithms by measuring the regret of chosen points with respect to the known global minimum.
Therefore, we created another search space with only one intermediate node for both normal and reduction cells, and 3 operation choices in each edge, namely 3 × 3 SepConv , SkipConnection , and 3 × 3 MaxPooling .
The total number of possible architectures in this space is 81, all of which we evaluated a-priori.
We dub this space S5 .
[Figure]
We ran DARTS on this search space three times for each dataset and compared its result to the baseline of Random Search with weight sharing (RS-ws) by Li & Talwalkar (2019).
Figure 2 shows the test regret of the architectures selected by DARTS (blue) and RS-ws (green) throughout the search.
DARTS manages to ﬁnd an architecture close to the global minimum, but around epoch 40 the test performance deteriorated.
Note that the search model validation error (dashed red line) did not deteriorate but rather converged, indicating that the architectural parameters are overﬁtting to the validation set.
In contrast, RS-ws stays relatively constant throughout the search; when evaluating only the ﬁnal architecture found, RS-ws indeed outperformed DARTS.
S6: encoder-decoder architecture for disparity estimation.
To study whether our ﬁndings generalize beyond image recognition, we also analyzed a search space for a very different problem: ﬁnding encoder-decoder architectures for the dense regression task of disparity estimation; please refer to Appendix E for details.
We base this search space on AutoDispNet (Saikia et al., 2019), which used DARTS for a space containing normal , downsampling and upsampling cells.
We again constructed a reduced space.
Similarly to the image classiﬁcation search spaces, we found the normal cell to be mainly composed of parameter-less operations (see Figure 23 in Appendix G).
As expected, this causes a large generalization error (see ﬁrst row in Table 2 of our later experiments).
[Figure] [Figure]
T HE R OLE OF D OMINANT E IGENVALUES OF ∇ 2 α L valid
We now analyze why DARTS fails in all these cases.
Motivated by Section 2.1, we will have a closer look at the largest eigenvalue λ αmax of the Hessian matrix of validation loss ∇ 2 α L valid w
[Figure]
r
[Figure]
t. the architectural parameters α .
4. L ARGE ARCHITECTURAL EIGENVALUES AND GENERALIZATION PERFORMANCE
One may hypothesize that DARTS performs poorly because its approximate solution of the bi-level optimization problem by iterative optimization fails, but we actually observe validation errors to progress nicely: Figure 3 (left) shows that the search model validation error converges in all cases, even though the cell structures selected here are the ones in Figure 1.
Rather, the architectures DARTS ﬁnds do not generalize well.
This can be seen in Figure 3 (middle).
There, every 5 epochs, we evaluated the architecture deemed by DARTS to be optimal according to the α values.
Note that whenever evaluating on the test set, we retrain from scratch the architecture obtained after applying the argmax to the architectural weights α .
As one can notice, the architectures start to degenerate after a certain number of search epochs, similarly to the results shown in Figure 2.
We hypothesized that this might be related to sharp local minima as discussed in
Section 2.1.
To test this hypothesis, we computed the full Hessian ∇ 2 α L valid of the validation loss w.r.t.
the architectural parameters on a randomly sampled mini-batch.
Figure 3 (right) shows that the dominant eigenvalue λ αmax (which serves as a proxy for the sharpness) indeed increases in standard DARTS, along with the test error (middle) of the ﬁnal architectures, while the validation error still decreases (left).
We also studied the correlation between λ αmax and test error more directly, by measuring these two quantities for 24 different architectures (obtained via standard DARTS and the regularized versions we discuss in Section 5).
For the example of space S1 on CIFAR-10, Figureshows that λ α max indeed strongly correlates with test error (with a Pearson correlation coefﬁcient of 0.867).
4. L ARGE ARCHITECTURAL EIGENVALUES AND PERFORMANCE DROP AFTER PRUNING
[Figure]
model vs. the sharpness of minima (by means of λ αmax )
One reason why DARTS performs poorly when the architectural eigenvalues are large (and thus the minimum is sharp) might be the pruning step at the end of DARTS: the optimal, continuous α ∗ from the search is pruned to obtain a discrete α ∗ d , somewhere in the neighbourhood of α ∗ .
In the case of a sharp minimum α ∗ , α ∗ d might have a loss function value signiﬁcantly higher than the minimum α ∗ , while in the case of a ﬂat minimum, α ∗ d is expected to have a similar loss function value.
To investigate this hypothesis, we measured the performance drop w.r.t.
to the search model weights incurred by this discretization step and correlated it with λ αmax .
The results in Figure 5 show that, indeed, low curvature never led to large performance drops.
Having identiﬁed this relationship, we now move on to avoid high curvature.
[Figure]
4. E ARLY STOPPING BASED ON LARGE EIGENVALUES OF ∇ 2 α L valid
We propose a simple early stopping methods to avoid large curvature and thus poor generalization.
We emphasize that simply stopping the search based on validation performance (as one would do in the case of training a network) does not apply here as NAS directly optimizes validation performance, which – as we have seen in Figure 2 – keeps on improving.
[Figure]
[Table]
[Figure]
[Table]
Instead, we propose to track λ αmax over the course of architecture search and stop whenever it increases too much.
To implement this idea, we use a simple heuristic that worked off-the-shelf without any tuning.
Let λ αmax ( i ) denote the value of λ αmax smoothed over k =epochs around i ; then, we stop if λ αmax ( i − k ) /λ αmax ( i ) < 0 .
75 and return the architecture from epoch
i − k .
By this early stopping heuristic, we do not only avoid exploding eigenvalues, which are correlated with poor generalization (see Figure 4), but also shorten the time of the search.
Table 1 shows the results for running DARTS with this early stopping criterion (DARTS-ES) across S1-S4 and all three image classiﬁcation datasets.
Note that we never use the test data when applying the early stopping mechanism.
Early stopping signiﬁcantly improved DARTS for all settings without ever harming it.
R EGULARIZATION OF INNER OBJECTIVE IMPROVES GENERALIZATION OF ARCHITECTURES
As we saw in Section 4.1, sharper minima (by means of large eigenvalues) of the validation loss lead to poor generalization performance.
In our bi-level optimization setting, the outer variables’ trajectory depends on the inner optimization procedure.
Therefore, we hypothesized that modifying the landscape of the inner objective L train could redirect the outer variables α to ﬂatter areas of the architectural space.
We study two ways of regularization (data augmentation in Section 5.1 and Lregularization in Section 5.2) and ﬁnd that both, along with the early stopping criterion from Section 4.3, make DARTS more robust in practice.
We emphasize that we do not alter the regularization of the ﬁnal training and evaluation phase, but solely that of the search phase.
The setting we use for all experiments in this paper to obtain the ﬁnal test performance is described in Appendix C.
5. R EGULARIZATION VIA DATA AUGMENTATION
We ﬁrst investigate the effect of regularizing via data augmentation, namely masking out parts of the input and intermediate feature maps via Cutout (CO, DeVries & Taylor (2017)) and ScheduledDropPath (DP, Zoph et al.
(ScheduledDropPath is a regularization technique, but we list it here since we apply it together with Cutout), respectively, during architecture search.
We ran
DARTS with CO and DP (with and without our early stopping criterion, DARTS-ES) with different maximum DP probabilities on all three image classiﬁcation datasets and search spaces S1-S4.
[Figure] [Figure]
Figure 6 summarizes the results: regularization improves the test performance of DARTS and DARTSES in all cases, sometimes very substantially, and at the same time kept the dominant eigenvalue relatively low (Figure 12).
This also directly results in smaller drops in accuracy after pruning, as discussed in Section 4.2; indeed, the search runs plotted in Figure 5 are the same as in this section.
Figure 16 in the appendix explicitly shows how regularization relates to the accuracy drops.
We also refer to further results in the appendix: Figure 10 (showing test vs. validation error) and Table 5 (showing that overﬁtting of the architectural parameters is reduced).
Similar observations hold for disparity estimation on
S6, where we vary the strength of standard data augmentation methods, such as shearing or brightness change, rather then masking parts of features, which is unreasonable for this task.
The augmentation strength is described by an “augmentation scaling factor” (Appendix E).
Table 2 summarizes the results.
We report the average end point error (EPE), which is the Euclidean distance between the predicted and ground truth disparity maps.
Data augmentation avoided the degenerate architectures and substantially improved results.
5. I NCREASED L 2 REGULARIZATION
As a second type of regularization, we also tested different L 2 regularization factors 3 i · 10 − 4 for i ∈ { 1 , 3 , 9 , 27 , 81 } .
Standard DARTS in fact does already include a small amount of L 2 regularization; i = 1 yields its default.
Figure 7 shows that DARTS’ test performance (solid lines) can be signiﬁcantly improved by higher L 2 factors across all datasets and spaces, while keeping the dominant eigenvalue low (Figure 13).
DARTS with early stopping (dashed lines) also beneﬁts from additional regularization.
Again, we observe the implicit regularization effect on the outer objective which reduces the overﬁtting of the architectural parameters.
We again refer to Table 2 for disparity estimation; Appendix F shows similar results for language modelling (Penn TreeBank).
5. P RACTICAL R OBUSTIFICATION OF DARTS BY R EGULARIZING THE I NNER O BJECTIVE
Based on the insights from the aforementioned analysis and empirical results, we now propose two alternative simple modiﬁcations to make DARTS more robust in practice without having to manually tune its regularization hyperparameters.
DARTS with adaptive regularization One option is to adapt DARTS’ regularization hyperparameters in an automated way, in order to keep the architectural weights in areas of the validation loss objective with smaller curvature.
The simplest off-the-shelf procedure towards this desiderata would be to increase the regularization strength whenever the dominant eigenvalue starts increasing rapidly.
Algorithm 1 (DARTS-ADA, Appendix D.1) shows such a procedure.
We use the same stopping criterion as in DARTS-ES (Section 4.3), roll back DARTS to the epoch when this criterion is met, and continue the search with a larger regularization value R for the remaining epochs (larger by a factor of η ).
This procedure is repeated whenever the criterion is met, unless the regularization value exceeds some maximum predeﬁned value R max .
Multiple DARTS runs with different regularization strength Liu et al.
(2019) already suggested to run the search phase of DARTS four times, resulting in four architectures, and to return the best of these four architectures w.r.t.
validation performance when retrained from scratch for a limited number of epochs.
We propose to use the same procedure, with the only difference that the four runs use different amounts of regularization.
The resulting RobustDARTS (R-DARTS) method is conceptually very simple, trivial to implement and likely to work well if any of the tried regularization strengths works well.
[Table]
Table 3 evaluates the performance of our practical robustiﬁcations of
DARTS,
DARTS-ADA and
R-DARTS (based on either L2 or
ScheduledDropPath regularization), by comparing them to the original
DARTS, DARTS-ES and Random
Search with weight sharing (RS-ws).
For each of these methods, as proposed in the DARTS paper (Liu et al., 2019), we ran the search four independent times with different random seeds and selected the architecture used for the ﬁnal evaluation based on a validation run as described above.
As the table shows, in accordance with Li & Talwalkar (2019), RS-ws often outperformed the original DARTS; however, with our robustiﬁcations, DARTS typically performs substantially better than RS-ws.
DARTS-ADA consistently improved over standard DARTS for all benchmarks, indicating that a gradual increase of regularization during search prevents ending up in the bad regions of the architectural space.
Finally, RobustDARTS yielded the best performance and since it is also easier to implement than DARTS-ES and DARTS-ADA, it is the method that we recommend to be used in practice.
search spaces.
We show mean ± stddev for 5 repetitions (based on 4 fresh subruns each as in Table 3); for the more expensive PTB we could only afford 1 such repetition
[Table]
Finally, since the evaluations in this paper have so far focussed on smaller subspaces of the original DARTS search space, the reader may wonder how well RobustDARTS works on the full search spaces.
As Table 4 shows, RobustDARTS performed similarly to DARTS for the two original benchmarks from the
DARTS paper (PTB and CIFAR-10), on which DARTS was developed and is well tuned; however, even when only changing the dataset to CIFAR-100 or SVHN, RobustDARTS already performed signiﬁcantly better than DARTS, underlining its robustness.
C ONCLUSIONS
We showed that the generalization performance of architectures found by DARTS is related to the eigenvalues of the Hessian matrix of the validation loss w.r.t.
the architectural parameters.
Standard DARTS often results in degenerate architectures with large eigenvalues and poor generalization.
Based on this observation, we proposed a simple early stopping criterion for DARTS based on tracking the largest eigenvalue.
Our empirical results also show that properly regularizing the inner objective helps controlling the eigenvalue and therefore improves generalization.
Our ﬁndings substantially improve our understanding of DARTS’ failure modes and lead to much more robust versions.
They are consistent across many different search spaces on image recognition tasks and also for the very different domains of language modelling and disparity estimation.
Our code is available for reproducibility.
A M ORE DETAIL ON DARTS
Here we present a detailed description of DARTS architectural update steps.
We ﬁrstly provide the general formalism which computes the gradient of the outer level problem in (1) by means of the implicit function theorem.
Afterwards, we present how DARTS computes the gradient used to update the architectural parameters α .
A. D ERIVATIVE WITH SMOOTHED NON - QUADRATIC LOWER LEVEL PROBLEM
Consider the general deﬁnition of the bi-level optimization problem as given by (1) and (2).
Given that f is twice continuously differentiable and that all stationary points are local minimas, one can make use of the implicit function theorem to ﬁnd the derivative of the solution map θ ∗ ( y ) w.r.t.
y (Bengio, 2000).
Under the smoothness assumption, the optimality condition of the lower level (2) is ∇ θ f
( y, θ ) = 0 , which deﬁnes an implicit function θ ∗ ( y ) .
With the assumption that min θ f ( y, θ ) has a solution, there exists a ( y, θ ∗ ) such that ∇ θ f ( y, θ ∗ ) = 0 .
Under the condition that ∇ θ f ( y, θ ∗ ) = 0 is continuously differentiable and that θ ∗ ( y ) is continuously differentiable at y , implicitly differentiating the last equality from both sides w.r.t.
y and applying the chain rule, yields:
[Equation]
Assuming that the Hessian ∇ 2 θ f
( y, θ ∗ ) is invertible, we can rewrite (3) as follows:
[Equation]
Applying the chain rule to (1) for computing the total derivative of F with respect to y yields:
[Equation]
where we have omitted the evaluation at ( y, θ ∗ ) .
Substituting (4) into (5) and reordering yields:
[Equation]
equation 6 computes the gradient of F , given the function θ ∗ ( y ) , which maps outer variables to the inner variables minimizing the inner problem.
However, in most of the cases obtaining such a mapping is computationally expensive, therefore different heuristics have been proposed to approximate dF/dy (Maclaurin et al., 2015; Pedregosa, 2016; Franceschi et al., 2017; 2018).
A. DARTS ARCHITECTURAL GRADIENT COMPUTATION
DARTS optimization procedure is deﬁned as a bi-level optimization problem where L valid is the outer objective (1) and L train is the inner objective (2):
[Equation]
where both losses are determined by both the architecture parameters α (outer variables) and the network weights w (inner variables).
Based on Appendix A.1, under some conditions, the total derivative of L valid w.r.t.
α evaluated on ( α, w ∗ ( α )) would be:
[Equation]
where ∇ α = ∂∂α , ∇ w = ∂∂w and ∇ 2 α,w = ∂ 2 ∂α∂w .
Computing the inverse of the Hessian is in general not possible considering the high dimensionality of the model parameters w , therefore resolving to gradient-based iterative algorithms for ﬁnding w ∗ is necessary.
However, this would also require to optimize the model parameters w till convergence each time α is updated.
If our model is a deep neural network it is clear that this computation is expensive, therefore Liu et al.
(2019) propose to approximate w ∗ ( α ) by updating the current model parameters w using a single gradient descent step:
[Equation]
where ξ is the learning rate for the virtual gradient step DARTS takes with respect to the model weights w .
From equation 10 the gradient of w ∗ ( α ) with respect to α is
[Equation]
By setting the evaluation point w ∗ =
w − ξ ∇ w L train ( α, w ) and following the same derivation as in Appendix A.1, we obtain the DARTS architectural gradient approximation:
[Equation]
where the inverse Hessian ∇ 2 w L − 1 train in (9) is replaced by the learning rate ξ .
This expression however contains again an expensive vector-matrix product.
Liu et al.
(2019) reduce the complexity by using the ﬁnite difference approximation around w ± =
w ± (cid:15) ∇ w L valid ( α, w ∗ ) for some small (cid:15)
= 0 .
(cid:107)∇ w L valid ( α, w ∗ ) (cid:107) 2 to compute the gradient of ∇ α L train ( α, w ∗ ) with respect to w as
[Equation]
In the end, combining equation 12 and equation 13 gives the gradient to compute the architectural updates in DARTS:
[Equation]
In all our experiments we always use ξ = η (also called second order approximation in Liu et al.
(2019)), where η is the learning rate used in SGD for updating the parameters w .
B C ONSTRUCTION OF S1 FROM S ECTION
We ran DARTS two times on the default search space to ﬁnd the two most important operations per mixed operation.
Initially, every mixed operation consists of 8 operations.
After the ﬁrst DARTS run, we drop the 4 (out of 8) least important ones.
In the second DARTS run, we drop the 2 (out of the remaining 4) least important ones.
S1 is then deﬁned to contain only the two remaining most important operations per mixed op.
Refer to Figure 8 for an illustration of this pre-optimized space.
C F INAL A RCHITECTURE E VALUATION
Similar to the original DARTS paper (Liu et al., 2019), the architecture found during the search are scaled up by increasing the number of ﬁlters and cells and retrained from scratch to obtain the ﬁnal test performance.
For CIFAR-100 and SVHN we use 16 number of initial ﬁlters and 8 cells when training architectures from scratch for all the experiments we conduct.
The rest of the settings is the same as in Liu et al.
On CIFAR-10, when scaling the ScheduledDropPath drop probability, we use the same settings for training from scratch the found architectures as in the original DARTS paper, i.e. 36 initial ﬁlters and 20 stacked cells.
However, for search space S2 and S4 we reduce the number of initial ﬁlters to 16 in order to avoid memory issues, since the cells found with more regularization usually are composed only with separable convolutions.
When scaling the L 2 factor on CIFAR-10 experiments we use 16 initial ﬁlters and 8 stacked cells, except the experiments on S1, where the settings are the same as in Liu et al.
i.e. 36 initial ﬁlters and 20 stacked cells.
[Figure]
Note that although altering the regularization factors during DARTS search
[Figure]
when training the ﬁnal architectures from scratch we always use the same values for them as in Liu et al.
(2019), i.e.
ScheduledDropPath maximum drop probability linearly increases from 0 towards 0.2 throughout training, Cutout is always enabled with cutout probability 1.0, and the L 2 regularization factor is set to 3 ·
10 − 4 .
D A DDITIONAL EMPIRICAL RESULTS
[Figure]
t r
[Figure]
g
[Figure]
t ( % )
[Figure]
[Table]
S 0. S 0.
[Table] [Table]
[Table] [Table]
S
[Table] [Table]
D. A DAPTIVE DARTS DETAILS
We evaluated DARTS-ADA (Section 5.3) with R = 3 · 10 − 4 (DARTS default), R max = 3 · 10 −and η
= 10 on all the search spaces and datasets we use for image classiﬁcation.
The results are shown in Table 3 (DARTS-ADA).
The function train and eval conducts the normal DARTS search for one epoch and returns the architecture at the end of that epoch’s updates and the stop value if a decision was made to stop the search and rollback to stop epoch .
/* E: epochs to search; R : initial regularization value; R max : maximal regularization value; stop criter: stopping criterion; η : regularization increase factor */ Input : E, R , R max , stop criter, η /* start search for E epochs */ for epoch in E do /* run DARTS for one epoch and return stop=
True together with the stop epoch */ /* and the architecture at stop epoch if the criterion is met */ stop, stop epoch, arch ← train and eval(stop criter) ; if stop & R ≤ R max then /* start DARTS from stop epoch with a larger R */ arch ← DARTS ADA(E - stop epoch, η · R , R max , stop criter, η ) ; break end end
Output: arch
[Figure] [Figure]
[Figure] [Figure]
[Figure] [Figure]
[Figure]
Max. Drop probability
[Figure]
V a li d a t i o n a cc u r a c y d r o p ( % )
[Figure]
[Table]
SVHN
S
[Table] [Table]
[Table] [Table]
[Table] [Table]
[Table] [Table]
S
[Table] [Table]
[Table] [Table]
[Table] [Table]
[Table] [Table]
[Table] [Table]
S
[Table]
89 ±
[Table] [Table]
00 ±
[Table] [Table]
78 ±
[Table] [Table]
58 ±
[Table]
S
[Table]
42 ±
[Table] [Table]
90 ±
[Table] [Table]
55 ±
[Table] [Table]
52 ±
[Table]
D. A CLOSER LOOK AT THE EIGENVALUES
Over the course of all experiments from the paper, we tracked the largest eigenvalue across all conﬁguration and datasets to see how they evolve during the search.
Figures 12 and 13 shows the results across all the settings for image classiﬁcation.
It can be clearly seen that increasing the inner objective regularization, both in terms of L 2 or data augmentation, helps controlling the largest eigenvalue and keeping it to a small value, which again helps explaining why the architectures found with stronger regularization generalize better.
The markers on each line highlight the epochs where DARTS is early stopped.
As one can see from Figure 4, there is indeed some correlation between the average dominant eigenvalue throughout the search and the test performance of the found architectures by DARTS.
Figures 14 and 15 (top 3 rows) show the full spectrum (sorted based on eigenvalue absolute values) at the end of search, whilst bottom 3 rows plot the distribution of eigenvalues in the eigenspectrum.
As one can see, not only the dominant eigenvalue is larger compared to the cases when the regularization is stronger and the generalization of architectures is better, but also the other eigenvalues in the spectrum have larger absolute value, indicating a sharper objective landscape towards many dimensions.
Furthermore, from the distribution plots note the presence of more negative eigenvalues whenever the architectures are degenerate (lower regularization value) indicating that DARTS gets stuck in a point with larger positive and negative curvature of the validation loss objective, associated with a more degenerate Hessian matrix.
E D ISPARITY E STIMATION E. D ATASETS
We use the FlyingThings3D dataset (Mayer et al., 2016) for training AutoDispNet.
It consists of rendered stereo image pairs and their ground truth disparity maps.
The dataset provides a training and testing split consisting of 21 , 818 and 4248 samples respectively with an image resolution of 960 × 540 .
We use the Sintel dataset ( Butler et al.
for testing our networks.
Sintel is another synthetic dataset from derived from an animated movie which also provides ground truth disparity maps ( 1064 samples) with a resolution of 1024 × 436 .
E. T RAINING
We use the AutoDispNet-C architecture as described in Saikia et al.
However, we use the smaller search which consists of three operations: MaxPool 3 × 3 , SepConv 3 × 3 , and
SkipConnect .
For training the search network, images are downsampled by a factor of two and trained for 300 k mini-batch iterations.
During search, we use SGD and ADAM to optimize the inner and outer objectives respectively.
Differently from the original AutoDispNet we do not warmstart the search model weights before starting the architectural parameter updates.
The extracted network is also trained for 300 k mini-batch iterations but full resolution images are used.
Here,
ADAM is used for optimization and the learning rate is annealed to 0 from 1 e − 4 , using a cosine decay schedule.
E. E FFECT OF REGULARIZATION ON THE INNER OBJECTIVE
To study the effect of regularization on the inner objective for AutoDispNet-C we use experiment with two types of regularization: data augmentation and of L 2 regularization on network weights.
We note that we could not test the early stopping method on AutoDispNet since AutoDispNet relies on custom operations to compute feature map correlation (Dosovitskiy et al., 2015) and resampling, for which second order derivatives are currently not available (which are required to compute the
Hessian).
Data augmentation.
Inspite of fairly large number of training samples in FlyingThings3D, data augmentation is crucial for good generalization performance.
Disparity estimation networks employ spatial transformations such as translation, cropping, shearing and scaling.
Additionally, appearance transformations such as additive Gaussian noise, changes in brightness, contrast, gamma and color are also applied.
Parameters for such transformations are sampled from a uniform or Gaussian distribution (parameterized by a mean and variance).
In our experiments, we vary the data augmentation strength by multiplying the variance of these parameter distributions by a ﬁxed factor, which we dub the augmentation scaling factor .
The extracted networks are evaluated with the same augmentation parameters.
The results of increasing the augmentation strength of the inner objective can be seen in Table 2.
We observe that as augmentation strength increases DARTS ﬁnds networks with more number of parameters and better test performance.
The best test performance is obtained for the network with maximum augmentation for the inner objective.
At the same time the search model validation error increases when scaling up the augmentation factor, which again enforces the argument that the overﬁtting of architectural parameters is reduced by this implicit regularizer.
L2 regularization.
We study the effect of increasing regularization strength on the weights of the network.
The results are shown in Table 2.
Also in this case best test performance is obtained with the maximum regularization strength.
F R ESULTS ON P ENN T REEBANK
Here we investigate the effect of more L 2 regularization on the inner objective for searching recurrent cells on Penn Treebank (PTB).
We again used a reduced search space with only ReLU and identity mapping as possible operations.
The rest of the settings is the same as in (Liu et al., 2019).
We run DARTS search four independent times with different random seeds, each with four Lregularization factors, namely 5 × 10 − 7 (DARTS default), 15 × 10 − 7 , 45 × 10 − 7 and 135 × 10 − 7 .
Figure 17 shows the test perplexity of the architectures found by DARTS with the aforementioned
L 2 regularization values.
As we can see, a stronger regularization factor on the inner objective makes the search procedure more robust.
The median perplexity of the discovered architectures gets better as we increase the L 2 factor from 5 × 10 − 7 to 45 × 10 − 7 , while the search model (one-shot) validation mean perplexity increases.
This observation is similar to the ones on image classiﬁcation shown in
Figure 9, showing again that properly regularizing the inner objective helps reduce overﬁtting the architectural parameters.
[Figure]
i
[Figure]
G D ISCOVERED CELLS ON SEARCH SPACES S1-S
[Figure] [Figure]
by DARTS on CIFAR-100 and SVHN when ran with its default hyperparameters on spaces S1-S4.
Notice the dominance of parameter-less operations such as skip connection and pooling ops.
c_{k-2} 0 max_pool_3x
[Table]
c_{k-1} max
[Table]
pool
[Table]
3x
[Table]
sep
[Table]
c_{k-2} 0 avg_pool
[Table] [Table]
{k
[Table] [Table]
[Figure]
[Table]
[Figure]
[Table]
[Figure] [Figure]
[Figure]
k
[Figure] [Figure]
[Table]
_3x 3 sep_conv_3x c_{k} sep_conv_3x3 sep_conv_3x (j) S2 (SVHN) c_{k-2}
sep_conv_3x 1 sep_conv_3x3 c_{k-1} sep_conv_3x 2 sep_conv_3x3 sep_conv_3x c_{k} sep_conv_3x 3 sep_conv_3x3 sep_conv_3x (k) S3 (SVHN) c_{k-2} sep_conv_3x 1 sep_conv_3x3sep_conv_3x3sep_conv_3x c_{k-1} sep_conv_3x sep_conv_3x3 sep_conv_3x sep_conv_3x c_{k} (l) S4 (SVHN)
[Figure]
[Figure]
_3x 1 skip_connect c_{k-1} skip_connect sep_conv_3x 2 sep_conv_3x c_{k} sep_conv_3x (
[Figure]
Upsampling cell c_{k-2} 0 sep_conv_3x 1 sep_conv_3x c_{k-1} sep_conv_3x sep_conv_3x 2 sep_conv_3x c_{k} sep_conv_3x (a) Normal cell c_{k-2} 0 sep_conv_3x 1 sep_conv_3x c_{k-1} skip_connect sep_conv_3x 2 sep_conv_3x c_{k} sep_conv_3x (b
[Figure]
Reduction cell c_skip 0 skip_connect 1 skip_connect 2 skip_connect c_pred_{k-1} sep_conv_3x3 c_{k-2} sep_conv_3x c_{k-1} sep_conv_3x sep_conv_3x3 sep_conv_3x c_{k} (
[Figure]
Upsampling cell c_{k-2} 0 sep_conv_3x  sep_conv_3x c_{k-1} sep_conv_3x 1 sep_conv_3x3 sep_conv_3x
c_{k} sep_conv_3x (a) Normal cell c_{k-2}
skip_connect c_{k-1} max_pool_3x 1 skip_connect sep_conv_3x 2 sep_conv_3x c_{k} sep_conv_3x
[Figure]
cell c_skip 0 sep_conv_3x 1 sep_conv_3x3 c_pred_{k-1} sep_conv_3x c_{k-2}
skip_connect c_{k-1} sep_conv_3x 2 sep_conv_3x
sep_conv_3x c_{k
} sep_conv_3x
[Figure]
[Table]
S
[Table]
[Figure]
[Figure]
2 max_pool_3x 3 max_pool_3x c_{k-1} max_pool_3x skip_connect
{+skip_connect c_{k} skip_connect (e) S c_{k-2} 
+}
{+skip_connect 1 skip_connect c_{k-1} sep_conv_3x skip_connect 2 skip_connect 3 skip_connect c_{k} skip_connect skip_connect (f) S c_{k-2} 0 skip_connect 3 skip_connect c_{k-1} sep_conv_3x 1 sep_conv_3x 2 sep_conv_3x3sep_conv_3x3skip_connect c_{k} skip_connect (g) S c_{k-2} 0 sep_conv_3x 1 sep_conv_3x 2 sep_conv_3x 3 sep_conv_3x c_{k-1} sep_conv_3x3 sep_conv_3x3 sep_conv_3x sep_conv_3x
+}
c_{k} (h) S
H M ORE REGULARIZATION IN THE INNER OBJECTIVE
[Figure]