Abstract:Computing the matrix square root and its inverse in a differentiable manner is important in a variety of computer vision tasks. Previous methods either adopt the Singular Value Decomposition (SVD) to explicitly factorize the matrix or use the Newton-Schulz iteration (NS iteration) to derive the approximate solution. However, both methods are not computationally efficient enough in either the forward pass or the backward pass. In this paper, we propose two more efficient variants to compute the differentiable matrix square root and the inverse square root. For the forward propagation, one method is to use Matrix Taylor Polynomial (MTP), and the other method is to use Matrix Pad\'e Approximants (MPA). The backward gradient is computed by iteratively solving the continuous-time Lyapunov equation using the matrix sign function. A series of numerical tests show that both methods yield considerable speed-up compared with the SVD or the NS iteration. Moreover, we validate the effectiveness of our methods in several real-world applications, including de-correlated batch normalization, second-order vision transformer, global covariance pooling for large-scale and fine-grained recognition, attentive covariance pooling for video recognition, and neural style transfer. The experimental results demonstrate that our methods can also achieve competitive and even slightly better performances. The Pytorch implementation is available at https://github.com/KingJamesSong/FastDifferentiableMatSqrt

Fast Compensated Algorithms for the Reciprocal Square Root, the Reciprocal Hypotenuse, and Givens Rotations

A Fast Compensated Algorithm for Computing Givens Rotations

A Correctly Rounded Newton Step for the Reciprocal Square Root

Low Cost Design for Elementary Function Approximation Based on Piecewise Quadratic Interpolation

Fast and accurate approximation algorithms for computing floating point square root

Fast approximate reciprocal approximations for iterative algorithms

Inexactness and Correction of Floating-Point Reciprocal, Division and Square Root

Effective Quadratic Error Bounds for Floating-Point Algorithms Computing the Hypotenuse Function

Fast Differentiable Matrix Square Root

Fast Algorithm For Composition Of The Bernstein Polynomials

Algorithms for Square Root of Semi-Infinite Quasi-Toeplitz M-Matrices

The LLL Algorithm Using Fast Givens

A General Methodology and Architecture for Arbitrary Complex Number Nth Root Computation

Efficient numerical algorithms for multi-precision and multi-accuracy calculation of the error functions and Dawson integral with complex arguments

Improving the numerical stability of fast matrix multiplication

Low-Power and Low-Latency Hardware Implementation of Approximate Hyperbolic and Exponential Functions for Embedded System Applications

A Low-Latency CORDIC Algorithm Based on Pre-Rotation and Its Application on Computation of Arctangent Function

Fast Differentiable Matrix Square Root and Inverse Square Root

Square-root RTS smoothing algorithms

GH CORDIC-Based Architecture for Computing $N$ Th Root of Single-Precision Floating-Point Number

Efficient scaling and squaring method for the matrix exponential