Related papers: Accurate complex Jacobi rotations
We introduce a class of doubly infinite complex Jacobi matrices determined by a simple convergence condition imposed on the diagonal and off-diagonal sequences. For each Jacobi matrix belonging to this class, an analytic function, called a…
We give an analogy of Jacobi's formula, which relates the hypergeometric function with parameters $(1/4,1/4,1)$ and theta constants. By using this analogy and twice formulas of theta constants, we obtain a transformation formula for this…
Asymptotic approximations to the zeros of Jacobi polynomials are given, with methods to obtain the coefficients in the expansions. These approximations can be used as standalone methods for the non-iterative computation of the nodes of…
Landen transformation formulas, which connect Jacobi elliptic functions with different modulus parameters, were first obtained over two hundred years ago by changing integration variables in elliptic integrals.We rediscover known results as…
The complete elliptic integral of the first and second kind, K(k) and E(k), appear in a multitude of physics and engineering applications. Because there is no known closed-form, the exact values have to be computed numerically. Here,…
Given a structured matrix $A$ we study the problem of finding the closest normal matrix with the same structure. The structures of our interest are: Hamiltonian, skew-Hamiltonian, per-Hermitian, and perskew-Hermitian. We develop a…
First-order perturbative calculation of the frequency-shifts caused by special relativity is performed for a charged particle confined in a Penning trap. The perturbed motion is approximated by the Jacobian elliptic functions which describe…
An improvement on precision of recursive function simulation in IEEE floating point standard is presented. It is shown that the average of rounding towards negative infinite and rounding towards positive infinite yields a better result than…
We use precise asymptotic expansions for Jacobi functions $\phi^{(\alpha,\beta)}_\lambda$ parameters $\alpha$, $\beta$ satisfying $\alpha>1/2$, $\alpha>\beta>-1/2$, to generalizing classical H\"ormander-type multiplier theorem for the…
We present a fast sweeping method for a class of Hamilton-Jacobi equations that arise from time-independent problems in optimal control theory. The basic method in two dimensions uses a four point stencil and is extremely simple to…
Elliptic functions are largely studied and standardized mathematical objects. The two usual approaches are due to Jacobi and Weierstrass. From a contour integral which allowed us to unify many summation formulae (Euler-MacLaurin, Poisson,…
Finite-precision floating point arithmetic unavoidably introduces rounding errors which are traditionally bounded using a worst-case analysis. However, worst-case analysis might be overly conservative because worst-case errors can be…
This work proposes a higher-order iterative framework for solving matrix equations, inspired by the structure and functionality of neural networks. A modification of the classical Jacobi iterative method is introduced to compute…
Asymptotic rates of convergence of optimal SOR applied to linear systems with consistently ordered 2-cyclic matrices have been extensively studied in the case where the Jacobi eigenvalues are are real and contained in an interval centered…
We describe two main classes of one-sided trigonometric and hyperbolic Jacobi-type algorithms for computing eigenvalues and eigenvectors of Hermitian matrices. These types of algorithms exhibit significant advantages over many other…
Generating 2-by-2 unitary matrices in floating-precision arithmetic is a delicate task. One way to reduce the accumulation error is to use less floating-point operations to compute each of the entries in the 2-by-2 unitary matrix. This…
In this document, some elements of the theory and algorithmics corresponding to the existence and computability of approximate joint eigenpairs for finite collections of matrices with applications to model order reduction, are presented.…
Guided by the Jacobi's work published the year before his death about the rotation of a rigid body, the behavior of the rotation matrix describing the dynamics of the free rigid body is studied. To illustrate this dynamics one draws on a…
In this note, we propose a symplectic algorithm for the stable manifolds of the Hamilton-Jacobi equations combined with an iterative procedure in [Sakamoto-van~der Schaft, IEEE Transactions on Automatic Control, 2008]. Our algorithm…
In this paper, we introduce and analyze arbitrarily high-order quadrature rules for evaluating the two-dimensional singular integrals of the forms \begin{align} I_{i,j} = \int_{\mathbb{R}^2}\phi(x)\frac{x_ix_j}{|x|^{2+\alpha}} \d x, \quad…