REVIEW 4 major objections 4 minor 2 cited by
An Adaptive Orthogonal Convolution Scheme for Efficient and Flexible CNN Architectures
T0 review · 4 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read A kernel-composition scheme builds orthogonal convolutions that natively support stride, transposition, groups, and dilation.
desk verdict AOC is a genuinely useful orthogonal-convolution construction with native stride and other modern features, proven under circular padding; the zero-padding accuracy setting is only quasi-orthogonal, and the overhead numbers are slightly oversold. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing mechanism is the block-convolution operator $\circledast$, which computes the equivalent kernel of two composed convolutions without executing either one: $(B \circledast A) \star x = B \star (A \star x)$, producing a kernel of size $k_A + k_B - 1$. AOC pairs a BCOP factor — orthogonal at stride one for any kernel size, built from stacked symmetric projectors $N$ and $I - N$ — with an RKO factor, orthogonal precisely when its kernel size equals its stride, a fact the paper proves by reordering the strided convolution into a permutation matrix followed by a matrix product. The internal channel dimension $c = \max(c_i, \lfloor c_o/s^2 \rfloor)$ is what aligns the two factors' orthogonality types, and an efficient implementation (batched block-convolution and a parallel associative scan) is what makes the construction cheap enough for large images.
What would settle it
Build the Toeplitz matrix of a trained AOC layer on a small image ($8\times8$) using zero padding and compute its singular values: any singular value departing from 1 by more than $10^{-4}$ would falsify exact orthogonality in deployed settings. A complementary check is to verify the reported 'AOC accurate' result directly — the paper itself reports 0% certified accuracy at radius $36/255$, which a padding-aware bound should reproduce.
Extended reading notes
Core claim
The paper's central claim is that the kernel $K_{\mathrm{AOC}} = K_{\mathrm{RKO}} \circledast K_{\mathrm{BCOP}}$, with internal channel count $c = \max(c_i, \lfloor c_o / s^2 \rfloor)$, is an orthogonal convolution with stride $s$ and kernel size $k$. The construction fuses a BCOP kernel — built by stacking symmetric projectors into $1\times2$ and $2\times1$ orthogonal convolutions and composing them — with an RKO kernel, a kernel reshaped from an orthogonal matrix, which this paper proves to be orthogonal exactly when its kernel size equals its stride. The choice of $c$ ensures that both factors are row-orthogonal or both column-orthogonal, so the composition inherits orthogonality from the corresponding property of orthogonal matrices. From this core result the paper derives native transposed orthogonal convolutions (the transpose of an orthogonal convolution is again orthogonal), grouped orthogonal convolutions (the Toeplitz matrix is block diagonal), and dilated orthogonal convolutions, and it notes that the covered parameter range $s \le k$ is exhaustive because no orthogonal kernel exists when the stride exceeds the kernel size.
Load-bearing premise
The proofs assume circular padding at the image boundary, so with the standard zero padding used in most architectures the layer is only quasi-orthogonal along the borders, and the exact norm preservation that tight robustness certificates rely on is not guaranteed.
Editorial extensions
If this is right
- Native strided orthogonal convolutions replace stride emulation through input reshaping, cutting the parameter count of strided layers by a factor of $s^2$ and removing the cubic-cost blowup of emulated methods.
- Transposed orthogonal convolutions become available for upsampling, so U-Nets, VAEs, and invertible normalizing-flow layers can be built from norm-preserving components.
- Grouped and dilated orthogonal convolutions let depthwise and large-receptive-field architectures stay inside the 1-Lipschitz framework where robustness certificates are tight.
- The reported scaling behaviour — the overhead falls as batch size grows, reaching about 13% on ImageNet-scale ResNet-34 — puts certifiably robust training on large datasets within a practical compute budget.
- Because the covered configurations include every case where an orthogonal kernel can exist (stride no larger than kernel size), the method is a full coverage of the feasible parameter space rather than a partial patch.
Reading between the lines
- One consequence the authors leave implicit: with the circular-padding assumption in the proofs, any zero-padding deployment inherits a border effect, so certificates should be computed with a padding-aware Lipschitz bound rather than trusting exact orthogonality; the paper's own high-accuracy run reports 0% provable accuracy at radius $36/255$.
- If the scalability numbers hold, the practical blocker for certified robustness at ImageNet scale shifts from layer cost to data and training budget, which are the constraints already identified by the robust-training literature.
- A natural follow-up would be to swap the BCOP factor for a complete parametrization of orthogonal convolutions (BCOP provably misses some of them); the same $\circledast$ machinery would then add striding to a strictly more expressive base class.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes AOC (Adaptive Orthogonal Convolution), a construction that composes a BCOP kernel with an RKO kernel through the block-convolution operator to obtain convolutional layers that are orthogonal while supporting native stride, transposition, grouping, and dilation. The main theoretical claims are Propositions 2.4, 2.5, 2.7, and 2.8, with proofs in Appendix G. The paper also contributes an efficient implementation (the Orthogonium package) and evaluates AOC on CIFAR-10 and ImageNet-1K for certified robustness and on a ResNet-34 scalability benchmark. The authors report that AOC achieves competitive robustness certificates in the circular-padding setting and much lower training overhead than prior orthogonal convolution methods at scale.
Significance. If the orthogonality claims hold, AOC is a practically valuable contribution: it is, to my knowledge, the first explicit construction that natively supports strided, transposed, grouped, and dilated orthogonal convolutions, and it comes with an open-source implementation and an unusually extensive unit-test suite (1442 tests reported in Appendix D). The proofs in Appendix G are mostly clear and are stated under explicit assumptions. The main caveat is that the orthogonality proofs assume circular padding, whereas one of the two experimental configurations uses zero padding and is only quasi-orthogonal; this boundary issue is acknowledged in Appendix H but is not reflected in the abstract and introduction, where the method is presented as delivering strict orthogonality for standard CNN use. The scalability results are encouraging but are reported without statistical variation and with a small numerical inconsistency against the abstract.
major comments (4)
- [Section 2.1, Appendix H, Table 2] The orthogonality proofs are explicitly conditioned on circular padding (Section 2.1: 'we assume circular padding in all proofs'), yet the 'AOC accurate' configuration in Appendix H uses zero padding and is described as '1-Lipschitz and quasi-orthogonal (orthogonal everywhere except for the images border)'. Under zero padding the strict orthogonality used in Proposition 2.3 does not hold, and Table 2 shows the consequence: 'AOC accurate' attains 91.5% clean accuracy but 0% provable accuracy at radius 36/255. The abstract and introduction nevertheless present AOC as providing strictly orthogonal layers for standard CNN architectures. This is a load-bearing scope mismatch that must be fixed by either restricting the orthogonality and certification claims to the circular-padding setting, adding a complete treatment of the zero-padding boundary, or clearly repositioning the zero-padding experiments as quasi-orthogonal rather than orthogonal.
- [Abstract and Section 3.2 / Table 3] The abstract states that AOC incurs 'only a 10% slowdown compared to unconstrained models in realistic Imagnet 1k training setup', but Section 3.2 reports 'a training time only 13% slower than its unconstrained counterpart' and Table 3 lists 1.75x at batch 128 and 1.13x at batch 512. The numbers should be made consistent, and the exact configuration (batch size, architecture, hardware) behind the advertised percentage should be stated in the abstract or immediately in Section 3.2.
- [Table 2 and Appendix H] The robustness comparison in Table 2 mixes numbers quoted from other papers that use different architectures, training recipes, loss parameters, and data augmentation (see Appendix H, Table 4), and the 'AOC robust*' row uses extra data while other rows are marked with an asterisk only in the footnote. The comparison is therefore indicative rather than controlled. I recommend reporting per-configuration details in the table caption, adding error bars or multiple seeds for the AOC rows, and marking which rows use extra data directly in the table.
- [Appendix G.5] The proof of Proposition 2.5 chooses the internal channel c = max(ci, floor(co/s^2)), but the inequalities in the proof are stated in terms of co/s^2 rather than floor(co/s^2). The argument is valid in the integer cases shown, but the proof should explicitly handle the non-integer case to justify that floor(co/s^2) always satisfies the required row- or column-orthogonality condition. This is a small but load-bearing gap in the central construction.
minor comments (4)
- [Throughout] There are several typos and infelicities, including 'Adaptative' in the title, 'Imagnet 1k' in the abstract, 'smae setup' in Appendix C.2, and 'we evaluate' beginning a sentence with a lowercase letter in the caption of Table 2. These should be corrected in a final pass.
- [Table 3] The scalability table reports single timing measurements without error bars or repeated trials. Since the central scalability claim depends on small overhead differences, at least a few repeated runs with standard deviations should be reported for the headline batch sizes.
- [Section 3.1 / Observation 1] The sentence 'This comes at the cost of 0% provable accuracy (certificates smaller than the radius ϵ = 36/255)' is accurate, but the text should explicitly connect this to the quasi-orthogonality under zero padding rather than presenting it only as an accuracy-robustness trade-off.
- [Appendix G.4] The RKO proof assumes that the input height and width are multiples of the stride and that the input is partitioned into non-overlapping blocks. This assumption should be stated in the proposition itself, not only in the proof, so readers know the exact domain of validity.
Circularity Check
No significant circularity: the AOC orthogonality result is proved from the Toeplitz/block-convolution algebra rather than assumed or fitted, and the self-citations are not load-bearing.
full rationale
The central claim (Proposition 2.5, Appendix G.5) derives orthogonality of KAOC = KRKO ? KBCOP by direct matrix algebra: after fixing c = max(ci, floor(co/s^2)), the two factors are shown to be both row- or both column-orthogonal, so Proposition 2.3 applies. RKO's orthogonality for k = s is proved in the paper via a permutation argument (Appendix G.4), and BCOP orthogonality is reproduced from Li et al. with proof in Appendix G.3; neither is assumed as a black box. Transposed, grouped, and dilated cases follow from the same Toeplitz/block-diagonal structure (G.6-G.8). The evaluation uses external benchmarks and comparisons, so no fitted parameter is renamed as a prediction. Self-citations exist (RKO originates in Serrurier et al. 2021, with overlapping authors, and the s > k impossibility is attributed to Achour et al. 2022, co-authored by Mamalet), but these are used as provenance or as external theorems, not as evidence for the key k = s orthogonality step, which is proved in this paper. Appendix H does acknowledge that the zero-padding 'AOC accurate' setting is only quasi-orthogonal and achieves 0% provable accuracy at radius 36/255; this is a clearly stated limitation of the experimental setup, not a circular derivation. The orthogonality proofs are conditional on circular padding, and that condition is explicit, so the derivation is self-contained rather than circular.
Assumptions & free parameters
assumptions (4)
- standard math The block-convolution operator ⊛ correctly represents composition of convolutions and is associative.
- domain assumption Convolution is represented as a Toeplitz matrix and orthogonality is defined via that matrix, assuming circular padding.
- domain assumption The BCOP parametrization is incomplete, and AOC inherits this limitation.
- domain assumption No orthogonal convolution exists when stride s > kernel size k.
Cite this review
Pith. "Pith review of An Adaptive Orthogonal Convolution Scheme for Efficient and Flexible CNN Architectures." pith.science (2026). https://pith.science/paper/DYKFZ5EJ
@misc{pith2026250107930,
author = {Pith},
title = {Pith review of: An Adaptive Orthogonal Convolution Scheme for Efficient and Flexible CNN Architectures},
year = {2026},
howpublished = {\url{https://pith.science/paper/DYKFZ5EJ}},
note = {Machine review of arXiv:2501.07930}
}
read the original abstract
Orthogonal convolutional layers are valuable components in multiple areas of machine learning, such as adversarial robustness, normalizing flows, GANs, and Lipschitz-constrained models. Their ability to preserve norms and ensure stable gradient propagation makes them valuable for a large range of problems. Despite their promise, the deployment of orthogonal convolution in large-scale applications is a significant challenge due to computational overhead and limited support for modern features like strides, dilations, group convolutions, and transposed convolutions. In this paper, we introduce AOC (Adaptative Orthogonal Convolution), a scalable method that extends a previous method (BCOP), effectively overcoming existing limitations in the construction of orthogonal convolutions. This advancement unlocks the construction of architectures that were previously considered impractical. We demonstrate through our experiments that our method produces expressive models that become increasingly efficient as they scale. To foster further advancement, we provide an open-source python package implementing this method, called Orthogonium ( https://github.com/deel-ai/orthogonium ) .
Figures
Figures from the paper (2 more)
Forward citations
Cited by 2 Pith papers
-
Distributed Retraction-Free and Communication-Efficient Optimization on the Stiefel Manifold
EF-Landing provably converges at O(1/sqrt(N K)) for distributed stochastic problems on the Stiefel manifold while using compressed communication and no retraction.
-
HOFT: Householder Orthogonal Fine-tuning
HOFT and SHOFT fine-tune foundation models with two Householder-built orthogonal matrices, matching or beating LoRA, DoRA, OFT, BOFT and HRA on reasoning, translation, image generation and math.
Reference graph
Works this paper leans on
-
[1]
Checking of numerical Instabilities: Issues arising from floating-point precision, such as those introduced by small epsilon values added to avoid division by zero
-
[2]
Checking the orthogonality of a layer under stride, group, transposition, and dilation conditions
Checking for implementation discrepancies: Differences between mathematical formalism and its translation to popular frameworks (e.g., SOC proofs assume circular padding, while its implementation uses zero padding). Checking the orthogonality of a layer under stride, group, transposition, and dilation conditions. The numerical stability and the convergenc...
-
[3]
increase the margin, increase the model size and number of epochs until it reaches 100% training accuracy
-
[4]
repeat steps 2. and 3. until the desired robustness is achieved. The target margin is m = 3 2 √ 2 to match the parameters used by (Araujo et al., 2023). In practical contexts, m should be set to m = ϵ √ 2 αL where ϵ is the targeted robustness radius, L the model’s Lipschitz constant and α the scaling factor applied on the data (for instance when standard ...
work page 2023
-
[7]
This method, though accurate, is computationally expensive for large input images
Explicit SVD on Toeplitz Matrices: Using the impulse response approach, we construct the Toeplitz matrix for any padding and stride, allowing direct computation of singular values. This method, though accurate, is computationally expensive for large input images
-
[8]
Acc.” refers to the high-accuracy setting, and “Rob
Product Bound for BCOP and RKO Kernels: The upper bound for the BCOP kernel is computed using standard methods, while the SVD of the reshaped RKO kernel is used for direct evaluation. Unit testing of the implementation. We used both of these two approaches in our unit tests. This enables us to ensure that the second method (which is faster and more scalab...
work page 2022
-
[9]
start with a given architecture and a low margin m = 0
-
[10]
The model reached its maximum generalization for this margin
increase data augmentation until train accuracy goes below 100%. The model reached its maximum generalization for this margin
Show all 13 references
-
[456]
Kiani, B., Balestriero, R., LeCun, Y ., and Lloyd, S
pmlr, 2015. Kiani, B., Balestriero, R., LeCun, Y ., and Lloyd, S. projunn: efficient method for training deep networks with unitary matrices. Advances in Neural Information Processing Systems, 35:14448–14463, 2022. Kingma, D. P. Auto-encoding variational bayes. arXiv preprint ...
2015 arXiv
-
[582]
B´ethune, L
PMLR, 2019. B´ethune, L. Deep learning with Lipschitz constraints. PhD thesis, Universit´e de Toulouse, 2024. 9 An Adaptive Orthogonal Convolution Scheme for Efficient and Flexible CNN Architectures B´ethune, L., Boissin, T., Serrurier, M., Mamalet, F., Friedrich, C., and Gonz...
2019 arXiv
-
[1820]
Li, Q., Haque, S., Anil, C., Lucas, J., Grosse, R
URL http://eudml.org/doc/203444. Li, Q., Haque, S., Anil, C., Lucas, J., Grosse, R. B., and Ja- cobsen, J.-H. Preventing gradient attenuation in lipschitz constrained convolutional networks. Advances in neural information processing systems, 32, 2019. Li, Y ., Zhang, X., and C...
2019 arXiv
-
[2022]
Although already proven by previous work, proof of this can be found in Appendix G.3
noted that additional constraints are needed, proposing a half-rank symmetric projector to construct a1 × 2 orthogonal convolution: from a column-orthogonal matrix M ∈ Rc× c 2 5, the matrix N = M MT ∈ Rc×c is a symmetric projector that satisfies: N = N 2 = N T and (I − N ) = (...
2022
-
[2023]
The residual connections were removed, and a constant value of 1.15 replaced the scalar factors
with the original convolutions replaced by AOC convolutions (which include downsampling convolutions). The residual connections were removed, and a constant value of 1.15 replaced the scalar factors. This resulted in a network with a tighter Lipschitz constant evaluation, whic...
2009
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.