We consider perturbation bounds and condition numbers for a complex indefinite linear algebraic system, which is of interest in science and engineering. Some existing results are improved, and illustrative numerical examples are provided.

High-dimensional two-sided space fractional diffusion equations with variable diffusion coefficients are discussed. The problems can be solved by an implicit finite difference scheme that is proven to be uniquely solvable, unconditionally stable and first-order convergent in the infinity norm. A nonsingular multilevel circulant pre-conditoner is proposed to accelerate the convergence rate of the Krylov subspace linear system solver efficiently. The preconditoned matrix for fast convergence is a sum of the identity matrix, a matrix with small norm, and a matrix with low rank under certain conditions. Moreover, the preconditioner is practical, with an O(N logN) operation cost and O(N) memory requirement. Illustrative numerical examples are also presented.

A multigrid method is proposed to compute the ground state solution of Bose-Einstein condensations by the finite element method based on the multilevel correction for eigenvalue problems and the multigrid method for linear boundary value problems. In this scheme, obtaining the optimal approximation for the ground state solution of Bose-Einstein condensates includes a sequence of solutions of the linear boundary value problems by the multigrid method on the multilevel meshes and some solutions of nonlinear eigenvalue problems some very low dimensional finite element space. The total computational work of this scheme can reach almost the same optimal order as solving the corresponding linear boundary value problem. Therefore, this type of multigrid scheme can improve the overall efficiency for the simulation of Bose-Einstein condensations. Some numerical experiments are provided to validate the efficiency of the proposed method.

We present a quasi-linear iterative method for solving a system of $m$ -nonlinear coupled differential equations. We provide an error analysis of the method to study its convergence criteria. In order to show the efficiency of the method, we consider some computational examples of this class of problem. These examples validate the accuracy of the method and show that it gives results which are convergent to the exact solutions. We prove that the method is accurate, fast and has a reasonable rate of convergence by computing some local and global error indicators.

The h-p version of the continuous Petrov-Galerkin time steppingmethod is analyzed for nonlinear initial value problems. AnL∞-error bound explicit with respect to the localdiscretization and regularity parameters is derived. Numerical examples areprovided to illustrate the theoretical results.

For any given matrix A ∈ℂnxn, a preconditioner tU(A) called thesuperoptimal preconditioner was proposed in 1992 by Tyrtyshnikov. It has beenshown that tU(A) is an efficient preconditioner forsolving various structured systems, for instance, Toeplitz-like systems. In thispaper, we construct the superoptimal preconditioners for different functions ofmatrices. Let f be a function of matrices from ℂnxn to ℂnxn. For any A ∈ ℂ nxn, one may construct two superoptimal preconditioners for f(A):tU(f(A)) and f(tU(A)).We establish basic properties of tU(f(A)) andf(tU(A)) for different functions ofmatrices. Some numerical tests demonstrate that the proposed preconditioners arevery efficient for solving the system f(A)x = b.

Conjugate symplectic eigenvalue problems arise in solving discrete linear-quadratic optimal control problems and discrete algebraic Riccati equations. In this article, backward errors of approximate pairs of conjugate symplectic matrices are obtained from their properties. Several numerical examples are given to illustrate the results.

In this paper, a multilevel correction scheme is proposed to solve the Steklov eigenvalue problem by nonconforming finite element methods. With this new scheme, the accuracy of eigenpair approximations can be improved after each correction step which only needs to solve a source problem on finer finite element space and an Steklov eigenvalue problem on the coarsest finite element space. This correction scheme can increase the overall efficiency of solving eigenvalue problems by the nonconforming finite element method. Furthermore, as same as the direct eigenvalue solving by nonconforming finite element methods, this multilevel correction method can also produce the lower-bound approximations of the eigenvalues.

This paper illuminates the derivation, applicability, and efficiency of the multiplicative Runge–Kutta method, derived in the framework of geometric multiplicative calculus. The removal of the restrictions of geometric multiplicative calculus on positive-valued functions of real variables and the fact that the multiplicative derivative does not exist at the roots of the function are presented explicitly to ensure that the proposed method is universally applicable. The error and stability analyses are also carried out explicitly in the framework of geometric multiplicative calculus. The method presented is applied to various problems and the results are compared to those obtained from the ordinary Runge–Kutta method. Moreover, for one example, a comparison of the computation time against relative error is worked out to illustrate the general advantage of the proposed method.

With the aim of solving in a four dimensional phase space a multi-scale Vlasov-Poisson system, we propose in a Particle-In-Cell framework a robust time-stepping method that works uniformly when the small parameter vanishes. As an exponential integrator, the scheme is able to use large time steps with respect to the typical size of the solution’s fast oscillations. In addition, we show numerically that the method has accurate long time behaviour and that it is asymptotic preserving with respect to the limiting Guiding Center system.

We propose a new Adomian decomposition method (ADM) using an integrating factor for the Emden–Fowler equation. With this method, we are able to solve certain Emden–Fowler equations for which the traditional ADM fails. Numerical results obtained from testing our linear and nonlinear models are far more reliable and efficient than those from existing methods. We also present a complete error analysis and a convergence criterion for this method. One drawback of the traditional ADM is that the interval of convergence of the Adomian truncated series is very small. Some techniques, such as Pade approximants, can enlarge this interval, but they are too complicated. Here, we use a continuation technique to extend our method to a larger interval.

We study the radius of absolute monotonicity $R$ of rational functions with numerator and denominator of degree $s$ that approximate the exponential function to order $p$ . Such functions arise in the application of implicit $s$ -stage, order $p$ Runge–Kutta methods for initial value problems, and the radius of absolute monotonicity governs the numerical preservation of properties like positivity and maximum-norm contractivity. We construct a function with $p=2$ and $R>2s$ , disproving a conjecture of van de Griend and Kraaijevanger. We determine the maximum attainable radius for functions in several one-parameter families of rational functions. Moreover, we prove earlier conjectured optimal radii in some families with two or three parameters via uniqueness arguments for systems of polynomial inequalities. Our results also prove the optimality of some strong stability preserving implicit and singly diagonally implicit Runge–Kutta methods. Whereas previous results in this area were primarily numerical, we give all constants as exact algebraic numbers.

In this study, Liénard equations in their general form are treated using the Adomian decomposition method. The special structure of the Liénard equation is exploited to obtain a numerically efficient algorithm suitable for solution by a computer program.

In this paper we present computer-assisted proofs of a number of results in theoretical fluid dynamics and in quantum mechanics. An algorithm based on interval arithmetic yields provably correct eigenvalue enclosures and exclosures for non-self-adjoint boundary eigenvalue problems, the eigenvalues of which are highly sensitive to perturbations. We apply the algorithm to: the Orr–Sommerfeld equation with Poiseuille profile to prove the existence of an eigenvalue in the classically unstable region for Reynolds number R=5772.221818; the Orr–Sommerfeld equation with Couette profile to prove upper bounds for the imaginary parts of all eigenvalues for fixed R and wave number α; the problem of natural oscillations of an incompressible inviscid fluid in the neighbourhood of an elliptical flow to obtain information about the unstable part of the spectrum off the imaginary axis; Squire’s problem from hydrodynamics; and resonances of one-dimensional Schrödinger operators.

Implicit Runge–Kutta methods have a special role in the numerical solution of stiff problems, such as those found by applying the method of lines to the partial differential equations arising in physical modelling. Of particular interest in this paper are the high-order methods based on Gaussian quadrature and the efficiently implementable singly implicit methods.

In this paper we look at the performance of trigonometric integrators applied to highly oscillatory differential equations. It is widely known that some of the trigonometric integrators suffer from low-order resonances for particular step sizes. We show here that, in general, trigonometric integrators also suffer from higher-order resonances which can lead to loss of nonlinear stability. We illustrate this with the Fermi–Pasta–Ulam problem, a highly oscillatory Hamiltonian system. We also show that in some cases trigonometric integrators preserve invariant or adiabatic quantities but at the wrong values. We use statistical properties such as time averages to further evaluate the performance of the trigonometric methods and compare the performance with that of the mid-point rule.

An exchangeable sequence of random variables is constructed with all finite-dimensional distribution functions having an Archimedean copula (as defined by Schweizer and Sklar (1983)). Through a monotone transformation of this exchangeable sequence, we obtain and characterize the class of exchangeable sequences possessing the max-stable property as defined by De Haan and Rachev (1989). Several parametric examples are given.

For Fredholm equations of the first kind with continuous kernels we investigate the uniform convergence of a general class of regularization methods. Applications are made to Tikhonov regularization and Landweber's iteration method.