大规模问题に対する精度保証付き数値计算に成功

掲载日2023.06.09

概要

岩手大学理工学部物理?材料理工学科数理?物理コース宫岛信也教授は、サイズのとても大きな行列の実数乗とベクトルとの積を含む区間をコンピュータ上で高速に求める手法を確立しました。ただし、扱う行列は、応用上よく現れるある構造をもっていることを仮定しています。宮島教授は、その構造を利用することで、必要となる演算回数と記憶容量を大幅に削減するための定理とその証明を与えました。
本成果により、精度保証付き数値計算法が抱える課題の１つである大規模問題への適用の困難性を解決するための新たな方向性が見いだされました。本研究成果は、令和５年９月１５日にエルゼビアが発刊する著名な国際学術誌Applied Mathematics and Computationより全世界へ公開されます（本研究論文は、エルゼビアが運営するページScience Directにおいて、先行公開されています）。

背景

コンピュータは计算を间违えることがあります。この原因の１つはコンピュータの中で数字を记忆するのに使用する桁数には限りがあることです。「それでは、桁数を増やせばいいのでは？」と思うかもしれませんが、必ずしもそれで解决できる訳ではありません。実际、次のような例が存在します。

「８桁记忆できる环境で、ある値を计算すると１.１７２６...となった。次に、１７桁记忆できる环境で计算しても结果は１.１７２６...となった。さらに、３４桁记忆できる环境で计算しても结果は１.１７２６...となった。しかし、正しい値は－０.８２７３...である。」

この例は「ある桁数で计算した结果とそれよりも长い桁数で计算した结果が一致していれば、その结果は正しい」とは限らないことを示しています。

この问题に対処するための１つのアプローチとして、正しい値を含む区间をコンピュータを使って计算するための手法が研究されてきました。このような区间が得られれば、その区间の中心は正しい値の近似、半径はその近似が取りうるズレ（误差）の上限とみなすことができます。よって、半径が小さければ、中心は高信頼な近似であると断言できます。また、このような区间を利用して、これまで証明できなかった数学上の命题が証明できるようになりました。このような区间を求める手法を精度保証付き数値计算法といいます。

精度保証付き数値计算が抱える课题の１つは大规模问题に适用するのが难しいことです。（区间を求めない）普通の计算では、演算回数や记忆容量を节约するための様々な工夫がこれまでに考案されてきました。しかし、これらの工夫を精度保証付き数値计算に流用できるかというと、いつでも流用できる訳ではありません。その结果、普通の计算で扱える问题の规模と精度保証付き数値计算で扱える问题の规模には、大きな隔たりがあります。

研究成果

本論文では、これまでの精度保証付き数値計算では扱うことができなかった、サイズのとても大きな行列の実数乗とベクトルとの積を含む区間を高速に求める手法を確立しました。ただし、扱う行列は、応用上よく現れるある構造をもっていることを仮定しています。この構造をもった行列の実数乗とベクトルとの積は、例えば非整数階偏微分方程式の解を求める際に必要となります。宫岛信也教授は、その構造を利用することで、必要となる演算回数と記憶容量を大幅に削減するための定理とその証明を与えました。その結果、あるパソコンにおいて７０分かかっていた計算が０.４９秒で済むようになりました。また、あるパソコンではこれまでメモリ不足のため計算できなかった区間が計算できるようになりました。これらの成果はエルゼビアが発刊する著名な国際学術誌Applied Mathematics and Computationに掲載されました。本論文の成果は上述の精度保証付き数値計算が解決すべき課題に光を当てています。

掲载论文

題　目： Fast verified computation for real powers of large matrices with Kronecker structure
著　者： Shinya Miyajima
誌　名： Applied Mathematics and Computation
公表日： 15 September 2023

本研究は、以下の研究事业の成果の一部として得られました。
?文部科学省科学研究費補助金?基盤研究(C)「非整数階微分方程式系の解に対する精度保証付き数値計算法の研究」研究代表者：宫岛信也

本件に関する问い合わせ先

理工学部物理?材料理工学科数理?物理コース宫岛信也
019-621-6303
miyajima@iwate-u.ac.jp

猫咪视频