
doi: 10.1137/0910046
The theme of this paper is that the primary computational bottleneck in the solution of stiff ordinary differential equations (ODEs) and the parallel solution of nonstiff ODEs is the implicitness of the ODE rather than the approximation of the integration process (or in conventional terminology, numerical stability rather than accuracy), and therefore it may be fruitful to apply (at least conceptually) the iterative techniques needed to overcome implicitness in continuous time, before discretization—to waveforms rather than values at a point in time. Several classical iterations, based on splitting, are discussed, but the emphasis is on those not based on a partitioning of the ODE system. The shifted Picard iteration is proposed as a compromise between the cheap but slow Picard iteration and the fast but expensive Newton iteration. By varying the shift parameter from one iteration to the next, a good rate of convergence seems possible. As an alternative, the author also examines the more classical acceleration technique applied to the Picard iteration. Some experimental results are given. However, the practical aspects of discretization are beyond the scope of this paper.
| selected citations These citations are derived from selected sources. This is an alternative to the "Influence" indicator, which also reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically). | 30 | |
| popularity This indicator reflects the "current" impact/attention (the "hype") of an article in the research community at large, based on the underlying citation network. | Average | |
| influence This indicator reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically). | Top 10% | |
| impulse This indicator reflects the initial momentum of an article directly after its publication, based on the underlying citation network. | Top 10% |
