Related papers: ADI iteration for Lyapunov equations: a tangential…
The basic idea of an extremely fast convergent iterative method, the Forth-and-Back Implicit Lambda Iteration (FBILI), is briefly described and the applications of the method to various RT problems are listed and discussed.
The tangent dynamics of the Lyapunov modes and their dynamics as generated numerically - {\it the numerical dynamics} - is considered. We present a new phenomenological description of the numerical dynamical structure that accurately…
We develop two new variants of alternating direction methods of multipliers (ADMM) and two parallel primal-dual decomposition algorithms to solve a wide range class of constrained convex optimization problems. Our approach relies on a novel…
The correct computation of orbits of discrete dynamical systems on the interval is considered. Therefore, an arbitrary-precision floating-point approach based on automatic error analysis is chosen and a general algorithm is presented. The…
Performing Bayesian inference on large spatio-temporal models requires extracting inverse elements of large sparse precision matrices for marginal variances, as well as estimating model hyperparameters. Although direct matrix factorizations…
The numerical solution of large-scale Lyapunov matrix equations with symmetric banded data has so far received little attention in the rich literature on Lyapunov equations. We aim to contribute to this open problem by introducing two…
The classical Ka\v{c}anov scheme for the solution of nonlinear variational problems can be interpreted as a fixed point iteration method that updates a given approximation by solving a linear problem in each step. Based on this observation,…
We introduce a method which enables a recurrent dynamics model to be temporally abstract. Our approach, which we call Adaptive Skip Intervals (ASI), is based on the observation that in many sequential prediction tasks, the exact time at…
The Asymptotic Iteration Method (AIM) is a technique for solving analytically and approximately the linear second-order differential equation, especially the eigenvalue problems that frequently appear in theoretical and mathematical…
We derive a formula that simplifies the original asymptotic iteration method formulation to find the energy eigenvalues for the analytically solvable cases. We then show that there is a connection between the asymptotic iteration and the…
Recently, a nonlinear Poisson equation has been introduced to model nonlinear and nonlocal hyperpolarization effects in electrostatic solute-solvent interaction for biomolecular solvation analysis. Due to a strong nonlinearity associated…
We present the Deep Picard Iteration (DPI) method, a new deep learning approach for solving high-dimensional partial differential equations (PDEs). The core innovation of DPI lies in its use of Picard iteration to reformulate the typically…
In this work, we study two-dimensional diffusion-wave equations with variable exponent, modeling mechanical diffusive wave propagation in viscoelastic media with spatially varying properties. We first transform the diffusion-wave model into…
We introduce differentiable indirection -- a novel learned primitive that employs differentiable multi-scale lookup tables as an effective substitute for traditional compute and data operations across the graphics pipeline. We demonstrate…
In this paper a new restarting method for Krylov subspace matrix exponential evaluations is proposed. Since our restarting technique essentially employs the residual, some convergence results for the residual are given. We also discuss how…
The multi-variable affine index polynomial was defined by the author in previous work. The aim of this short note is to update the definition so it is generalizable to virtual tangles and to show it is compatible with tangle decomposition.…
This paper presents a new algorithm KIOPS for computing linear combinations of $\varphi$-functions that appear in exponential integrators. This algorithm is suitable for large-scale problems in computational physics where little or no…
We propose a new relative-error inexact version of the alternating direction method of multipliers (ADMM) for convex optimization. We prove the asymptotic convergence of our main algorithm as well as pointwise and ergodic…
We present here a new splitting method to solve Lyapunov equations of the type $AP + PA^T=-BB^T$ in a Kronecker product form. Although that resulting matrix is of order $n^2$, each iteration of the method demands only two operations with…
By means of a linear scaling of the variables we convert a singular bifurcation equation in $\R^n$ into an equivalent equation to which the classical implicit function theorem can be directly applied. This allows to deduce the existence of…