Isospectral schemes for exact algorithms for stochastic processes and PDE’s
摘要
Many problems in life sciences require the understanding of averages of a time-continuum evolution of a random phenomena or solutions to a parabolic equation defined on a continuum. The nature of both continua necessitates to implement approximation schemes in both time-space variables. We propose a procedure that identifies exactly these averages as averages of random phenomena defined on a discrete set whose time evolution is perfectly randomized. Consequently, both time and state sets are amenable to precise computation, effectively reducing the overall error to that inherent in simple Monte Carlo or integration methods. Our approach also extends to parabolic problems and enables to identify exact discretization schemes for solving the associated PDEs. These identifications are based on isospectral schemes called gateway relations as they relate two seemingly disconnected worlds, the continuum and the lattice, and their recent refinement, the interweaving relations. We show promising numerical examples leveraging the gateway relationship between the semigroups of squared Bessel and CIR processes, in