Abstract:An additive Runge-Kutta method is used for the time stepping, which integrates the linear stiff terms by an explicit singly diagonally implicit Runge-Kutta (ESDIRK) method and the nonlinear terms by an explicit Runge-Kutta (ERK) method. In each time step, the implicit solve is performed by the recently developed Hierarchical Poincaré-Steklov (HPS) method. This is a fast direct solver for elliptic equations that decomposes the space domain into a hierarchical tree of subdomains and builds spectral collocation solvers locally on the subdomains. These ideas are naturally combined in the presented method since the singly diagonal coefficient in ESDIRK and a fixed time-step ensures that the coefficient matrix in the implicit solve of HPS remains the same for all time stages. This means that the precomputed inverse can be efficiently reused, leading to a scheme with complexity (in two dimensions) $\mathcal{O}(N^{1.5})$ for the precomputation where the solution operator to the elliptic problems is built, and then $\mathcal{O}(N \log N)$ for the solve in each time step. The stability of the method is proved for first order in time and any order in space, and numerical evidence substantiates a claim of stability for a much broader class of time discretization methods. Numerical experiments supporting the accuracy of efficiency of the method in one and two dimensions are presented.

Extrapolating Solution Paths of Polynomial Homotopies towards Singularities with PHCpack and phcpy

A Family of Spacial PH Curves Represented by Algebraic-Trigonometric Functions and Their Applications

Geometric Hermite Interpolation by a Family of Spatial Algebraic-Trigonometric PH Curves

Extrapolating on Taylor Series Solutions of Homotopies with Nearby Poles

Chiral Glass Phase in Ceramic Superconductors

Certified homotopy tracking using the Krawczyk method

GPU Accelerated Newton for Taylor Series Solutions of Polynomial Homotopies in Multiple Double Precision

A Special Homotopy Continuation Method For A Class of Polynomial Systems

Polyhedral Homotopies in Cox Coordinates

Early Ending In Homotopy Path-Tracking For Real Roots

A stochastic homotopy tracking algorithm for parametric systems of nonlinear equations

On G 2 Approximation of Planar Algebraic Curves under Certified Error Control by Quintic Pythagorean-Hodograph Splines

Exponential node clustering at singularities for rational approximation, quadrature, and PDEs

A stratified polyhedral homotopy method for sampling positive dimensional zero sets of polynomial systems

Double Saddle-Point Preconditioning for Krylov Methods in the Inexact Sequential Homotopy Method

PolyCD: Optimization via Cycling through the Vertices of a Polytope

Hermite Spectral Method with Hyperbolic Cross Approximations to High-Dimensional Parabolic Pdes

Weight calculation and convergence analysis of polyharmonic spline (PHS) with polynomials for different stencils

Fast and high-order approximation of parabolic equations using hierarchical direct solvers and implicit Runge-Kutta methods

Continuation Path Learning for Homotopy Optimization

A Holomorphic Embedding Based Continuation Method for Identifying Multiple Power Flow Solutions