High performance verified computing using C-XSC

被引:0
作者
Walter Krämer
机构
[1] University of Wuppertal,Scientific Computing/Software Engineering, Faculty of Mathematics and Natural Sciences
来源
Computational and Applied Mathematics | 2013年 / 32卷
关键词
Verified computing; Self-validating methods; High performance computing; Parallelization; Thread-safety; Sparse methods; C-XSC; Primary 65G20; Secondary 65G30;
D O I
暂无
中图分类号
学科分类号
摘要
So called self-validating or self-verifying numerical methods allow to prove mathematical statements (existence of a fixed point, of a solution of an ODE, of a zero of a continuous function, of a global minimum within a given range, etc.) using a digital computer. To validate the assertions of the underlying mathematical theorems only fast finite precision machine arithmetic is used. The results are absolutely rigorous. We report on the accuracy as well as on the efficiency of the C++ class library C-XSC, our well known open source software tool designed to facilitate self-verifying numerical calculations. We focus mainly on solvers for dense and sparse interval linear systems. In recent years, these solvers have been improved significantly with respect to high performance computing within our bilateral Probral project HPVC (see Acknowledgments). As a motivating nontrivial example, where we need in an intermediate step an efficient solver for large dense interval linear systems, the computation of a verified functional enclosure for the solution of an integral equation is briefly discussed. The newest version C-XSC 2.5.1 released on June 9, 2011 allows using C-XSC in multi-threaded environments. The library as well as some further packages not mentioned in this paper are open source and freely available from the web site of the author’s research group Scientific Computing/Software Engineering at the University of Wuppertal: http://www2.math.uni-wuppertal.de/org/WRST/index_de.html.
引用
收藏
页码:385 / 400
页数:15
相关论文
共 32 条
[1]  
Choi J(1996)ScaLAPACK: a portable linear algebra library for distributed memory computers. Design issues and performance Comput Phys Comm 97 1-15
[2]  
Demmel J(2004)Algorithm 832: UMFPACK, an unsymmetric-pattern multifrontal method ACM Trans Math Software 30 196-199
[3]  
Dhillon I(1992)Sparse matrices in MATLAB: design and implementation SIAM J Matrix Anal Appl 13 333-356
[4]  
Dongarra J(2008)Dense linear system: a parallel self-verified solver Int J Parallel Program 36 412-425
[5]  
Ostrouchov S(2011)Efficient parallel solvers for large dense systems of linear interval equations Reliab Comput 15 193-206
[6]  
Petitet A(1986)The arithmetic of the digital computer: a new approach SIAM Rev 28 1-40
[7]  
Stanley K(2005)Accurate sum and dot product SIAM J Sci Comput 26 6-132
[8]  
Walker D(2010)A solver for complex-valued parametric linear systems Serdica J Comput 4 123-S442
[9]  
Whaley R(1995)Verified computation of the solution of large sparse linear systems Z Angew Math Mech (ZAMM) 75 S439-60
[10]  
Davis TA(1999)Fast and parallel interval arithmetic BIT Numer Math 39 53-452