Thanks to Sergey Fomel for writing much of the core nonstationary PEF code used for this algorithm, and also for insightful discussions.