Related papers: Viskovatov algorithm for Hermite-Pad\'e polynomial…
This paper presents an algorithm for solving multiobjective optimization problems involving composite functions, where we minimize a quadratic model that approximates $F(x) - F(x^k)$ and that can be derivative-free. We establish theoretical…
Let $\mathbf{f} = (f_1, \ldots, f_R)$ be a system of polynomials with integer coefficients in which the degrees need not all be the same. We provide sufficient conditions for which the system of equations $f_j (x_1, \ldots, x_n) = 0 \ (1…
We propose a method for interpolating divergence-free continuous magnetic fields via vector potential reconstruction using Hermite interpolation, which ensures high-order continuity for applications requiring adaptive, high-order ordinary…
We study higher-dimensional interlacing Fibonacci sequences, generated via both Chebyshev type functions and $m$-dimensional recurrence relations. For each integer $m$, there exist both rational and integer versions of these sequences,…
We give reconstruction algorithms for subclasses of depth-3 arithmetic circuits. In particular, we obtain the first efficient algorithm for finding tensor rank, and an optimal tensor decomposition as a sum of rank-one tensors, when given…
Quantum signal processing is a powerful framework in quantum algorithms, playing a central role in Hamiltonian simulation and related applications. The sequence of polynomials implemented at each step of this protocol provides a polynomial…
We consider two basic algorithmic problems concerning tuples of (skew-)symmetric matrices. The first problem asks to decide, given two tuples of (skew-)symmetric matrices $(B_1, \dots, B_m)$ and $(C_1, \dots, C_m)$, whether there exists an…
We obtain extensions of the Poincar\'e and Perron theorems for higher order recurrence relations and apply them to obtain an inverse type theorem for row sequences of (type II) Hermite-Pad\'e approximation of a vector of formal power…
Zeilberger's algorithm provides a method to compute recurrence and differential equations from given hypergeometric series representations, and an adaption of Almquist and Zeilberger computes recurrence and differential equations for…
Present day electromagnetic field calculations have limitations that are due to techniques employing edge-based discretization methods. While these vector finite element methods solve the issues of tangential continuity of fields and the…
We prove new linear independence results for the values of generalized hypergeometric functions ${}_pF_q$ at several distinct algebraic points, over arbitrary algebraic number fields. Our approach combines constructions of type II Pad\'{e}…
A unified theory of orthogonal polynomials of a discrete variable is presented through the eigenvalue problem of hermitian matrices of finite or infinite dimensions. It can be considered as a matrix version of exactly solvable Schr\"odinger…
We introduce a bi-graded polynomial that encodes the cohomology groups of the wild Hitchin system of type~$A_{n-1}$, constructed using an irregular singularity (determined by an integer~$m$) and an arbitrary regular singularity~$f$. When…
In this work, we investigate the sequence of monic q-Hermite I-Sobolev type orthogonal polynomials of higher-order, denoted as $\{\mathbb{H}_{n}(x;q)\}_{n\geq 0}$, which are orthogonal with respect to the following non-standard inner…
We present a new closed form for the interpolating polynomial of the general univariate Hermite interpolation that requires only calculation of polynomial derivatives, instead of derivatives of rational functions. This result is used to…
We present an algorithm to compute the primary decomposition of a submodule $\mathcal{N}$ of the free module $\Z[x_1, \ldots, x_n]^m$. For this purpose we use algorithms for primary decomposition of ideals in the polynomial ring over the…
We study multipoint Pad\'e approximants of type $(n,n)$ for the Hurwitz zeta function $f(a)=\zeta(s,a)$ with $\Re s>1$, constructed at quantile nodes $a_{n,j}=n\alpha_{n,j}$ generated by a real-analytic density $\kappa$ on…
We prove that for |x|,|t|<1, -1 <q \leq1 and n\geq0: \Sigma_{i\geq0}((t^{i})/((q)_{i}))h_{n+i}(x|q) = h_{n}(x|t,q) \Sigma_{i\geq0}((t^{i})/((q)_{i}))h_{i}(x|q), where h_{n}(x|q) and h_{n}(x|t,q) are respectively the so called q-Hermite and…
Lanczos methods for solving $\textit{A}\textbf{x}=\textbf{b}$ consist in constructing a sequence of vectors $(\textbf{x}_k), k=1,...$ such that $\textbf{r}_{k}=\textbf{b}-\textit{A}\textbf{x}_{k}=\textit{P}_{k}(\textit{A})\textbf{r}_{0}$,,…
We consider the computation of two normal forms for matrices over the univariate polynomials: the Popov form and the Hermite form. For matrices which are square and nonsingular, deterministic algorithms with satisfactory cost bounds are…