Abstract:We show that assuming the availability of the processor with variable precision arithmetic, we can compute matrix-by-matrix multiplications in $O(N^2log_2N)$ computational complexity. We replace the standard matrix-by-matrix multiplications algorithm $\begin{bmatrix}A_{11}&A_{12}\\A_{21}&A_{22}\end{bmatrix}\begin{bmatrix}B_{11}&B_{12}\\B_{21}&B_{22}\end{bmatrix}=\begin{bmatrix}A_{11}B_{11}+A_{12}B_{21}&A_{11}B_{12}+A_{12}B_{22}\\A_{21}B_{11}+A_{22}B_{21}&A_{21}B_{12}+A_{22}B_{22}\end{bmatrix}$ by $\begin{bmatrix}A_{11}&A_{12}\\A_{21}&A_{22}\end{bmatrix}\begin{bmatrix}B_{11}&B_{12}\\B_{21}&B_{22}\end{bmatrix}=\Bigl\lfloor\begin{bmatrix} (A_{11}+\epsilon A_{12})(B_{11}+1/{\epsilon}B_{21})&(A_{11}+\epsilon A_{12})(B_{12}+1/{\epsilon}B_{22})\\(A_{21}+\epsilon A_{22})(B_{11}+1/{\epsilon}B_{21})&(A_{21}+\epsilon A_{22})(B_{12}+1/{\epsilon}B_{22})\end{bmatrix}\Bigr\rfloor \mod \frac{1}{\epsilon}$. The resulting computational complexity for $N\times N$ matrices can be estimated from recursive equation $T(N)=4(N/2)^2$ (multiplication of a matrix by number)+$4(N/2)^2$ (additions of matrices)+$2N^2$ (floor and modulo)+$4T(N/2)$ (recursive calls) as $O(N^2log_2N)$. The novelty of the method lies in the observation, somehow ignored by other matrix-by-matrix multiplication algorithms, that we can multiply matrix entries by non-integer numbers to improve computational complexity. In other words, while having a processor that can compute multiplications, additions, modulo and floor operations with variable precision arithmetic in $O(1)$, we can obtain a matrix-by-matrix multiplication algorithm with $O(N^2log_2N)$ computational complexity. We also present a MATLAB code using VPA variable precision arithmetic emulator that can multiply matrices of size $N\times N$ using $(4log_2N+1)N^2$ variable precision arithmetic operations. This emulator uses $O(N)$ digits to run our algorithm.

On the complexity of computing prime tables

Computing $π(N)$: An elementary approach in $\tilde{O}(\sqrt{N})$ time

More Comprehensive Computational Number Theory

Small gaps between almost-twin primes

On algorithms to find p-ordering

An Algebraic-Geometry Approach to Prime Factorization

On Primes $p$ such that $p-b$ Has a Large Power Factor and Few Other Divisors

Prime polynomials in short intervals and in arithmetic progressions

Matrix-by-matrix multiplication algorithm with $O(N^2log_2N)$ computational complexity for variable precision arithmetic

The largest prime factor of $n^2+1$ and improvements on subexponential $ABC$

Improving MATLAB's isprime performance without arbitrary-precision arithmetic

The largest prime factor of and improvements on subexponential

A New Theorem on the Prime-Counting Function

On the theory of prime producing sieves

On Faster Integer Calculations using Non-Arithmetic Primitives

From explicit estimates for the primes to explicit estimates for the Möbius function II

Linear-Time and Constant-Space Algorithms to compute Multi-Sequences that arise in Enumerative Combinatorics (and Elsewhere)

A Theorem on Prime Numbers

p-composites of the main sequence of odd numbers as building blocks for $π(x)$

A Note on the Complexity of Comparing Succinctly Represented Integers, with an Application to Maximum Probability Parsing

The sum-product problem for integers with few prime factors