Convergence analysis and remedies for higher-mode smooth Brownian bridges
Establish the convergence behavior of SDC-SBB methods using additional terms in the Karhunen–Loève smooth Brownian bridge expansion, and determine effective remedies for the resulting high-frequency quadrature and substepping difficulties as the time step decreases.
References
Additional terms in the expansion not only require specialized quadrature for the computation of $I_S$, but pose restrictions on the number of collocation nodes due to substepping by an Euler method in the computation of the correction sweep, see eq.~eq:correctionsweep. As a rule of thumb, for an explicit Euler method, the stepsize, i.e., here the distance between two quadrature nodes, should be less than $1/4$-th of the oscillation period. For a term involving $\cos(j\pi t/\Delta T)$ with $t\in [0,\Delta T]$ this would result in $M>2j$ collocation nodes. Options to avoid the need for fine discretization include introducing terms to the expansion that reduce the variance error without high-frequency oscillations, or using higher-order substeppers. A thorough analysis of these convergence issues and potential remedies is out of scope for this paper, but will be a topic for future research.