Related papers: Computing the Fr\'echet Derivative of the Polar De…
We study several variants of decomposing a symmetric matrix into a sum of a low-rank positive semidefinite matrix and a diagonal matrix. Such decompositions have applications in factor analysis and they have been studied for many decades.…
The Riemann-Liouville formula for fractional derivatives and integrals (differintegration) is used to derive formulae for matrix order derivatives and integrals. That is, the parameter for integration and differentiation is allowed to…
Many applications in computational science require computing the elements of a function of a large matrix. A commonly used approach is based on the the evaluation of the eigenvalue decomposition, a task that, in general, involves a…
Computing the polar decomposition and the related matrix sign function has been a well-studied problem in numerical analysis for decades. Recently, it has emerged as an important subroutine within the Muon optimizer for training deep neural…
Let A and E be Hermitian self-adjoint matrices, where A is fixed and E a small perturbation. We study how the eigenvalues and eigenvectors of A+E depend on E, with the aim of obtaining first order formulas (and when possible also second…
We present a derivation of the gradients of feedforward neural networks using Fr\'echet calculus which is arguably more compact than the ones usually presented in the literature. We first derive the gradients for ordinary neural networks…
A novel construction of polar codes with dynamic frozen symbols is proposed. The proposed codes are subcodes of extended BCH codes, which ensure sufficiently high minimum distance. Furthermore, a decoding algorithm is proposed, which…
Given a square, nonsingular matrix of univariate polynomials $\mathbf{F} \in \mathbb{K}[x]^{n \times n}$ over a field $\mathbb{K}$, we give a fast, deterministic algorithm for finding the Hermite normal form of $\mathbf{F}$ with complexity…
We establish pointwise formulas for the shape derivative of solutions to the Dirichlet problem associated with the fractional Laplacian. Specifically, we consider the equation $(-\Delta)^s u = h$ in $\Omega$ and $u=0$ in $\Omega^c$, where…
In the present paper we consider the functional dissipativity of the Dirichlet problem for systems of partial differential operators of the form $\partial_{h} ({\mathop{\mathscr A}\nolimits}^{hk}(x)\partial_{k})$ (${\mathop{\mathscr…
A new method is presented for Fourier decomposition of the Helmholtz Green Function in cylindrical coordinates, which is equivalent to obtaining the solution of the Helmholtz equation for a general ring source. The Fourier coefficients of…
The unitary polar factor of a matrix F is the unitary matrix Q realizing the minimum of the norm of F-Q over all unitary matrices Q. Tracing back the development on the optimality of the polar factor to its presumable roots, in this paper…
QR factorisation plays an important role in matrix computations. Within the context of optimisation and of automatic differentiation of such computations, we need to compute the derivative of this factorisation. For tall matrices, however,…
There are various definitions of right and left polar decompositions of an $m\times n$ matrix $F \in \mathbb{K}^{m\times n}$ (where $\mathbb{K}=\mathbb{C}$ or $\mathbb{R}$) with respect to bilinear or sesquilinear products defined by…
Fractional derivative can be defined as a fractional power of derivative. The commutator (i/h)[H, ], which is used in the Heisenberg equation, is a derivation on a set of observables. A derivation is a map that satisfies the Leibnitz rule.…
We give an effective procedure to explicitly find the decomposition of a polarized abelian variety into its simple factors if a period matrix is known. Since finding this datum is not easy, we also provide two methods to compute the period…
Let $T$ be a square matrix with a real spectrum, and let $f$ be an analytic function. The problem of the approximate calculation of $f(T)$ is discussed. Applying the Schur triangular decomposition and the reordering, one can assume that $T$…
A procedure to obtain differentiation matrices is extended straightforwardly to yield new differentiation matrices useful to obtain derivatives of complex rational functions. Such matrices can be used to obtain numerical solutions of some…
We show that the unitary factor U in the polar decomposition of a nonsingular matrix Z = U H is the minimizer for both ||Log(Q^* Z)|| and ||sym (Log(Q^*Z))|| over unitary Q, for any given invertible complex n-times-n matrix Z, for any…
One deals with degenerations by coordinate sections of the square generic Hankel matrix over a field $k$ of characteristic zero, along with its main related structures, such as the determinant of the matrix, the ideal generated by its…