In this paper, preconditioned iterative methods for solving two-dimensional space-fractional diffusion equations are considered. The fractional diffusion equation is discretized by a second-order finite difference scheme, namely, the Crank-Nicolson weighted and shifted Grünwald difference (CN-WSGD) scheme proposed in [W. Tian, H. Zhou and W. Deng, A class of second order difference approximation for solving space fractional diffusion equations, Math. Comp., 84 (2015) 1703-1727]. For the discretized linear systems, we first propose preconditioned iterative methods to solve them. Then we apply the D’Yakonov ADI scheme to split the linear systems and solve the obtained splitting systems by iterative methods. Two preconditioned iterative methods, the preconditioned generalized minimal residual (preconditioned GMRES) method and the preconditioned conjugate gradient normal residual (preconditioned CGNR) method, are proposed to solve relevant linear systems. By fully exploiting the structure of the coefficient matrix, we design two special kinds of preconditioners, which are easily constructed and are able to accelerate convergence of iterative solvers. Numerical results show the efficiency of our preconditioners.

Scaling of proposals for Metropolis algorithms is an important practical problem in Markov chain Monte Carlo implementation. Analyses of the random walk Metropolis for high-dimensional targets with specific functional forms have shown that in many cases the optimal scaling is achieved when the acceptance rate is approximately 0.234, but that there are exceptions. We present a general set of sufficient conditions which are invariant to orthonormal transformation of the coordinate axes and which ensure that the limiting optimal acceptance rate is 0.234. The criteria are shown to hold for the joint distribution of successive elements of a stationary pth-order multivariate Markov process.

In this paper we consider the splitting method first introduced in rare event analysis. In this technique, the sample paths are split into R multiple copies at various stages to speed up the simulation. Given the cost, the optimization of the algorithm suggests taking all the transition probabilities to be equal; nevertheless, in practice, these quantities are unknown. We address this problem by presenting an algorithm in two phases.

This paper proposes and analyzes discrete-time approximations to a class of diffusions, with an emphasis on preserving certain important features of the continuous-time processes in the approximations. We start with multivariate diffusions having three features in particular: they are martingales, each of their components evolves within the unit interval, and the components are almost surely ordered. In the models of the term structure of interest rates that motivate our investigation, these properties have the important implications that the model is arbitrage-free and that interest rates remain positive. In practice, numerical work with such models often requires Monte Carlo simulation and thus entails replacing the original continuous-time model with a discrete-time approximation. It is desirable that the approximating processes preserve the three features of the original model just noted, though standard discretization methods do not. We introduce new discretizations based on first applying nonlinear transformations from the unit interval to the real line (in particular, the inverse normal and inverse logit), then using an Euler discretization, and finally applying a small adjustment to the drift in the Euler scheme. We verify that these methods enforce important features in the discretization with no loss in the order of convergence (weak or strong). Numerical results suggest that these methods can also yield a better approximation to the law of the continuous-time process than does a more standard discretization.

We study the present value Z∞ = ∫0∞ e-Xt-dYt where (X,Y) is an integrable Lévy process. This random variable appears in various applications, and several examples are known where the distribution of Z∞ is calculated explicitly. Here sufficient conditions for Z∞ to exist are given, and the possibility of finding the distribution of Z∞ by Markov chain Monte Carlo simulation is investigated in detail. Then the same ideas are applied to the present value Z-∞ = ∫0∞ exp{-∫0tRsds}dYt where Y is an integrable Lévy process and R is an ergodic strong Markov process. Numerical examples are given in both cases to show the efficiency of the Monte Carlo methods.

Quasi-birth-death processes are commonly used Markov chain models in queueing theory, computer performance, teletraffic modeling and other areas. We provide a new, simple algorithm for the matrix-geometric rate matrix. We demonstrate that it has quadratic convergence. We show theoretically and through numerical examples that it converges very fast and provides extremely accurate results even for almost unstable models.

This paper presents a state space and time discretization for the general average impulse control of piecewise deterministic Markov processes (PDPs). By combining several previous results we show that under some continuity, boundedness and compactness conditions on the parameters of the process, boundedness of the discretizations, and compactness of the state space, the discretized problem will converge uniformly to the original one. An application to optimal capacity expansion under uncertainty is given.

Past work relating to the computation of time-dependent state probabilities in M/M/1 queueing systems is reviewed, with emphasis on methods that avoid Bessel functions. A new series formula of Sharma [13] is discussed and its connection with traditional Bessel function series is established. An alternative new series is developed which isolates the steady-state component for all values of traffic intensity and which turns out to be computationally superior. A brief comparison of our formula, Sharma's formula, and a classical Bessel function formula is given for the computation time of the probability that an initially empty system is empty at time t later.