Related papers: An extension of the Chudnovsky algorithm
The class of input-output systems representable as Chen-Fliess series arises often in control theory. One well known drawback of this representation, however, is that the iterated integrals which appear in these series are algebraically…
In this paper, we analyze several variants of a simple method for generating prime numbers with fewer random bits. To generate a prime $p$ less than $x$, the basic idea is to fix a constant $q\propto x^{1-\varepsilon}$, pick a uniformly…
We design a class of Chudnovsky-type algorithms multiplying k elements of a finite extension of order n a finite field K. We prove that these algorithms give a tensor decomposition of the k-multiplication for which the rank is linear in n…
We present a new recursive generation algorithm for prefix normal words. These are binary strings with the property that no substring has more 1s than the prefix of the same length. The new algorithm uses two operations on binary strings,…
The time evolution of Markovian open quantum systems is governed by Lindblad master equations, whose solution can be formally written as the Lindbladian exponential acting on the initial density matrix. By expanding this Lindbladian…
By solving a control problem and using Malliavin calculus, explicit derivative formula is derived for the semigroup $P_t$ generated by the Gruschin type operator on $\R^{m}\times \R^{d}:$ $$L (x,y)=\ff 1 2 \bigg\{\sum_{i=1}^m \pp_{x_i}^2…
Topological invariants of a dataset, such as the number of holes that survive from one length scale to another (persistent Betti numbers) can be used to analyze and classify data in machine learning applications. We present an improved…
We significantly enhance the simulation accuracy of initial Trotter circuits for Hamiltonian simulation of quantum systems by integrating first-order Riemannian optimization with tensor network methods. Unlike previous approaches, our…
R. L. Graham and H. O. Pollak observed that the sequence $$u_1=1,\qquad u_{n+1}=\lfloor \sqrt{2} (u_n+1/2)\rfloor, \quad n\geq 1,$$ has the curious property that the sequence of numbers $(u_{2n+1}-2u_{2n-1})_{n\geq 1}$ denotes the binary…
Various types of expansions in series of Chebyshev-Hermite polynomials currently used in astrophysics for weakly non-normal distributions are compared, namely the Gram-Charlier, Gauss-Hermite and Edgeworth expansions. It is shown that the…
We present a complexity reduction algorithm for a family of parameter-dependent linear systems when the system parameters belong to a compact semi-algebraic set. This algorithm potentially describes the underlying dynamical system with…
A single-index model (SIM) provides for parsimonious multi-dimensional nonlinear regression by combining parametric (linear) projection with univariate nonparametric (non-linear) regression models. We show that a particular Gaussian process…
Digital System Research has pioneered the mathematics and design for a new class of computing machine using residue numbers. Unlike prior art, the new breakthrough provides methods and apparatus for general purpose computation using several…
In this work we present a Gaussian process that arise from the iteration of p fractional Ornstein-Uhlenbeck processes generated by the same fractional Brownian motion. This iteration results, when the values of lambdas are pairwise…
Nowadays, parallel computing is ubiquitous in several application fields, both in engineering and science. The computations rely on the floating-point arithmetic specified by the IEEE754 Standard. In this context, an elementary brick of…
We have calculated on the computer the sum $\bar{\BB}_M$ of reciprocals of all 47 known Mersenne primes with the accuracy of over 12000000 decimal digits. Next we developed $\bar{\BB}_M$ into the continued fraction and calculated…
We propose a new yet natural algorithm for learning the graph structure of general discrete graphical models (a.k.a. Markov random fields) from samples. Our algorithm finds the neighborhood of a node by sequentially adding nodes that…
We revisit the problem of learning from untrusted batches introduced by Qiao and Valiant [QV17]. Recently, Jain and Orlitsky [JO19] gave a simple semidefinite programming approach based on the cut-norm that achieves essentially…
A general quantum circuit can be simulated classically in exponential time. If it has a planar layout, then a tensor-network contraction algorithm due to Markov and Shi has a runtime exponential in the square root of its size, or more…
It is known that computing the permanent of the matrix $1+A$, where $A$ is a finite-rank matrix, requires a number of operations polynomial in the matrix size. Motivated by the boson-sampling proposal of restricted quantum computation, I…