clsqr1

所属分类:超算/并行计算
开发工具:Visual C++
文件大小:19KB
下载次数:1
上传日期:2010-12-22 17:31:09
上 传 者aabbas7
说明:  Hight Performance Computing Demos

文件列表:
clsqr1 (0, 2007-09-03)
clsqr1\lsqr.c (33857, 2007-09-03)
clsqr1\lsqr.h (16052, 2007-09-03)
clsqr1\test_lsqr.c (8341, 2007-09-03)
clsqr1\test_lsqr.h (582, 2000-10-07)
clsqr1\test_lstp.c (3788, 2007-09-03)
clsqr1\test_mult.c (2619, 2007-09-03)
clsqr1\test_prog.c (679, 2007-09-03)
clsqr1\makefile (837, 2000-10-07)

lsqr_c.README The software for LSQR (C version) is provided by SOL, Stanford University under the terms of the OSI Common Public License (CPL): http://www.opensource.org/licenses/cpl1.0.php 08 Sep 1999: First set of files provided by James W. Howse . 11 Oct 2000: Files unpacked by Michael Saunders. Lines longer than 80 characters had extraneous carriage returns. Tried to repair them, but not sure if I/O formats will be correct. 09 Dec 2005: Beware that this C version of LSQR allows x0 to be input. I RECOMMEND THAT YOU ALWAYS USE x0 = 0 (an n-vector of zeros). Otherwise, you must be conscious of the fact that LSQR will be solving A dx = r0 = b - A x0 and adding parts of dx to x0 AS IT GOES. If x0 is very close to a solution (that is, if dx is very small compared to x0), many digits of accuracy in dx will be lost. It is safer to solve A dx = r0 (with x0 = 0) and then compute x = x0 + dx yourself at the end. In either case, providing x0 will generally not save work (will not reduce the number of LSQR iterations) unless you use larger tolerances atol and btol than you would for Ax = b. 31 Aug 2007: Some bugs fixed (see lsqr.c). Note that some of the rectangular test problems at the end of test_prog.c are really too ill-conditioned for any least-squares method to give an accurate solution. Ignore the "LSQR seems to have failed" messages if ||A'r|| is quite small. 31 Aug 2007: Simplified C version of LSQR contributed by Michael Friedlander. Added as clsqr2 to http://www.stanford.edu/group/SOL/lsqr.html 02 Aug 2007: Dima Sorkin recommends changes to improve portability. Also the following compile command: gcc -ansi -pedantic -Wall -c lsqr.c Please send comments to Michael Saunders, SOL, Stanford University saunders@stanford.edu 650-723-1875 ----------------------------------------------------------------------------- The C version of LSQR involves the following files: lsqr.c lsqr.h makefile test_lsqr.c test_lsqr.h test_lstp.c test_mult.c test_prog.c test_lstp.c contains lstp, a least-square test problem generator. It constructs data for a matrix A and a vector b. test_mult.c computes products Ax and A'y. test_lsqr.c calls lstp to generate a problem, then calls lsqr to solve it. test_prog.c calls test_lsqr several times to generate and solve various problems.

近期下载者

相关文件


收藏者