1998 Fiscal Year Annual Research Report
Project/Area Number |
09680327
|
Research Institution | The University of Tokyo |
Principal Investigator |
小柳 義夫 東京大学, 大学院・理学系研究科, 教授 (60011673)
|
Co-Investigator(Kenkyū-buntansha) |
西田 晃 東京大学, 大学院・理学系研究科, 助手 (60302808)
須田 礼仁 名古屋大学, 大学院工学研究科, 講師 (40251392)
|
Keywords | 偏微分方程式 / 離散化 / 並列処理 / 前処理 / 共役勾配法 / ICCG法 / Jocobi法 / 多重積分 |
Research Abstract |
本研究では、偏微分方程式を不規則な構造をもった格子で離散化した線形方程式の解法を研究し、併せてその並列計算機上への実装を研究した。 移流項のない偏微分方程式を離散化して得られる,疎で対称な大規模連立一次方程式に対しては,前処理つき共役勾配法(Preconditioned CG法PCG法)が有効な解法として知られている.逐次計算では,問題に依存した性質を利用しない方法の中で最も速い前処理の一つとして,不完全コレスキー分解を使ったPCG法(ICCG法)があり,広く用いられているが,並列計算における汎用の前処理について,逐次計算でのICCG法ほど決定的なものはないといってよい.本研究では,逐次計算では有効ではないが並列性の高い,前処理としての定常的な反復法に注目した.そして,線Jacobi法を変形した定常的な反復法から導かれる,行列形式のみを利用する並列性の高い前処理を提案,実装した.さらに,数値実験により,点Jacobi前処理と比較して収束が速いことを示した. 領域を不規則な構造をもった格子で離散化する際の性質を分析するために、与えられた節点についての種々の三角形分割による多重積分の誤差について分析した。単項式の線形結合である関数族全体に対する総合的な誤差評価として、BarnhillとDavisの方法によって誤差評価をする実験を行なった結果,誤差を最小にするような最適な三角形分割法は,実用上よく使われるDelaunayの三角形分割法に必ず一致するとは限らないが,あまり積分公式に影響されない,何らかの良い性質をもったものであることが示唆された.
|
Research Products
(6 results)
-
[Publications] R.Suda,C.Y.Hsieh,Y.Oyauagi: "The new self test Mont eCarlo method" Proc.High Performance Compating Conference. vol.2. 1277-1284 (1998)
-
[Publications] A.Moriyama and Y.Oyanagi: "Triangulation and error of numerical integration" Proc.of the 27th Numerical Analpis symposium. 83-86 (1998)
-
[Publications] A.Nishida: "Least Squares Arnoldi for Large Nonspmmetric Eigenproblems" Proc.of 5th Copper Mountain Corberence on gterative metholy. vol.2. (1998)
-
[Publications] T.Sekine and Y.Oyanagi: "Quari-Monte Carlo method on Ising model simulation" Proc.of the 27th Numecical Snalysis Lepnpoium. 2-5 (1998)
-
[Publications] 野田敏達・小柳義夫: "領域分割法とその並列化" 京都大学数理解析研究所講究録. Vol.-1040. 196-203 (1998)
-
[Publications] Y.Hamana and Y.Oyanagi: "短距離分子動力学法の並列化" IPSJ SIG Notes. 98(72). 25-30 (1998)