Efficient and Parallel Solution of High-order Continuous Time Galerkin for Dissipative and Wave Propagation Problems

03/09/2023

∙

by   Zhiming Chen, et al.

∙

0

∙

share

We propose efficient and parallel algorithms for the implementation of the high-order continuous time Galerkin method for dissipative and wave propagation problems. By using Legendre polynomials as shape functions, we obtain a special structure of the stiffness matrix which allows us to extend the diagonal Padé approximation to solve ordinary differential equations with source terms. The unconditional stability, hp error estimates, and hp superconvergence at the nodes of the continuous time Galerkin method are proved. Numerical examples confirm our theoretical results.