cd /news/artificial-intelligence/the-end-of-matrix-multiplications-we… · home topics artificial-intelligence article
[ARTICLE · art-86730] src=pub.towardsai.net ↗ pub= topic=artificial-intelligence verified=true sentiment=· neutral

The End of Matrix Multiplications: Welcome to Addition

A new wave of AI hardware and software co-design is replacing heavy floating-point matrix multiplications with addition-only, ternary weight sets {-1, 0, 1}, enabling up to a 32-fold reduction in weight memory and a 5.3-fold reduction in KV cache footprints, according to research by Ma et al. (2024) and Su et al. (2026). This shift allows trillion-parameter models to run at human-reading speeds on commodity CPU-only edge devices, bypassing the memory bandwidth bottleneck. The industry is moving toward native training of 1-bit architectures, though challenges such as dynamic activation outliers and representational collapse remain.

read18 min views1 publishedAug 4, 2026

Visualizing the foundational architectural transition from heavy floating-point matrix multiplications to lightweight, addition-only hardware pipelines.

Think of the modern silicon wafer as a highway system in a city where everyone has a private jet but the roads are permanently jammed because of a bridge toll. A few years ago, I stood in a liquid-cooled hyper-scale data center, listening to the deafening, jet-engine roar of massive server racks running transformer workloads. We were burning megawatts of electricity not to perform complex computations, but to play a glorified, expensive game of digital hot potato. High-precision weights and activations were being shuffled endlessly back and forth across tiny copper buses, hitting the physical ‘memory wall’ at every turn (Wang et al., 2024). It was in that moment of deafening hum that I realized our current paradigm of brute-forcing larger, floating-point matrices is structurally bankrupt (Wang et al., 2025).

📊 Executive Summary:The memory bandwidth bottleneck of modern autoregressive decoding can be bypassed by replacing high-precision floating-point operations (FP16/BF16) with ternary weight sets{-1, 0, 1} (Ma et al., 2024). Recent hardware-software co-designs, such as the W1A4KV2 inference standard, demonstrate up to a 32-fold reduction in weight memory (Ma et al., 2024) and a 5.3-fold reduction in KV cache footprints (Su et al., 2026). These approaches achieve lossless, multiplier-free processing at human-reading speeds natively on commodity edge CPU architectures (Wang et al., 2025; Ma et al., 2024).

The silicon straw that is breaking the back of our most beloved deep learning models is not a lack of compute cores (Wang et al., 2025). The ultimate future of deep learning is a multiplier-free architecture where matrix multiplications are entirely replaced by simple addition and subtraction (Ma et al., 2024). By compiling neural networks into extreme low-bit discrete mathematical representations, such as 1-bit or 1.58-bit ternary weights, we can unlock up to a 32-fold reduction in memory requirements (Ma et al., 2024). This allows massive trillion-parameter models to run at human-reading speed entirely on cheap, consumer-grade, CPU-only edge devices (Wang et al., 2025; Ma et al., 2024).

For years, the industry has relied on post-training quantization, shaving 16-bit weights down to 8-bit or 4-bit representations to squeeze models into smaller devices (Wang et al., 2025; Dettmers et al., 2022). While this post-hoc compression works as a temporary patch, it introduces a structural compromise, degrading the model’s accuracy as the math is aggressively flattened (Huang et al., 2026; Ma et al., 2024). The true vanguard of AI infrastructure is moving toward native training of 1-bit architectures, where the networks are designed from birth to operate within highly restricted mathematical regimes (Ma et al., 2024). This transition, however, is not a simple walk in the park, as it forces us to confront massive, systemic structural anomalies that threaten the very stability of our models (Huang et al., 2026; Ma et al., 2024). To successfully build this addition-only future, we must conquer several brutal algorithmic demons. We have to wrestle with dynamic activation outliers that spike wildly during runtime, causing uniform quantization grids to completely collapse (Dettmers et al., 2022; Huang et al., 2026). We must also prevent representational collapse, where gradients homogenize during training and turn a brilliant transformer into a pile of random noise (Huang et al., 2026; Oberländer et al., 2026). Finally, we must defuse the sprawling key-value (KV) cache bottleneck, which threatens to balloon our memory footprint out of control as sequence context windows scale to millions of tokens (Su et al., 2026; Liu et al., 2024).

This post will serve as your deep-dive manual to this architectural revolution. We will map out how the industry is systematically moving away from heavy floating-point operations toward elegant, hardware-aligned ternary arithmetic (Ma et al., 2024). We will dissect the mathematics of orthogonal rotations, the clever packing of sparse weights, and the asymmetric compression of dynamic memories (Huang et al., 2026; Ashkboos et al., 2024; Liu et al., 2024). Ultimately, we will see why the future of AI belongs not to the massive silicon foundries churning out ultra-complex multipliers, but to the clean simplicity of integer addition.

“Multiplication is the cost of computation; addition is the currency of intelligence.”— Mohit Sewak, Ph.D.

To understand why naive quantization fails so spectacularly, we must look at a fascinating phenomenon that emerges spontaneously as transformers scale: activation outliers (Dettmers et al., 2022). In the early days of deep learning, we assumed that activation values flowing through our networks were well-behaved, normally distributed tensors (Xiao et al., 2023). However, research has revealed that a tiny fraction of dimensions within the latent representation space spontaneously develop massive, highly concentrated magnitudes (Dettmers et al., 2022). In models such as T5–11B, these rare activation dimensions can reach staggering values exceeding 10⁻⁻, which are over ten thousand times larger than the median value of the surrounding tensor (Su et al., 2026; Zhao et al., 2025). Crucially, these spikes are not random artifacts; they are deeply tied to layer depth rather than parameter count alone, escalating in intensity as information travels deeper into the network architecture (Zhao et al., 2025; Dettmers et al., 2022).

Demonstrating how extreme activation outliers crush uniform quantization grids, forcing normal distribution features into zero.

🔍 Fact Check:Empirical analyses of emergent transformer anomalies show that activation outliers in deep networks like T5–11B can exceed magnitudes of 10⁻⁻, representing a value over 10,000 times larger than the median neighboring tensor (Su et al., 2026; Zhao et al., 2025). This extreme variance causes uniform quantization grids to collapse instantly.

Many early engineers looked at these massive activation spikes and assumed they were simply high-frequency noise that could be safely clipped or normalized away. This was a catastrophic mistake. Distillation and pruning studies have proved that these emergent outliers actually encode critical syntactic and semantic representations of human language (Zhao et al., 2025). They are the high-altitude peaks of the model’s cognitive landscape, housing the essential logical pathways that drive zero-shot reasoning. If you aggressively truncate or squash these anomalies, the model’s performance immediately takes a nosedive, leading to severe linguistic degradation and structural incoherence (Dettmers et al., 2022; Huang et al., 2026).

This brings us to the mathematical tragedy of the “squashing” phenomenon (Huang et al., 2026; Ma et al., 2024). When applying standard absolute maximum (absmax) scaling to map high-precision floats to an integer grid, the quantization scale factor is bound by the largest value in the tensor (Dettmers et al., 2022). Mathematically, we express this integer mapping as:

Quantized Value = round(Real Value × (Target Integer Range ÷ Tensor Max)) When the Tensor Max is a staggering 10⁻⁻ and the normal, informative features hover around a quiet 0.1, the scaling factor becomes so small that the entire non-outlier distribution is rounded directly to absolute zero (Dettmers et al., 2022). This uniform compression effectively lobotomizes the network, turning highly nuanced semantic features into a barren, silent desert (Huang et al., 2026; Ma et al., 2024).

The challenge is further compounded by a pronounced asymmetric shift budget between weights and activations. Weight outliers are inherently static; because they are fixed once the model is trained, we can easily preprocess them offline using out-of-band extraction or advanced profiling (Dettmers et al., 2022). Activations, however, are dynamic beasts produced at runtime in direct response to the specific, unpredictable sequence of incoming tokens (Dettmers et al., 2022; Wang et al., 2024). This dynamic volatility makes activation outliers the most stubborn bottleneck in the entire low-bit computing landscape (Wang et al., 2024; Oberländer et al., 2026). We cannot simply hardcode a fixed scaling factor, as the network’s internal values will fluctuate wildly depending on whether it is reading a recipe or a line of Python code.

Early attempts to bypass this problem relied on complex, dual-path mixed-precision architectures. In systems like LLM.int8(), the network dynamically identifies outlier channels at runtime and extracts them into a high-precision 16-bit matrix multiplication path, while quantizing the remaining 99.9% of normal features to 8-bit (Dettmers et al., 2022). Other strategies use Hessian-based calculations to isolate and protect the dimensions most sensitive to reconstruction error (Dettmers et al., 2022). While these hybrid approaches successfully protect model accuracy, they introduce highly irregular memory access patterns and non-uniform data formats (Huang et al., 2026; Dettmers et al., 2022). The constant branching and pipeline divergence inside the GPU significantly dilute our hardware gains, leaving us with complex code and disappointing real-world speedups (Huang et al., 2026).

Illustrating the BitNet b1.58 ternary weight architecture and Sherry 3:4 hardware packing framework.

Frustrated by the overhead of patching pre-trained networks, the research community shifted from post-training quantization to a bolder frontier: native, quantization-aware training. This paradigm shift is epitomized by the BitNet b1.58 architecture, which replaces conventional full-precision linear layers with custom BitLinear modules (Ma et al., 2024). Instead of treating quantization as an afterthought, BitLinear designs the network from the ground up to operate in extreme low-bit regimes (Ma et al., 2024). Every single parameter in the weight matrices is dynamically restricted to a ternary set of {-1, 0, 1} during the forward pass (Ma et al., 2024). This strict mathematical constraint requires exactly log₂ 3 ≈ 1.58 bits of storage per weight, giving the architecture its elegant name (Ma et al., 2024).

The structural inclusion of the 0 state is a quiet stroke of genius that distinguishes ternary networks from binary ones. In a pure binary network restricted to {-1, 1}, the model is forced to keep every connection active, which creates immense representational noise (Ma et al., 2024). The introduction of the 0 state acts as an implicit, continuous pruning agent that allows the network to natively filter out background noise (Ma et al., 2024). This simple addition dramatically preserves the model’s representational capacity without requiring complex, sparse-pruning algorithms (Ma et al., 2024). By allowing weights to choose silence, the model can focus its limited bit-width strictly on the channels that convey genuine, high-value information.

The operational consequence of this ternary restriction is nothing short of revolutionary. Because any multiplication by the set {-1, 0, 1} is mathematically equivalent to simple addition, subtraction, or pass-through operations, the hardware completely eliminates the need for power-hungry multipliers (Ma et al., 2024). Heavy matrix multiplications, which historically consumed the lion’s share of a chip’s thermal and energy budget, are transformed into element-wise additions (Ma et al., 2024). This mathematical simplification allows massive, trillion-parameter models to run at blazing-fast human-reading speeds entirely on consumer-grade CPUs (Wang et al., 2025; Ma et al., 2024). Utilizing custom runtimes like bitnet.cpp, commodity devices can comfortably decode 5 to 7 tokens per second, democratizing frontier-class AI execution (Wang et al., 2025; Ma et al., 2024).

To prove that this was not merely an academic toy, researchers trained the open-source BitNet b1.58 2B4T model on a staggering 4 trillion tokens (Ma et al., 2024). The results were a revelation to the machine learning community: the ternary model matched the perplexity and reasoning capacity of SOTA, full-precision baselines of equivalent scale (Ma et al., 2024). However, a physical roadblock quickly emerged when mapping these models to standard commodity chips. Our modern CPUs and GPUs are hardwired to process information in strict power-of-two formats, such as 2-bit, 4-bit, and 8-bit registers (Huang et al., 2026; Ashkboos et al., 2024). Forcing a 1.58-bit model into a standard 2-bit register causes a painful 0.42-bit wastage per weight, while irregular 1.67-bit packing leads to misaligned memory access that stalls compute pipelines (Huang et al., 2026; Ashkboos et al., 2024).

To resolve this tension and achieve true hardware alignment, developers have introduced highly optimized sparsity-enabled hardware packing frameworks. The Sherry framework elegantly resolves this issue by introducing a fine-grained 3:4 sparsity pattern that packs blocks of 4 ternary weights into exactly 5 bits of space, yielding a regularized 1.25-bit overall width (Huang et al., 2026). This structural adjustment completely eliminates bit wastage, restoring perfect power-of-two alignment on standard CPU architectures while providing a 10% speedup (Huang et al., 2026).

Mapping how the Arenas residual module injects gradient variance to prevent weight trapping and representational collapse.

🔍 Fact Check:The Sherry framework eliminates the 0.42-bit packing waste of ternary weights by enforcing a 3:4 fine-grained sparsity pattern. This packs four weights into exactly 5 bits to achieve a regularized 1.25-bit width, yielding a 10% speedup on commodity CPUs (Huang et al., 2026).

Meanwhile, frameworks like STBLLM explore sub-1-bit boundaries by applying N:M structural sparsity directly to binary networks (Dong et al., 2024). They utilize a Standardized Importance (SI) routing metric based on parameter and input feature norms, ensuring that only the most critical pathways are preserved (Dong et al., 2024).

When you compress a massive neural network’s parameters into a highly restricted 1-bit or ternary space, you are playing with fire. The mathematical landscape of the model becomes incredibly fragile, exposing it to the constant threat of representational collapse (Huang et al., 2026; Oberländer et al., 2026). This catastrophe occurs when the network’s parameters lose their geometric diversity and expressive capacity, causing the model’s perplexity to explode to infinity or degrade to random chance (Huang et al., 2026; Oberländer et al., 2026). It is the mathematical equivalent of a black hole, where all gradient paths converge into a singular, lifeless state. This collapse manifests differently depending on whether the system is trained from scratch or quantized post-hoc.

During Quantization-Aware Training (QAT), the primary culprit of collapse is a phenomenon known as weight trapping (Huang et al., 2026). Because we must backpropagate continuous gradients through discrete ternary parameters, we rely on a Straight-Through Estimator (STE) to bridge the mathematical gap (Huang et al., 2026). However, this approximation can cause gradients to homogenize over time, locking continuous weights into localized, sub-optimal regions (Huang et al., 2026). The continuous weights accumulate in these valleys and fail to cross the discrete quantization thresholds, permanently halting parameter updates and killing representational diversity. This issue is further exacerbated if you treat weights and dynamic activations symmetrically in early training epochs, as highly volatile activation ranges instantly trigger collapse before the network can stabilize (Oberländer et al., 2026).

For Post-Training Quantization (PTQ), the threat of collapse stems from cascading, non-independent errors. When you force a pre-trained layer into extreme low-bit representations, the resulting quantization error does not remain isolated (Oberländer et al., 2026). It propagates forward, amplifying and multiplying through every subsequent deep layer in the transformer stack (Oberländer et al., 2026). Without explicit mathematical correction, this cumulative bias builds like an avalanche, ultimately destroying the final layer’s representational coherence. Ablation studies on frameworks like GRINQH reveal that uncorrected 1-bit models suffer catastrophic perplexity explosions exceeding 10⁻⁻ on standard language benchmarks, rendering the model completely inoperable (Oberländer et al., 2026). To immunize networks against these catastrophic failures, researchers have developed a suite of pioneering stabilization techniques. The Sherry framework tackles weight trapping by introducing a specialized “Arenas” module, which functions as an Annealing residual synapse (Huang et al., 2026). This module actively injects heterogeneous gradients during the backward pass, breaking up gradient homogenization and allowing continuous weights to traverse quantization boundaries (Huang et al., 2026). By keeping the parameter distribution dynamic and diverse, the Arenas module prevents latent weights from settling into lifeless, frozen states (Huang et al., 2026). This active injection of variance ensures that the model continues to learn and adapt, even when constrained to a microscopic mathematical footprint.

Visualizing orthogonal Walsh-Hadamard rotation smoothing spiky activation outliers for uniform 4-bit quantization.

To address cumulative errors in post-training scenarios, other frameworks focus on hierarchical precision distribution. The GRINQH framework leverages hierarchical bit-slicing and dynamic post-training bias correction to mitigate cascading errors across deep layers (Oberländer et al., 2026). It carefully measures the sensitivity of each layer, dynamically allocating more precision to fragile boundary layers while aggressively compressing more robust intermediate states (Oberländer et al., 2026). Simultaneously, frameworks like NMP-QAT prevent early activation collapse by initializing dynamic activations at a safe 4-bit width (Oberländer et al., 2026). This wider range allows stable gradients to flow during the critical early epochs of training, only progressively tightening the activation precision as the network converges (Oberländer et al., 2026).

When these stabilization techniques are correctly applied during native training, an incredibly beautiful topological phenomenon occurs. Over several epochs, the ternary network naturally organizes itself into what researchers call a “quantization-valley” structure (Ma et al., 2024). The latent weights systematically migrate away from ambiguous, near-zero boundaries, polarizing themselves into distinct, highly decisive tri-modal clusters (Ma et al., 2024). This deep polarization acts as an organic, structural shield against representation collapse (Ma et al., 2024). It allows the model to withstand aggressive hardware pruning and sparse masking patterns that would immediately shatter a conventional, full-precision network.

Even if we successfully compress our model weights to 1.58 bits, we are only solving half of the puzzle. Autoregressive inference is a two-sided equation, requiring both static weights and dynamic activation states to process each incoming token (Dettmers et al., 2022; Wang et al., 2024). If our activations remain locked in full-precision formats like FP16 or INT8, our hardware registers are forced to perform expensive, mixed-precision conversions (Dettmers et al., 2022; Wang et al., 2024). The compute cores cannot execute uniform low-bit pipelines, and our theoretical hardware speedups vanish (Wang et al., 2024). To build a truly fast, efficient addition-only engine, we must conquer the challenge of low-bit activation quantization.

As we saw earlier, the primary roadblock to compressing activations is the presence of spiky, heavy-tailed outlier dimensions (Dettmers et al., 2022; Huang et al., 2026). Historically, enforcing a 4-bit activation limit caused immediate, severe accuracy drops because the tight integer bins could not accommodate these massive spikes without squashing normal features (Huang et al., 2026; Ma et al., 2024). The definitive mathematical solution to this bottleneck lies in the concept of orthogonal rotation, specifically leveraging the Walsh-Hadamard Transform (Ashkboos et al., 2024; Wang et al., 2024; Wang et al., 2025). This technique, pioneered by frameworks like SmoothQuant,

Dong, P., Li, L., Du, D., Chen, Y., Tang, Z., Wang, Q., Xue, W., Luo, W., Liu, Q.-f., Guo, Y.-T., & Chu, X. (2024). STBLLM: Breaking the 1-bit barrier with structured binary LLMs. arXiv preprint arXiv:2408.01803. https://arxiv.org/abs/2408.01803

Huang, H., Wu, D., Hu, Q., Yu, G., Yang, J., Zhu, J., Liu, X., & Wu, D. (2026). Sherry: Hardware-efficient 1.25-bit ternary quantization via fine-grained sparsification. arXiv preprint arXiv:2601.07892. https://arxiv.org/abs/2601.07892

Ma, S., Wang, H., Ma, L., Wang, L., Wang, W., Huang, S., Dong, L., Wang, R., Xue, J., & Wei, F. (2024). The Era of 1-bit LLMs: All Large Language Models are in 1.58 Bits. arXiv preprint arXiv:2402.17764. https://arxiv.org/abs/2402.17764

Ashkboos, S., Mohtashami, A., Croci, M. L., Li, B., Cameron, P., Jaggi, M., Hoefler, T., & Hensman, J. (2024). QuaRot: Outlier-free 4-bit inference in rotated LLMs. arXiv preprint arXiv:2404.00456. https://arxiv.org/abs/2404.00456

Dettmers, T., Lewis, M., Belkada, Y., & Zettlemoyer, L. (2022). LLM.int8(): 8-bit matrix multiplication for transformers at scale. In Advances in Neural Information Processing Systems (NeurIPS 2022), 35, 30318–30332. https://arxiv.org/abs/2208.07339

Oberländer, J., Finkbeiner, J., Schöfmann, C. M., & Neftci, E. (2026). GRINQH: Graded input-based quantization hierarchy for efficient LLM generation. arXiv preprint arXiv:2606.23419. https://arxiv.org/abs/2606.23419

Wang, H., Ma, S., & Wei, F. (2024). BitNet a4.8: 4-bit activations for 1-bit LLMs. arXiv preprint arXiv:2411.04965. https://arxiv.org/abs/2411.04965 Wang, H., Ma, S., & Wei, F. (2025). BitNet v2: Native 4-bit activations with Hadamard transformation for 1-bit LLMs. arXiv preprint arXiv:2504.18415. https://arxiv.org/abs/2504.18415

Xiao, G., Lin, J., Seznec, M., Wu, H., Demouth, J., & Han, S. (2023). SmoothQuant: Accurate and efficient post-training quantization for large language models. In Proceedings of the 40th International Conference on Machine Learning (ICML 2023), 38087–38099. https://arxiv.org/abs/2211.10438

Zhao, T., Singh, K. Y., Appalaraju, S., Tang, P., Wu, Y. N., & Li, L. E. (2025). On the analysis and distillation of emergent outlier properties in pre-trained language models. In Proceedings of the 2025 Conference of the Americas Chapter of the Association for Computational Linguistics: Human Language Technologies (NAACL 2025), 8475–8507. https://doi.org/10.18653/v1/2025.naacl-long.430

Liu, Z., Yuan, J., Jin, H., Zhong, S., Xu, Z., Braverman, V., Chen, B., & Hu, X. (2024). KIVI: A tuning-free asymmetric 2-bit quantization for KV cache. In Proceedings of the 41st International Conference on Machine Learning (ICML 2024). https://arxiv.org/abs/2402.02750

Su, Z., Yang, R., Zhang, C., Liu, Y., Zhang, Y., Wu, W., Xiong, J., Du, D., Zhuang, X., Qian, Y., Xie, Y., Wu, Y.-C., Yang, H., & Wong, N. (2026). OScaR: The Occam’s razor for extreme KV cache quantization in LLMs and beyond. arXiv preprint arXiv:2605.19660. https://arxiv.org/abs/2605.19660

Zhang, Y., Huang, Y., Yang, B., Venkitesh, B., Locatelli, A., Ye, H., Cai, T., Lewis, P., & Chen, D. (2024). KV cache is 1 bit per channel: Efficient large language model inference with coupled quantization. In Advances in Neural Information Processing Systems (NeurIPS 2024), 37, 3304–3331. https://arxiv.org/abs/2405.04503

Zhu, R., Zhao, Q., Zhang, Z., Ni, Z., Eshraghian, J. K., & UCSC (2023). SpikeGPT: Generative pre-trained language model with spiking neural networks. arXiv preprint arXiv:2302.13939. https://arxiv.org/abs/2302.13939

Disclaimer: The views and opinions expressed in this article are personal and do not necessarily reflect the official policy or position of any associated agencies, organizations, or the India AI Mission. AI assistance was utilized in the research, drafting, and ideation of this article. Licensed under CC BY-ND 4.0.

The End of Matrix Multiplications: Welcome to Addition was originally published in Towards AI on Medium, where people are continuing the conversation by highlighting and responding to this story.

── more in #artificial-intelligence 4 stories · sorted by recency
── more on @ma et al. 3 stories trending now
sponsored brought to you by zahid.host 4,200+ EU-deployed projects
reading about agents? ship yours in a single git push.

Run your AI side-project on zahid.host

EU-based hosting, git-push deploys, automatic HTTPS, no cold starts. Free tier with a custom domain — perfect for shipping the agent you just read about.

$git push zahid main
Live at https://your-agent.zahid.host
Get free account → Pricing
from €0/mo · no card required
LIVE [news/the-end-of-matrix-mu…] indexed:0 read:18min 2026-08-04 ·