LAPACK-3.2.2软件包
于2010年6月21日星期一更新
由田纳西州J.Langou和科罗拉多州丹佛市J.Languo创建的文件LAPACK-3.2.2是一个错误修复版本。LAPACK-3.2.2现在包括CMAKE构建系统。错误更正
==============

(*)bug0014::CSYEQUB/ZSYEQUB/SSYEQUB/DSYEQUB/CHEEQUB:EXTERNAL和INTRINSIC列表不准确(或甚至不存在) 
(*)svn::r660o Alexander V.Kobotov(Intel)于2009年4月6日星期一向“lapack@cs.utk.edu".c(he/sy)eqb:外部列表和lsame中没有slamch,内部没有描述D和S版本以及WORK数组也存在同样的问题。o朱利安·兰古(2009年5月10日星期日SVN第660版)承诺(*)bug0015::当CTGSEN/DTGSEN/STGSEN/ZTGSEN中的IJOB=0与标头描述冲突时,引用WORK(1)和IWORK
(*)svn::r658659o Alexander V.Kobotov(Intel)于2009年4月6日星期一向“lapack@cs.utk.edu".“(d/s)tgsen.f:iwork(1)总是被引用:第455行:iwork如果IJOB=0,则不应该,因此NULL指针会导致segfault。"C和Z版本以及WORK数组也存在同样的问题。o由Julien Langou承诺(SVN版本658和659,2009年5月10日星期日)(*)bug0016::CGEESX/DGEESX/SGEESX/ZGEESX的工作区查询中出现问题
(*)svn::r657o Alexander V.Kobotov(Intel)于2009年4月6日星期一向“lapack@cs.utk.edu".“(c/z)geesx:no lwork=-1分支,执行lquery时返回info=-15。”在WORKSPACE查询后,S和D版本也没有退出。o由Julien Langou承诺(SVN第657版,2009年5月10日星期日)(*)bug0017::在CHBGST/ZHBGST中,引用DO 240和DO 470循环中RWORK的未初始化元素
(*)svn::r656o Lawrence Mulholland(NAG)于2009年5月6日星期三发布补丁o参见论坛主题1408o由Julien Langou承诺(SVN第656版,2009年5月10日星期日)(*)bug0019::'make-install'太容易生成“***错误:操作时间为零”
(*)bug0019::“make-install”可以生成负触发器计数
(*)svn::r693o Sergei Steshenko于2008年11月20日星期四报道o参见:电子邮件o由Julie Langou承诺(SVN第693版,2009年9月11日星期五)(*)bug0021::在DGESDD中,工作区查询提供的值小于标头中给定的最小值,以运行例程
(*)svn::r649o由Guy Bencteux于2008年12月6日星期六报道o参见论坛主题846o朱莉·兰古承诺(SVN第649版,2009年3月17日星期二)(*)bug0023::SROTMG和DROTMG使用了弃用的Fortran ASSIGN语句和赋值的GOTO语句
(*)svn::r691o由Marco Atzeri于2009年3月6日星期五报道o参见论坛主题1420o由Julie Langou承诺(SVN第691版,2009年9月10日星期四)(*)bug0031::注释中的空白修复
(*)svn::r692o由“Jan”于2009年6月2日星期二报道o参见论坛主题1449o由Julie Langou承诺(SVN第692版,2009年9月11日星期五)(*)bug0038::使用DGELSS时DLASET中超出范围的元素访问
(*)svn::r695o Tim Hopkins于2009年7月31日星期五发送的错误报告o看论坛主题1582o由Julie Langou承诺(SVN第695版,2009年9月25日星期五)(*)bug0041::中打印格式错误/测试/EIG/[s,c,d,z]drvgbx。(f)
(*)svn::r696o Michael Chuvelev(Intel)于2009年7月29日星期二向“lapack@cs.utk.edu".“lapack-3.2.1/TESTING/EIG/zdrvgbx.f中的错误”o打印格式错误,如果zgbsvxx中出现错误,这可能会导致分段错误。o由Julien Langou承诺(SVN第696版,2009年9月29日星期二)(*)bug0042::中打印格式错误/测试/EIG/[s,c,d,z]drvgbx。(f)
(*)svn::r697o Kevin Wadleigh于2009年9月29日星期二向“lapack@cs.utk.edu".“LAPACK例程ieeeck.f和编译器优化”o应修改LAPACK函数IEEECK,以避免编译器优化代码行NAN6=NAN5*0.0并将其替换为NAN6=0.0。我看到过在默认优化级别使用两个不同的编译器时会发生这种情况。但是,使用NAN6=NAN5*零如附件中所示,避免发生这种情况。o Julie Langou承诺(SVN修订版697,2009年9月30日星期三)(*)bug0043::VL,VU在lapack-3.2.1 dsyevr中引用,即使范围为.ne.'V'
(*)svn::r699o Michael Chuvelev(Intel)于2009年10月14日星期三向“lapack@cs.utk.edu".“VL、VU在lapack-3.2.1 dsyevr中引用,即使范围为.ne.'V'”o由Julien Langou承诺(SVN第699版,2009年10月20日星期二)(*)bug0044::N=1时dsyevr.f中的bug
(*)svn::r704o Ashutosh Mahajan(阿贡国家实验室数学和计算机科学部)于2009年12月11日星期五向“lapack@cs.utk.edu".o由Julien Langou承诺(SVN第704版,2009年12月15日星期二)(*)bug0045::N=2时dstemr.f中的bug
(*)svn::r705o Ashutosh Mahajan(阿贡国家实验室数学和计算机科学部)于2009年12月10日(星期四)向“lapack@cs.utk.edu".o由Julien Langou承诺(SVN第705版,2009年12月15日星期二)(*)bug0038::DGELSD中没有LIWORK查询
(*)svn::r698o 2009年8月18日星期二Mathew Cross的错误报告o看论坛主题1601o由Julie Langou承诺(SVN第698版,2009年10月20日星期四)(*)bug0046::ZGESDD中记录错误的RWORK工作区
(*)svn::r729o 2010年3月15日星期一Zbigniew的错误报告o参见论坛主题1779o朱莉·兰古承诺(SVN 729版,2010年3月18日星期四)(*)bug0047::*gerfsx中的R/C参数描述不正确
(*)svn::r726o Michael Chuvelev(Intel)于2010年3月1日星期三向“lapack@cs.utk.edu".o James Demmel发送的错误修复“R/C仅按用途输入”“这些例程中没有FACT参数”o朱莉·兰古承诺(SVN修订版r726,2010年3月9日星期二)(*)bug0048::如果NaN位于输入矩阵中,则*gebal中可能会发生挂起。
(*)svn::r730o Alexander Kobotov(Intel)于2010年2月3日星期三向“lapack@cs.utk.edu".o朱莉·兰古承诺(SVN修订版r730,2010年3月9日星期二)(*)bug0049::在快速退出(n.eq.0)或在*laed8中等级1修饰符足够小的情况下,输出参数GIVPTR保持未初始化状态。
(*)序列号:r731o Alexander Kobotov(Intel)于2010年2月3日星期三向“lapack@cs.utk.edu".o朱莉·兰古承诺(SVN修订版r731,2010年3月9日星期二)(*)bug0050::STFSM.f中矩阵A第348-249行后面的额外参数导致SGEMM调用崩溃。
(*)svn::r732o Alexander Kobotov(Intel)于2010年2月3日星期三向“lapack@cs.utk.edu".o朱莉·兰古承诺(SVN修订版r732,2010年3月9日星期二)(*)由于xTGEX2中的阈值太严格,xDGGES中的错误0051::INFO=N+2。
(*)svn::r735o Ondra Kamenik于2010年3月16日星期二发送的错误报告o参见论坛主题1783o我们已在xGEQR2中将xLARFG更改为xLARFP。xGEQR2由xDGGES使用。这个因此,矩阵是不同的(尽管它们都是正确的)对于xLARFP生成的矩阵,xTGEX2超出了THRES值10%,导致错误标志。我们决定将阈值增加2倍。尽管该错误是由xLARFG更改为xLARFP、xLARFP触发的并不是真的感到内疚。o朱莉·兰古于2010年3月18日星期四确认o Guillaume Revy和Jim Demmel提供的修复o朱利安·兰古承诺(SVN修订版r735,2010年4月2日星期四)(*)bug0052::调用中的NRHS参数不正确?SYTRS来自?拉丁美洲_{po,sy,he}rfsx扩展。
(*)svn::r736o Michael Chuvelev(Intel)于2010年4月2日星期五向“lapack@cs.utk.edu".o朱莉·兰古承诺(SVN修订版r736,2010年4月2日星期五)(*)bug0053::修复测试套件文件cebchvxx中的外部声明问题。(f)
(*)序列号:r738o Michael Chuvelev(Intel)于2010年4月5日星期五向“lapack@cs.utk.edu".o朱莉·兰古承诺(SVN修订版r738,2010年4月5日星期五)(*)bug0054::修复注释某些参数描述的类型错误。
(*)svn::r739o Paul Roberts(NAG)于2010年4月12日星期五向“lapack@cs.utk.edu".o朱莉·兰古承诺(SVN修订版r739,2010年4月12日星期五)(*)bug0055::[csz]syequb中缺少工作和uplo的参数描述
(*)svn::r740o Paul Roberts(NAG)于2010年4月13日星期五向“lapack@cs.utk.edu".o朱莉·兰古承诺(SVN修订版r740,2010年4月13日星期五)(*)bug0058::xLASQ3中的变量N0是INPUT/OUTPUT。(评论说是INPUT。)
(*)svn::r750o 2009年10月16日,William Gandler(NIH/CIT)提供了错误报告和错误修复;4月23日,2010; 2010年5月1日和5月3日。至“lapack@cs.utk.edu".o朱莉·兰古承诺(SVN修订版r750,2010年5月27日星期五)(*)bug0059::如果实际工作空间的分配与工作空间查询中报告的完全相同,则ZGELSD在大型矩阵上崩溃(eq M=N=NHRS>=88)。
(*)svn::r765o Alexander Kobotov(Intel)于2010年6月23日星期三向“lapack@cs.utk.edu".o根本原因是错误估计了ZLALSA中的实际工作空间大小(在路径zgelsd->zlalsd->ZLALSA上)。它假设RWORK数组长度应为max(N,(SMLSZ+1)*NRHS*3),其中SMLSIZ=25。但事实上,该函数调用zlals0,它要求RWORK为(K*(1+NRHS)+2*NRHS),其中K可能与矩阵A中的列数相同。CGELSD也是如此。o朱莉·兰古承诺(SVN修订版r765,2010年6月24日星期四)
    
(*)bug0060::dtgsy2中的Bug,缺少IWORK调用
(*)svn::r767o Vasile Sima于2009年2月24日向“lapack@cs.utk.edu".o错误是“IWORK”必须出现在第640行。(它包含在1999年版本中。)注意,语句IE=(I+1)-1很奇怪;此外,在其他类似的地方,IWORK的用法如下。o朱莉·兰古承诺(SVN修订版r767,2010年6月24日星期四)(*)cpot01和zpot01中的工作空间类型错误
(*)svn::r655o Julie Langou于2009年4月30日报告的Bugo CPOT01和ZPOT01需要的工作区类型,它们需要是REAL(或DOUBLE PRECISION)而不是复合体(或双复合体)。o由Julie Langou承诺(SVN修订版r655,于2009年5月3日星期日)(*)将DLARFG恢复到3.1.0,将DLARAFG移动到DLARFP并改进DLARFP
(*)svn::(r706、r707、r708、r710、r712、r713)+(r755)o Pat Quillen(MathWorks)报告的Bugo Householder反射器生成例程dlarfg(用于双精度)已由dlarfp在LAPACK-3.2中更改为非负对角线元素(在QR的R系数的对角线上因式分解)。中的dlarfp子例程中存在问题LAPACK-3.2。错误报告来自MathWorks的Pat Quillen参与了错误修复:Pat Quillen、Jim Demmel、Sven Hammarling、,马克·霍曼(Mark Hoemmen)、纪尧姆·雷维(Guillaume Revy)、威廉·卡汉(William Kahan)、朱莉·兰古(Julie Langou)和朱利安(Julien)朗古。o Jim Demmel、Guillaume Revy和William Kahan的大量工作发现新的Householder反射器保证了QR中R的对角线与旧反射器:1) 计算出的反射器的正交度可以减小约4倍,(吉姆、纪尧姆和威廉不知道为什么。)2) 这一新程序本质上对过流/欠流更为敏感,这一点现在已经很清楚了,这里省略了细节。传统的Householder反射器没有出现这些问题。o所以决定1) 尽可能多地修复xLARFP例程。此问题已由更正朱利安(参见(r706、r707、r708、r710、r712、r713))基于吉姆的想法。2) 退出LAPACK中的更改以恢复使用旧反射器(xlarfg)无处不在(例如在xGEQRF、xGEQR2等中)。3) 当你有很多尾随的零时,保留可以避免失败的更改,因为这是一个独立的好主意。4) 将名称从xLARFP更改为xLARFGP,并引入新的QR例程(xGEQRFP和xGEQR2P),使用xlarfp保证新的正极R的对角线,适合需要它的人。5) 这需要调整测试套件6) 注意:只有QR例程具有xGEQRFP。虽然看起来很不错一些用户关心带非负对角线R.Noone的QR分解似乎关心LQ、QL或RQ上的此功能。所以我们省略了这些例行程序。7) 步骤2至6由Julie在r755期间完成。(*)bug0025::DGESDD中存在严重问题:DGESD返回负奇异值,INFO=0
(*)svn::r760o由Guy Bencteux于2008年12月6日星期六报道o参见论坛主题847o当时,朱莉无法复制该病毒。该错误高度依赖于BLAS。o 2010年6月19日星期六,ESIAL的Bruno Pinçon报告了另一个具有相同问题的矩阵o参见论坛主题1868o我们能够使用gfortran参考LAPACK重现MacOS的问题和vecLib BLAS,我们无法使用参考LAPACK和参考BLAS重现此问题。o Julien Langou提出并承诺了一个补丁。(见论坛主题1868。)o注:有强有力的证据表明盖伊·本克特克斯和这只臭虫都是一样的,尽管我们从来没有重现过盖伊的问题。(*)注释修复o第662版:SRC/clarrv.f、SRC/dlaeda.f、FRC/dlarrv.f、SLC/slaeda.f、SRC/slarrv.f、SLC/zlarrv.fo修订版664:SRC/chesvxx.f、SRC/cposvxx.f、SRC/csysvxx.f、SRC/dposvxx.f、SRC/dsysvxx.f、SRC/sposvxx.f,SRC/ssysvxx.f、SRC/zhesvxx.f、SRC/zposvxx.f/SRC/zsysvxx.fo第692版:SRC/dgbrfsx.f、SRC/dgbsvxx.f、SRC/dgejsv.f、SRC/dgerfsx.f、SRC/dgesc2.f、SRC/dgesvj.f、,SRC/dgesvxx.f、SRC/dggbal.f、SRC/dgsvj0.f、SRC/dgsvj1.f、SRC/disnan.f、SRS/dla_lin_berr.f、,SRC/dla_porcond.f、SRC/dlaein.f、SRC/dlaexc.f、SRC/dlaisnan.f、SRC/dlaneg.f、SLC/dlaqr.f、,SRC/dlaqr3.f、SRC/dlarra.f、SRC/dlarre.f、FRC/dlarrf.f、SLC/dlarrj.f、SRC/dlarrv.f、,SRC/dlarscl2.f、SRC/dlatdf.f、SRC/dpbequ.f、SRC/dpftri.f、SLC/dpftrs.f、FRC/dporfsx.f、,SRC/dpstf2.f、SRC/dpstrf.f、SRC/dsgesv.f、FRC/dspcon.f、RSC/dsposv.f、SLC/dsycon.f和SRC/dsyequb.f,SRC/dsygvx.f、SRC/dsyrfsx.f、SRC/dsysvxx.f、SRC/dtftri.f、SRC/dtfttp.f、FRC/dtfttr.f、SLC/dtgex2.f、,SRC/dtgexc.f、SRC/dtpttf.f、SRC/dtrttf.fo第714版:SRC/chfrk.f、SRC/chla_transtype.f、SRC/cla_gbamv.f、SRC/cla_gbrpvgrw.f、,SRC/cla_heamv.f、SRC/cla_ herpvgrw.f、SRC/cla_lin_berr.f、SRC/cla_porpvgrw.f、SRC/cla_rpvgrw,SRC/cla_syamv.f、SRC/cla_syrpvgrw.f、SRC/cla_wwaddw.f,SRC/dla_gbrpvgrw.f、SRC/dla _geamv.f、SRC/dla _ lin_berr.f、SRC/dla _ orpvgrw.f、SRC/dla _ rpvgrw.f、,SRC/dla_syamv.f、SRC/dla_syrpvgrw.f、SRC/dla_waddw.f、SRC/dlaisnan.f、SRC/dlarrf.f、SRC/dlasdt.f,SRC/dlasq3.f、SRC/dsfrk.f、SRC/dtfsm.f、SRC/ilaclc.f、SRC/ilaclr.f、SRC:iladlc.f、SRC-iladlc/f、SRC/iladlr.f、,SRC/ilaslc.f、SRC/ilassr.f、SRC/ilazlc.f,SRC/sla_gbrpvgrw.f、SRC/sla_geamv.f、SRC/sla_ lin_berr.f、SRC/sla_porpvgrw.f、SRC/sla_rpvgrw.f、,SRC/sla_syamv.f、SRC/sla_syrpvgrw.f、SRC/sla_wwaddw.f,SRC/slaisna.f、SRC/slaeg.f、SRC/slarrf.f、,SRC/slasdt.f、SRC/slassq3.f、SRC/ssfrk.f、SRC/stfsm.f、SRC/xerbla_array.f、SRC/zhfrk.f和SRC/zla_gbamv.f,SRC/zla_gbrpvgrw.f、SRC/zla _geamv.f、SRC/zla _ heamv.f,SRC/zla_rpvgrw.f、SRC/zla_rpvgrw.f、SRC/zla_syamv.f、SRC/zla_syrpvgrw.f、SRC/zla_waddw.f、SRC/ztfsm.f0第718版:SRC/disnan.f、SRC/dlaisnan.f、SRC/dlarrf.f、SRC/sisnan.f、SRC/slaisnan.h、SRC/srarrb.f、,SRC/slarrc.f、SRC/slARR.f、SRC/slarrf.f、SRC/slarrj.f、SLC/slarrrr.f、SRC/slarrv.fo第719版:SRC/cgeqpf.f、SRC/chseqr.f、SRC/dgeqpf.f、SRC/dhseqr.f、SRC/dlanv2.f、SRC/sgeqpf.f、SRC/shseqr.f,SRC/slav2.f、SRC/zgeqpf.f、SRC/zhseqr.fo第726版:SRC/cgerfsx.f、SRC/dgerfsx.f、SRC/sgerfsx.f、SRC/zgerfsx.fo第733版:SRC/cpstf2.f、SRC/cpstrf.f、SRC/ctfsm.f、SRC/dsgesv.f、SLC/dsposv.f、,SRC/spstrf.f、SRC/stfsm.f、SRC/zcgesv.f、SRC/zcposv.f、SLC/zpstf2.f、SRC/zpstrf.fo第737版:SRC/ilaprec.fo第741版:TESTING/MATGEN/clahilb.f、TESTING/MATGEN/clarror.f、ETSTING/MATGEN/clatm1.f、TEST/MATGEN/clatm2.f,TESTING/MATGEN/clatm3.f、TESTING/MATGEN/clatm6.f、TESTING/MATGEN/clatme.f、TEST/MATGEN/clatmr.f、,TESTING/MATGEN/clatms.f、TESTING/MATGEN/clatmt.f、TEST/MATGEN/dlahilb.f、TETSING/MATGEN/dlatm1.f、,TESTING/MATGEN/dlatm2.f、TESTING/MATGEN/dlatm3.f、TEST/MATGEN/dlatme.f、TETSING/MAT GEN/dlat mr.f、,TESTING/MATGEN/dlatms.f、TESTING/MATGEN/dlatmt.f、STING/MATGEN/slahilb.f、ETING/MATGEN/slatm1.f、,TESTING/MATGEN/slatm2.f、TESTING/MATGEN/slatm3.f、TETSING/MAT GEN/slame.f、TEST/MAT GEN/slatemr.f、,TESTING/MATGEN/slams.f、TESTINS/MATGEN/slatmt.f、TEST/MATGEN/zlahilb.f、STESTING/MATGEN/zlaror.f、,TESTING/MATGEN/zlatm1.f、TESTING/MATGEN/zlatm2.f、TEST/MATGEN/zlatm3.f、TETSING/MAT GEN/zlatem.f、,测试/MATGEN/zlatmr.f、测试/MATGEN/zlatms.f、测试/MATGEN/zlatmt.fo第748版,SRC/dtfttp.f、SRC/dtfttr.f、SRC/dtpttf.f、FRC/dtrttf.f,SRC/slansf.f、SLC/spftrf.f、,SRC/spftri.f、SRC/spfctrs.f、SRC/stfsm.f、SRC-stftri.f、SRC/stfttp.f、SRC/stfttr.f、,SRC/stpttf.f、SRC/strttf.fo第749版:SRC/dla_gbrcond.f、SRC/sla_gbrbond.fo第756版:SRC/dlasd0.f、SRC/slasd0.fo第757版:SRC/dbdsdc.f、SRC/sbdsdc.fo第758版:SRC/clalsd.f、SRC/dlalsd.f、SRC/dlasd1.f、SRC/dlasd3.f、SRC-dlasd6.f、SRC/dlasd8.f、SRC/dlasda.f、,SRC/slalsd.f、SRC/slasd1.f、SRC/slasd3.f、SRC-slasd6.f、SRC/slasd8.f、SRC/slasda.f、SRC/zlalsd。(f)以前的发行说明
======================o个LAPACK-3.2.1软件包o个LAPACK-3.2软件包