(8 intermediate revisions by the same user not shown)
Line 1: Line 1:
 
==Abstract==
 
==Abstract==
  
We derive an algorithm for computing the wave-kernel functions cosh<math>\surd A</math>
+
We derive an algorithm for computing the wave-kernel functions <math>\cos h \sqrt A</math> and <math>\sin h \sqrt A</math> for an arbitrary square matrix <math>A</math>, where <math>\sin h c z = \sin h\frac{(z)}z</math>. The algorithm is based on Padé approximation and the use of double angle formulas. We show that the backward error of any approximation to <math>\cos h\sqrt{A}</math> can be explicitly expressed in terms of a hypergeometric function. To bound the backward error we derive and exploit a new bound for <math>\|A^k\|^{1/k}</math> that is sharper than one previously obtained by Al-Mohy and Higham [SIAM J. Matrix Anal. Appl., 31 (2009), pp. 970--989]. The amount of scaling and the degree of the Padé approximant are chosen to minimize the computational cost subject to achieving backward stability for <math>\cosh\sqrt{A}</math> in exact arithmetic. Numerical experiments show that the algorithm behaves in a forward stable manner in floating-point arithmetic and is superior in this respect to the general purpose Schur--Parlett algorithm applied to these functions.
  
==2 The main text<!-- You can enter and format the text of this document by selecting the ‘Edit’ option in the menu at the top of this frame or next to the title of every section of the document. This will give access to the visual editor. Alternatively, you can edit the source of this document (Wiki markup format) by selecting the ‘Edit source’ option.
+
==Full Document==
  
Most of the documents in Scipedia are written in English (write your manuscript in American or British English, but not a mixture of these). Anyhow, specific publications in other languages can be published in Scipedia. In any case, the documents published in other languages must have an abstract written in English.
+
<pdf>Media:Draft_Samper_112556199_4125_18m1170352.pdf</pdf>
 
+
 
+
2.1 Subsections
+
 
+
Divide your article into clearly defined and numbered sections. Subsections should be numbered 1.1, 1.2, etc. and then 1.1.1, 1.1.2, ... Use this numbering also for internal cross-referencing: do not just refer to 'the text'. Any subsection may be given a brief heading. Capitalize the first word of the headings.
+
 
+
 
+
2.2 General guidelines
+
 
+
Some general guidelines that should be followed in your manuscripts are:
+
 
+
*  Avoid hyphenation at the end of a line.
+
 
+
*  Symbols denoting vectors and matrices should be indicated in bold type. Scalar variable names should normally be expressed using italics.
+
 
+
*  Use decimal points (not commas); use a space for thousands (10 000 and above).
+
 
+
*  Follow internationally accepted rules and conventions. In particular use the international system of units (SI). If other quantities are mentioned, give their equivalent in SI.
+
 
+
 
+
2.3 Tables, figures, lists and equations
+
 
+
Please insert tables as editable text and not as images. Tables should be placed next to the relevant text in the article. Number tables consecutively in accordance with their appearance in the text and place any table notes below the table body. Be sparing in the use of tables and ensure that the data presented in them do not duplicate results described elsewhere in the article.
+
 
+
Graphics may be inserted directly in the document and positioned as they should appear in the final manuscript.
+
 
+
Number the figures according to their sequence in the text. Ensure that each illustration has a caption. A caption should comprise a brief title. Keep text in the illustrations themselves to a minimum but explain all symbols and abbreviations used. Try to keep the resolution of the figures to a minimum of 300 dpi. If a finer resolution is required, the figure can be inserted as supplementary material
+
 
+
For tabular summations that do not deserve to be presented as a table, lists are often used. Lists may be either numbered or bulleted. Below you see examples of both.
+
 
+
1. The first entry in this list
+
 
+
2. The second entry
+
 
+
2.1. A subentry
+
 
+
3. The last entry
+
 
+
* A bulleted list item
+
 
+
* Another one
+
 
+
You may choose to number equations for easy referencing. In that case they must be numbered consecutively with Arabic numerals in parentheses on the right hand side of the page. Below is an example of formulae that should be referenced as eq. (1].
+
 
+
 
+
2.4 Supplementary material
+
 
+
Supplementary material can be inserted to support and enhance your article. This includes video material, animation sequences, background datasets, computational models, sound clips and more. In order to ensure that your material is directly usable, please provide the files with a preferred maximum size of 50 MB. Please supply a concise and descriptive caption for each file. -->==
+
 
+
 
+
 
+
 
+
==3 Bibliography<!--
+
Citations in text will follow a citation-sequence system (i.e. sources are numbered by order of reference so that the first reference cited in the document is [1], the second [2], and so on) with the number of the reference in square brackets. Once a source has been cited, the same number is used in all subsequent references. If the numbers are not in a continuous sequence, use commas (with no spaces) between numbers. If you have more than two numbers in a continuous sequence, use the first and last number of the sequence joined by a hyphen
+
 
+
You should ensure that all references are cited in the text and that the reference list. References should preferably refer to documents published in Scipedia. Unpublished results should not be included in the reference list, but can be mentioned in the text. The reference data must be updated once publication is ready. Complete bibliographic information for all cited references must be given following the standards in the field (IEEE and ISO 690 standards are recommended). If possible, a hyperlink to the referenced publication should be given. See examples for Scipedia’s articles [1], other publication articles [2], books [3], book chapter [4], conference proceedings [5], and online documents [6], shown in references section below. -->==
+
 
+
 
+
 
+
 
+
==4 Acknowledgments<!-- Acknowledgments should be inserted at the end of the document, before the references section. -->==
+
 
+
 
+
 
+
 
+
==5 References<!--[1] Author, A. and Author, B. (Year) Title of the article. Title of the Publication. Article code. Available: http://www.scipedia.com/ucode.
+
 
+
[2] Author, A. and Author, B. (Year) Title of the article. Title of the Publication. Volume number, first page-last page.
+
 
+
[3] Author, C. (Year). Title of work: Subtitle (edition.). Volume(s). Place of publication: Publisher.
+
 
+
[4] Author of Part, D. (Year). Title of chapter or part. In A. Editor & B. Editor (Eds.), Title: Subtitle of book (edition, inclusive page numbers). Place of publication: Publisher.
+
 
+
[5] Author, E. (Year, Month date). Title of the article. In A. Editor, B. Editor, and C. Editor. Title of published proceedings. Paper presented at title of conference, Volume number, first page-last page. Place of publication.
+
 
+
[6] Institution or author. Title of the document. Year. [Online] (Date consulted: day, month and year). Available: http://www.scipedia.com/document.pdf.
+
-->==
+

Latest revision as of 16:31, 5 May 2020

Abstract

We derive an algorithm for computing the wave-kernel functions and for an arbitrary square matrix , where . The algorithm is based on Padé approximation and the use of double angle formulas. We show that the backward error of any approximation to can be explicitly expressed in terms of a hypergeometric function. To bound the backward error we derive and exploit a new bound for that is sharper than one previously obtained by Al-Mohy and Higham [SIAM J. Matrix Anal. Appl., 31 (2009), pp. 970--989]. The amount of scaling and the degree of the Padé approximant are chosen to minimize the computational cost subject to achieving backward stability for in exact arithmetic. Numerical experiments show that the algorithm behaves in a forward stable manner in floating-point arithmetic and is superior in this respect to the general purpose Schur--Parlett algorithm applied to these functions.

Full Document

The PDF file did not load properly or your web browser does not support viewing PDF files. Download directly to your device: Download PDF document
Back to Top

Document information

Published on 01/01/2018

DOI: 10.1137/18M1170352
Licence: CC BY-NC-SA license

Document Score

0

Views 2
Recommendations 0

Share this document

claim authorship

Are you one of the authors of this document?