• Search Research Projects
  • Search Researchers
  • How to Use
  1. Back to previous page

Integrable Algorithms: Basis of Higher Accurate Computations with Positivity

Research Project

Project/Area Number 17H02858
Research Category

Grant-in-Aid for Scientific Research (B)

Allocation TypeSingle-year Grants
Section一般
Research Field Foundations of mathematics/Applied mathematics
Research InstitutionOsaka Seikei University (2019, 2021)
Kyoto University (2017-2018)

Principal Investigator

Nakamura Yoshimasa  大阪成蹊大学, その他部局等, 教授 (50172458)

Co-Investigator(Kenkyū-buntansha) 木村 欣司  福井大学, 学術研究院工学系部門, 准教授 (10447899)
高田 雅美  奈良女子大学, 生活環境科学系, 講師 (20397574)
關戸 啓人  京都大学, 国際高等教育院, 特定講師 (40718235)
前田 一貴  福知山公立大学, 情報学部, 講師 (80732982)
Project Period (FY) 2017-04-01 – 2020-03-31
Project Status Completed (Fiscal Year 2021)
Budget Amount *help
¥17,810,000 (Direct Cost: ¥13,700,000、Indirect Cost: ¥4,110,000)
Fiscal Year 2019: ¥5,330,000 (Direct Cost: ¥4,100,000、Indirect Cost: ¥1,230,000)
Fiscal Year 2018: ¥5,980,000 (Direct Cost: ¥4,600,000、Indirect Cost: ¥1,380,000)
Fiscal Year 2017: ¥6,500,000 (Direct Cost: ¥5,000,000、Indirect Cost: ¥1,500,000)
Keywords特異値分解 / 高精度線形計算 / 可積分アルゴリズム / 片側ヤコビ法 / 両側ヤコビ法 / 直交QD法 / ランチョス法 / 並列計算 / 原点シフト / 櫻井・杉浦法 / 並列計算の実装 / 相対精度 / 部分特異値分解 / 大規模スパース行列 / ランチョスアルゴリズム / QR分解 / 応用可積分系 / 固有値計算
Outline of Final Research Achievements

We studied integrable algorithms from various points of view. The Jacobi method, which is known as the singular value decomposition method for dense matrices, is more accurate but slower than the standard solution method that combines the Householder transformation and the QR method. The Jacobi method has one side and both sides, and only the one-sided Jacobi method is implemented in LAPACK in the United States. We implemented the two-sided Jacobi method and confirmed that it is faster and more accurate than the LAPACK version of the one-sided Jacobi method. On the other hand, the LAPACK version of the one-sided Jacobi method has the disadvantage that the eigenvectors are not sufficiently orthogonal, although the calculation time for small and medium-sized matrices is not so large. It was shown that high-precision implementation of the one-sided Jacobi method is possible by using high-precision Givens rotation and fused multiply-add (FMA) operation.

Academic Significance and Societal Importance of the Research Achievements

米国の標準数値計算ソフトウェアライブラリのLAPACKでは、特異値分解のための片側ヤコビ法の実装コードが公開され、高精度な特異値分解に資するものとなっているが、両側ヤコビ法の実装コードは見当たらない。本研究では、可積分アルゴリズムとみなすことのできる様々なアルゴリズムの高精度な実装を研究してきたが、最終年度は連続極限が可積分系となることが知られているヤコビ法を取り上げ、両側ヤコビ法の実装コードを開発してLAPACK版の片側ヤコビ法よりも高速かつ高精度であることを確かめている。また、高精度ギブンズ回転と積和演算(FMA)の利用によって片側ヤコビ法の高精度な実装を実現している。

Report

(4 results)
  • 2021 Final Research Report ( PDF )
  • 2019 Annual Research Report
  • 2018 Annual Research Report
  • 2017 Annual Research Report
  • Research Products

    (25 results)

All 2022 2021 2020 2019 2018 2017 Other

All Journal Article (11 results) (of which Peer Reviewed: 11 results,  Open Access: 8 results) Presentation (10 results) (of which Int'l Joint Research: 4 results,  Invited: 1 results) Book (3 results) Remarks (1 results)

  • [Journal Article] Accelerating the numerical computation of positive roots of polynomials using suitable combination of lower bounds2022

    • Author(s)
      Masami Takata, Takuto Akiyama, Sho Araki, Hiroyuki Ishigami, Kinji Kimura and Yoshimasa Nakamura
    • Journal Title

      IPSJ Trans. Math. Modeling Appl.

      Volume: 15 Pages: 18-30

    • Related Report
      2019 Annual Research Report
    • Peer Reviewed / Open Access
  • [Journal Article] On an implementation of the one-sided Jacobi method with high accuracy2021

    • Author(s)
      Masami Takata, Sho Araki, Takahiro Miyamae, Kinji Kimura and Yoshimasa Nakamura
    • Journal Title

      IPSJ Trans. Math. Modeling Appl.

      Volume: 14 Pages: 68-75

    • Related Report
      2019 Annual Research Report
    • Peer Reviewed / Open Access
  • [Journal Article] On an implementation of the two-sided Jacobi method2021

    • Author(s)
      Sho Araki, Masana Aoki, Masami Takata, Kinji Kimura and Yoshimasa Nakamura
    • Journal Title

      IPSJ Trans. Math. Modeling Appl.

      Volume: 14 Pages: 12-20

    • Related Report
      2019 Annual Research Report
    • Peer Reviewed / Open Access
  • [Journal Article] Extension of the Lanczos-Phillips algorithm with Laurent biorthogonal polynomials and its application to the Thron continued fractions2020

    • Author(s)
      Kaoru Sawa and Yoshimasa Nakamura
    • Journal Title

      JSIAM Letters

      Volume: 12 Issue: 0 Pages: 1-4

    • DOI

      10.14495/jsiaml.12.1

    • NAID

      130007786409

    • ISSN
      1883-0609, 1883-0617
    • Related Report
      2019 Annual Research Report
    • Peer Reviewed / Open Access
  • [Journal Article] A discrete three-wave system of Kahan-Hirota-Kimura type and the QRT mapping2019

    • Author(s)
      Yuko Takae, Masami Takata, Kinji Kimura, and Yoshimasa Nakamura
    • Journal Title

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

      Volume: 1 Pages: 142-148

    • Related Report
      2018 Annual Research Report
    • Peer Reviewed
  • [Journal Article] Improvement of the thick-restart Lanczos method in single precision floating point arithmetic using Givens rotations2019

    • Author(s)
      Masana Aoki, Masami Takata, Kinji Kimura, and Yoshimasa Nakamura
    • Journal Title

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

      Volume: 1 Pages: 149-155

    • Related Report
      2018 Annual Research Report
    • Peer Reviewed
  • [Journal Article] On an implementation of the two-sded Jacobi method2019

    • Author(s)
      Sho Araki, Masana Aoki, Masami Takata, Kinji Kimura and Yoshimasa Nakamura
    • Journal Title

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

      Volume: 1

    • Related Report
      2018 Annual Research Report
    • Peer Reviewed
  • [Journal Article] Another generalization of the box ball system with many kinds of balls2018

    • Author(s)
      Maeda Kazuki
    • Journal Title

      Journal of Integrable Systems

      Volume: 3 Issue: 1 Pages: 5033954-5033954

    • DOI

      10.1093/integr/xyy007

    • Related Report
      2017 Annual Research Report
    • Peer Reviewed / Open Access
  • [Journal Article] Implementation of computing partial singular value decomposition for proncipal component analysis using ARPACK2018

    • Author(s)
      Masami Takata, Sho Araki, Kinji Kimura, Yuki Fujii, Yoshimasa Nakamura
    • Journal Title

      IPSJ Trans. Math. Modeling Appl.

      Volume: 11-1 Pages: 37-44

    • Related Report
      2017 Annual Research Report
    • Peer Reviewed / Open Access
  • [Journal Article] Improvement of the augmented implicitly restarted Lanczos bidiagonalization method in single precision floating point arithmetic2018

    • Author(s)
      Yuya Ishida, Masami Takata, Kinji Kimura, Yoshimasa Nakamura
    • Journal Title

      IPSJ Trans. Math. Modeling Appl.

      Volume: 11-3 Pages: 19-25

    • Related Report
      2017 Annual Research Report
    • Peer Reviewed / Open Access
  • [Journal Article] Application of the Lanczos-Phillips algorithm to continued fractions and its extension with orthogonal polynomials2018

    • Author(s)
      Kaoru Sawa and Yoshimasa Nakamura,
    • Journal Title

      JSIAM Letters

      Volume: 10 Issue: 0 Pages: 57-60

    • DOI

      10.14495/jsiaml.10.57

    • ISSN
      1883-0609, 1883-0617
    • Related Report
      2017 Annual Research Report
    • Peer Reviewed / Open Access
  • [Presentation] 融合積和演算を利用した特異値分解のための両側ヤコビ法の実装について2019

    • Author(s)
      荒木翔, 高田雅美, 木村欣司, 中村佳正
    • Organizer
      令和元年度応用数理学会年会
    • Related Report
      2018 Annual Research Report
  • [Presentation] 割線法を利用した特異値分解のための両側ヤコビ法の実装について2019

    • Author(s)
      荒木翔,高田雅美,木村欣司,中村佳正
    • Organizer
      第48 回数値解析シンポジウム
    • Related Report
      2018 Annual Research Report
  • [Presentation] Kahan-廣田-木村型離散three wave system とQRT写像2018

    • Author(s)
      高江宥光, 木村欣司, 中村佳正
    • Organizer
      平成30 年度応用数理学会年会
    • Related Report
      2017 Annual Research Report
  • [Presentation] Fast computation method of column space by using the DQDS method and the OQDS method2018

    • Author(s)
      Sho Araki, Hiroki Tanaka, Masami Takata, Kinji Kimura and Yoshimasa Nakamura
    • Organizer
      The 2018 International Conference on Parallel and Distributed Processing Techniques and Applications (PDPTA2018)
    • Related Report
      2017 Annual Research Report
    • Int'l Joint Research
  • [Presentation] Improvement of computing partial singular value decomposition for dense matrix using thick-restart-Lanczos method2018

    • Author(s)
      Masana Aoki, Masami Takata, Kinji Kimura and Yoshimasa Nakamura
    • Organizer
      The 2018 International Conference on Parallel and Distributed Processing Techniques and Applications (PDPTA2018)
    • Related Report
      2017 Annual Research Report
    • Int'l Joint Research
  • [Presentation] ランチョス・フィリップスアルゴリズムの可積分性とその直交多項式への応用2018

    • Author(s)
      澤薫, 中村佳正
    • Organizer
      第14 回応用数理学会研究部会連合発表会
    • Related Report
      2017 Annual Research Report
  • [Presentation] 単精度浮動小数点数環境にも有効なAugmented Implicitly Restarted Lanczos Bidiagonalization 法の改良2018

    • Author(s)
      石田遊也, 高田雅美, 木村欣司, 中村佳正
    • Organizer
      第14 回応用数理学会研究部会連合発表会,
    • Related Report
      2017 Annual Research Report
  • [Presentation] Parallel symmetric tridiagonal eigensolver based on bisection and inverse iteration algorithms with the blocked classical Gram-Schmidt reorthogonalization and the shifted Cholesky QR decomposition2018

    • Author(s)
      Kinji Kimura, Masayuki Osawa, Sho Araki, Hiroyuki Ishigami, Yoshimasa Nakamura
    • Organizer
      SIAM Conference on Parallel Processing for Scientific Computing
    • Related Report
      2017 Annual Research Report
    • Int'l Joint Research
  • [Presentation] Acceleration of inverse iteration method using Cholesky QR decomposition with shifts2018

    • Author(s)
      Masayuki Osawa, Sho Araki, Hiroyuki Ishigami, Kinji Kimura, Yoshimasa Nakamura
    • Organizer
      International Workshop on Eigenvalue Problems: Algorithms; Software and Applications, in Petascale Computing (EPASA2018)
    • Related Report
      2017 Annual Research Report
    • Int'l Joint Research
  • [Presentation] 応用可積分系分野の基礎方程式としての戸田格子方程式,2017

    • Author(s)
      中村佳正
    • Organizer
      平成29 年度応用数理学会年会
    • Related Report
      2017 Annual Research Report
    • Invited
  • [Book] 計算力学レクチャーコース 固有値計算と特異値計算2019

    • Author(s)
      長谷川,今村,山田,櫻井,荻田,相島,木村,中村共著
    • Total Pages
      201
    • Publisher
      丸善出版,
    • ISBN
      9784621304730
    • Related Report
      2018 Annual Research Report
  • [Book] 解析学百科:可積分系の数理2018

    • Author(s)
      中村,高崎,辻本,尾角,井ノ口(第5章応用可積分系/執筆分担)
    • Total Pages
      432
    • Publisher
      朝倉書店
    • ISBN
      9784254117271
    • Related Report
      2017 Annual Research Report
  • [Book] 数値線形代数の数理とHPC2018

    • Author(s)
      櫻井, 松尾, 片桐編(2.2.2 mdLVs法/執筆分担)
    • Total Pages
      314
    • Publisher
      共立出版
    • ISBN
      9784320019553
    • Related Report
      2017 Annual Research Report
  • [Remarks] LAPROGNC (Lin. Alg. PROGrams in Nume. Comput.)

    • URL

      http://syskiso.fuee.u-fukui.ac.jp/~kkimur/LAPROGNC/LAPROGNC-j.html

    • Related Report
      2019 Annual Research Report

URL: 

Published: 2017-04-28   Modified: 2023-01-30  

Information User Guide FAQ News Terms of Use Attribution of KAKENHI

Powered by NII kakenhi