Improving the accuracy of the Chebyshev rational approximation method using substeps

Aarno E. Isotalo, Maria Pusa

Research output: Contribution to journalArticleScientificpeer-review

17 Citations (Scopus)


The Chebyshev rational approximation method (CRAM) for solving the decay and depletion of nuclides is shown to have a remarkable decrease in error when advancing the system with the same time step and microscopic reaction rates as the previous step. This property is exploited here to achieve high accuracy in any end-of-step solution by dividing a step into equidistant substeps. The computational cost of identical substeps can be reduced significantly below that of an equal number of regular steps, as the lower-upper decompositions for the linear solutions required in CRAM need to be formed only on the first substep. The improved accuracy provided by substeps is most relevant in decay calculations, where there have previously been concerns about the accuracy and generality of CRAM. With substeps, CRAM can solve any decay or depletion problem with constant microscopic reaction rates to an extremely high accuracy for all nuclides with concentrations above an arbitrary limit.
Original languageEnglish
Pages (from-to)65-77
JournalNuclear Science and Engineering
Issue number1
Publication statusPublished - 2016
MoE publication typeA1 Journal article-refereed


  • CRAM
  • decay calculations
  • substeps


Dive into the research topics of 'Improving the accuracy of the Chebyshev rational approximation method using substeps'. Together they form a unique fingerprint.

Cite this