ANALYSIS OF THE TRUNCATED SPIKE ALGORITHM

Download

index.pdf

Date

2008-01-01

Author

MİKKELSEN, CCK
Manguoğlu, Murat

Metadata

Show full item record

This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

Item Usage Stats

164
views

0
downloads

The truncated SPIKE algorithm is a parallel solver for linear systems which are banded and strictly diagonally dominant by rows. There are machines for which the current implementation of the algorithm is faster and scales better than the corresponding solver in ScaLAPACK (PDDBTRF/PDDBTRS). In this paper we prove that the SPIKE matrix is strictly diagonally dominant by rows with a degree no less than the original matrix. We establish tight upper bounds on the decay rate of the spikes as well as the truncation error. We analyze the error of the method and present the results of some numerical experiments which show that the accuracy of the truncated SPIKE algorithm is comparable to LAPACK and ScaLAPACK.

Subject Keywords

Analysis

URI

https://hdl.handle.net/11511/40174

Journal

SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS

DOI

https://doi.org/10.1137/080719571

Collections

Department of Computer Engineering, Article

Suggestions

OpenMETU
Core

NORMAL SOLVABILITY OF ELLIPTIC BOUNDARY-VALUE-PROBLEMS ON ASYMPTOTICALLY FLAT MANIFOLDS ERKIP, AK; SCHROHE, E (Elsevier BV, 1992-10-01) Normal solvability is shown for a class of boundary value problems on Riemannian manifolds with noncompact boundary using a concept of weighted pseudodifferential operators and weighted Sobolev spaces together with Lopatinski-Shapiro type boundary conditions. An essential step is to show that the standard normal derivative defined in terms of the Riemannian metric is in fact a weighted pseudodifferential operator of the considered class provided the metric is compatible with the symbols.
Local operator spaces, unbounded operators and multinormed C*-algebras Dosiev, Anar (Elsevier BV, 2008-10-01) In this paper we propose a representation theorem for local operator spaces which extends Ruan's representation theorem for operator spaces. Based upon this result, we introduce local operator systems which are locally convex versions of the operator systems and prove Stinespring theorem for local operator systems. A local operator C*-algebra is an example of a local operator system. Finally, we investigate the injectivity in both local operator space and local operator system senses, and prove locally conv...
AN OPTIMAL-CONTROL PROBLEM WITH NONLINEAR ELLIPTIC STATE-EQUATIONS Leblebicioğlu, Mehmet Kemal (Elsevier BV, 1992-02-01) In this article some of the results for optimal control of linear systems have been generalized to a nonlinear case. This is achieved by employing standard techniques of the nonlinear theory. After demonstrating the existence of optimal controls, finite element method is used to discretize the problem. The resulting finite dimensional problem is solved by a special algorithm. The theoretical discussions are completed by proving that approximate solutions are reduced to exact solutions as the element size te...
Development of an incompressible navier-stokes solver with alternating cell direction implicit method on structured and unstructured quadrilateral grids Baş, Onur; Tuncer, İsmail Hakkı; Department of Aerospace Engineering (2007) In this research, the Alternating Cell Direction Implicit method is used in temporal discretisation of the incompressible Navier-Stokes equations and compared with the well known and widely used Point Gauss Seidel scheme on structured and quadrilateral unstructured meshes. A two dimensional, laminar and incompressible Navier-Stokes solver is developed for this purpose using the artificial compressibility formulation. The developed solver is used to obtain steady-state solutions with implicit time stepping m...
Inverse Sturm-Liouville Systems over the whole Real Line Altundağ, Hüseyin; Taşeli, Hasan; Department of Mathematics (2010) In this thesis we present a numerical algorithm to solve the singular Inverse Sturm-Liouville problems with symmetric potential functions. The singularity, which comes from the unbounded domain of the problem, is treated by considering the limiting case of the associated problem on the symmetric finite interval. In contrast to regular problems which are considered on a finite interval the singular inverse problem has an ill-conditioned structure despite of the limiting treatment. We use the regularization t...

Citation Formats

C. MİKKELSEN and M. Manguoğlu, “ANALYSIS OF THE TRUNCATED SPIKE ALGORITHM,” SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, pp. 1500–1519, 2008, Accessed: 00, 2020. [Online]. Available: https://hdl.handle.net/11511/40174.