• 研究課題をさがす
  • 研究者をさがす
  • KAKENの使い方
  1. 課題ページに戻る

2011 年度 実施状況報告書

正則化法による逆問題の高精度近似理論と次世代計算環境による数値的実現

研究課題

研究課題/領域番号 23740075
研究機関京都大学

研究代表者

藤原 宏志  京都大学, 情報学研究科, 助教 (00362583)

研究期間 (年度) 2011-04-28 – 2014-03-31
キーワード多倍長計算 / 正則化法 / 数値的不安定性 / 非適切問題 / 高精度シミュレーション
研究概要

本研究は、非適切問題の離散化で現れる数値的に不安定な問題に対し、多倍長数値計算を適用することにより高精度数値計算を目指すものである。特に、正則化法に代表される安定化手法の近似誤差、離散化の誤差、および丸め誤差のそれぞれの影響に適切に対処することによる高精度数値計算手法の確立を目指すものである。初年度は、非適切問題の数値計算における数値的不安定性を定量的に調べるための指標を提案した。この指標は、数値的に不安定であるほど丸め誤差の増大が大きくなることに着目し、多倍長計算において計算精度が可変であるという特徴をもちいて事後誤差解析として算出するものである。この指標が実際に利用される数値計算アルゴリズムに適用可能であることを示すために、複素逆 Laplace 変換の数値計算手法のうち細野の方法と Sheen の方法に適用し、それぞれの数値計算における不安定性を調べた。その結果、いずれも数学的には Bromwich の方法に対する数値計算アルゴリズムでありながら、Sheen の方法においては逆変換像を求める計算点の値が増大するごとに不安定性が増大することが示されたのに対し、細野の方法においては、不安定性は計算点の値に依存せず、ほぼ一定であることが示された。これは細野の方法で採用されている安定化手法である核近似の特徴であると考えられ、引き続き、この核近似の意味について調べることが重要であると考えられる。また、眞鍋秀悟氏(申請者が研究指導を務る修士大学院生)との共同研究により、近年並列計算の分野で利用される画面描画用途のプロセサ GPU をもちい、その並列演算性能の高さに着目した多倍長精度のベクトル演算ライブラリの構築をおこなった。

現在までの達成度 (区分)
現在までの達成度 (区分)

2: おおむね順調に進展している

理由

本研究の目的は、多倍長計算環境を利用し、非適切問題の高精度数値計算手法を確立するものである。そのために理論的な側面および多倍長計算環境の構築という計算科学の側面を合わせた研究である。初年度は、当初の計画に従って理論的な側面の研究として、典型的な非適切問題である複素逆 Laplace 変換に対し、安定化手法を適用する高精度数値計算を実現した。また、この複素逆 Laplace 変換を例として、数値的不安定性を定量的に評価するための指標を提案し、具体的に2つの複素 Laplace逆変換手法の特性を明らかにしたが、その原因の究明には至らなかったため、引き続き検討が必要であると考えられる。一方、これと合わせて、研究2年目に予定していた多倍長計算の高速化についても研究を開始し、大学院生との共同研究により GPU を利用した並列多倍長計算環境の構築をおこなった。この多倍長計算環境は固定小数点演算に限られているものの、科学技術計算に頻繁に現れるベクトル演算の高速化を実現するものである。得られた環境を複素逆 Laplace 変換に適用した結果、簡単な例においても既存の多倍長計算環境の10倍程度の高速化が達成され、当初の予定を上回る進展があった。以上を合わせて、全体としておおむね順調に進展していると判断した。

今後の研究の推進方策

研究2年目にあたる24年度は、多倍長計算環境の高速化を目的として、最新のプロセサを活用した高速多倍長計算環境の構築をおこなう。23年度において、申請者を指導教員とする大学院生との共同研究により、画面描画用途のプロセサ GPU をもちいた高速多倍長計算環境を構築した。これは固定小数点演算のベクトル計算に特化したものであり、上述の複素逆Laplace変換においては既存の多倍長計算環境の10倍程度の高速化を達成した一方で、一般的な科学技術計算で必要となる浮動小数点演算へは対応しなかった。これは既存のGPUが分岐命令に不向きである特性 (divergence branch) を考慮してのことである。24年度には GPU に限定することなく、最新プロセサを利用した多倍長計算の高速化に取り組む。特に Intel 社が近年発表した SIMD演算命令セット AVX は、命令レベル並列での高速多倍長計算環境の実装に向くと考えており、その利用を試みる。ただし、これが多倍長計算に不向きな場合には、MPI などのプロセスレベルでの並列計算をもちいた高速化を試みる。また、理論的な側面においては、初年度に提案した数値的安定性の指標の特徴の検討も続ける。一般的に函数方程式で記述される非適切問題の離散問題の数値的不安定性は、誤差の高周波成分が急激に増大することに起因すると考えられている。一方、本年度に提案した指標は桁落ち現象に起因する誤差の増大を捉えていることも特徴であり、これら両者の関係を明らかにすることは、非適切問題の数値計算を阻害する原因について新たな知見を与え、高精度数値計算に向けて重要なステップであると考えられる。

次年度の研究費の使用計画

23年度中に新たな計算機を購入するための予算を計上していたが、既存の計算環境のみでの複素逆Laplace変換の研究が進展したため、新規の計算機の購入を見送った。これに対して、研究の2年度は多倍長計算環境の高速化が主たる研究の対象となるため、新たな計算機の購入を予定している。特に24年度初めに Intel 社から新たなプロセサ Ivy Bridge (The Third Generation Intel Core Processor) が発表されている。このプロセサでは SIMD 演算のための新しい命令セットが提供されており、加えて GPU が統合されていることから、科学技術計算を目的とする本研究の多倍長計算環境の高速化を試みるのに最適なプロセサであると判断し、その購入を計画している。また、6月には台湾で応用数学の東アジアでの国際研究集会 EASIM が開催され、11月には台湾大学において逆問題に関するワークショップが開催される予定であり、これらに参加して研究発表および討論をおこなうための海外出張を予定している。また、従来より共同研究を進めている斎藤三郎教授(アヴェイロ大)との打ち合わせのため一週間程度の海外出張も予定している。

  • 研究成果

    (12件)

すべて 2012 2011 その他

すべて 雑誌論文 (6件) (うち査読あり 2件) 学会発表 (5件) 備考 (1件)

  • [雑誌論文] 非適切問題の数値計算における丸め誤差の評価について2012

    • 著者名/発表者名
      藤原宏志
    • 雑誌名

      第61回理論応用力学講演会講演論文集

  • [雑誌論文] 3次元Voigt型粘弾性方程式の数値計算の試み2012

    • 著者名/発表者名
      藤原宏志
    • 雑誌名

      第61回理論応用力学講演会講演論文集

      巻: 巻号なし ページ: 電子版のため頁番号なし

  • [雑誌論文] 多倍長計算による高精度計算の実現に向けて2011

    • 著者名/発表者名
      藤原宏志
    • 雑誌名

      応用数理

      巻: Vol.21 ページ: 247-258

  • [雑誌論文] GPU上の多倍長数値計算環境による複素逆Laplace変換の高速化2011

    • 著者名/発表者名
      眞鍋秀悟、藤原宏志
    • 雑誌名

      計算数理工学論文集

      巻: 11 ページ: 31-34

    • 査読あり
  • [雑誌論文] Approximate Steady State Models for Magnetic Resonance Elastography2011

    • 著者名/発表者名
      Y. Jiang, H. Fujiwara, and G. Nakamura
    • 雑誌名

      SIAM J. on Applied Mathematics

      巻: 71 ページ: 1965-1989

    • 査読あり
  • [雑誌論文] High-Accurate Computation of One-Loop Integrals by Several Hundred Digits Multiple-Precision Arithmetic2011

    • 著者名/発表者名
      H. Fujiwara
    • 雑誌名

      Proceedings of 3rd Computational Particle Physics Workshop

      巻: PoS(CPP2010)013 ページ: オンラインのためページなし

  • [学会発表] 非適切問題の数値計算における丸め誤差の評価について2012

    • 著者名/発表者名
      藤原宏志
    • 学会等名
      第61回理論応用力学講演会
    • 発表場所
      東京大学生産技術研究所
    • 年月日
      2012年3月8日
  • [学会発表] 3次元Voigt型粘弾性方程式の数値計算の試み2012

    • 著者名/発表者名
      藤原宏志
    • 学会等名
      第61回理論応用力学講演会
    • 発表場所
      東京大学生産技術研究所
    • 年月日
      2012年3月7日
  • [学会発表] Remarks on Numerical Instability of Complex Inverse Laplace Transforms Using Multiple-Precision Arithmetic2011

    • 著者名/発表者名
      H. Fujiwara
    • 学会等名
      The 7th EASIAM Conference 2011
    • 発表場所
      Waseda University, Kitakyushu Campus
    • 年月日
      Jun. 29, 2011
  • [学会発表] 光トモグラフィーの基礎 - 次代を拓くためのアプローチ2011

    • 著者名/発表者名
      藤原宏志
    • 学会等名
      日本ヒト脳機能マッピング学会(招待講演)
    • 発表場所
      国立京都国際会館
    • 年月日
      2011年9月1日
  • [学会発表] Direct treatment of numerical instability of inverse Laplace transforms using multiple-precision arithmetic2011

    • 著者名/発表者名
      H. Fujiwara
    • 学会等名
      Laplace Transform Methods and Thier Applications(招待講演)
    • 発表場所
      National Institute for Mathematical Science, Daejeon
    • 年月日
      2011, Nov. 4
  • [備考] 実装を進めている多倍長計算環境は、以下の URL にて公開している。

    • URL

      http://www-an.acs.i.kyoto-u.ac.jp/~fujiwara/exflib

URL: 

公開日: 2013-07-10  

サービス概要 検索マニュアル よくある質問 お知らせ 利用規程 科研費による研究の帰属

Powered by NII kakenhi