皮埃尔·鲁克斯;元朗沃罗宁;桑卡拉纳拉亚南，斯里拉姆 验证多项式不变量的数值半定规划求解器。 （英语） Zbl 1425.68081号 形式方法系统。设计。 53，第2期，286-312（2018）. 摘要：在许多程序验证任务中，半定规划（SDP）解算器越来越多地用作原语，以合成和验证各种系统（包括程序、混合系统和随机模型）的多项式不变量。一方面，它们为半代数约束的推理提供了一种易于处理的替代方法。然而，由于“数值问题”，结果往往不可靠，其中包括许多原因，如浮点错误、病态问题、严格可行性失败，以及更普遍的，用于求解SDP的算法的细节。这些问题影响最终数值结果是否可信。在本文中，我们简要介绍了SDP求解器在静态分析社区中的新兴应用。我们报告了将SDP解算器用于常见不变合成任务的风险，描述了可能导致不可靠答案的常见故障。接下来，我们将演示保半定的现有工具通常证明不足以满足我们需求的编程。最后，我们提出了一个验证半定规划的解决方案，该方案可用于检查求解器输出解的可靠性，以及一个填充过程，该过程可检查求解器所输出解的可行附近解的存在性。我们报道了一些成功的涉及填充程序的初步实验。 引用于4文件 理学硕士： 68号30 软件工程的数学方面（规范、验证、度量、需求等） 65千5 数值数学规划方法 90C22型 半定规划 关键词：静态分析;多项式不变量;数值优化;平方和松弛;半定规划;浮点运算;形式证明 软件：塞杜米;VSDP公司;红线批注;DIMACS公司;CSDP公司;d真实;SDPT3系统;QEPCAD公司;轨道库;标准差分法;莫塞克;优化软件的决策树;github;frlib公司;无SPOT;CalCS公司 PDF格式BibTeX公司 XML格式引用 \textit{P.Roux}等人，《形式方法系统》。设计。53，第2号，286--312（2018；Zbl 1425.68081） 全文： 内政部 哈尔 参考文献： [1] AdjéA，Garoche P-L，Magron，V（2015）使用平方和优化的基于属性的多项式不变量生成。包含：SAS，第235-251页 [2] AdjéA，Gaubert S，Goubaulté（2010）将策略迭代与半定松弛耦合，以计算静态分析中的精确数值不变量。In:员工持股计划·Zbl 1260.68082号 [3] Ahmadi AA、Majumdar A（2014）DSOS和SDSOS优化：基于LP和SOCP的平方和优化替代方案。参加：信息科学和系统年会（CISS） [4] Alipanahi，B。；Krislock，N。；Ghodsi，A。；Wolkowicz，H。；唐纳森，L。；Li，M.，通过半定规划从NOESY距离约束确定蛋白质结构，计算生物学杂志，20，296-310，（2013）·doi:10.1089/cmb.2012.0089 [5] Alizadeh，F。；海伯里，J-PA；Overton，ML，半定规划的原对偶内点方法：收敛速度、稳定性和数值结果，SIAM J Optim，8746-768，（1998）·Zbl 0911.65047号 ·doi:10.1137/S1052623496304700 [6] Allamigeon X，Gaubert S，Goubault E，Putot S，Stott N（2015）一种可扩展的代数方法，用于推断切换系统的二次不变量。收件人：EMSOFT [7] Anjos MF，Lasserre JB（2012）半定、二次曲线和多项式优化简介。收录：半定、二次曲线和多项式优化手册。施普林格·Zbl 1235.90002号 [8] Bagnara R、Rodríguez-Carbonell E、Zafanella E（2005）使用凸多面体生成基本半代数不变量。输入：SAS·兹比尔1141.68446 [9] Basu S、Pollock R、Roy M-F（2006）《实代数几何中的算法》。施普林格，柏林·兹伯利1102.14041 [10] Ben-Tal A，Ghaoui LE（2009）稳健优化。普林斯顿应用数学系列。普林斯顿大学出版社·Zbl 1221.90001号 ·doi:10.1515/9781400831050 [11] Ben Sassi MA，Sankaranarayanan S，Chen X，Abraham E（2015）多项式lyapunov函数合成的多项式正性线性松弛。IMA J数学控制信息·Zbl 1397.93189号 [12] Bernstein SN（1912）《Weierstrass-fonde e sur le calcul des probabiles》。哈尔科夫数学社区2 [13] Borchers B（1999）CSDP，半定编程的C库。Optim方法软件 [14] Borwein，JM，有限维范围抽象凸规划的最优性表征，J.Aust。数学。Soc.序列号。A、 30390-411（1980）·Zbl 0469.90088号 ·doi:10.1017/S1446788700017882 [15] Borwein JM，Wolkowicz H（1980/81）con-convex编程问题的面部简化。J.奥斯特。数学。Soc.序列号。A、， [16] JM Borwein；Wolkowicz，H.，正则化抽象凸规划，数学分析应用杂志，83，495-530，（1981）·Zbl 0467.90076号 ·doi:10.1016/0022-247X（81）90138-4 [17] Boyd S，Vandenberghe L（2004）凸优化。剑桥大学出版社·Zbl 1058.90049号 ·doi:10.1017/CBO9780511804441 [18] Burkowski，F。；张，Y-L；Wolkowicz，H.，《蛋白质构象中选择旋转异构体的半定规划的有效使用》，Inform J Compute，26748-766，（2014）·Zbl 1304.90224号 ·doi:10.1287/ijoc.2014.0596 [19] Chakarov A，Voronin Y-L，Sankaranarayanan S（2016）几乎确定持久性和递归性的演绎证明。输入：TACAS [20] 张，Y-L；舒尔，S。；Wolkowicz，H。；Bailey，DH（编辑）；Bauschke，HH（编辑）；Borwein，P.（编辑）；Garvan，F.（编辑）；Thera，M.（编辑）；Vanderwerff，J.（编辑）；Wolkowicz，H.（编辑），退化半定程序的预处理和正则化，（2013），柏林 [21] Cheung Y-L（2013）通过面约简实现半定规划的预处理和约简：理论与实践。滑铁卢大学博士论文 [22] Collins GE（1975）通过柱面代数分解消除实闭场的量词。In：自动机理论和形式语言 [23] Collins GE，Hong H（1991）量词消去的部分柱面代数分解。J Symb计算机·Zbl 0754.68063号 [24] Cousot P（2005）通过参数抽象证明程序不变性和终止性。拉格朗日松弛和半定规划。输入：VMCAI·Zbl 1111.68503号 [25] Cousot P，Cousot R（1977）抽象解释：通过构造或近似不动点对程序进行静态分析的统一格模型。输入：POPL·Zbl 1149.68389号 [26] Dang T，Gawlitza TM（2011）基于模板的仿射混合自动机无界时间验证。输入：APLAS·Zbl 1348.68100号 [27] Demmel J（1989）关于Cholesky中的浮点错误。LAPACK工作说明14 CS-89-87。美国田纳西州诺克斯维尔市田纳西大学计算机科学系 [28] Dolzmann A，Sturm T（1997）REDLOG：计算机代数与计算机逻辑相结合。ACM SIGSAM公牛 [29] Dür M，Jargalsaikhan B，Still G（2012）Slater条件在线性圆锥规划中是通用的 [30] Ghaoui，L。；Oustry，F。；Lebret，H.，不确定半定规划的鲁棒解，SIAM J Optim，9，33-52，（1998）·Zbl 0960.93007号 ·doi:10.1137/S1052623496305717 [31] Farouki RT（2012）《伯恩斯坦多项式基础：百年回顾》。计算机辅助Geomet-Des·Zbl 1252.65039号 [32] 费伦埃（2010）《从控制系统到控制软件》。控制系统，IEEE [33] Fränzle M，Herde C，Teige T，Ratschan S，Schubert T（2007）具有复杂布尔结构的大型非线性算术约束系统的高效求解。SAT/CP积分的J满意布尔模型计算专题·Zbl 1144.68371号 [34] Gao S，Kong S，Clarke EM（2013）Dreal:实际非线性理论的SMT求解器。In：自动扣减国际会议（CADE），第208-214页·Zbl 1381.68268号 [35] Gaubert S，Goubault E，Taly A，Zennou S（2007）通过关系域上的策略迭代进行静态分析。In:员工持股计划·Zbl 1187.68151号 [36] Gawlitza T，Seidl H（2007）通过策略迭代精确计算不动点。In:员工持股计划·Zbl 1187.68152号 [37] Gawlitza TM，Monniaux D（2011）通过SMT解决改进策略。在：员工持股计划·兹比尔1326.68093 [38] Gawlitza TM，Seidl H（2010）《精确计算二次区域的放松抽象语义》。输入：SAS·Zbl 1306.68025号 [39] Gruber，G。；Rendl，F.，半定规划中不适定问题的计算经验，计算优化应用，21，201-212，（2002）·Zbl 0988.90024号 ·doi:10.1023/A:1013716917710 [40] Handelman D（1988）用紧凸多面体上的正线性函数表示多项式。太平洋数学杂志·Zbl 0659.52002年 [41] Harrison J（2007）通过平方和验证非线性实公式。输入：TPHOL·Zbl 1144.68357号 [42] Härter V，Jansson C，Lange M（2016）VSDP：验证半定规划。http://www.ti3.tuhh.de/jansson/vsdp/。3月28日访问 [43] Helmberg C（2012）半定规划。https://www-user.tu-chemnitz.de/helmberg/semidef.html。上次更新时间：·Zbl 1008.90044号 [44] Henrion D，Naldi S，Din MS El（2015）线性矩阵不等式的精确算法。arXiv预打印arXiv:1508.03715 [45] IEEE计算机学会。浮点运算IEEE标准。IEEE标准754-20082008 [46] Jansson C，Chaykin D，Keil C（2007）半定规划中最优值的严格误差界。SIAM J数字分析·Zbl 1167.90009号 [47] Kaltoffen E，Li B，Yang Z，Zhi L（2012）通过有理系数有理函数的平方和在全局多项式优化中的精确证明。J Symb计算机·Zbl 1229.90115号 [48] Krislock，N。；Wolkowicz，H.，使用半定表示和面部还原的显式传感器网络定位，SIAM J Optim，20，2679-2708，（2010）·Zbl 1229.90250号 ·doi:10.1137/090759392 [49] Lasserre JB（2001）《多项式全局优化与矩问题》。SIAM J优化·Zbl 1010.90061号 [50] Löfberg J（2009）实践中的前置和后置平方和程序。In:IEEE自动控制事务·Zbl 1367.90002号 [51] Maréchal A，FouilhéA，King T，Monniaux D，ël Périn M（2016）使用Handelman定理的多元多项式的多面体逼近。输入：VMCAI [52] Martin-Dorelé，Roux P（2017）使用数值解算器和浮点计算实现多项式正性的自反策略。参见：Yves B，Viktor V（eds）第六届ACM SIGPLAN认证程序和证明会议记录，2017年CPP，2017年1月16日至17日，法国巴黎，第90-99页。ACM、， [53] Mittelmann HD（2016）优化软件的决策树：半定规划。http://plato.asu.edu/sub/nlores.html#semidef。3月28日访问 [54] Monniaux D，Corbineau P（2011）《关于在退化案件中产生积极的tellensatz证人》。输入：ITP·Zbl 1342.68296号 [55] MOSEK ApS（2015）MOSEK C优化器API手册7.1版（第40版） [56] Nakata M（2010）半定规划求解器的高精度多精度算术版本的数值评估：SDPA-GMP、-QD和-DD。In：计算机辅助控制系统设计 [57] Nesterov Y，Nemirovskii A（1994）凸规划中的内点多项式算法。Soc工业应用数学·Zbl 0824.90112号 [58] 聂，J。；Ranestad，K。；Sturmfels，B.，《半定规划的代数度》，《数学程序》，122，379-405，（2008）·Zbl 1184.90119号 ·doi:10.1007/s10107-008-0253-6 [59] Nuzzo P、Puggelli A、Seshia SA、Sangiovanni-Vincentelli AL（2010）计算：非线性凸约束的SMT求解。摘自：Bloem R，Sharygina N（编辑）第十届计算机辅助设计形式方法国际会议论文集。瑞士，FMCAD 2010，卢加诺，10月20日至23日，第71-79页。电气与电子工程师协会 [60] Oulamara M，Venet AJ（2015）用高维椭球体和圆锥外推进行抽象解释。输入：CAV·Zbl 1381.68178号 [61] Parrilo P（2000）结构化半定程序和鲁棒性和优化中的半代数几何方法。加州理工学院博士论文 [62] Parrilo PA（2003）半代数问题的半定规划松弛。数学课程·Zbl 1043.14018号 [63] Pataki，G。；Bailey，DH（编辑）；Bauschke，HH（编辑）；Borwein，P.（编辑）；Garvan，F.（编辑）；Thera，M.（编辑）；Vanderwerff，J.（编辑）；Wolkowicz，H.（编辑），二次曲线线性规划中的强对偶：面约简和扩展对偶，613-634，（2013），柏林·兹比尔1282.90231 [64] Percenter F，Parrilo P（2014）《面部部分缩小：通过psd锥体近似值的简化等效sdps》。arXiv预印本arXiv：1408.4685·Zbl 1405.90098号 [65] Percenter F，Parrilo PA（2015）《SDP面部整容工具》。https://github.com/frankpercenter/frlib。2015年4月10日访问 [66] Peyrl H，Parrilo PA（2008）用有理系数计算平方和分解。Theor计算机科学·Zbl 1156.65062号 [67] Platzer A、Quesel J-D、Rümmer P（2009）《真实世界验证》。输入：CADE [68] Prajna S，Jadbabaie A（2004）使用屏障证书进行安全验证。输入：HSCC·兹比尔1135.93317 [69] Putinar M（1993）紧半代数集上的正多项式。印第安纳大学数学J·Zbl 0796.12002号 [70] 拉马纳，MV；Tunçel，L。；Wolkowicz，H.，半定规划的强对偶性，SIAM J Optim，7641-662，（1997）·Zbl 0891.90129号 ·doi:10.1137/S1052623495288350 [71] Reid G，Wang F，Wolkowicz H，Wu W（2016）arxiv:1504.00931。3月31日访问 [72] Roux P（2016）舍入误差界的形式证明——应用于自动正定性检查。J自动原因·Zbl 1409.68263号 [73] Roux P，Voronin Y-L，Sankaranarayanan S（2016）验证多项式不变量的数值半定规划解。在Xavier R（ed）静态分析第23届国际研讨会上，SAS 2016，英国爱丁堡，2016年9月8日至10日，计算机科学讲稿第9837卷，第424-446页。施普林格·Zbl 1394.68085号 [74] Rump SM（2006）《积极确定性的验证》。BIT数字数学·Zbl 1101.65039号 [75] Sankaranarayanan S，Sipma H，Manna Z（2008）构建混合系统的不变量。形式化方法系统设计·Zbl 1133.68365号 [76] Sankaranarayanan S，Sipma HB，Manna Z（2005）《使用数学规划对线性系统进行可缩放分析》。输入：VMCAI·Zbl 1111.68514号 [77] Schmieta SH，Pataki G（2016）混合半定二次线性规划的DIMACS库的解决方案质量报告。http://dimacs.rutgers.edu/Challenges/Seventh/Instances/error_report.html。[在线；2016年3月23日访问] [78] Sherali HD，Cihan H，Tuncbilek CH（1991）使用重新表述线性化技术的多项式规划的全局优化算法。J全球优化·Zbl 0787.90088号 [79] Shor NZ（1987）多项式函数的全局最小界类。控制论1987。原文为俄语：Kibernetika [80] Shoukry Y、Nuzzo P、Sangiovanni-Vincentelli AL、Seshia SA、Pappas GJ、Tabuada P（2017）SMC：可满足性模凸优化。收录于：Goran F，Sayan M（eds）《第20届混合系统国际会议论文集：计算与控制》，HSCC 2017，匹兹堡，宾夕法尼亚州，美国，2017年4月18日至20日，第19-28页。ACM公司·Zbl 1366.68102号 [81] Sturm，JF，线性矩阵不等式的误差界，SIAM J Optim，101228-1248，（2000）·Zbl 0999.90027号 ·doi:10.1137/S1052623498338606 [82] Sturm JF（1999）使用SeDuMi 1.02，一个用于对称锥体优化的MATLAB工具箱。Optim方法软件 [83] Sturm，JF，混合半定和二阶锥优化问题的内点方法的实现，Optim methods Softw，17，1105-1154，（2002）·Zbl 1032.90021号 ·网址：10.1080/1055678021000045123 [84] Tarski A（1951）初等代数和几何的决策方法。加州大学伯克利分校出版社，技术报告·Zbl 0044.25102号 [85] Tuncel L（2010）组合优化中的多面体和半定规划方法。美国数学学会·Zbl 1207.90005号 [86] TüTüncüRH，Toh KC，Todd MJ（2003）使用SDPT3求解半定二次线性程序。数学编程·Zbl 1030.90082号 [87] 范登伯格，L。；Boyd，S.，《半定规划》，SIAM Rev，38，49-95，（1996）·Zbl 0845.65023号 ·数字对象标识代码：10.1137/1038003 [88] Waki，H。；Muramatsu，M.，《一种用于发现稀疏SOS表示的面部简化算法》，Oper Res Lett，38，361-365，（2010）·Zbl 1205.90219号 ·doi:10.1016/j.orl.2010.05.011 [89] Waki，H。；Muramatsu，M.，二次曲线优化问题的简化算法，《优化理论应用杂志》，158188-215，（2013）·Zbl 1272.90048号 ·doi:10.1007/s10957-012-0219-y [90] Waki H，Nakata M，Muramatsu M（2011）解决多项式优化中半定规划问题的内点方法的奇怪行为。计算优化应用程序·Zbl 1264.90162号 [91] Weispfenning V（1997）实代数的量化器消去——二次型及其后的情形。In：应用代数和纠错码（AAECC）·Zbl 0867.03003号 [92] Wolkowicz，H。；Zhao，Q.，图划分问题的半定规划松弛，离散应用。数学。，96/97, 461-479, (1999) ·Zbl 0932.90030号 ·doi:10.1016/S0166-218X（99）00102-X [93] Wolkowicz H，Saigal R，Vandenberghe L（2000）半定规划手册。多德雷赫特Kluwer学术出版社·Zbl 0962.90001号 ·doi:10.1007/978-1-4615-4381-7 [94] Yamashita M、Fujisawa K、Nakata K、Naka M、Fukuda M、Kobayashi K、Goto K（2010）半定程序的高性能软件包：SDPA 7。技术报告B-460，东京理工大学 [95] Zhao Q，Karisch SE，Rendl F，Wolkowicz H（1998）二次分配问题的半定规划松弛。J Comb Optim 2（1）：71-109。组合优化问题的半定规划和内点方法（多伦多，ON，1996）·Zbl 0904.90145号 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。