Source author record

Simon J. A. Malham

Simon J. A. Malham appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.PR math.NA math-ph math.AP math.CO math.MP nlin.SI q-fin.PR math.DS math.OC math.SG math.SP math.ST q-fin.CP quant-ph Statistics Theory

Catalog footprint

What is connected

11works

16topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Applications of Grassmannian flows to integrable systems

We show how many classes of partial differential systems with local and nonlocal nonlinearities are linearisable in the sense that they are realisable as Fredholm Grassmannian flows. In other words, time-evolutionary solutions to such systems can be constructed from solutions to the corresponding underlying linear partial differential system, by solving a linear Fredholm equation. For example, it is well-known that solutions to classical integrable partial differential systems can be generated by solving a corresponding linear partial differential system for the scattering data and then solving the linear Fredholm (or Volterra) integral equation known as the Gel'fand-Levitan-Marchenko equation. In this paper and in a companion paper, Doikou et al. [DMSW:graphflows], we both, survey the classes of nonlinear systems that are realisable as Fredholm Grassmannian flows, and present new example applications of such flows. We also demonstrate the usefulness of such a representation. Herein we extend the work of Poppe and demonstrate how solution flows of the non-commutative potential Korteweg de Vries and nonlinear Schrodinger systems are examples of such Grassmannian flows. In the companion paper we use this Grassmannian flow approach as well as an extension to nonlinear graph flows, to solve Smoluchowski coagulation and related equations.

preprint2021arXiv

Series expansions and direct inversion for the Heston model

Efficient sampling for the conditional time integrated variance process in the Heston stochastic volatility model is key to the simulation of the stock price based on its exact distribution. We construct a new series expansion for this integral in terms of double infinite weighted sums of particular independent random variables through a change of measure and the decomposition of squared Bessel bridges. When approximated by series truncations, this representation has exponentially decaying truncation errors. We propose feasible strategies to largely reduce the implementation of the new series to simulations of simple random variables that are independent of any model parameters. We further develop direct inversion algorithms to generate samples for such random variables based on Chebyshev polynomial approximations for their inverse distribution functions. These approximations can be used under any market conditions. Thus, we establish a strong, efficient and almost exact sampling scheme for the Heston model.

preprint2021arXiv

The non-commutative Korteweg--de Vries hierarchy and combinatorial Poppe algebra

We give a constructive proof, to all orders, that each member of the non-commutative potential Korteweg-de Vries hierarchy is a Fredholm Grassmannian flow and is therefore linearisable. Indeed we prove this for any linear combination of fields from this hierarchy. That each member of the hierarchy is linearisable, and integrable in this sense, means that the time evolving solution can be generated from the solution to the corresponding linear dispersion equation in the hierarchy, combined with solving an associated linear Fredholm equation representing the Marchenko equation. Further, we show that within the class of polynomial partial differential fields, at every order, each member of the non-commutative potential Korteweg--de Vries hierarchy is unique. Indeed, we prove to all orders, that each such member matches the non-commutative Lax hierarchy field, which is therefore a polynomial partial differential field. We achieve this by constructing the abstract combinatorial algebra that underlies the non-commutative potential Korteweg-de Vries hierarchy. This algebra is the non-commutative polynomial algebra over the real line generated by the set of all compositions endowed with the Poppe product. This product is the abstract representation of the product rule for Hankel operators pioneered by Ch. Poppe for integrable equations such as the Sine-Gordon and Korteweg-de Vries equations. Integrability of the hierarchy members translates, in the combinatorial algebra, to proving the existence of a `Poppe polynomial' expansion for basic compositions in terms of `linear signature expansions'. Proving the existence of such Poppe polynomial expansions boils down to solving a linear algebraic problem for the expansion coefficients, which we solve constructively to all orders.

preprint2015arXiv

Flows and stochastic Taylor series in Ito calculus

For stochastic systems driven by continuous semimartingales an explicit formula for the logarithm of the Ito flow map is given. A similar formula is also obtained for solutions of linear matrix-valued SDEs driven by arbitrary semimartingales. The computation relies on the lift to quasi-shuffle algebras of formulas involving products of Ito integrals of semimartingales. Whereas the Chen-Strichartz formula computing the logarithm of the Stratonovich flow map is classically expanded as a formal sum indexed by permutations, the analogous formula in Ito calculus is naturally indexed by surjections. This reflects the change of algebraic background involved in the transition between the two integration theories.

preprint2014arXiv

Efficient almost-exact Levy area sampling

We present a new method for sampling the Levy area for a two-dimensional Wiener process conditioned on its endpoints. An efficient sampler for the Levy area is required to implement a strong Milstein numerical scheme to approximate the solution of a stochastic differential equation driven by a two-dimensional Wiener process whose diffusion vector fields do not commute. Our method is simple and complementary to those of Gaines-Lyons and Wiktorsson, and amenable to quasi-Monte--Carlo implementation. It is based on representing the Levy area by an infinite weighted sum of independent Logistic random variables. We use Chebychev polynomials to approximate the inverse distribution function of sums of independent Logistic random variables in three characteristic regimes. The error is controlled by the degree of the polynomials, we set the error to be uniformly 10^(-12). We thus establish a strong almost-exact Levy area sampling method. The complexity of our method is square logarithmic. We indicate how our method can contribute to efficient sampling in higher dimensions.

preprint2013arXiv

Computing the Maslov index for large systems

We address the problem of computing the Maslov index for large linear symplectic systems on the real line. The Maslov index measures the signed intersections (with a given reference plane) of a path of Lagrangian planes. The natural chart parameterization for the Grassmannian of Lagrangian planes is the space of real symmetric matrices. Linear system evolution induces a Riccati evolution in the chart. For large order systems this is a practical approach as the computational complexity is quadratic in the order. The Riccati solutions, however, also exhibit singularites (which are traversed by changing charts). Our new results involve characterizing these Riccati singularities and two trace formulae for the Maslov index as follows. First, we show that the number of singular eigenvalues of the symmetric chart representation equals the dimension of intersection with the reference plane. Second, the Cayley map is a diffeomorphism from the space of real symmetric matrices to the manifold of unitary symmetric matrices. We show the logarithm of the Cayley map equals the arctan map (modulo 2i) and its trace measures the angle of the Langrangian plane to the reference plane. Third, the Riccati flow under the Cayley map induces a flow in the manifold of unitary symmetric matrices. Using the natural unitary action on this manifold, we pullback the flow to the unitary Lie algebra and monitor its trace. This avoids singularities, and is a natural robust procedure. We demonstrate the effectiveness of these approaches by applying them to a large eigenvalue problem. We also discuss the extension of the Maslov index to the infinite dimensional case.

preprint2013arXiv

Levy Processes and Quasi-Shuffle Algebras

We investigate the algebra of repeated integrals of semimartingales. We prove that a minimal family of semimartingales generates a quasi-shuffle algebra. In essence, to fulfill the minimality criterion, first, the family must be a minimal generator of the algebra of repeated integrals generated by its elements and by quadratic covariation processes recursively constructed from the elements of the family. Second, recursively constructed quadratic covariation processes may lie in the linear span of previously constructed ones and of the family, but may not lie in the linear span of repeated integrals of these. We prove that a finite family of independent Levy processes that have finite moments generates a minimal family. Key to the proof are the Teugels martingales and a strong orthogonalization of them. We conclude that a finite family of independent Levy processes form a quasi-shuffle algebra. We discuss important potential applications to constructing efficient numerical methods for the strong approximation of stochastic differential equations driven by Levy processes.

preprint2012arXiv

Algebraic structure of stochastic expansions and efficient simulation

We investigate the algebraic structure underlying the stochastic Taylor solution expansion for stochastic differential systems.Our motivation is to construct efficient integrators. These are approximations that generate strong numerical integration schemes that are more accurate than the corresponding stochastic Taylor approximation, independent of the governing vector fields and to all orders. The sinhlog integrator introduced by Malham & Wiese (2009) is one example. Herein we: show that the natural context to study stochastic integrators and their properties is the convolution shuffle algebra of endomorphisms; establish a new whole class of efficient integrators; and then prove that, within this class, the sinhlog integrator generates the optimal efficient stochastic integrator at all orders.

preprint2012arXiv

Chi-square simulation of the CIR process and the Heston model

The transition probability of a Cox-Ingersoll-Ross process can be represented by a non-central chi-square density. First we prove a new representation for the central chi-square density based on sums of powers of generalized Gaussian random variables. Second we prove Marsaglia's polar method extends to this distribution, providing a simple, exact, robust and efficient acceptance-rejection method for generalized Gaussian sampling and thus central chi-square sampling. Third we derive a simple, high-accuracy, robust and efficient direct inversion method for generalized Gaussian sampling based on the Beasley-Springer-Moro method. Indeed the accuracy of the approximation to the inverse cumulative distribution function is to the tenth decimal place. We then apply our methods to non-central chi-square variance sampling in the Heston model. We focus on the case when the number of degrees of freedom is small and the zero boundary is attracting and attainable, typical in foreign exchange markets. Using the additivity property of the chi-square distribution, our methods apply in all parameter regimes.

preprint2010arXiv

An introduction to SDE simulation

We outline the basic ideas and techniques underpinning the simulation of stochastic differential equations. In particular we focus on strong simulation and its context. We also provide illustratory examples and sample matlab algorithms for the reader to use and follow. Our target audience is advanced undergraduate and graduate students interested in learning about simulating stochastic differential equations. We try to address the FAQs we have encountered.

preprint2009arXiv

Stochastic expansions and Hopf algebras

We study solutions to nonlinear stochastic differential systems driven by a multi-dimensional Wiener process. A useful algorithm for strongly simulating such stochastic systems is the Castell--Gaines method, which is based on the exponential Lie series. When the diffusion vector fields commute, it has been proved that at low orders this method is more accurate in the mean-square error than corresponding stochastic Taylor methods. However it has also been shown that when the diffusion vector fields do not commute, this is not true for strong order one methods. Here we prove that when there is no drift, and the diffusion vector fields do not commute, the exponential Lie series is usurped by the sinh-log series. In other words, the mean-square error associated with a numerical method based on the sinh-log series, is always smaller than the corresponding stochastic Taylor error, in fact to all orders. Our proof utilizes the underlying Hopf algebra structure of these series, and a two-alphabet associative algebra of shuffle and concatenation operations. We illustrate the benefits of the proposed series in numerical studies.