summaries: The paper studies the non-convex optimization problem of training CNNs with ReLU activations under different choices for the CNN architecture, and shows how these can be framed as convex problems with a poly time complexity w.r.t. relevant variables. The derived convex problems provide valuable insights on how the CNN architecture induces different weight regularizers by giving them in explicit form -- these show a rich connection between the architecture and regularizer. [Summary] This paper focuses on training convolutional neural networks (CNNs) by using convex optimization techniques. By taking the dual of the nonconvex training problems, (and the dual of its dual), the main contribution of the paper is to show the strong duality between the convex problem and its original nonconvex training problems. This result has been proved for multi-layer CNNs with one ReLU layer and three-layer CNNs with two ReLU layers. The paper considers several types of CNN and proposes convex reformulations for non-convex problems of training these networks. As a result a polynomial complexity is shown for the training problem. The results are also interpreted as implicit regularization induced by the choice of the architecture. Finally, numerical experiments are made to support the theoretical findings and show that in the predicted regime, SGD for the original problem converges to the global minimizer given by the convex reformulation.
abstractText: We study training of Convolutional Neural Networks (CNNs) with ReLU activations and introduce exact convex optimization formulations with a polynomial complexity with respect to the number of data samples, the number of neurons, and data dimension. More specifically, we develop a convex analytic framework utilizing semi-infinite duality to obtain equivalent convex optimization problems for several twoand three-layer CNN architectures. We first prove that two-layer CNNs can be globally optimized via an `2 norm regularized convex program. We then show that multi-layer circular CNN training problems with a single ReLU layer are equivalent to an `1 regularized convex program that encourages sparsity in the spectral domain. We also extend these results to three-layer CNNs with two ReLU layers. Furthermore, we present extensions of our approach to different pooling methods, which elucidates the implicit architectural bias as convex regularizers.
authors.name: POLYNOMIAL TIME
authors.name: Tolga Ergen
authors.name: Mert Pilanci
references.authors: Akshay Agrawal Robin Verschueren Steven Diamond Stephen Boyd
references.title: A rewriting system for convex optimization problems
references.venue: Journal of Control and Decision,
references.authors: Francis Bach
references.title: Breaking the curse of dimensionality with convex neural networks
references.venue: The Journal of Machine Learning Research,
references.authors: Eugene Belilovsky Michael Eickenberg Edouard Oyallon
references.title: Greedy layerwise learning can scale to imagenet
references.venue: In International Conference on Machine Learning, pp. 583–593,
references.authors: Yoshua Bengio Nicolas L Roux Pascal Vincent Olivier Delalleau Patrice Marcotte
references.title: Convex neural networks. In Advances in neural information processing
references.authors: Guy Blanc Neha Gupta Gregory Valiant Paul Valiant
references.title: Implicit regularization for deep neural networks driven by an ornstein-uhlenbeck like process
references.venue: CoRR, abs/1904.09080,
references.authors: Digvijay Boob Santanu S Dey Guanghui Lan
references.title: Complexity of training relu neural network
references.venue: arXiv preprint arXiv:1809.10787,
references.authors: Thomas M Cover
references.title: Geometrical and statistical properties of systems of linear inequalities with applications in pattern recognition
references.venue: IEEE transactions on electronic computers,
references.authors: Steven Diamond Stephen Boyd
references.title: CVXPY: A Python-embedded modeling language for convex optimization
references.venue: Journal of Machine Learning Research,
references.authors: Herbert Edelsbrunner Joseph O’Rourke Raimund Seidel
references.title: Constructing arrangements of lines and hyperplanes with applications
references.venue: SIAM Journal on Computing,
references.authors: Tolga Ergen Mert Pilanci
references.title: Convex duality and cutting plane methods for over-parameterized neural networks
references.venue: In OPT-ML workshop,
references.authors: Tolga Ergen Mert Pilanci
references.title: Convex programs for global optimization of convolutional neural networks in polynomial-time
references.venue: In OPT-ML workshop,
references.authors: Tolga Ergen Mert Pilanci
references.title: Revealing the structure of deep neural networks via convex duality
references.venue: arXiv preprint arXiv:2002.11219,
references.authors: Tolga Ergen Mert Pilanci
references.title: Convex geometry and duality of over-parameterized neural networks
references.venue: arXiv preprint arXiv:2002.11219,
references.authors: Tolga Ergen Mert Pilanci
references.title: Convex geometry of two-layer relu networks: Implicit autoencoding and interpretable models
references.venue: Proceedings of the Twenty Third International Conference on Artificial Intelligence and Statistics,
references.authors: Tolga Ergen Arda Sahiner Batu Ozturkler John Pauly Morteza Mardani Mert Pilanci
references.title: Demystifying batch normalization in relu networks: Equivalent convex optimization models and implicit regularization
references.authors: Robert Geirhos Patricia Rubisch Claudio Michaelis Matthias Bethge Felix A Wichmann Wieland Brendel
references.title: Imagenet-trained cnns are biased towards texture; increasing shape bias improves accuracy and robustness
references.venue: In International Conference on Learning Representations,
references.authors: Michael Grant Stephen Boyd
references.title: CVX: Matlab software for disciplined convex programming, version 2.1
references.authors: Suriya Gunasekar Jason D Lee Daniel Soudry Nati Srebro
references.title: Implicit bias of gradient descent on linear convolutional networks
references.venue: Advances in Neural Information Processing Systems
references.authors: Hartmut Maennel Olivier Bousquet Sylvain Gelly
references.title: Gradient descent quantizes relu network features
references.venue: arXiv preprint arXiv:1803.08367,
references.authors: L. Meier S. van de Geer P. Bühlmann
references.title: The group lasso for logistic regression
references.venue: Journal of the Royal Statistical Society, Series B,
references.authors: Behnam Neyshabur Ryota Tomioka Nathan Srebro
references.title: In search of the real inductive bias: On the role of implicit regularization in deep learning
references.venue: arXiv preprint arXiv:1412.6614,
references.authors: Piyush C Ojha
references.title: Enumeration of linear threshold functions from the lattice of hyperplane intersections
references.venue: IEEE Transactions on Neural Networks,
references.authors: Rahul Parhi Robert D. Nowak
references.title: Minimum "norm" neural networks are splines, 2019
references.authors: Mert Pilanci Tolga Ergen
references.title: Neural networks are convex regularizers: Exact polynomial-time convex optimization formulations for two-layer networks, 2020
references.authors: Nasim Rahaman Aristide Baratin Devansh Arpit Felix Draxler Min Lin Fred Hamprecht Yoshua Bengio Aaron Courville
references.title: On the spectral bias of neural networks
references.venue: In International Conference on Machine Learning,
references.authors: Benjamin Recht Maryam Fazel Pablo A Parrilo
references.title: Guaranteed minimum-rank solutions of linear matrix equations via nuclear norm minimization
references.venue: SIAM review,
references.authors: Saharon Rosset Grzegorz Swirszcz Nathan Srebro Ji Zhu
references.title: L1 regularization in infinite dimensional feature spaces
references.venue: In International Conference on Computational Learning Theory,
references.authors: Walter Rudin
references.title: Principles of Mathematical Analysis
references.authors: Arda Sahiner Tolga Ergen John M. Pauly Mert Pilanci
references.title: Vector-output relu neural network problems are copositive programs: Convex analysis of two layer networks and polynomial-time algorithms
references.venue: In International Conference on Learning Representations,
references.authors: Pedro Savarese Itay Evron Daniel Soudry Nathan Srebro
references.title: How do infinite width bounded norm networks look in function space
references.venue: CoRR, abs/1902.05040,
references.authors: Alexander Shapiro
references.title: Semi-infinite programming, duality, discretization and optimality conditions
references.authors: Karen Simonyan Andrew Zisserman
references.title: Very deep convolutional networks for large-scale image recognition
references.venue: arXiv preprint arXiv:1409.1556,
references.authors: Maurice Sion
references.title: On general minimax theorems
references.venue: Pacific J. Math.,
references.authors: Richard P Stanley
references.title: An introduction to hyperplane arrangements
references.venue: Geometric combinatorics,
references.authors: RH Tütüncü KC Toh MJ Todd
references.title: Sdpt3—a matlab software package for semidefinite-quadraticlinear programming, version 3.0
references.venue: Web page http://www. math. nus. edu. sg/mattohkc/sdpt3. html,
references.authors: RO Winder
references.title: Partitions of n-space by hyperplanes
references.venue: SIAM Journal on Applied Mathematics,
references.authors: Ming Yuan Yi Lin
references.title: Model selection and estimation in regression with grouped variables
references.venue: Journal of the Royal Statistical Society: Series B (Statistical Methodology),
references.authors: Chiyuan Zhang Samy Bengio Moritz Hardt Benjamin Recht Oriol Vinyals
references.title: Understanding deep learning requires rethinking generalization
references.venue: arXiv preprint arXiv:1611.03530,
sections.heading: 1 INTRODUCTION
sections.text: Convolutional Neural Networks (CNNs) have shown a remarkable success across various machine learning problems (LeCun et al., 2015). However, our theoretical understanding of CNNs still remains restricted, where the main challenge arises from the highly non-convex and nonlinear structure of CNNs with nonlinear activations such as ReLU. Hence, we study the training problem for various CNN architectures with ReLU activations and introduce equivalent finite dimensional convex formulations that can be used to globally optimize these architectures. Our results characterize the role of network architecture in terms of equivalent convex regularizers. Remarkably, we prove that the proposed methods are polynomial time with respect to all problem parameters.
  Convex neural network training was previously considered in Bengio et al. (2006); Bach (2017). However, these studies are restricted to two-layer fully connected networks with infinite width, thus, the optimization problem involves infinite dimensional variables. Moreover, it has been shown that even adding a single neuron to a neural network leads to a non-convex optimization problem which cannot be solved efficiently (Bach, 2017). Another line of research in Parhi & Nowak (2019); Ergen & Pilanci (2019; 2020a;b;c;d); Pilanci & Ergen (2020); Savarese et al. (2019); Gunasekar et al. (2018); Maennel et al. (2018); Blanc et al. (2019); Zhang et al. (2016) focuses on the effect of implicit and explicit regularization in neural network training and aims to explain why the resulting network generalizes well. Among these studies, Parhi & Nowak (2019); Ergen & Pilanci (2020b;c;d); Savarese et al. (2019) proved that the minimum `2 norm two-layer network that perfectly fits a one dimensional dataset outputs the linear spline interpolation. Moreover, Gunasekar et al. (2018) studied certain linear convolutional networks and revealed an implicit non-convex quasi-norm regularization. However, as the number of layers increases, the regularization approaches to `0 quasi-norm, which is not computationally tractable. Recently, Pilanci & Ergen (2020) showed that two-layer CNNs with linear activations can be equivalently optimized as nuclear and `1 norm regularized convex problems. Although all the norm characterizations provided by these studies are insightful for future research, existing results are quite restricted due to linear activations, simple settings or intractable problems.
  Shallow CNNs and their representational power: As opposed to their relatively simple and shallow architecture, CNNs with two/three layers are very powerful and efficient models. Belilovsky et al. (2019) show that greedy training of two/three layer CNNs can achieve comparable performance to deeper models, e.g., VGG-11(Simonyan & Zisserman, 2014). However, a full theoretical understanding and interpretable description of CNNs even with a single hidden layer is lacking in the literature.
  Our contributions: Our contributions can be summarized as follows:
  • We develop convex programs that are polynomial time with respect to all input parameters: the number of samples, data dimension, and the number of neurons to globally train CNNs. To the best of our knowledge, this is the first work characterizing polynomial time trainability of nonconvex CNN models. More importantly, we achieve this complexity with explicit and interpretable convex optimization problems. Consequently, training CNNs, especially in practice, can be further accelerated by leveraging extensive tools available from convex optimization theory.
  • Our work reveals a hidden regularization mechanism behind CNNs and characterizes how the architecture and pooling strategies, e.g., max-pooling, average pooling, and flattening, dramatically alter the regularizer. As we show, ranging from `1 and `2 norm to nuclear norm (see Table 1 for details), ReLU CNNs exhibit an extremely rich and elegant regularization structure which is implicitly enforced by architectural choices. In convex optimization and signal processing, `1, `2 and nuclear norm regularizations are well studied, where these structures have been applied in compressed sensing, inverse problems, and matrix completion. Our results bring light to unexplored and promising connections of ReLU CNNs with these established disciplines.
  Notation and preliminaries: We denote matrices/vectors as uppercase/lowercase bold letters, for which a subscript indicates a certain element/column. We use Ik for the identity matrix of size k. We denote the set of integers from 1 to n as [n]. Moreover, ‖ · ‖F and ‖ · ‖∗ are Frobenius and nuclear norms and Bp := {u ∈ Cd : ‖u‖p ≤ 1} is the unit `p ball. We also use 1[x ≥ 0] as an indicator. To keep the presentation simple, we will use a regression framework with scalar outputs and squared loss. However, we also note that all of our results can be extended to vector outputs and arbitrary convex regression and classification loss functions. We present these extensions in Appendix. In our regression framework, we denote the input data matrix and the corresponding label vector as X ∈ Rn×d and y ∈ Rn, respectively. Moreover, we represent the patch matrices, i.e., subsets of columns, extracted from X as Xk ∈ Rn×h, k ∈ [K], where h denotes the filter size. With this notation, {Xku}Kk=1 describes a convolution operation between the filter u ∈ Rh and the data matrix X. Throughout the paper, we will use the ReLU activation function defined as (x)+ = max{0, x}. However, since CNN training problems with ReLUs are not convex in their conventional form, below we introduce an alternative formulation for this activation, which will be crucial for our derivations.
  Prior Work (Pilanci & Ergen, 2020): Recently, Pilanci & Ergen (2020) introduced an exact convex formulation for training two-layer fully connected ReLU networks in polynomial time for training data X ∈ Rn×d of constant rank, where the model is a standard two-layer scalar output network fθ(X) := ∑m j=1 (Xuj)+ αj . However, this model has three main limitations. First, as noted by the authors, even though the algorithm is polynomial time, i.e., O(nr), provided that r := rank(X), the complexity is exponential in r = d, i.e., O(nd), if X is full rank. Additionally, as a direct consequence of their model, the analysis is limited to fully connected architectures. Although they briefly analyzed some CNN architectures in Section 4, as emphasized by the authors, these are either fully linear (without ReLU) or separable over the patch index k as fully connected models, which do not correspond to weight sharing in classical CNN architectures in practice. Finally, their analysis does not extend to three-layer architectures with two ReLU layers since the analysis of two ReLU layers is significantly more challenging. On the contrary, we prove that classical CNN architectures can be globally optimized by standard convex solvers in polynomial time independent of the rank
  1The results on two-layer CNNs are presented in Appendix A.4. 2This refers to an L-layer network with only one ReLU layer and circular convolutions.
sections.heading: 1.1 HYPERPLANE ARRANGEMENTS
sections.text: LetH be the set of all hyperplane arrangement patterns of X, defined as the following set
  H := ⋃{ {sign(Xw)} : w ∈ Rd } ,
  which has finitely many elements, i.e., |H| ≤ NH < ∞, NH ∈ N. We now define a collection of sets that correspond to positive signs for each element in H, by S := { {∪hi=1{i}} : h ∈ H } . We first note that ReLU is an elementwise function that masks the negative entries of a vector or matrix. Hence, given a set S ∈ S, we define a diagonal mask matrix D(S) ∈ Rn×n defined as D(S)ii := 1[i ∈ S]. Then, we have an alternative representation for ReLU as (Xw)+ = D(S)Xw given D(S)Xw ≥ 0 and (In −D(S)) Xw ≤ 0. Note that these constraints can be compactly defined as (2D(S)− In) Xw ≥ 0. If we denote the cardinality of S as P , i.e., the number of regions in a partition of Rd by hyperplanes passing through the origin and are perpendicular to the rows of the data matrix X with r := rank(X) ≤ min(n, d), then P can be upper-bounded as follow
  P ≤ 2 r−1∑ k=0 ( n− 1 k ) ≤ 2r ( e(n− 1) r )r (Ojha, 2000; Stanley et al., 2004; Winder, 1966; Cover, 1965) (see Appendix A.2 for details).
sections.heading: 1.2 CONVOLUTIONAL HYPERPLANE ARRANGEMENTS
sections.text: We now define a notion of hyperplane arrangements for CNNs, where we introduce the patch matrices {Xk}Kk=1 instead of directly operating on X. We first construct a new data matrix as M = [X1; X2; . . . XK ] ∈ RnK×h. We then define convolutional hyperplane arrangements as the hyperplane arrangements for M and denote the cardinality of this set as Pconv . Then, we have
  Pconv ≤ 2 rc−1∑ k=0 ( nK − 1 k ) ≤ 2rc ( e(nK − 1) rc )rc where rc := rank(M) ≤ h and K = ⌊ d−h stride ⌋ + 1. Note that when the filter size h is fixed, Pconv is polynomial in n and d. Similarly, we consider hyperplane arrangements for circular CNNs followed by a linear pooling layer, i.e., XUw, where U ∈ Rd×d is a circulant matrix generated by the elements u ∈ Rh. Then, we define circular convolutional hyperplane arrangements and denote the cardinality of this set as Pcconv , which is exponential in the rank of the circular patch matrices, i.e., rcc.
  Remark 1.1. There exist P hyperplane arrangements of X where P is exponential in r. Thus, if X is full rank, r = d, then P can be exponentially large in the dimension d. As we will show, this makes the training problem for fully connected networks challenging. On the other hand, for CNNs, the number of relevant hyperplane arrangements Pconv is exponential in rc. If M is full rank, then rc = h d and accordingly Pconv P . This shows that the parameter sharing structure in CNNs enables a significant reduction in the number of possible hyperplane arrangements. Consequently, as shown in the sequel and Table 2, our results imply that the complexity of training problem is significantly lower compared to fully connected networks.
sections.heading: 2 TWO-LAYER CNNS
sections.text: In this section, we present exact convex formulation for two-layer CNN architectures.
sections.heading: 2.1 TWO-LAYER CNNS WITH AVERAGE POOLING
sections.text: We first consider an architecture with m filters, average pooling3, i.e., is defined as fθ(X) :=∑
  j ∑ k (Xkuj)+ wj with parameters θ := {uj , wj}, and standard weight decay regularization,
  which can be trained via the following problem
  p∗1 = min{uj ,wj}mj=1
  1
  2 ∥∥∥∥∥∥ m∑ j=1 K∑ k=1 (Xkuj)+ wj − y ∥∥∥∥∥∥ 2
  2
  + β
  2 m∑ j=1 ( ‖uj‖22 + w2j ) , (1)
  where uj ∈ Rh and w ∈ Rm are the filter and output weights, respectively, and β > 0 is a regularization parameter. After a rescaling (see Appendix A.3), we obtain the following problem
  p∗1 = min{uj ,wj}mi=1 uj∈B2,∀j
  1
  2 ∥∥∥∥∥∥ m∑ j=1 K∑ k=1 (Xkuj)+ wj − y ∥∥∥∥∥∥ 2
  2
  + β‖w‖1. (2)
  Then, taking dual with respect to w and changing the order of min-max yields the weak dual
  p∗1 ≥ d∗1 = max v −1 2 ‖v − y‖22 + 1 2 ‖y‖22 s.t. max u∈B2 ∣∣∣∣∣ K∑ k=1 vT (Xku)+ ∣∣∣∣∣ ≤ β, (3) which is a semi-infinite optimization problem and the dual can be obtained as a finite dimensional convex program using semi-infinite optimization theory (Goberna & López-Cerdá, 1998). The same dual also corresponds to the bidual of equation 1. Surprisingly, strong duality holds when m exceeds a threshold. Then, using strong duality, we characterize a set of optimal filter weights as the extreme point of the constraint in equation 3. Below, we use this characterization to derive an exact convex formulation for equation 1. Theorem 2.1. Let m be a number such that m ≥ m∗ for some m∗ ∈ N,m∗ ≤ n + 1, then strong duality holds for equation 3, i.e., p∗1 = d ∗ 1, and the equivalent convex program for equation 1 is
  min {ci,c′i} Pconv i=1 ci,c ′ i∈R h,∀i
  1
  2 ∥∥∥∥∥ Pconv∑ i=1 K∑ k=1 D(Ski )Xk (c ′ i − ci)− y ∥∥∥∥∥ 2
  2
  + β Pconv∑ i=1 (‖ci‖2 + ‖c′i‖2) (4)
  s.t. (2D(Ski )− In)Xkci ≥ 0, (2D(Ski )− In)Xkc′i ≥ 0,∀i, k. Moreover, an optimal solution to equation 1 with m∗ filters can be constructed as follows
  (u∗j1i , w ∗ j1i) =
  ( c′ ∗
  i√ ‖c′∗i ‖2
  , √ ‖c′∗i ‖2 ) if ‖c′ ∗ i ‖2 > 0
  (u∗j2i , w ∗ j2i) = ( c∗i√ ‖c∗i ‖2 , − √ ‖c∗i ‖2 ) if ‖c∗i ‖2 > 0 ,
  where {c′∗i , c∗i } Pconv i=1 are optimal, m ∗ := ∑Pconv i=1 1[‖c∗i ‖2 6= 0] + ∑Pconv i=1 1[‖c′ ∗
  i ‖2 6= 0], and jsi ∈ [|Js|] given the definitions J1 := {i1 : ‖c′i1‖ > 0} and J2 := {i2 : ‖ci2‖ > 0}. 4
  Therefore, we obtain a finite dimensional convex formulation with 2hPconv variables and 2nPconvK constraints for the non-convex problem in equation 1. Since Pconv is polynomial in n and d given a fixed rc ≤ h, equation 4 can be solved by a standard convex optimization solver in polynomial time. Remark 2.1. Table 2 shows that for fixed rank rc, or fixed filter size h, the complexity is polynomial in all problem parameters: n (number of samples), m (number of filters, i.e., neurons), and d (dimension). The filter size h is typically a small constant, e.g., h = 9 for 3× 3 filters. We also note that for fixed n and rank(X) = d, the complexity of fully connected networks is exponential in d, which cannot be improved unless P = NP even for m = 2 (Boob et al., 2018; Pilanci & Ergen, 2020). However, this result shows that CNNs can be trained to global optimality with polynomial complexity as a convex program.
  3We define the average pooling operation as ∑K
  k=1 (Xkuj)+, which is also known as global average pooling. 4Since our proof technique is similar for different CNNs, we present only the proof of Theorem 2.1 in Section
  5. The rest of the proofs can be found in Appendix (including the strong duality results in A.7).
  Interpreting non-convex CNNs as convex variable selection models: Interestingly, we have the sum of the squared `2 norms of the weights (i.e., weight decay regularization) in the non-convex problem equation 1 as the regularizer, however, the equivalent convex program in equation 4 is regularized by the sum of the `2 norms of the weights. This particular regularizer is known as group `1 norm, and is well-studied in the context of sparse recovery and variable selection (Yuan & Lin, 2006; Meier et al., 2008). Hence, our convex program reveals an implicit variable selection mechanism in the original non-convex problem. More specifically, the original features in X are mapped to higher dimensions via convolutional hyperplane arrangements as {D(Ski )Xk} Pconv i=1 and followed by a convex variable selection strategy using the group `1 norm. Below, we show that this implicit regularization changes significantly with the CNN architecture and pooling strategies and can range from `1 and `2 norms to nuclear norm.
sections.heading: 2.2 TWO-LAYER CNNS WITH MAX POOLING
sections.text: Here, we consider the architecture with max pooling, i.e., fθ(X) = ∑ j maxpool ( {(Xkuj)+}k ) wj , which is trained as follows
  p∗1 = min{uj ,wj}mi=1 uj∈B2,∀j
  1
  2 ∥∥∥∥∥∥ m∑ j=1 maxpool ( {(Xkuj)+} K k=1 ) wj − y ∥∥∥∥∥∥ 2
  2
  + β‖w‖1, (5)
  where maxpool(·) is an elementwise max function over the patch index k. Then, taking dual with respect to w and changing the order of min-max yields
  p∗1 ≥ d∗1 = max v −1 2 ‖v − y‖22 + 1 2 ‖y‖22 s.t. max u∈B2 ∣∣vTmaxpool ({(Xku)+}Kk=1)∣∣ ≤ β. (6) Theorem 2.2. Let m be a number such that m ≥ m∗ for some m∗ ∈ N,m∗ ≤ n + 1, then strong duality holds for equation 6, i.e., p∗1 = d ∗ 1, and the equivalent convex program for equation 5 is
  min {ci,c′i} Pconv i=1 ci,c ′ i∈R h,∀i
  1
  2 ∥∥∥∥∥ Pconv∑ i=1 K∑ k=1 D(Ski )Xk (c ′ i − ci)− y ∥∥∥∥∥ 2
  2
  + β Pconv∑ i=1 (‖ci‖2 + ‖c′i‖2) (7)
  s.t. (2D(Ski )− In)Xkci ≥ 0, (2D(Ski )− In)Xkc′i ≥ 0,∀i, k, D(Ski )Xkci ≥ D(Ski )Xjci, D(Ski )Xkc′i ≥ D(Ski )Xjc′i,∀i, j, k.
  Moreover, an optimal solution to equation 5 can be constructed from equation 7 as in Theorem 2.1.
  We note that max pooling corresponds to the last two linear constraints of the above program. Hence, max pooling can be interpreted as additional regularization, which constraints the parameters further.
sections.heading: 3 MULTI-LAYER CIRCULAR CNNS
sections.text: In this section, we first consider L-layer circular CNNs with L− 2 pooling layers before ReLU, i.e., fθ(X) = ∑ j (X ∏ l Uljw1j)+ w2j , which is trained via the following non-convex problem
  p∗2 = min {{ulj}L−2l=1 ,w1j ,w2j} m j=1
  ulj∈UL,∀l,j
  1
  2 ∥∥∥∥∥∥ m∑ j=1 ( X L−2∏ l=1 Uljw1j ) + w2j − y ∥∥∥∥∥∥ 2
  2
  + β
  2 m∑ j=1 ( ‖w1j‖22 + w22j ) , (8)
  where Ulj ∈ Rd×d is a circulant matrix generated using ulj ∈ Rhl and UL := {(u1, . . . ,uL−2) : ul ∈ Rhl ,∀l ∈ [L− 2]; ∥∥∥∏L−2l=1 Ul∥∥∥2
  F ≤ 1} and we include unit norm constraints w.l.o.g.
  Theorem 3.1. Let m be a number such that m ≥ m∗ for some m∗ ∈ N,m∗ ≤ n + 1, then strong duality holds for equation 8, i.e., p∗2 = d ∗ 2, and the equivalent convex problem is
  min {ci,c′i} Pcconv i=1
  ci,c ′ i∈C d,∀i
  1
  2 ∥∥∥∥∥ Pcconv∑ i=1 D(Si)X̃ (c ′ i − ci)− y ∥∥∥∥∥ 2
  2
  + β
  d L−2 2 Pcconv∑ i=1 (‖ci‖1 + ‖c′i‖1) (9)
  s.t. (2D(Si)− In)X̃ci ≥ 0, (2D(Si)− In)X̃c′i ≥ 0,∀i,
  where X̃ = XF and F ∈ Cd×d is the DFT matrix. Additionally, as in Theorem 2.1, we can construct an optimal solution to equation 8 from equation 9.5
  Remarkably, although the sum of the squared `2 norms in the non-convex problem in equation 8 stand for the standard weight decay regularizer, the equivalent convex program in equation 9 is regularized by the sum of the `1 norms which encourages sparsity in the spectral domain X̃. Thus, even with the simple choice of the weight decay in the non-convex problem, the architectural choice for a CNN implicitly employs a more sophisticated regularizer that is revealed by our convex optimization approach. We further note that in the above problem D(Si)X̃ are the spectral features of a subset of data points which are seperated by a hyperplane from all the other spectral features. While such spectral features can be very predictive for images in many applications, we believe that our convex program also sheds light into the undesirable bias of CNNs, e.g., towards certain textures and low frequencies (Geirhos et al., 2018; Rahaman et al., 2019).
sections.heading: 4 THREE-LAYER CNNS WITH TWO RELU LAYERS
sections.text: Here, we consider three-layer CNNs with two ReLU layers, which has the following primal problem
  p∗3 = min{uj ,w1j ,w2j}mj=1 uj∈B2
  1
  2 ∥∥∥∥∥∥ m∑ j=1 ( K∑ k=1 (Xkuj)+ w1jk ) + w2j − y ∥∥∥∥∥∥ 2
  2
  + β
  2 m∑ j=1 ( ‖w1j‖22 + w22j ) (10)
  with fθ(X) = ∑ j (∑ k (Xkuj)+ w1jk ) + w2j and the following convex equivalent problem. Theorem 4.1. Let m be a number such that m ≥ m∗ for some m∗ ∈ N,m∗ ≤ n + 1, then strong duality holds for equation 8, i.e., p∗3 = d ∗ 3, and the equivalent convex problem is
  min {cijk,c′ijk}ijk cijk,c ′ ijk∈R h
  1
  2 ∥∥∥∥∥∥ P2∑ j=1 DS2j P1∑ i=1 K∑ k=1 IijkD(Sk1i)Xk ( c′ijk − cijk ) − y ∥∥∥∥∥∥ 2
  2
  + β P1∑ i=1 P2∑ j=1 ( ‖Cij‖F + ‖C′ij‖F )
  s.t. (2D(S2j)− In) P1∑ i=1 K∑ k=1 IijkD(Sk1i)Xkcijk ≥ 0, (2D(Sk1i)− In)Xkcijk ≥ 0,∀i, j, k (11)
  (2D(S2j)− In) P1∑ i=1 K∑ k=1 IijkD(Sk1i)Xkc′ijk ≥ 0, (2D(Sk1i)− In)Xkc′ijk ≥ 0,∀i, j, k.
  where P1 and P2 are the number hyperplane arrangements for the first and second layers, Iijk ∈ {±1} are sign patterns to enumerate all possible sign patterns of the second layer weights, and Cij = [cij1 . . . cijK ] (see Appendix A.10 for further details).
  It is interesting to note that, although the sum of the squared `2 norms in the non-convex problem equation 10 is the standard weight decay regularizer, the equivalent convex program equation 11 is regularized by the sum of the Frobenius norms that promote matrix group sparsity, where the groups are over the patch indices. Note that this is similar to equation 4 except an extra summation due to having one more ReLU layer. Therefore, we observe that adding more convolutional layers with ReLU implicitly regularizes for group sparsity over a richer hierarchical representation of the data via two consecutive hyperplane arrangements.
sections.heading: 5 PROOF OF THE MAIN RESULT (THEOREM 2.1)
sections.text: Here, we provide our proof technique for Theorem 2.1. We first focus on the single-sided constraint
  max u∈B2 K∑ k=1 vT (Xku)+ ≤ β, (12)
  5The details are presented in Appendix A.9
  where the maximization problem can be written as
  max Sk⊆[n] Sk∈S max u∈B2
  K∑ k=1 vTD(Sk)Xku s.t. (2D(Sk)− In)Xku ≥ 0, ∀k. (13)
  Since the maximization is convex and strictly feasible for fixed D(Sk), equation 13 can be written as
  max Sk⊆[n] Sk∈S min αk≥0 max u∈B2
  K∑ k=1 ( vTD(Sk)Xk + α T k (2D(S k)− In)Xk ) u
  = max Sk⊆[n] Sk∈S min αk≥0
  ∥∥∥∥∥ K∑ k=1 vTD(Sk)Xk + α T k (2D(S k)− In)Xk ∥∥∥∥∥ 2 .
  We now enumerate all hyperplane arrangements and index them in an arbitrary order, i.e., denoted as( S1i , . . . , S K i ) , where i ∈ [Pconv], Pconv = |SK |, SK := {(S1i , . . . , SKi ) : Ski ∈ S,∀k, i}. Then,
  equation 12 ⇐⇒ ∀i ∈ [Pconv], min αk≥0 ∥∥∥∥∥ K∑ k=1 vTD(Ski )Xk + α T k (2D(S k i )− In)Xk ∥∥∥∥∥ 2 ≤ β
  ⇐⇒ ∀i ∈ [Pconv], ∃αik ≥ 0 s.t. ∥∥∥∥∥ K∑ k=1 vTD(Ski )Xk + α T ik(2D(S k i )− In)Xk ∥∥∥∥∥ 2 ≤ β.
  We now use the same approach for the two-sided constraint in equation 3 to obtain the following
  max v αik,α ′ ik≥0 −1 2 ‖v − y‖22 + 1 2 ‖y‖22 s.t. ∥∥∥∥∥ K∑ k=1 vTD(Ski )Xk + α T ik(2D(S k i )− In)Xk ∥∥∥∥∥ 2 ≤ β (14)
  ∥∥∥∥∥ K∑ k=1 −vTD(Ski )Xk + α′ T ik (2D(S k i )− In)Xk ∥∥∥∥∥ 2 ≤ β, ∀i.
  Note that this problem is convex and strictly feasible for v = αik = α′ik = 0. Therefore, Slater’s conditions and consequently strong duality holds, and equation 14 can be written as
  min λi,λ′i≥0 max v αik,α ′ ik≥0 −1 2 ‖v − y‖22 + 1 2 ‖y‖22 + Pconv∑ i=1 λi
  ( β − ∥∥∥∥∥ K∑ k=1 vTD(Ski )Xk + α T ik(2D(S k i )− In)Xk ∥∥∥∥∥ 2 )
  + Pconv∑ i=1 λ′i
  ( β − ∥∥∥∥∥ K∑ k=1 −vTD(Ski )Xk + α′ T ik (2D(S k i )− In)Xk ∥∥∥∥∥ 2 ) . (15)
  Next, we first introduce new variables zi, z′i ∈ Rh. Then, by recalling Sion’s minimax theorem (Sion, 1958), we change the order of the inner max-min as follows
  min λi,λ′i≥0 min zi∈B2 z′i∈B2 max v αik,α ′ ik≥0 −1 2 ‖v − y‖22 + 1 2 ‖y‖22 + Pconv∑ i=1 λi
  ( β + ( K∑ k=1 vTD(Ski )Xk + α T ik(2D(S k i )− In)Xk ) zi )
  + Pconv∑ i=1 λ′i
  ( β + ( K∑ k=1 −vTD(Ski )Xk + α′ T ik (2D(S k i )− In)Xk ) z′i ) . (16)
  We now compute the maximum with respect to v,αik,α′ik analytically to obtain the following
  min λi,λ′i≥0 min zi∈B2 z′i∈B2
  1
  2 ∥∥∥∥∥ Pconv∑ i=1 K∑ k=1 D(Ski )Xk (λ ′ iz ′ i − λizi)− y ∥∥∥∥∥ 2
  2
  + β Pconv∑ i=1 (λi + λ ′ i)
  s.t. (2D(Ski )− In)Xkzi ≥ 0, (2D(Ski )− In)Xkz′i ≥ 0,∀i, k. (17)
  Then, we apply a change of variables and define ci = λizi and c′i = λ ′ iz ′ i. Thus, we obtain
  min ci,c′i∈Rh
  1
  2 ∥∥∥∥∥ Pconv∑ i=1 K∑ k=1 D(Ski )Xk (c ′ i − ci)− y ∥∥∥∥∥ 2
  2
  + β Pconv∑ i=1 (‖ci‖2 + ‖c′i‖2)
  s.t. (2D(Ski )− In)Xkci ≥ 0, (2D(Ski )− In)Xkc′i ≥ 0,∀i, k, (18)
  since λi = ‖ci‖2 and λ′i = ‖c′i‖2 are feasible and optimal. Then, using the prescribed {u∗j , w∗j }m ∗
  j=1, we evaluate the non-convex objective in equation 1 as follows
  p∗1 ≤ 1
  2 ∥∥∥∥∥∥ m∗∑ j=1 K∑ k=1 (Xku ∗ j )+w ∗ j − y ∥∥∥∥∥∥ 2
  2
  + β
  2 Pconv∑ i=1,c′
  ∗ i 6=0
  ∥∥∥∥∥ c′ ∗ i√ ‖c′∗i ‖2 ∥∥∥∥∥ 2
  2
  + ∥∥∥∥√‖c′∗i ‖2∥∥∥∥2 2  + β
  2 Pconv∑ i=1,c∗i 6=0 ∥∥∥∥∥ c∗i√‖c∗i ‖2 ∥∥∥∥∥ 2 2 + ∥∥∥∥√‖c∗i ‖2∥∥∥∥2 2  which has the same objective value with equation 18. Since strong duality holds for the convex program, p∗1 = d ∗ 1, which is equal to the value of equation 18 achieved by the prescribed parameters.
sections.heading: 6 NUMERICAL EXPERIMENTS
sections.text: In this section6,7, we present numerical experiments to verify our claims. We first consider a synthetic dataset, where (n, d) = (6, 20), X ∈ R6×20 is generated using a multivariate normal distribution with zero mean and identity covariance, and y = [1 − 1 1 − 1 − 1 1]T . We then train the threelayer circular CNN model in equation 8 using SGD and the convex program equation 9. In Figure 1, we plot the regularized objective value with respect to the computation time with 5 different independent realizations for SGD. We also plot both the non-convex objective in equation 8 and the convex objective in equation 9 for our convex program, where optimal prescribed parameters are used to convert the solution of the convex program to the original non-convex CNN architecture (see Appendix A.9). In Figure 1a, we use 5 filters with h = 3 and stride 1, where only one trial converges to the optimal objective value achieved by both our convex program and feasible network. As m increases, all the trials are able to converge to the optimal objective value in Figure 1b. We also evaluate the same model on a subset of MNIST (LeCun) and CIFAR10 (Krizhevsky et al., 2014) for binary classification. Here, we first randomly sample the dataset and then select (n, d,m, h, stride) = (99, 50, 20, 3, 1) and a batch size of 10 for SGD. Similarly for CIFAR10, we select (n, d,m, h, stride) = (99, 50, 40, 3, 1) and use a batch size of 10 for SGD. In Figure 2, we plot both the regularized objective values in equation 8 and equation 9, and the corresponding test accuracies with the computation time. Since the number of filters is large enough, all the SGD trials converge the optimal value provided by our convex program.
sections.heading: 7 CONCLUDING REMARKS
sections.text: We studied various non-convex CNN training problems and introduced exact finite dimensional convex programs. Particularly, we provide equivalent convex characterizations for ReLU CNN architectures in a higher dimensional space. Unlike the previous studies, we prove that these equivalent characterizations have polynomial complexity in all input parameters and can be globally optimized via convex optimization solvers. Furthermore, we show that depending on the type of a CNN architecture, equivalent convex programs might exhibit different norm regularization structure, e.g., `1, `2, and nuclear norm. Thus, we claim that the implicit regularization phenomenon in modern neural networks architectures can be precisely characterized as convex regularizers. Therefore, extending our results to deeper networks is a promising direction. We also conjecture that the proposed convex approach can also be used to analyze popular heuristic techniques to train modern deep learning architectures. For example, after our work, Ergen et al. (2021) studied batch normalization through our convex framework and revealed an implicit patchwise whitening effect. Similarly, Sahiner et al. (2021) extended our model to vector outputs. More importantly, in the light of our results, efficient optimization algorithms can be developed to exactly (or approximately) optimize deep CNN architectures for large scale experiments in practice, which is left for future research.
sections.heading: ACKNOWLEDGEMENTS
sections.text: This work was partially supported by the National Science Foundation under grants IIS-1838179 and ECCS-2037304, Facebook Research, Adobe Research and Stanford SystemX Alliance.
sections.heading: Table of Contents
sections.heading: A Appendix 12
sections.text: A.1 Additional numerical results . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12 A.2 Constructing hyperplane arrangements in polynomial time . . . . . . . . . . . . 12 A.3 Equivalence of the `1 penalized objectives . . . . . . . . . . . . . . . . . . . . 14 A.4 Two-layer linear CNNs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 A.5 Extensions to vector outputs . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17 A.6 Extensions to arbitrary convex loss functions . . . . . . . . . . . . . . . . . . . 17 A.7 Strong duality results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18 A.8 Proof of Theorem 2.2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19 A.9 Proof of Theorem 3.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 A.10 Proof of Theorem 4.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
sections.heading: A APPENDIX
sections.text: In this section, we present additional materials and proofs of the main results that are not included in the main paper due to the page limit.
sections.heading: A.1 ADDITIONAL NUMERICAL RESULTS
sections.text: Here, we present additional numerical experiments to further verify our theory. We first perform an experiment with another synthetic dataset, where X ∈ R6×15 is generated using a multivariate normal distribution with zero mean and identity covariance, and y = [1 − 1 1 1 1 − 1]T . In this case, we use the two-layer CNN model in equation 1 and the corresponding convex program in equation 4. In Figure 3, we perform the experiment using m = 5, 8, 15 filters of size h = 10 and stride 5, where we observe that as the number of filters increases, the ratio of the trials converging to the optimal objective value increases as well.
  In order to apply our convex approach in Theorem 2.1 to larger scale experiments, we now introduce an unconstrained version of the convex program in equation 4 as follows
  min {ci,c′i} Pconv i=1 ci,c ′ i∈R h,∀i
  1
  2 ∥∥∥∥∥ Pconv∑ i=1 K∑ k=1 D(Ski )Xk (c ′ i − ci)− y ∥∥∥∥∥ 2
  2
  + β Pconv∑ i=1 (‖ci‖2 + ‖c′i‖2) (19)
  + ρ1T Pconv∑ i=1 K∑ k=1 (( −(2D(Ski )− In)Xkci ) + + ( −(2D(Ski )− In)Xkc′i ) + ) ,
  where ρ > 0 is a trade-off parameter. Since the problem in equation 19 is in an unconstrained form, we can directly optimize its parameters using conventional algorithms such as SGD. Hence, we use PyTorch to optimize the parameters of a two-layer CNN architecture using both the non-convex objective in equation 1 and the convex objective in equation 19, where we use the full CIFAR-10 dataset for binary classification, i.e., (n, d) = (10000, 3072). In Figure 4, we provide the training objective and the test accuracy of each approach with respect to the number of epochs. Here, we observe that the optimization on the convex formulation achieves lower training objective and higher test accuracy compared to the classical optimization on the non-convex problem.
sections.heading: A.2 CONSTRUCTING HYPERPLANE ARRANGEMENTS IN POLYNOMIAL TIME
sections.text: In this section, we discuss the number of distinct hyperplane arrangements, i.e., P , and present algorithm that enumerates all the distinct arrangements in polynomial time.
  We first consider the number of all distinct sign patterns sign(Xw) for all w ∈ Rd. This number corresponds to the number of regions in a partition of Rd by hyperplanes passing through the origin,
  and are perpendicular to the rows of X. Here, one can replace the dimensionality d with the rank of the data matrix X, i.e., denoted as r, without loss of generality. Let us first introduce the Singular Value Decomposition of X in a compact form as X = UΣVT , where U ∈ Rn×r, Σ ∈ Rr×r, and V ∈ Rr×d. Then, for a given vector w ∈ Rd, Xw = Uw′, where w′ = ΣVTw, w′ ∈ Rr. Hence, the number of distinct sign patterns sign(Xw) for all possible w ∈ Rd is equal to the number of sign patterns sign(Uw′) for all possible w′ ∈ Rr. Consider an arrangement of n hyperplanes in Rr, where n ≥ r. Let us denote the number of regions in this arrangement by Pn,r. In Ojha (2000); Cover (1965), it is shown that this number satisfies
  Pn,r ≤ 2 r−1∑ k=0 ( n− 1 k ) .
  For hyperplanes in general position, the above inequality is in fact an equality. In Edelsbrunner et al. (1986), the authors present an algorithm that enumerates all possible hyperplane arrangements O(nr) time, which can be used to construct the data for the convex programs we present throughout the paper.
sections.heading: A.3 EQUIVALENCE OF THE `1 PENALIZED OBJECTIVES
sections.text: In this section, we prove the equivalence between the original problems with `2 regularization and their `1 penalized versions. We also note that similar equivalence results were also presented in Savarese et al. (2019); Neyshabur et al. (2014); Ergen & Pilanci (2019; 2020c;d). We start with the equivalence between equation 1 and equation 2. Lemma A.1. The following two problems are equivalent:
  min {uj ,wj}mj=1
  1
  2 ∥∥∥∥∥∥ m∑ j=1 K∑ k=1 (Xkuj)+ wj − y ∥∥∥∥∥∥ 2
  2
  + β
  2 m∑ j=1 ( ‖uj‖22 + w2j )
  = min {uj ,wj}mj=1 uj∈B2,∀j
  1
  2 ∥∥∥∥∥∥ m∑ j=1 K∑ k=1 (Xkuj)+ wj − y ∥∥∥∥∥∥ 2
  2
  + β m∑ j=1 ‖w‖1.
  Proof of Lemma A.1. We rescale the parameters as ūj = γjuj and w̄j = wj/γj , for any γj > 0. Then, the output becomes
  m∑ j=1 K∑ k=1 (Xkūj)+w̄j = m∑ j=1 K∑ k=1 (Xkujγj)+ wj γj = m∑ j=1 K∑ k=1 (Xuj)+wj ,
  which proves that the scaling does not change the network output. In addition to this, we have the following basic inequality
  1
  2 m∑ j=1 (‖uj‖22 + w2j ) ≥ m∑ j=1 (|wj | ‖uj‖2),
  where the equality is achieved with the scaling choice γj = ( |wj | ‖uj‖2 ) 1 2
  is used. Since the scaling operation does not change the right-hand side of the inequality, we can set ‖uj‖2 = 1,∀j. Therefore, the right-hand side becomes ‖w‖1. Now, let us consider a modified version of the problem, where the unit norm equality constraint is relaxed as ‖uj‖2 ≤ 1. Let us also assume that for a certain index j, we obtain ‖uj‖2 < 1 with wj 6= 0 as an optimal solution. This shows that the unit norm inequality constraint is not active for uj , and hence removing the constraint for uj will not change the optimal solution. However, when we remove the constraint, ‖uj‖2 →∞ reduces the objective value since it yields wj = 0. Therefore, we have a contradiction, which proves that all the constraints that correspond to a nonzero wj must be active for an optimal solution. This also shows that replacing ‖uj‖2 = 1 with ‖uj‖2 ≤ 1 does not change the solution to the problem.
  Next, we prove the equivalence between equation 8 for L = 3 and equation 30. Lemma A.2. The following two problems are equivalent:
  min {uj ,w1j ,w2j}mj=1
  uj∈B2,∀j
  1
  2 ∥∥∥∥∥∥ m∑ j=1 (XUjw1j)+ w2j − y ∥∥∥∥∥∥ 2
  2
  + β
  2 m∑ j=1 ( ‖w1j‖22 + w22j )
  = min {uj ,w1j ,w2j}mj=1 uj ,w1j∈B2,∀j
  1
  2 ∥∥∥∥∥∥ m∑ j=1 (XUjw1j)+ w2j − y ∥∥∥∥∥∥ 2
  2
  + β‖w2‖1.
  Proof of Lemma A.2. We rescale the parameters as w̄1j = γjw1j and w̄2j = w2j/γj , for any γj > 0. Then, the output becomes
  m∑ j=1 (XUjw̄1j)+w̄2j = m∑ j=1 (XUjw1jγj)+ w2j γj = m∑ j=1 (XUjw1j)+w2j ,
  which proves that the scaling does not change the network output. In addition to this, we have the following basic inequality
  1
  2 m∑ j=1 (‖w1j‖22 + w22j) ≥ m∑ j=1 (‖w1j‖2 |w2j |),
  where the equality is achieved with the scaling choice γj = ( |w2j | ‖w1j‖2 ) 1 2
  is used. Since the scaling operation does not change the right-hand side of the inequality, we can set ‖w1j‖2 = 1,∀j. Therefore, the right-hand side becomes ‖w2‖1. The rest of the proof directly follows from the proof of Lemma A.1.
sections.heading: A.4 TWO-LAYER LINEAR CNNS
sections.text: We now consider two-layer linear CNNs, for which the training problem is
  min {uj ,wj}mj=1
  1
  2 ∥∥∥∥∥∥ K∑ k=1 m∑ j=1 Xkujwjk − y ∥∥∥∥∥∥ 2
  2
  + β
  2 m∑ j=1 ( ‖uj‖22 + ‖wj‖22 ) . (20)
  Theorem A.1. (Pilanci & Ergen, 2020) The equivalent convex program for equation 20 is
  min {zk}Kk=1,zk∈Rh
  1
  2 ∥∥∥∥∥ K∑ k=1 Xkzk − y ∥∥∥∥∥ 2
  2
  + β ‖[z1, . . . , zK ]‖∗ . (21)
  Proof of Theorem A.1. We first apply a rescaling (as in Lemma A.1) to the primal problem in equation 20 as follows
  min {uj ,wj}mj=1
  uj∈B2
  1
  2 ∥∥∥∥∥∥ K∑ k=1 m∑ j=1 Xkujwjk − y ∥∥∥∥∥∥ 2
  2
  + β m∑ j=1 ‖wj‖2.
  Then, taking the dual with respect to the output layer weights wj yields
  max v −1 2 ‖v − y‖22 + 1 2 ‖y‖22 s.t. max u∈B2 √∑ k ( vTXku )2 ≤ β. (22) Let us then reparameterize the problem above as follows
  max M,v −1 2 ‖v − y‖22 + 1 2 ‖y‖22 s.t. σmax (M) ≤ β, M = [XT1 v . . . XTKv],
  where σmax(M) represent the maximum singular value of M. Then the Lagrangian is as follows
  L(λ,Z,M,v) = −1 2 ‖v − y‖22 + 1 2 ‖y‖22 + λ (β − σmax (M)) + trace(ZTM)− trace(ZT [XT1 v . . . XTKv])
  = −1 2 ‖v − y‖22 + 1 2 ‖y‖22 + λ (β − σmax (M)) + trace(ZTM)− vT K∑ k=1 Xkzk
  where λ ≥ 0. Then maximizing over M and v yields the following dual form
  min {zk}Kk=1,zk∈Rh
  1
  2 ∥∥∥∥∥ K∑ k=1 Xkzk − y ∥∥∥∥∥ 2
  2
  + β ‖[z1 . . . zK ]‖∗ ,
  where ‖[z1 . . . zK ]‖∗ = ‖Z‖∗ = ∑ i σi(Z) is the `1 norm of singular values, i.e., nuclear norm (Recht et al., 2010).
  The regularized training problem for two-layer circular CNNs as follows
  min {uj ,wj}mi=1
  1
  2 ∥∥∥∥∥∥ m∑ j=1 XUjwj − y ∥∥∥∥∥∥ 2
  2
  + β
  2 m∑ j=1 ( ‖uj‖22 + ‖wj‖22 ) (23)
  where Uj ∈ Rd×d is a circulant matrix generated by a circular shift modulo d using uj ∈ Rh. Theorem A.2. (Pilanci & Ergen, 2020) The equivalent convex program for equation 23 is
  min z∈Cd
  1
  2 ∥∥∥X̃z− y∥∥∥2 2 + β√ d ‖z‖1, (24)
  where X̃ = XF and F ∈ Cd×d is the DFT matrix.
  Proof of Theorem A.2. We first apply a rescaling (as in Lemma A.1) to the primal problem in equation 23 as follows
  min {uj ,wj}mi=1
  uj∈B2
  1
  2 ∥∥∥∥∥∥ m∑ j=1 XUjwj − y ∥∥∥∥∥∥ 2
  2
  + β m∑ j=1 ‖wj‖2
  and then taking the dual with respect to the output layer weights wj yields
  max v −1 2 ‖v − y‖22 + 1 2 ‖y‖22 s.t. max D∈D ‖vTXFDFH‖2 ≤ β,
  where D := {D : ‖D‖2F ≤ d}. In the problem above, we use the eigenvalue decomposition U = FDFH , where F ∈ Cd×d is the DFT matrix and D ∈ Cd×d is a diagonal matrix defined as D := diag( √ dFu). We also note that the unit norm constraint in the primal problem, i.e., uj ∈ B2, is equivalent to Dj ∈ D since Dj = diag( √ dFuj) and ‖Dj‖2F = d‖uj‖22 due the properties of circulant matrices.
  Now let us first define a variable change as X̃ = XF. Then, the problem above can be equivalently written as
  max v −1 2 ‖v − y‖22 + 1 2 ‖y‖22 s.t. max D∈D ‖vT X̃D‖2 ≤ β.
  Since D is a diagonal matrix with a norm constraint on its diagonal entries, for an arbitrary vector s ∈ Cn, we have
  ‖sTD‖2 = √√√√ n∑ i=1 |si|2|Dii|2 ≤ smax √√√√ n∑ i=1 |Dii|2 = smax √ d,
  where smax := maxi |si|. If we denote the maximum index as imax := arg maxi |si|, then the upper-bound is achieved when
  Dii =
  {√ d if i = imax
  0, otherwise .
  Using this observation, the problem above can be further simplified as
  max v −1 2 ‖v − y‖22 + 1 2 ‖y‖22 s.t. ‖vT X̃‖∞ ≤ β√ d .
  Then, taking the dual of this problem gives the following
  min z∈Cd
  1
  2 ∥∥∥X̃z− y∥∥∥2 2 + β√ d ‖z‖1.
sections.heading: A.5 EXTENSIONS TO VECTOR OUTPUTS
sections.text: Here, we present the extensions of our approach to vector output. To keep the notation and presentation simple, we consider the vector output version of the two-layer linear CNN model in Section A.4. The training problem is as follows
  min {uj ,{wjk}Kk=1} m j=1
  1
  2 ∥∥∥∥∥∥ K∑ k=1 m∑ j=1 Xkujw T jk −Y ∥∥∥∥∥∥ 2
  F
  + β
  2 m∑ j=1 ( ‖uj‖22 + K∑ k=1 ‖wjk‖22 ) .
  The corresponding dual problem is given by
  max V −1 2 ‖V −Y‖2F + 1 2 ‖Y‖2F s.t. max u∈B2 √√√√ K∑ k=1 ‖VTXku‖22 ≤ β.
  The maximizers of the dual are the maximal eigenvectors of ∑K k=1 X T kVV
  TXk, which are optimal filters.
  We now focus on the dual constraint as in Proof of Theorem 2.1.
  max u∈B2 √√√√ K∑ k=1 ‖VTXku‖22 = maxu,s,gk∈B2 K∑ k=1 skg T k V TXku = max u,s,gk∈B2 K∑ k=1 sk 〈 V,Xkug T k
  〉 = max ‖Gk‖∗≤1
  s∈B2
  K∑ k=1 sk 〈V,XkGk〉 = max ‖Gk‖∗≤sk
  s∈B2
  K∑ k=1 〈V,XkGk〉
  Then, the rest of the derivations directly follow Section A.4.
sections.heading: A.6 EXTENSIONS TO ARBITRARY CONVEX LOSS FUNCTIONS
sections.text: In this section, we first show the procedure to create an optimal standard CNN architecture using the optimal weights provided by the convex program in. Then, we extend our derivations to arbitrary convex loss functions.
  In order to keep our derivations simple and clear, we use the regularized two-layer architecture in equation 1. For a given convex loss function `(·,y), the regularized training problem can be stated as follows
  p∗1 = min{uj ,wj}mj=1 `  m∑ j=1 K∑ k=1 (Xkuj)+wj ,y + β 2 m∑ j=1 (‖uj‖22 + w2j ) . (25)
  Then, the corresponding finite dimensional convex equivalent is
  min {ci,c′i} Pconv i=1 ci,c ′ i∈R h,∀i `
  ( Pconv∑ i=1 K∑ k=1 D(Ski )Xk(c ′ i − ci),y ) + β Pconv∑ i=1 K∑ k=1 (‖ci‖2 + ‖c′i‖2) (26)
  s.t. (2D(Ski )− In)Xkci ≥ 0, (2D(Ski )− In)Xkci ≥ 0, ∀i, k. We now define m∗ := ∑Pconv i=1 1[‖c∗i ‖2 6= 0] + ∑Pconv i=1 1[‖c′ ∗ i ‖2 6= 0], where {c∗i , c′ ∗ i } Pconv i=1 are the optimal weights in equation 26. Theorem A.3. The convex program equation 26 and the non-convex problem equation 25, where m ≥ m∗ has identical optimal values. Moreover, an optimal solution to equation 25 can be constructed from an optimal solution to equation 26 as follows
  (u∗j1i , w ∗ j1i) =
  ( c′ ∗
  i√ ‖c′∗i ‖2
  , √ ‖c′∗i ‖2 ) if ‖c′ ∗ i ‖2 > 0
  (u∗j2i , w ∗ j2i) = ( c∗i√ ‖c∗i ‖2 , − √ ‖c∗i ‖2 ) if ‖c∗i ‖2 > 0 ,
  where {c′∗i , c∗i } Pconv i=1 are the optimal solutions to equation 26.
  Proof of Theorem A.3. We first note that there will be m∗ vectors {c′∗i , c∗i }. Constructing {u∗j , w∗j }m ∗
  j=1 as stated in the theorem, and plugging in the non-convex objective equation 25, we obtain the value
  p∗1 ≤ ` m∗∑ j=1 K∑ k=1 (Xku ∗ j )+w ∗ j ,y + β 2 Pconv∑ i=1,c′
  ∗ i 6=0
  ∥∥∥∥∥ c′ ∗ i√ ‖c′∗i ‖2 ∥∥∥∥∥ 2
  2
  + ∥∥∥∥√‖c′∗i ‖2∥∥∥∥2 2  + β
  2 Pconv∑ i=1,c∗i 6=0 ∥∥∥∥∥ c∗i√‖c∗i ‖2 ∥∥∥∥∥ 2 2 + ∥∥∥∥√‖c∗i ‖2∥∥∥∥2 2  which is identical to the objective value of the convex program equation 26. Since the value of the convex program is equal to the value of it’s dual d∗1 in the dual, we conclude that p ∗ 1 = d ∗ 1, which is equal to the value of the convex program equation 26 achieved by the prescribed parameters.
  We also show that our dual characterization holds for arbitrary convex loss functions
  min {uj ,wj}mj=1 uj∈B2,∀j `
   m∑ j=1 K∑ k=1 (Xkuj)+wj ,y + β‖w‖1, (27) where `(·,y) is a convex loss function.
  Theorem A.4. The dual of equation 27 is given by
  max v −`∗(v) s.t. ∣∣∣∣∣ K∑ k=1 vT (Xku)+ ∣∣∣∣∣ ≤ β, ∀u ∈ B2 , where `∗ is the Fenchel conjugate function defined as
  `∗(v) = max z zTv − `(z,y) .
  Proof of Theorem A.4. The proof follows from classical Fenchel duality (Boyd & Vandenberghe, 2004). We first describe equation 27 in an equivalent form as follows
  min {uj ,wj}mj=1,z uj∈B2,∀j
  `(z,y) + β‖w‖1 s.t. z = m∑ j=1 K∑ k=1 (Xkuj)+wj , .
  Then the dual function is
  g(v) = min {uj ,wj}mj=1,z uj∈B2,∀j
  `(z,y)− vT z + vT m∑ j=1 K∑ k=1 (Xkuj)+wj + β‖w‖1.
  Therefore, using the classical Fenchel duality (Boyd & Vandenberghe, 2004) yields the claimed dual form.
sections.heading: A.7 STRONG DUALITY RESULTS
sections.text: Proposition A.1. Given m ≥ m∗, strong duality holds for equation 3, i.e., p∗1 = d∗1.
  We first review the basic properties of infinite size neural networks and introduce technical details to derive the dual of equation 3. We refer the reader to Rosset et al. (2007); Bach (2017) for further details. Let us first consider a measurable input space X with a set of continuous basis functions (i.e., neurons or filters in our context) ψu : X → R, which are parameterized by u ∈ B2. Next, we use real-valued Radon measures with the uniform norms (Rudin, 1964). Let us consider a signed Radon measure denoted as µ. Now, we can use µ to formulate an infinite size neural network as f(x) = ∫ u∈B2 ψu(x)dµ(u), where x ∈ X is the input. The norm for µ is usually defined as its
  total variation norm, which is the supremum of ∫ u∈B2 g(u)dµ(u) over all continuous functions g(u) that satisfy |g(u)| ≤ 1. Now, we consider the case where the basis functions are ReLUs, i.e., ψu = ( xTu ) +
  . Then, the output of a network with finitely many neurons, say m neurons, can be written as
  f(x) = m∑ j=1 ψujwj
  which can be obtained by selecting µ as a weighted sum of Dirac delta functions, i.e., µ =∑m j=1 wjδ(u − uj). In this case, the total variation norm, denoted as ‖µ‖TV , corresponds to the `1 norm ‖w‖1. Now, we ready to derive the dual of equation 3, which can be stated as follows (see Section 8.6 of Goberna & López-Cerdá (1998) and Section 2 of Shapiro (2009) for further details)
  d∗1 ≤ p1,∞ = min µ
  1
  2 ∥∥∥∥∥ ∫ u∈B2 K∑ k=1 (Xku)+ dµ(u)− y ∥∥∥∥∥ 2
  2
  + β‖µ‖TV . (28)
  Although equation 28 involves an infinite dimensional integral form, by Caratheodory’s theorem, we know that the integral can be represented as a finite summation, to be more precise, a summation of at most n+ 1 Dirac delta functions (Rosset et al., 2007). If we denote the number of Dirac delta functions as m∗, where m∗ ≤ n+ 1, then we have
  p1,∞ = min {uj ,wj}m ∗ j=1
  uj∈B2,∀j
  1
  2 ∥∥∥∥∥∥ m∗∑ j=1 K∑ k=1 (Xkuj)+ wj − y ∥∥∥∥∥∥ 2
  2
  + β‖w‖1
  = p∗1
  provided that m ≥ m∗. We now need to show that strong duality holds, i.e., p∗1 = d∗1. We first note that the semi-infinite problem equation 3 is convex. Then, we prove that the optimal value is finite. Since β > 0, we know that v = 0 is strictly feasible, and achieves 0 objective value. Moreover, since −‖y − v‖22 ≤ 0, the optimal objective value p∗1 is finite. Therefore, by Theorem 2.2 of Shapiro (2009), strong duality holds, i.e., p∗1,∞ = d ∗ 1 provided that the solution set of equation 3 is nonempty and bounded. We also note that the solution set of equation 3 is the Euclidean projection of y onto a convex, closed and bounded set since (Xku)+ can be expressed as the union of finitely many convex closed and bounded sets.
sections.heading: A.8 PROOF OF THEOREM 2.2
sections.text: The proof follows the proof of Proposition A.1. The dual of equation 6 is as follows
  d∗1 ≤ p1,∞ = min µ
  1
  2 ∥∥∥∥∫ u∈B2 maxpool ( {(Xku)+} K k=1 ) dµ(u)− y ∥∥∥∥2 2 + β‖µ‖TV ,
  which has the following finite equivalent
  p1,∞ = min {uj ,wj}m ∗ j=1
  uj∈B2,∀j
  1
  2 ∥∥∥∥∥∥ m∗∑ j=1 maxpool ( {(Xkuj)+} K k=1 ) wj − y ∥∥∥∥∥∥ 2
  2
  + β‖w‖1
  = p∗1
  provided that m ≥ m∗. We now need to show that strong duality holds, i.e., p∗1 = d∗1. Since maxpool(·) can be expressed as the union of finitely many convex, closed and bounded sets, the rest of the strong duality results directly follow from the proof of Proposition A.1.
  We now focus on the single-sided dual constraint
  max u∈B2
  vTmaxpool( {
  (Xku)+} K k=1
  ) ≤ β,
  which can be written as
  max Sk⊆[n] Sk∈S max u∈B2
  K∑ k=1 vTD(Sk)Xku s.t. (2D(Sk)− In)Xku ≥ 0, ∀k,
  D(Sk)Xku ≥ D(Sk)Xju,∀j, k ∈ [K], K∑ k=1 D(Sk) = In.
  We again enumerate all hyperplane arrangements and index them in an arbitrary order, where we define the overall set as SK := {(S1i , . . . , SKi ) : Ski ∈ S,∀k, i; ∑K k=1 D(S k i ) = In,∀i} and Pconv = |SK |. Then, following the same steps in (13)–(17) gives the following convex problem
  min wi,w′i∈Rh
  1
  2 ∥∥∥∥∥ Pconv∑ i=1 K∑ k=1 D(Ski )Xk (w ′ i −wi)− y ∥∥∥∥∥ 2
  2
  + β Pconv∑ i=1 (‖wi‖2 + ‖w′i‖2) (29)
  s.t. (2D(Ski )− In)Xkwi ≥ 0, (2D(Ski )− In)Xkw′i ≥ 0,∀i, k, D(Ski )Xkwi ≥ D(Ski )Xjwi, D(Ski )Xkw′i ≥ D(Ski )Xjw′i,∀i, j, k.
  We now note that there will be m∗ pairs {w′∗i ,w∗i }. Then, we can construct a set of weights {u∗j , w∗j }m ∗
  j=1 as defined in the theorem and evaluate the non-convex objective in equation 5 using these weights as follows
  p∗1 ≤ 1
  2 ∥∥∥∥∥∥ m∗∑ j=1 maxpool ( { ( Xku ∗ j ) + }Kk=1 ) w∗j − y ∥∥∥∥∥∥ 2
  2
  + β
  2 Pconv∑ i=1,w′
  ∗ i 6=0
  ∥∥∥∥∥ w′ ∗ i√ ‖w′∗i ‖2 ∥∥∥∥∥ 2
  2
  + ∥∥∥∥√‖w′∗i ‖2∥∥∥∥2 2  + β
  2 Pconv∑ i=1,w∗i 6=0 ∥∥∥∥∥ w∗i√‖w∗i ‖2 ∥∥∥∥∥ 2 2 + ∥∥∥∥√‖w∗i ‖2∥∥∥∥2 2  which has the same objective value with equation 29. Since strong duality holds for the convex program, we have p∗1 = d ∗ 1, which is equal to the value of the convex program equation 29 achieved by the prescribed parameters above.
sections.heading: A.9 PROOF OF THEOREM 3.1
sections.text: By using a rescaling for each w1j and w2j , equation 8 can be equivalently stated as
  p∗2 = min {{ulj}L−2l=1 ,w1j ,w2j} m j=1
  w1j∈B2,ulj∈UL,∀l,j
  1
  2 ∥∥∥∥∥∥ m∑ j=1 ( X L−2∏ l=1 Uljw1j ) + w2j − y ∥∥∥∥∥∥ 2
  2
  + β‖w2‖1. (30)
  Let us denote the eigenvalue decomposition of Ulj as Ulj = FDljFH , where F ∈ Cd×d is the DFT matrix and Dlj ∈ Cd×d is a diagonal matrix. Then, we again take the dual with respect to w2 and change the order of min-max as follows
  p∗2 ≥ d∗2 = max v −1 2 ‖v − y‖22 + 1 2 ‖y‖22 s.t. max
  Dlj∈DL w1j∈B2
  ∣∣∣∣∣∣vT ( XF L−2∏ l=1 DljF Hw1j ) + ∣∣∣∣∣∣ ≤ β, ∀j, (31)
  where DL := {(D1, . . . ,DL−2 : Dl ∈ Cd×d,∀l ∈ [L − 2]; ∥∥∥∏L−2l=1 Dl∥∥∥2
  F ≤ dL−2}. Below we
  prove that strong duality holds for equation 31.
  In order to obtain the dual of the semi-infinite problem in equation 31, we again take dual with respect to v (see Appendix A.7 and Goberna & López-Cerdá (1998); Shapiro (2009) for further details), which yields
  d∗2 ≤ p2,∞ = min µ
  1
  2 ∥∥∥∥∥∥ ∫ θL∈ΘL ( XF L−2∏ l=1 DlF Hw1 ) + dµ(θL)− y ∥∥∥∥∥∥ 2
  2
  + β‖µ‖TV
  where ΘL := {(D1, . . . ,DL−2,w1) : Dl ∈ DL,∀l ∈ [L − 2]; w1 ∈ B2}. Then, selecting µ = ∑m∗ j=1 w2j δ(θL − θLj), where m∗ ≤ n+ 1, gives
  p2,∞ = min {{Dlj}L−2l=1 ,w1j ,w2j} m∗ j=1
  Dlj∈DL,w1j∈B2,∀j,l
  1
  2 ∥∥∥∥∥∥ m∗∑ j=1 ( XF L−2∏ l=1 DljF Hw1j ) + w2j − y ∥∥∥∥∥∥ 2
  2
  + β‖w2‖1
  = p∗2 provided that m ≥ m∗ holds. Then, the rest of the strong duality proof directly follows from Proof of Proposition A.1.
  We now focus on the single-sided dual constraint
  max Dl∈DL w̃1∈B2
  vT ( X̃
  L−2∏ l=1 Dlw̃1 ) + ≤ β,
  which can be written as
  max S⊆[n] S∈S max Dl∈DL w̃1∈B2
  vTD(S1)X̃ L−2∏ l=1 Dlw̃1 s.t. (2D(S)− In)X̃ L−2∏ l=1 Dlw̃1 ≥ 0. (32)
  Since the inner maximization is convex (after a variable change as q = ∏L−2 l=1 Dlw̃1) and there exists a strictly feasible solution for a fixed D(S) matrix, equation 32 can also be written as
  max S⊆[n] S∈S min α≥0 max Dl∈DL w̃1∈B2
  vTD(Si)X̃z + α T (2D(Si)− In)X̃ L−2∏ l=1 Dlw̃1
  = max S⊆[n] S∈S min α≥0 ‖vTD(S)X̃ + αT (2D(S)− In)X̃‖∞d
  L−2 2 .
  We now enumerate all hyperplane arrangements and index them in an arbitrary order, which are denoted as D(Si), where i ∈ [Pcconv]. Then, we have
  equation 32 ⇐⇒ ∀i ∈ [Pcconv], min α≥0 ‖vTD(Si)X̃ + αT (2D(Si)− In)X̃‖∞d
  L−2 2 ≤ β
  ⇐⇒ ∀i ∈ [Pcconv], ∃αi ≥ 0 s.t. ‖vTD(Si)X̃ + αTi (2D(Si)− In)X̃‖∞d L−2
  2 ≤ β. We now use the same approach for the two-sided constraint in equation 31 to represent equation 31 as a finite dimensional convex problem as follows
  max v αi,α ′ i≥0 −1 2 ‖v − y‖22 + 1 2 ‖y‖22 (33) s.t. ‖vTD(Si)X̃ + αTi (2D(Si)− In)X̃‖∞d L−2 2 ≤ β, ‖ − vTD(Si)X̃ + α′ T i (2D(Si)− In)X̃‖∞d L−2
  2 ≤ β, ∀i. We note that the above problem is convex and strictly feasible for v = αi = α′i = 0. Therefore, equation 33 can be written as
  min λi,λ′i≥0 max v αi,α ′ i≥0 −1 2 ‖v − y‖22 + 1 2 ‖y‖22 + Pcconv∑ i=1 λi ( β − ‖vTD(Si)X̃ + αTi (2D(Si)− In)X̃‖∞d L−2 2 )
  + Pcconv∑ i=1 λ′i ( β − ‖ − vTD(Si)X̃ + α′ T i (2D(Si)− In)X̃‖∞d L−2 2 ) . (34)
  Next, we introduce new variables zi, z′i ∈ Cd to represent equation 34 as
  min λi,λ′i≥0 max v
  αi,α ′ i≥0 min zi∈B1 z′i∈B1
  −1 2 ‖v − y‖22 + 1 2 ‖y‖22 + Pcconv∑ i=1 λi ( β + d L−2 2 ( vTD(Si)X̃ + α T i (2D(Si)− In)X̃ ) zi )
  + Pcconv∑ i=1 λ′i ( β + d L−2 2 ( −vTD(Si)X̃ + α′ T i (2D(Si)− In)X̃ ) z′i ) . (35)
  We note that the objective is concave in v,αi,α′i and convex in zi, z ′ i. Moreover the set B1 is convex and compact. We recall Sion’s minimax theorem (Sion, 1958) for the inner max-min problem and express the strong dual of the problem equation 35 as
  min λi,λ′i≥0 min zi∈B1 z′i∈B1 max v αi,α ′ i≥0 −1 2 ‖v − y‖22 + 1 2 ‖y‖22 + Pcconv∑ i=1 λi ( β + d L−2 2 ( vTD(Si)X̃ + α T i (2D(Si)− In)X̃ ) zi )
  + Pcconv∑ i=1 λ′i ( β + d L−2 2 ( −vTD(Si)X̃ + α′ T i (2D(Si)− In)X̃ ) z′i ) . (36)
  Now, we can compute the maximum with respect to v,αi,α′i analytically to obtain the following problem
  min λi,λ′i≥0 min zi∈B1 z′i∈B1
  1
  2 ∥∥∥∥∥dL−22 Pcconv∑ i=1 D(Si)X̃ (λ ′ iz ′ i − λizi)− y ∥∥∥∥∥ 2
  2
  + β Pcconv∑ i=1 (λi + λ ′ i) (37)
  s.t. (2D(Si)− In)X̃zi ≥ 0, (2D(Si)− In)X̃z′i ≥ 0,∀i.
  Now we apply a change of variables and define ci = d L−2 2 λizi and c′i = d L−2 2 λ′iz ′ i. Thus, we obtain
  min ci,c′i
  1
  2 ∥∥∥∥∥ Pcconv∑ i=1 D(Si)X̃ (c ′ i − ci)− y ∥∥∥∥∥ 2
  2
  + β
  d L−2 2 Pcconv∑ i=1 (‖ci‖1 + ‖c′i‖1) (38)
  s.t. (2D(Si)− In)X̃ci ≥ 0, (2D(Si)− In)X̃c′i ≥ 0,∀i,
  where we eliminate the variables λi, λ′i, since λi = ‖ci‖1/d L−2 2 and λ′i = ‖c′i‖1/d L−2
  2 are feasible and optimal.
  Optimal weight construction for equation 8: Given the optimal weights for the convex program in equation 38, i.e., denoted as {c∗i , c′ ∗ i } Pcconv i=1 , we use the following relation
  X̃c∗i = X̃diag
  ( d L−2 2 √ |c∗i | ‖c∗i ‖1 ) diag (√ |c∗i | d L−2 2 ) ejφi √ ‖c∗i ‖1 d L−2 2
  = X̃ ( L−2∏ l=1 diag (( |c∗i | ‖c∗i ‖1 ) 1 2(L−2) )) diag (√ |c∗i | d L−2 2 ) ejφi √ ‖c∗i ‖1 d L−2 2 ,
  where φi is defined such that c∗i = diag (|c∗i |) ejφi . Thus, we can directly set the parameters as follows
  D∗li = diag
  ( d 1 2 √ |c∗i | ‖c∗i ‖1 ) , w̃∗1i = diag (√ |c∗i | d L−2 2 ) ejφi , w∗2i = √ ‖c∗i ‖1 d L−2 2 ,
  which can be equivalently written as
  U∗li = Fdiag
  ( d 1 2 √ |c∗i | ‖c∗i ‖1 ) FH , w∗1i = √ |c∗i | d L−2 2 , w∗2i = √ ‖c∗i ‖1 d L−2 2
  to exactly match with the problem formulation in equation 8. We first note that ‖ ∏L−2 l=1 U ∗ li‖2F =
  ‖ ∏L−2 l=1 D ∗ li‖2F = dL−2,∀i, l, therefore, this set of parameters is feasible for equation 8. Now, we prove the optimality by showing that these parameters have the same regularization cost with the convex program in equation 38 as follows
  β
  2 Pccov∑ i=1 ( ‖w∗1i‖22 + w∗ 2 2i ) = β d L−2 2 Pcconv∑ i=1 ‖c∗i ‖1.
  The same steps can also be applied to c′ ∗
  i . Then, the rest of the proof directly follows from Theorem 2.1. Therefore, we prove that a set of optimal layer weights for equation 8, denoted as {{U∗lj} L−2 l=1 ,w ∗ 1j , w ∗ 2j}m ∗
  j=1, can be obtained from the optimal solution to equation 38, denoted as {c∗i , c′ ∗ i } Pcconv i=1 .
sections.heading: A.10 PROOF OF THEOREM 4.1
sections.text: We now provide the proof for the three-layer CNN architecture with two ReLU layers, which has the following primal optimization problem
  p∗3 = minuj ,w1j ,w2j uj∈B2
  1
  2 ∥∥∥∥∥∥ m∑ j=1 ( K∑ k=1 (Xkuj)+ w1jk ) + w2j − y ∥∥∥∥∥∥ 2
  2
  + β
  2 m∑ j=1 ( ‖w1j‖22 + w22j ) Then, the dual is
  p∗3 ≥ d∗3 = max v −1 2 ‖v − y‖22 + 1 2 ‖y‖22 s.t. max u,w1∈B2 ∣∣∣∣∣∣vT ( K∑ k=1 (Xku)+ w1k ) + ∣∣∣∣∣∣ ≤ β, (39) Dual of equation 39
  d∗3 ≤ p3,∞ = min µ
  1
  2 ∥∥∥∥∥∥ ∫ u,w1∈B2 ( K∑ k=1 (Xku)+ w1k ) + dµ(u,w1)− y ∥∥∥∥∥∥ 2
  2
  + β‖µ‖TV , (40)
  where strong duality holds, i.e., d∗3 = p3,∞, by Proof of Proposition A.1. Then, the finite equivalent is as follows
  p3,∞ = p3 = min uj ,w1j∈B2
  w2
  1
  2 ∥∥∥∥∥∥ m∗∑ j=1 ( K∑ k=1 (Xkuj)+ w1jk ) + w2j − y ∥∥∥∥∥∥ 2
  2
  + β‖w2‖1
  We now focus on a single-sided dual constraint
  max u,w1∈B2
  vT ( K∑ k=1 (Xku)+ w1k ) + ≤ β (41)
  which can be written as
  max S2,S k 1⊆[n]
  S2,S k 1∈S
  max u,w1∈B2
  vTD(S2) K∑ k=1 D(Sk1 )Xkuw1k s.t. (2D(S2)− In) K∑ k=1 D(Sk1 )Xkuw1k ≥ 0, (2D(Sk1 )− In)Xkuw1k ≥ 0
  = max Ik∈{±1} max S2,S k 1⊆[n]
  S2,S k 1∈S
  max w1∈B2 max ‖qk‖2≤|w1k|
  vTD(S2) K∑ k=1 IkD(Sk1 )Xkqk (42)
  s.t. (2D(S2)− In) K∑
  k=1
  IkD(Sk1 )Xkqk ≥ 0, (2D(Sk1 )− In)Xkqk ≥ 0,
  where we introduce the notation Ik ∈ {±1} to enumerate all possible sign patterns for w1k. Since the inner maximization is convex and there exists a strictly feasible solution for fixed D(Sk1 ),D(S2), and Ik equation 42 can also be written as
  max Ik∈{±1} max S2,S k 1⊆[n]
  S2,S k 1∈S
  min αk,γ≥0 max w1∈B2 max ‖qk‖2≤w1k K∑ k=1 vTD(Sk1 )Xkqk +α T k (2D(S k 1 )− In)Xkqk + IkγT (2D(S2)− In)D(Sk1 )Xkqk
  = max Ik∈{±1} max S2,S k 1⊆[n]
  S2,S k 1∈S
  min αk,γ≥0 max w1∈B2 K∑ k=1 ∥∥∥vTD(Sk1 )Xkqk +αTk (2D(Sk1 )− In)Xkqk + IkγT (2D(S2)− In)D(Sk1 )Xk∥∥∥ 2 |w1k|
  = max Ik∈{±1} max S2,S k 1⊆[n]
  S2,S k 1∈S
  min αk,γ≥0
  ( K∑
  k=1
  ∥∥∥vTD(Sk1 )Xkqk +αTk (2D(Sk1 )− In)Xkqk + IkγT (2D(S2)− In)D(Sk1 )Xk∥∥∥2 2
  ) 1 2
  ,
  Then, we have equation 41 ⇐⇒
  ∀i, j, min αk,γ≥0
  ( K∑
  k=1
  ∥∥∥vTD(Sk1i)Xkqk +αTk (2D(Sk1i)− In)Xkqk + IijkγT (2D(S2j)− In)D(Sk1i)Xk∥∥∥2 2
  ) 1 2
  ≤ β
  ⇐⇒
  ∀i, j, ∃αijk,γij ≥ 0 s.t.
  ( K∑
  k=1
  ∥∥∥vTD(Sk1i)Xkqk +αTijk(2D(Sk1i)− In)Xkqk + IijkγTij(2D(S2j)− In)D(Sk1i)Xk∥∥∥2 2
  ) 1 2
  ≤ β.
  We now use the same approach for the two-sided constraint as follows
  max v
  αijk,α ′ ijk≥0
  γij ,γ ′ ij≥0
  −1 2 ‖v − y‖22 + 1 2 ‖y‖22 (43)
  s.t.
  ( K∑
  k=1
  ∥∥∥vTD(Sk1i)Xkqk +αTijk(2D(Sk1i)− In)Xkqk + IijkγTij(2D(S2j)− In)D(Sk1i)Xk∥∥∥2 2
  ) 1 2
  ≤ β,
  ( K∑
  k=1
  ∥∥∥−vTD(Sk1i)Xkqk +αTijk(2D(Sk1i)− In)Xkqk + IijkγTij(2D(S2j)− In)D(Sk1i)Xk∥∥∥2 2
  ) 1 2
  ≤ β, ∀i, j.
  We note that the above problem is convex and strictly feasible for v = αijk = α′ijk = γij = γ ′ ij = 0.
  Therefore, Slater’s conditions and consequently strong duality holds (Boyd & Vandenberghe, 2004), and equation 43 can be written as
  min λij ,λ′ij≥0 max v
  αijk,α ′ ijk≥0
  γij ,γ ′ ij≥0
  −1 2 ‖v − y‖22 + 1 2 ‖y‖22 (44)
  + P1∑ i=1 P2∑ j=1 λij β −( K∑ k=1 ∥∥vTD(Sk1i)Xkqk + αTijk(2D(Sk1i)− In)Xkqk + IijkγTij(2D(S2j)− In)D(Sk1i)Xk∥∥22 ) 1 2  +
  P1∑ i=1 P2∑ j=1 λ′ij β −( K∑ k=1 ∥∥−vTD(Sk1i)Xkqk + αTijk(2D(Sk1i)− In)Xkqk + IijkγTij(2D(S2j)− In)D(Sk1i)Xk∥∥22 ) 1 2  . Next, we introduce new variables zijk, z′ijk ∈ Rh to represent equation 44 as
  min λij ,λ′ij≥0 max v
  αijk,α ′ ijk≥0
  γij ,γ ′ ij≥0
  min zijk:‖zijk‖2≤qijk z′ijk:‖z ′ ijk‖2≤q ′ ijk
  qij ,q ′ ij∈B2
  −1 2 ‖v − y‖22 + 1 2 ‖y‖22 (45)
  + P1∑ i=1 P2∑ j=1 λij
  ( β +
  K∑ k=1 ( vTD(Sk1i)Xkqk + α T ijk(2D(S k 1i)− In)Xkqk + IijkγTij(2D(S2j)− In)D(Sk1i)Xk ) zijk
  )
  + P1∑ i=1 P2∑ j=1 λ′ij
  ( β +
  K∑ k=1 ( −vTD(Sk1i)Xkqk + αTijk(2D(Sk1i)− In)Xkqk + IijkγTij(2D(S2j)− In)D(Sk1i)Xk ) z′ijk
  ) .
  Then, the strong dual of the problem equation 45 as
  min λij ,λ′ij≥0 min zijk:‖zijk‖2≤qijk z′ijk:‖z ′ ijk‖2≤q ′ ijk
  qij ,q ′ ij∈B2
  max v
  αijk,α ′ ijk≥0
  γij ,γ ′ ij≥0
  −1 2 ‖v − y‖22 + 1 2 ‖y‖22 (46)
  + P1∑ i=1 P2∑ j=1 λij
  ( β +
  K∑ k=1 ( vTD(Sk1i)Xkqk + α T ijk(2D(S k 1i)− In)Xkqk + IijkγTij(2D(S2j)− In)D(Sk1i)Xk ) zijk
  )
  + P1∑ i=1 P2∑ j=1 λ′ij
  ( β +
  K∑ k=1 ( −vTD(Sk1i)Xkqk + αTijk(2D(Sk1i)− In)Xkqk + IijkγTij(2D(S2j)− In)D(Sk1i)Xk ) z′ijk
  ) .
  Now, we can compute the maximum with respect to v,αijk,α′ijk,γij ,γ ′ ij analytically to obtain the following problem
  min λij ,λ′ij≥0 min zijk:‖zijk‖2≤qijk z′ijk:‖z ′ ijk‖2≤q ′ ijk
  qij ,q ′ ij∈B2
  1
  2 ∥∥∥∥∥∥ P2∑ j=1 DS2j P1∑ i=1 K∑ k=1 IijkD(Sk1i)Xk ( λ′ijz ′ ijk − λijzijk ) − y ∥∥∥∥∥∥ 2
  2
  + β P1∑ i=1 P2∑ j=1 ( λij + λ ′ ij ) (47)
  s.t. (2D(S2j)− In) P1∑ i=1 K∑ k=1 IijkD(Sk1i)Xkzijk ≥ 0, (2D(Sk1i)− In)Xkzijk ≥ 0,∀i, j, k
  (2D(S2j)− In) P1∑ i=1 K∑ k=1 IijkD(Sk1i)Xkz′ijk ≥ 0, (2D(Sk1i)− In)Xkz′ijk ≥ 0,∀i, j, k.
  Now we apply a change of variables and define cijk = λijzijk and c′ijk = λ ′ ijz ′ ijk. Thus, we obtain
  min cijk,c′ijk
  1
  2 ∥∥∥∥∥∥ P2∑ j=1 DS2j P1∑ i=1 K∑ k=1 IijkD(Sk1i)Xk ( c′ijk − cijk ) − y ∥∥∥∥∥∥ 2
  2
  + β P1∑ i=1 P2∑ j=1 (‖Cij‖F + ‖Cij‖F )
  s.t. (2D(S2j)− In) P1∑ i=1 K∑ k=1 IijkD(Sk1i)Xkcijk ≥ 0, (2D(Sk1i)− In)Xkcijk ≥ 0,∀i, j, k
  (2D(S2j)− In) P1∑ i=1 K∑ k=1 IijkD(Sk1i)Xkc′ijk ≥ 0, (2D(Sk1i)− In)Xkc′ijk ≥ 0,∀i, j, k.