BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Simon J. A. Malham

Simon J. A. Malham contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math-phmath.APmath.MPmath.NAmath.PRnlin.SIq-fin.PRmath.COmath.DSmath.SGmath.SPmath.STq-fin.CPquant-phStatistics Theory

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
15topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access