On the Leverrier-Faddeev algorithm for computing the Moore-Penrose inverse

We present an improvement in the implementation of the Leverrier-Faddeev algorithm for symbolic computation of the Moore-Penrose inverse of one-variable polynomial matrices, introduced in Linear Algebra Appl. 252, 35–60 (1997). Complexity analysis of the original algorithm and its improvement is presented. Algorithm and its improvement are implemented and compared in the symbolic computational package MATHEMATICA. We compare CPU time required for computation of some test matrices by means of the original algorithm and its improvement.     

The convergence of a new parallel algorithm for the Navier–Stokes equations

This paper deals with the convergence and stability of a new parallel algorithm and the error estimates for a particular case of the new parallel algorithm, which is used to solve the incompressible nonstationary Navier–Stokes equations. The theoretical results show that the scheme is (at least) conditionally stable and convergent.     

New results on condensed Cramer’s rule for the general solution to some restricted quaternion matrix equations

In this paper, we derive some condensed Cramer’s rules for the general solution, the least squares solution and the least norm solution to some restricted quaternion matrix equations, respectively. The findings of this paper extend some known results in the literature.     

Polynomial solution of the most general linear Fredholm integrodifferential–difference equations by means of Taylor matrix method

In this article, a Taylor matrix method is developed to find an approximate solution of the most general linear Fredholm integrodifferential–difference equations with variable coefficients under the mixed conditions in terms of Taylor polynomials. Also numerical examples are presented, which illustrate the pertinent features of the method. In some numerical examples, MAPLE modules are designed for the purpose of testing and using the method.     

Partitioned time discretization for parallel solution of coupled ODE systems

One decoupling method for multiphysics, multiscale, multidomain applications involves partitioning the problem via explicit time discretizations in the coupling terms. Specialized, problem-specific techniques are needed for the resulting partitioned methods to avoid time step restrictions which make long time calculations costly. This report studies unconditionally stable, uncoupled time stepping methods for a model problem sharing mathematical structure akin to the coupled atmosphere-ocean system. Three decoupled time stepping algorithms are introduced and their stability and consistency are rigorously examined. Numerical experiments are performed that study their stability and convergence properties.     

A matrix-less and parallel interpolation–extrapolation algorithm for computing the eigenvalues of preconditioned banded symmetric Toeplitz matrices

Numerical Algorithms - In the past few years, Bogoya, Böttcher, Grudsky, and Maximenko obtained the precise asymptotic expansion for the eigenvalues of a Toeplitz matrix Tn(f), under suitable...     

A new bound on Hrushovski’s algorithm for computing the Galois group of a linear differential equation

The complexity of computing the Galois group of a linear differential equation is of general interest. In a recent work, Feng gave the first degree bound on Hrushovski’s algorithm for computing the Galois group of a linear differential equation. This bound is the degree bound of the polynomials used in the first step of the algorithm for finding a proto-Galois group (see Definition 2.7) and is sextuply exponential in the order of the differential equation. In this paper, we use Szántó’s algorithm of triangular representation for algebraic sets to analyze the complexity of computing the Galois group of a linear differential equation and we give a new bound which is triple exponential in the order of the given differential equation.     

A method for the numerical solution of the Painlevé equations

A numerical method for solving the Cauchy problem for all the six Painlevé equations is proposed. The difficulty of solving these equations is that the unknown functions can have movable (that is, dependent on the initial data) singular points of the pole type. Moreover, the Painlevé III–VI equations may have singularities at points where the solution takes certain finite values. The positions of all these singularities are not a priori known and are determined in the process of solving the equation. The proposed method is based on the transition to auxiliary systems of differential equations in neighborhoods of the indicated points. The equations in these systems and their solutions have no singularities at the corresponding point and its neighborhood. Such auxiliary equations are derived for all Painlevé equations and for all types of singularities. Efficient criteria for transition to auxiliary systems are formulated, and numerical results illustrating the potentials of the method are presented.     

An interior-point path-following algorithm for computing a Leontief economy equilibrium

In this paper, we present an interior-point path-following algorithm for computing a Leontief economy equilibrium, that is, an exchange market equilibrium with Leontief utility functions, which is known to be in the complexity class of PPAD-complete. It is known that an equilibrium corresponds to a solution of a system of complementarities, so we construct a smooth homotopy interior-point path to tackle this system. We prove that there always exists a continuously differentiable path leading to a complementary solution of the nonlinear system and at the same time to a Leontief economy equilibrium associated with the solution. We also report preliminary computational results to show effectiveness of the path-following Newton method.     

On the solution of matrix‐valued linear stochastic differential equations driven by semimartingales

Given any semimartingale‐driven matrix‐valued linear stochastic differential equation, it is shown that the underlying homogeneous equation has a solution with a semigroup property. Under some commutativity assumptions it is the martix analogue of Doléans‐Dade's exponential. Some variation-of-constants formulas are given.     

Representation for the W-weighted Drazin inverse of linear operators

In this paper we study the W-weighted Drazin inverse of the bounded linear operators between Banach spaces and its representation theorem. Based on this representation, utilizing the spectral theory of Banach space operators, we derive an approximating expression of the W-weighted Drazin inverse and an error bound. Also, a perturbation theorem for the W-weighted Drazin inverse is uniformly obtained from the representation theorem.     

Hölder estimate for the solution of degenerate parabolic equations

Consider the degenerate parabolic equations of the type $$u_t = div A(x,t,u,Du) + b(x,t,u,Du)$$ which is of the same nature as $$u_t = div|Du|^p Du + |Du|^{p + 2} (p > 2).$$ This paper is to study the \(C^{1 + \alpha } ,\frac{{1 + \alpha }}{2}\) Hölder continuity of a class of degenerate parabolic equations and the existence and uniqueness of the initial boundary value problem.     

Using dual techniques to derive componentwise and mixed condition numbers for a linear function of a linear least squares solution

We prove duality results for adjoint operators and product norms in the framework of Euclidean spaces. We show how these results can be used to derive condition numbers especially when perturbations on data are measured componentwise relatively to the original data. We apply this technique to obtain formulas for componentwise and mixed condition numbers for a linear function of a linear least squares solution. These expressions are closed when perturbations of the solution are measured using a componentwise norm or the infinity norm and we get an upper bound for the Euclidean norm.      

Existence of a periodic solution for continuous and discrete periodic second-order equations with variable potentials

Green’s functions for new second-order periodic differential and difference equations with variable potentials are found, then used as kernels in integral operators to guarantee the existence of a positive periodic solution to continuous and discrete second-order periodic boundary value problems with periodic coefficient functions. A new version of the Leggett-Williams fixed point theorem is employed.