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

2018 年度 実績報告書

可積分アルゴリズム:正値性をもつ高精度計算基盤

研究課題

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

研究代表者

中村 佳正  京都大学, 情報学研究科, 教授 (50172458)

研究分担者 木村 欣司  サレジオ工業高等専門学校, その他部局等, 講師 (10447899)
關戸 啓人  京都大学, 国際高等教育院, 特定講師 (40718235)
前田 一貴  関西学院大学, 理工学部, 助教 (80732982)
高田 雅美  奈良女子大学, 生活環境科学系, 講師 (20397574)
研究期間 (年度) 2017-04-01 – 2020-03-31
キーワード可積分アルゴリズム / 片側ヤコビ法 / 両側ヤコビ法 / 直交QD法 / 櫻井・杉浦法
研究実績の概要

櫻井・杉浦法は左特異ベクトルのみを用いて計算を行うため直交QD法と相性が良いが,その一方で, 今倉らにより改良された櫻井・杉浦法に,LAPACKの片側ヤコビ法ではなく,我々の実装した片側ヤコビ法と両側ヤコビ法を搭載した場合の性能評価をGPGPUを利用して行うという研究計画に沿って研究を行った.
ハウスフォルダー2重対角化で前処理しQR法で特異ベクトルを計算するのが従来法である.複数の理由により,NVIDIA Jetson NanoではなくRaspberry pi 4を25台(実験24台+開発1台)購入した.また,ハウスフォルダー2重対角化で前処理し,AIRLB法の下位ルーチンとして,可積分アルゴリズムの一種である直交QD法で特異ベクトルを計算したところ,QR法よりも優れた解法であることを確認できた.AIRLB法では高精度な2重対角化ができるため,直交QD法を生かせた実装が可能である.
一方,2重対角化を経由しない片側ヤコビ法をSST法の下位ルーチンに使ったところ,最大特異値*EPS(計算機イプシロン、10^{-16})より小さい特異値の値であっても信頼できる値を計算できることが多いことがわかった.このため,櫻井・杉浦法には本研究で開発した片側ヤコビ法が最も適していることとなる.一方,両側ヤコビ法は高速性の観点で片側ヤコビ法よりも優れているため,高速性を必要とする特異スペクトル変換法に向いた計算法であることが判明した.
以上の研究で開発した実装コードはライブラリhttp://www-is.amp.i.kyoto-u.ac.jp/LAPROGNC/にて公開している.

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

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

理由

小規模行列の特異値分解を得意とするヤコビ法がある.ヤコビ法には,片側ヤコビ法と両側ヤコビ法があるが,両側ヤコビ法の適切な実装はまだない.本研究では,この問題を解決するために,両側ヤコビ法の3種類の実装方法を開発し提案している.この提案方法を用いる ことで,既存のLAPACKの片側ヤコビ法と比べ,より高精度に特異値分解できることを確認している.開発の過程で,当初の予定と違うアルゴリズムを使う必要が生じ,その分,数値実験が遅れたが,既存のQR法による特異ベクトル計算に対して結果的にはより良好な精度を実現することができた.

今後の研究の推進方策

工学では,最大特異値*EPS(計算機イプシロン、10^{-16})以上の特異値のみが重要で,これはAIRLB法で高精度に計算できる.しかし,これよりも小さい特異値は,
ハウスホルダー変換と同じく,AIRLB法でも高精度計算は困難である.一方,片側ヤコビ法は条件数の大きな行列には強いものの,特に単精度では近接特異値や密集特異値の場合の特異値分解の精度が良くない.そこで,直交QD法の出番となる.このように可積分アルゴリズムである直交QD法のより徹底した研究が待たれる.

  • 研究成果

    (6件)

すべて 2019

すべて 雑誌論文 (3件) (うち査読あり 3件) 学会発表 (2件) 図書 (1件)

  • [雑誌論文] A discrete three-wave system of Kahan-Hirota-Kimura type and the QRT mapping2019

    • 著者名/発表者名
      Yuko Takae, Masami Takata, Kinji Kimura, and Yoshimasa Nakamura
    • 雑誌名

      Proceedings of 2019 International Conference on Parallel and Distributed Processing Techniques and Applications

      巻: 1 ページ: 142-148

    • 査読あり
  • [雑誌論文] Improvement of the thick-restart Lanczos method in single precision floating point arithmetic using Givens rotations2019

    • 著者名/発表者名
      Masana Aoki, Masami Takata, Kinji Kimura, and Yoshimasa Nakamura
    • 雑誌名

      Proceedings of 2019 International Conference on Parallel and Distributed Processing Techniques and Applications

      巻: 1 ページ: 149-155

    • 査読あり
  • [雑誌論文] On an implementation of the two-sded Jacobi method2019

    • 著者名/発表者名
      Sho Araki, Masana Aoki, Masami Takata, Kinji Kimura and Yoshimasa Nakamura
    • 雑誌名

      Proceedings of 2019 International Conference on Parallel and Distributed Processing Techniques and Applications

      巻: 1 ページ: 156‐162

    • 査読あり
  • [学会発表] 融合積和演算を利用した特異値分解のための両側ヤコビ法の実装について2019

    • 著者名/発表者名
      荒木翔, 高田雅美, 木村欣司, 中村佳正
    • 学会等名
      令和元年度応用数理学会年会
  • [学会発表] 割線法を利用した特異値分解のための両側ヤコビ法の実装について2019

    • 著者名/発表者名
      荒木翔,高田雅美,木村欣司,中村佳正
    • 学会等名
      第48 回数値解析シンポジウム
  • [図書] 計算力学レクチャーコース 固有値計算と特異値計算2019

    • 著者名/発表者名
      長谷川,今村,山田,櫻井,荻田,相島,木村,中村共著
    • 総ページ数
      201
    • 出版者
      丸善出版,
    • ISBN
      978-4-621-30473-0

URL: 

公開日: 2021-01-27  

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

Powered by NII kakenhi