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.

Improving MATLAB's isprime performance without arbitrary-precision arithmetic

From Worst to Average Case to Incremental Search Bounds of the Strong Lucas Test

A polytime proof of correctness of the Rabin-Miller algorithm from Fermat's little theorem

On the complexity of computing prime tables

Improving the numerical stability of fast matrix multiplication

Four primality testing algorithms

Primes Representable by Polynomials and the Lower Bound of the Least Primes in Arithmetic Progressions

Lots and Lots of Perrin-Type Primality Tests and Their Pseudo-Primes

Sound Mixed-Precision Optimization with Rewriting

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

MATLAB Simulator of Level-Index Arithmetic

Faster arbitrary-precision dot product and matrix multiplication

Efficient implementation of elementary functions in the medium-precision range

Popular differences in primes along fractional powers

On a vectorized basic linear algebra package for prototyping codes in MATLAB

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

Refinements to the prime number theorem for arithmetic progressions

The use of the invariant's properties in the primality test and prime search

Novel performant primality test on a Pell's cubic

FastECPP over MPI

Numerical calculations of primes in short intervals