Alek Sobczyk

We have a new preprint with Anastasios Zouzias on prefix sum (scan) circuits. More details and link to the paper on linkedin: linkedin.com/posts/aleksand…

Congratulations to CNN for congratulating Bennett and Brassard for the Turin award! edition.cnn.com/2026/03/18/sci…

Am I the omly one who has not received any @icmlconf 2026 papers to review yet..?

To #ICML26 submitters, out of curiosity, did you choose:

@kfountou Or to remind them to stop slacking off and review the paper properly instead of using AI to do their job :D

AI conferences should use AI to warn reviewers when their questions have already been answered in the paper.

GPT-5.2 solves our COLT 2022 open problem: “Running Time Complexity of Accelerated L1-Regularized PageRank” using a standard accelerated gradient algorithm and a complementarity margin assumption. Link to the open problem: proceedings.mlr.press/v178/open-prob… All proofs were generated by GPT-5.2 Pro. The key bounds on the algorithm’s total work (in the COLT’22 open-problem setting) have been auto-formalized using a combination of GPT-5.2 Pro, @HarmonicMath's Aristotle, and Gemini 3 Pro (High) on Antigravity. Link to the proof: github.com/kfoynt/acceler… Link to the Lean code: github.com/kfoynt/acceler… Link to the informalization of the Lean code: github.com/kfoynt/acceler… Link to my GPT-5.2 prompts: chatgpt.com/share/693e3ce6… In addition to the formalization of the main result, I checked the proof myself twice. I hope I didn’t miss anything, but if I did, please let me know and I will try to fix it. Story behind the paper and relevant work In 2016, I worked on the convergence rate of the Iterative Soft-Thresholding Algorithm (ISTA) for l1-regularized PageRank. Link to the corresponding paper: link.springer.com/article/10.100… Surprisingly, the running time of the algorithm depends only on the number of non-zero nodes at optimality. It was only natural to ask the same question for accelerated methods, such as FISTA. However, we quickly realized that FISTA activates more nodes than the number of non-zeros at optimality, even though it eventually converges to the same active set. In practice, we would still observe that FISTA is fast. Link to empirical work: uwspace.uwaterloo.ca/items/693b002d… I tried for about three months to bound the total work of FISTA and other accelerated algorithms, and from time to time I would come back to the problem while I was a postdoctoral fellow. Eventually, I gave up. I gave it another try around 2021, and I failed again. I asked my excellent former student, Shenghao Yang, and he also failed, unfortunately. I asked a couple of prominent researchers if they think the problem is solvable, they quickly mentioned that it seemed hard. We ended up publishing it as an open problem at COLT 2022. In 2023, David Martínez-Rubio et al. provided the first successful solution. Their solution is “orthogonal” to what was proved by GPT-5.2. Link to their paper: proceedings.mlr.press/v195/martinez-… I loved their work btw, I also met David in person at ICML 2024, one of the few ML conferences I ever attended. Their proposed accelerated algorithm is not necessarily faster than ISTA; however, it does offer a new trade-off between the teleportation parameter of PageRank and the total work per iteration. More importantly, the proposed method isn’t necessarily practical, since it involves solving an expensive subproblem. To be fair, in the COLT 2022 problem, we didn’t impose the additional hard constraint of using standard accelerated methods. The problem was posed as a theoretical problem. The solution proved by GPT-5.2 establishes acceleration for the standard FISTA algorithm, which performs only one gradient computation per iteration. It also offers a clean parameterization of the total work with respect to a complementarity margin, which, for certain graph structures, shows a clear speed-up compared to ISTA. In 2024, Zhou et al. (dl.acm.org/doi/10.5555/37…) gave it another go. However, in my view, their work has important drawbacks. In particular, their guarantees for accelerated localized methods (e.g., localized Chebyshev / Heavy-Ball) assume a condition on the geometric mean of certain active-ratio factors (described as Θ(\sqrt{α})) in order to obtain an accelerated bound. Two distinctions matter for our setting: First, their accelerated runtime bounds are parameterized by evolving-set quantities and a residual-ratio assumption, which can be evaluated during a run but is not typically interpretable or verifiable a priori from graph structure alone. The solution by GPT-5.2 instead provides an explicit transient-phase bound in terms of a standard optimization-structure condition, and converts this directly into a total work bound. Second, they explicitly note that FISTA-style acceleration violates the monotonicity property needed to bound the per-iteration accessed volume, and emphasize that guaranteeing intermediate sparsity in accelerated frameworks is challenging. The margin-based analysis by GPT-5.2 directly targets this gap: even without any monotonicity of intermediate supports, GPT-5.2 bounded how much spurious activation can occur before the iterates enter a neighborhood of the unique minimizer, thereby yielding a concrete locality certificate for the accelerated proximal-gradient trajectory. Since 2024, every time OpenAI or Google released a new major model, I would give it a go. This time, with GPT-5.2, it seems to have worked.

@JannisBorn wrote a very nice teaser about our Neurips paper Quantum Doubly Stochastic Transformers (spotlight). Check it out in San Diego, and in EurIPS. You can find links to the paper, video, and poster below. #quantum #LLMs #Transformers linkedin.com/posts/jannis-b…

We have a new preprint with Almudena Carrera Vazquez @anedumlACV, one of the last (and favorite) works that I contributed to before leaving IBM. It is about quantum algorithms for approximating spectral gaps between consecutive pairs of eigenvalues. Link: arxiv.org/abs/2508.21002

Dear everyone, I have been using X less and less lately. Let's stay connected on Bluesky and/or Linkedin: bsky.app/profile/asobcz… linkedin.com/in/aleksandros…

An Exponential Separation Between Quantum and Quantum-Inspired Classical Algorithms for Machine Learning. arxiv.org/abs/2411.02087

Short career update: linkedin.com/pulse/%25CE%25… via @LinkedIn

The chance that this will occur in a real application is probably extremely rare, but it should be handled with care.

3) A numerical bug of the cuSPARSE gtsv solver that we reported in the 2015 paper is still there..! By planting a small adversarial non-invertible block of size 32*32, gtsv2 struggles to solve the linear system, while both gspike and scipy solve it to 16 digits of accuracy.

A little less than ~10 years ago we developed a cuda kernel to solve tridiagonal linear systems of equations. I was recently trying out WSL on my personal laptop, which has a dedicated GPU, and I thought to try to run that rusty code and see what happens. Some findings:

I will never complain about slow review process again.

I want to also thank the anonymous reviewers for quality feedback, and all the people who provided helpful discussions (listed in the preprint), especially Ryan Schneider and Daniel Kressner who promptly (and independently) pointed a technical detail in the initial draft.

A preprint can be found here: arxiv.org/abs/2311.10459 Appendix F can potentially be of interest to #compchem people as it is related to DFT.

I'm happy to share that a revised version of our paper has been accepted at #NeurIPS2024. I want to thank my co-authors, Marko Mladenović and my phd supervisor Mathieu Luisier for the excellent collaboration, and my manager @teodorolaino for all the support.

@kfountou Congrats Kimon and the team!

This paper was just accepted at NeurIPS. I am particularly happy about this because it originated as a course project by Robert in our Graph Neural Networks course.