__ Summary and Contributions__: ***Update after rebuttal***
I have read authors response. My opinion remains that this is a good paper and I recommend acceptance. Even though Lemma 1 can be derived without much problem from ref [5], It provides crucial theoretical grounding to the paper. As far as I know, all applications of such bounds rely precisely on their tightness, so I don't think comparison of the methods in robustness certification (or other) applications would necessarily complement the paper. This type of "verification" problems are much harder than training, so in my opinion it is ok to see experiments with smaller networks than used in industry.
***
The papers presents a hierarchical optimization approach, based on polynomial optimization, to compute an upper bound on the Lipschtiz constant of ReLU networks. The method is computationally expensive and so heuristics are employed to provide algorithms that are able to work with larger instances of networks, while always providing a valid upper bound on the Lipschtiz constant.

__ Strengths__: The semialgebraic description of the ReLU and its derivative is a tighter description than the ones appearing in previous work, thus they can potentially provide tighter bounds. Precisely, one of the main strengths of this paper is that it builds upon previous work, which introduced a polynomial optimization approach to provide upper bounds on Lipschitz constants. Thus, it has strong theoretical roots and many more improvements might be derived in the future.
It is also remarkable that the non-smoothness is properly handled, and in particular, the issues with the incorrectness of backpropagation on ReLU networks is correctly addressed, thus filling a theoretical hole usually present and overlooked in many other works.
As another strength, on the practical side the implementation of the methods seem to be more efficient in their use of memory and computing power, as the reported times and memory consumption, compared to the baselines. are much better. It would be a good idea to publicly share the code for the community to use.

__ Weaknesses__: One of the main weaknesses is the lack of comparison with the local lispchitz constant estimation baseline. If the code for the baseline was not available perhaps it could have been implemented easily?

__ Correctness__: The claims are properly supported with strong theory, I have checked proof of lemma 1 in appendix. The empirical methodology seems ok.

__ Clarity__: The paper writing is ok but some parts might benefit from improving readability. In particular the introduction seems a bit lacking in style. Also the abstract and title might need slight changes in wording.
Also there is some debatable claim "However this method is restricted to the L2-norm whereas most
deep learning applications are rather concerned with the Lā-norm." : the prevalence of the L_infinity norm is found in the robust-deep-learning literature. I dont think that it is ok to claim that such norm is 'globally' more important than the l2 norm in deep learning applications.

__ Relation to Prior Work__: The paper compares to previous work that follows the same polynomial optimization-based approach, while also mentioning and referencing other type of algorithms for this task. The differences are clearly stated.

__ Reproducibility__: Yes

__ Additional Feedback__:

__ Summary and Contributions__: This paper writes the ReLU and its generalized derivative in terms of polynomial inequalities, using this to write down a polynomial optimization problem whose value exactly equals the Lipschitz constant of a ReLU network, and gives a way of approximately solving it using a relaxation of Lasserre's hierarchy.

__ Strengths__: The underlying idea is very nice (expressing the ReLU and its generalized derivative with polynomial inequalities, and utilizing this in the optimization problem calculating the network's Lipschitz constant). I could see this being an important technique for analyzing ReLU networks, and as a core contribution is very deserving of praise. Also, it is noteworthy that the authors showed this formulation could be exploited by via a relaxation of the Lasserre/SDP hierarchy. The results compare quite favorably to reference 21, an LP relaxation for computing the Lipschitz constant.

__ Weaknesses__: It's difficult to see how this would work for modern neural nets, because the authors don't really go into detail on what happens with more layers, and how well the heuristic relaxation fares in that case (it looks like it might not fare that well because the relaxation appears to be dropping a lot of things with each layer).
Though, to be fair, it seems the previous LP method (which was accepted to ICLR) seems to suffer the same issue of scaling deeper, but because of computational infeasibility.
Also, from what I can tell the relaxation doesn't necessarily provide a valid upper bound (a problem from the SDP hierarchy might not, and a relaxation of it is even less likely to) -- wish this was commented on more.
*Edit* -- my above remark is wrong, yes their method provides a valid upper bound because it is a relaxation of a maximization problem.

__ Correctness__: The application of the reference generalizing differentiation to Lemma 1 was nice. Overall, the application of the Lasserre hierarchy and the relaxation used make sense.

__ Clarity__: As someone not too familiar with the area, I initially found the paper hard to understand, and had to read some of the textbook in the references to get a basic grasp on what was going on. Perhaps the paper is quite clear to someone with more experience with polynomial optimization, though. Actually, the writing is quite efficient, and I don't think there are any major omissions of information, but perhaps the paper could be written a bit friendlier.
Suggestions for clarity:
- When describing table 1 in lines 89-94, make it clear whose methods are what.
- To make the paper more self-contained, if there is not space in the main paper, I would suggest explaining the features of MomOpt-d i.e. what are the moment matrices, how is this a relaxation, etc. Or, to be more specific in what parts of the polynomial optimization textbook are relevant.
- line 194, I think x0 should be x_0?
- For MomLCEP-2, why not just have the constraint Ax+b=z?
- I think more context-providing statements would be good in general, e.g. when introducing Lasserre hierarchy, say you will use a relaxation of it to approximately solve LCEP, and are just giving an overview; more commentary on why you are considering 1-hidden layer in the main paper (in particular extending to multiple layers seems hacky, which would seem to make applying this to W-GAN difficult)

__ Relation to Prior Work__: The inclusion of the semialgebraic constraints determining the ReLU operation and its generalized derivative is a clear improvement upon the optimization problem posed by the paper this paper compares itself to. It's great that this paper provides a unique way of specifically addressing networks that use the ReLU, which is an important activation function. Furthermore, it makes sense why their method (with the heuristic relaxation) would be much faster than the one with the LP.

__ Reproducibility__: Yes

__ Additional Feedback__: * I reduced my score from 8 to 7 because of issues other reviewers brought up about lacking practicality; I agree with them that it would be good to show an application of the improved Lipschitz constant estimation on a downstream task. Also, perhaps the relaxation of the sparse structural constraints could be explained more clearly, and perhaps writing out a full example, or comparing numbers of variables/constraints of this reference and others, etc., could aid clarity. However, I still think the exact formulation of the Lipschitz constant is very nice.

__ Summary and Contributions__: This paper proposes a semidefinite programming hierarchy heuristic approach to compute the upper bounds on the Lipschitz constant of neural networks. In deriving this heuristic, the paper also proposes the semialgebraic expression of the generalized derivative of ReLU function. Empirical results show that the proposed method provides a trade-off between better bounds and the solver running time compared to the state-of-the-art methods.
This is a mediocre paper, and I am leaning to reject it because (1) there is a lack of experimental results on robustness certification or other applications that need the estimation of the Lipschitz constant of neural networks, and (2) even the experimental results provided in the paper do not show a clear advantage of the proposed method compared to other baseline methods.

__ Strengths__: The paper addresses an interesting topic of estimating the Lipschitz constant of a neural network, which is then useful for robustness certification. The background work and the proposed method are clearly presented.

__ Weaknesses__: First, other works in this area usually show the advantage of the estimated Lipschitz constant on robustness certification tasks or other downstream tasks, empirically or theoretically (see [1], [2], and [3]). Without those experiments and only looking at the upper bounds of the Lipschitz constant of neural networks, it is hard to evaluate the significance of the work.
Second, the upper bounds of the Lipschitz constant and the corresponding computational times also do not clearly show the advantage of the proposed method. For example, in table 2, when computing upper bounds of global Lipschitz constant and the solver running time on the trained network, while the proposed method achieves better upper bounds, it requires much more computational time. Again, this further emphasizes the need of additional experiments on downstream tasks such as robustness certification to justify the significance of the method.
[1] Aditi Raghunathan, Jacob Steinhardt, and Percy Liang. Semidefinite relaxations for certifying robustness to adversarial examples. In Advances in Neural Information Processing Systems, pages 10877ā10887, 2018.
[2] Yusuke Tsuzuku, Issei Sato, and Masashi Sugiyama. Lipschitz-margin training: Scalable certification of perturbation invariance for deep neural networks. In Advances in Neural Information Processing Systems, pages 6541ā6550, 2018.
[3] Krishnamurthy Dvijotham, Robert Stanforth, Sven Gowal, Timothy A Mann, and Pushmeet Kohli. A dual approach to scalable verification of deep networks. In UAI, pages 550ā559, 2018.

__ Correctness__: The claims and method are correct. Additional experiments are needed.

__ Clarity__: The paper is well-written.

__ Relation to Prior Work__: It is clearly discussed how the proposed method differs from previous works.

__ Reproducibility__: Yes

__ Additional Feedback__: I am willing to increase my score if the authors provide experiments on robustness certification or other downstream tasks.
Post Rebuttal Comments: This paper is not a theory paper. Also, it shows no practical relevance. After reading the paper, the reviews, and the discussion carefully, I decided to decrease my score for this paper from 4 to 3. My main concerns are: 1) the lack of experiments on downstream tasks such as robustness certification, and 2) the upper bounds of the Lipschitz constant and the corresponding computational times in Table 1 do not clearly show the advantage of the proposed method. The authors did not convince me on these two points in the rebuttal.

__ Summary and Contributions__: The manuscript proposed a polynomial optimization formulation to estimate the global and local Lipschitz constant of a multi-layer neural network. The polynomial optimization is then approximated using semidefinite programs based on both dense and sparse Lasserre's Hierarchy. Numerical experiments are performed to show the proposed methods' superior performance compare with LP based relaxations studied in [21]

__ Strengths__: This is a solid work with contributions mostly in the theoretical side. The proposed optimization can be used to numerically estimate the local and global Lipschitz constants (which are important for various neural network applications) for small- to medium- sided networks with ReLU activation.

__ Weaknesses__: The proposed methods and algorithms can only be applied to small networks. There remain tremendous challenges scaling the methods to networks of practical sizes. Most importantly, as the authors commented in the conclusion, the heuristic relation is designed for QCQP that applies to 1-hidden layer networks. (Though extensions are considered in Appendix E). All the numerical simulations are also performed for neural networks with at most one-hidden layer.

__ Correctness__: The claims and method are technically sound and correct to my knowledge.

__ Clarity__: Yes, the paper is well and clearly written.

__ Relation to Prior Work__: Comparison to LP based methods in [21] is extensively performed. Other relevant work are reviewed.

__ Reproducibility__: Yes

__ Additional Feedback__: Some minor comments:
- It's good to define the abbreviations of the algorithms shown in Table 1 in the main text.
- What is the unit (seconds I guess) for running time reported in Table 1.
- Define LBS before it's used in Table 1.
- Is HR-1 the same as SHOR?
- v in ||v||_* of equation (2) should be bold
- Vector-vector multiplications are used in several places, e.g., equation (4). I think they are component wise, but it's good to precisely state them somewhere.
- Perhaps include a reference for SDP in footnote 1 on page 4 if you think it's necessary to define SDP for the readers.
- Line 132 on page 4, what is y_i?
- Define LCEP in the main text when it's first used in line 130 on page 4 (even though you defined it in the section title).

__ Summary and Contributions__: [ QUESTION-ONLY NON-REVIEW FROM AC. SCORE IS FAKE. BUT PLEASE ANSWER! ]
1. Can the authors please clarify the theoretical contribution? The exact method appears direct from the definition of gradient, and the heuristic is based on the SoS hierarchy, but oddly, does not carry any approximation guarantees. Is there a hope of an approximation guarantee?
2. Can the authors please clarify the empirical contribution, specifically, a way this could be applied to standard networks? The present demonstration is on 80x80; can the authors predict what happens, including a time estimate, on something like AlexNet, VGGNet, ResNet? If the answer is large, can the others say how the methods already given in the paper are a significant step in lowering these numbers? Does any of the competing work handle these large practical networks?

__ Strengths__: .

__ Weaknesses__: .

__ Correctness__: .

__ Clarity__: .

__ Relation to Prior Work__: .

__ Reproducibility__: Yes

__ Additional Feedback__: