A micro/macro parallel-in-time (parareal) algorithm applied to a climate model with discontinuous non-monotone coefficients and oscillatory forcing

We present the application of a micro/macro parareal algorithm for a 1-D energy balance climate model with discontinuous and non-monotone coefficients and forcing terms. The micro/macro parareal method uses a coarse propagator, based on a (macroscopic) 0-D approximation of the underlying (microscopic) 1-D model. We compare the performance of the method using different versions of the macro model, as well as different numerical schemes for the micro propagator, namely an explicit Euler method with constant stepsize and an adaptive library routine. We study convergence of the method and the theoretical gain in computational time in a realization on parallel processors. We show that, in this example and for all settings, the micro/macro parareal method converges in fewer iterations than the number of used parareal subintervals, and that a theoretical gain in performance of up to 10 is possible

Noise reduction in coarse bifurcation analysis of stochastic agent-based models: an example of consumer lock-in

We investigate coarse equilibrium states of a fine-scale, stochastic agent-based model of consumer lock-in in a duopolistic market. In the model, agents decide on their next purchase based on a combination of their personal preference and their neighbours' opinions. For agents with independent identically-distributed parameters and all-to-all coupling, we derive an analytic approximate coarse evolution-map for the expected average purchase. We then study the emergence of coarse fronts when spatial segregation is present in the relative perceived quality of products. We develop a novel Newton-Krylov method that is able to compute accurately and efficiently coarse fixed points when the underlying fine-scale dynamics is stochastic. The main novelty of the algorithm is in the elimination of the noise that is generated when estimating Jacobian-vector products using time-integration of perturbed initial conditions. We present numerical results that demonstrate the convergence properties of the numerical method, and use the method to show that macroscopic fronts in this model destabilise at a coarse symmetry-breaking bifurcation.Comment: This version of the manuscript was accepted for publication on SIAD