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

実空間密度汎関数計算における超並列計算環境向け固有値解法の研究

Research Project

Project/Area Number 12J02480
Research Category

Grant-in-Aid for JSPS Fellows

Allocation TypeSingle-year Grants
Section国内
Research Field Engineering fundamentals
Research InstitutionUniversity of Tsukuba

Principal Investigator

二村 保徳  筑波大学, システム情報工学研究科, 特別研究員(DC2)

Project Period (FY) 2012 – 2013
Project Status Completed (Fiscal Year 2013)
Budget Amount *help
¥1,800,000 (Direct Cost: ¥1,800,000)
Fiscal Year 2013: ¥900,000 (Direct Cost: ¥900,000)
Fiscal Year 2012: ¥900,000 (Direct Cost: ¥900,000)
Keywords固有値問題 / 電子状態計算 / 密度汎関数法
Research Abstract

密度汎関数による電子状態計算におけるSelf consistent field (SCF)計算やバンド曲線計算において, 櫻井・杉浦法で複数の閉曲線を設置して固有値計算を行う際, 事前に固有値分布がわかっているとパラメータを効果的に設定することができる. 当特別研究員と共同研究者らはそのような状況に適した固有値分布の確率的推定法を2010年に提案したが, 第一年度に発表したshifted block conjugate gradient rQ法と本報告書記載の第1の雑誌論文でまとめられたBlock bilinear form近似法の考え方を組み合わせることで, この固有値分布確率的推定法の計算を高速に実行できるアルゴリズムを新たに開発した. このアルゴリズムを利用することで, 電子状態計算で現れる行列に対する固有値分布の確率的推定法の計算が従来の固有値を高精度に求める解法に比べ約1/27の計算時間で行うことができることを示した. この結果より提案手法が固有値解法の前処理として十分小さいコストで用いることができるといえる. この研究の成果は当特別研究員の博士論文において発表された. また, 当特別研究員らは櫻井・杉浦法の1万原子規模のシリコンナノワイヤのSCF計算で現れる標準固有値問題に対する性能を確かあるため, 京コンピュータを用いて数値実験を行った. この計算は約900万次元の疎行列の約2万の固有対を求める大規模な問題である. この実験の結果, 櫻井・杉浦法を並列性能は高いものの, 内部で用いたSBC GrQ法による線形方程式の求解に多くに反復回数を要するため, 従来解法である固有値問題向け共役勾配法に対し, 計算時間が多くなることが確認された.

Current Status of Research Progress
Current Status of Research Progress

1: Research has progressed more than it was originally planned.

Reason

当特別研究員は研究テーマである実空間密度汎関数計算の固有値問題の並列解法について着実に研究成果を出し, 計画以上の進展があった. 研究は交付申請書の研究の目的に記載されている通りに順調に進んだ. また, 研究活動において得られた成果をまとめ, 論文を投稿し, 厳正な査読を通じて学術雑誌に採録された. また研究成果を国際会議および国内会議において積極的に発表した. 特に国際会議3PGCIC2013においては研究成果が高く評価されたことによりBest Paper Awardを受賞した.

Strategy for Future Research Activity

研究実績の概要の最後で述べた計算においては並列資源を大量に用いた場合, 櫻井・杉浦法が共役勾配法より速く計算できる可能性があるが, そのために必要な並列資源の量が現実的でないため, さらなる高速化手法の開発が必要であるということが確認できた. このような結果を踏まえ, 今後当特別研究員は櫻井・杉浦法のさらなる高速化のため道筋を探っていく予定である.

Report

(2 results)
  • 2013 Annual Research Report
  • 2012 Annual Research Report
  • Research Products

    (13 results)

All 2014 2013 2012

All Journal Article (4 results) (of which Peer Reviewed: 4 results) Presentation (9 results)

  • [Journal Article] Block conjugate gradient type methods for the approximation of bilinear form $C^{H} A^{1} B$2014

    • Author(s)
      L. Du, Y. Futamura and T. Sak urai
    • Journal Title

      Computers & Mathematics with Applications

      Volume: 66(12) Issue: 12 Pages: 2446-2455

    • DOI

      10.1016/j.camwa.2013.09.023

    • Related Report
      2013 Annual Research Report
    • Peer Reviewed
  • [Journal Article] Numerical Construction of a Low-Energy Effective Hamiltonian in a Self-Consistent Bogohubov-de Gennes Approach of Superconduictivity2013

    • Author(s)
      Y. Nagai, Y. Shinohara, Y. Futamura. Y. Ota and T. Sakurai
    • Journal Title

      J. Phys. Soc. Jpn.

      Volume: 82(9)

    • Related Report
      2013 Annual Research Report
    • Peer Reviewed
  • [Journal Article] Multi-GPU scalable implementation of a contour-integral-based eigensoler for real symmetric dense generalized eigenvalue problems2013

    • Author(s)
      Y. Yano, Y. Futamura and T. Sakurai
    • Journal Title

      In Proc. 8^<TH> International Conference on P2P, Par allel, Cloud and Internet Computing

      Pages: 121-127

    • DOI

      10.1109/3pgcic.2013.24

    • Related Report
      2013 Annual Research Report
    • Peer Reviewed
  • [Journal Article] 非線形固有値問題の固有値密度推定法における適応的並列アルゴリズム2012

    • Author(s)
      山本和磨, 前田恭行, 二村保徳, 櫻井鉄也
    • Journal Title

      情報処理学会論文誌コンピューティングシステム

      Volume: 5(3) Pages: 22-29

    • NAID

      40019469241

    • Related Report
      2012 Annual Research Report
    • Peer Reviewed
  • [Presentation] A Hierarchical Parallel Software Package of a Complex Moment Based Eigensolver2014

    • Author(s)
      Yasunori Futamura and Tetsuya Sakurai
    • Organizer
      SIAM PP14
    • Place of Presentation
      Marriott Portland Downtown Waterfront Portland
    • Year and Date
      2014-02-20
    • Related Report
      2013 Annual Research Report
  • [Presentation] 周回積分を用いた固有値解法に基づく並列固有値解析エンジン2013

    • Author(s)
      二村保徳, 櫻井鉄也
    • Organizer
      日本応用数理学会
    • Place of Presentation
      アクロス福岡
    • Year and Date
      2013-09-09
    • Related Report
      2013 Annual Research Report
  • [Presentation] 周回積分を用いた固有値解法に基づく並列固有値解析エンジンの実装2013

    • Author(s)
      二村保徳, 櫻井鉄也
    • Organizer
      SWoPP北九州2013
    • Place of Presentation
      北九州国際会議場
    • Year and Date
      2013-08-01
    • Related Report
      2013 Annual Research Report
  • [Presentation] Conjugate gradient type method for shifted linear systems with multiple right hand sides2013

    • Author(s)
      Yasunori Futamura and Tetsuya Sakurai
    • Organizer
      EASIAM 2013
    • Place of Presentation
      The Newton hotel Bandung
    • Year and Date
      2013-06-18
    • Related Report
      2013 Annual Research Report
  • [Presentation] Density Functional Electronic Band Structure Calculations with a Complex Moment Based Eigensolver2013

    • Author(s)
      Yasunori Futamura, Tetsuya Sakurai, Shinnosuke Furuya, Jun-Ichi Iwata
    • Organizer
      SIAM Conference on Computational Science and Engineering
    • Place of Presentation
      Boston, MA, USA
    • Year and Date
      2013-02-25
    • Related Report
      2012 Annual Research Report
  • [Presentation] 櫻井・杉浦法におけるshifted Krylov部分空間反復法の応用2012

    • Author(s)
      二村保徳, 櫻井鉄也
    • Organizer
      日本応用数理学会「行列・固有値問題の解法とその応用」研究部会 第14回研究会
    • Place of Presentation
      筑波大学
    • Year and Date
      2012-11-20
    • Related Report
      2012 Annual Research Report
  • [Presentation] 周回積分を用いた固有値解法における線形方程式求解アルゴリズムとその電子状態計算への応用2012

    • Author(s)
      二村保徳, 櫻井鉄也, 古家真之介, 岩田潤一
    • Organizer
      日本応用数理学会年会
    • Place of Presentation
      稚内全日空ホテル
    • Year and Date
      2012-08-30
    • Related Report
      2012 Annual Research Report
  • [Presentation] Efficient Algorithm for Linear Systems Arising in Solutions of Eigenprobleis and its Application to Electronic-Structure Calculations2012

    • Author(s)
      Yasunori Futamura, Tetsuya Sakurai, Shinnosuke Furuya, Jun-Ichi Iwata
    • Organizer
      10th International Meeting on High-Perforiance Computing for Computational Science (VECPAR2012)
    • Place of Presentation
      RIKEN AICS, Kobe, Japan
    • Year and Date
      2012-07-20
    • Related Report
      2012 Annual Research Report
  • [Presentation] Efficient algorithm for solving linear systems arising from a sparse eigensolver2012

    • Author(s)
      Yasunori Futamura, Tetsuya Sakurai, Shinnosuke Furuya, Jun-Ichi Iwata
    • Organizer
      7th International Workshop on Parallel Matrix Algorithms and Applications
    • Place of Presentation
      Birkbeck University of London
    • Year and Date
      2012-06-30
    • Related Report
      2012 Annual Research Report

URL: 

Published: 2013-04-25   Modified: 2024-03-26  

Information User Guide FAQ News Terms of Use Attribution of KAKENHI

Powered by NII kakenhi