簡易檢索 / 詳目顯示

研究生: 陳俞廷
Yu-Ting Chen
論文名稱: 以周道積分解三維光子晶體之廣義特徵值問題
Solve three dimensional photonic crystal general eigenvalue problem by contour integral
指導教授: 黃聰明
Huang, Tsung-Ming
學位類別: 碩士
Master
系所名稱: 數學系
Department of Mathematics
論文出版年: 2014
畢業學年度: 102
語文別: 中文
論文頁數: 36
中文關鍵詞: 特徵值問題週道積分三維光子晶體平行運算
英文關鍵詞: general eigenvalue problem, contour integral, three dimensional photonic crystal, parallel computing
論文種類: 學術論文
相關次數: 點閱:101下載:8
分享至:
查詢本校圖書館目錄 查詢臺灣博碩士論文知識加值系統 勘誤回報
  • 在本文中,主要測試周道積分法(Contour integral)\cite {S2}的效能,測試並歸納出一些較好的使用原則。我們以此方法針對三維光子晶體的馬克斯威爾方程(Maxwell equation)取時諧波(time harmonic)後經過K. Yee的離散化程序\cite{T1}得到的廣義特徵值問題作效能的測試,觀察各參數間的交互作用,並且將其與E. Polizzi的FEAST\cite{E1,E2,E3}結合而提出了混合型的演算法。最後對這混合型的演算法與傳統的Lanczos作效能上的比較。文中的三維光子晶體馬克斯威爾方程包含兩種情形分別為簡易立方晶格(SC,Simple Cube)和面心立方晶格(FCC,Face Centered cube),並分為是否除去零空間(null space free)的兩種情況,於理論部分僅提供較簡易的簡易立方晶格的介紹,而我們的數值結果則著重於應用較廣的面心立方晶格。

    In this papper, we consider to solve general eigenvalue problem for three dimensional photonic crystal by contour integral, and focus on the solver's efficacy. At first, we take the time harmonic for three dimensional photonic crystal's Maxwell equation , and Discrete by Yee's scheme,then test the parameter for the sovler. We explain the implications of parameter for CIRR,and compare it with FEAST.After all, We propose a hybrid solver MLCIRR, it Combine CIRR and FEAST.

    導論 1 2 周道積分應用至瑞雷 -瑞茲投影特徵值解法 2 2.1 CIRR 的核心定理2 2.2 FEAST 演算法 4 2.3 CIRR 各項參數的意涵與作用6 2.3.1 M 7 2.3.2 L 8 2.3.3 N 9 2.3.4 積分路徑相關參數 9 2.4 混合型的 CIRR 10 3 馬克斯威爾方程與建構的矩陣 12 3.1 從馬克斯威爾方程到離散出的矩陣 14 3.2 以快速傅立葉變換得到的預處理矩陣 16 3.3 轉化為沒有零空間的問題 19 4 周道積分法與三為光子晶體的馬克斯威爾方程 22 4.1 Multi-level 的引入 23 4.2 線性系統與平移值之關聯 24 5 數值結果 25 5.1 線性系統迭代代價關係 25 5.2 移去零空間的必要性 27 5.3 參數 M 與 FEAST 的探討 27 5.4 參數 r 與 M 的交互影響 29 5.5 N 與 L 的比較 30 5.6 MLCIRR 與 Lancozs 的比較 31 6 結論 32 Reference 34

    J. Asakura, T. Sakurai, H. Tadano, T. Ikegami and K. Kimura, A numerical method for nonlinear eigenvalue problems using contour integrals, JSIAM Letters, Vol.1,52-55 (2009).

    Y. Futamura, H. Tadano, and T. Sakurai, Parallel stochastic estimation method of eigenvalue distribution, JSIAM Letters, Vol. 2, 127-130, (2010).

    T.-M. Huang, H.-E. Hsieh, W.-W. Lin and W. Wang, Eigendecomposition of the discrete double-curl operator with application to fast eigensolver for three dimensional photonic crystals, SIAM J. Matrix Anal. Appl., Vol. 34, 369-391(2013)

    T.-M. Huang, H.-E. Hsieh, W.-W. Lin, and W. Wang.Fast lanczos eigenvalue solvers for band structures of three dimensional photonic Crystals with face-centered cubic lattice(preprint)

    T.-M. HUANG, Y.-L. Huang.W.-W. Lin, and W.-C. Wang, A Null space free Jacobi-Davidson iteration for maxwell's operator(preprint)

    T.-M. Huang, W.-J. Chang, Y.-L. Huang, W.-W. Lin, W.-C. Wang and W. Wang, Preconditioning bandgap eigenvalue problems in three dimensional photonic crystals simulations, Journal of Computational Physics, Vol. 229,8684-8703(2010)

    T.-M. Huang, Y.-C. Kuo and W. Wang, Computing extremal eigenvalues for three-dimensional photonic crystals with wave vectors near the Brillouin zone center, J. Sci. Comput., Vol. 55, 529-551(2013)

    T. Ikegami and T. Sakurai, Contour integral eigensolver for non-Hermitian systems: a Rayleigh-Ritz-type approach, RANMEP2008Taiwanese J. Math., Vol. 14, pp. 825--837 (2010).

    T. Ikegami, T. Sakurai, U. Nagashima, A filter diagonalization for generalized eigenvalue problems based on the Sakurai-Sugiura projection method,CS-TR-08-13, Tsukuba (2008)

    A. Levin, D. Zhang, E. Polizzi, FEAST fundamental framework for electronic structure calculations: reformulation and solution of the muffintin problem computer physics communications, V183, I11, 2370-2375 (2012)

    E. Polizzi, Density-matrix-based algorithms for solving eingenvalue problems Phys. Rev. B., Vol. 79, 115112 (2009)

    E. Polizzi,A high-performance numerical library for solving eigenvalue problems: FEAST solver user's guide v2.1s,(2013)

    T. Sakurai, J. Asakura, H. Tadano, T. Ikegami and K. Kimura, A method for finding zeros of polynomial equations using a contour integral based eigensolver, Proc. Symbolic Numeric Computations 2009, Kyoto. 143-147 (2009)

    T. Sakurai, J. Asakura, H. Tadano and T. Ikegami, Error analysis for a matrix pencil of Hankel matrices with perturbed complex moments, JSIAM Letters, Vol. 1,76-79 (2009).

    T. Sakurai, H. Sugiura, A projection method for generalized eigenvalue problems, ISE-TR-02-189, Tsukuba,(2002)

    T. Sakurai, H. Tadano, CIRR: a Rayleigh-Ritz type method with contour integral for generalized eigenvalue problems, Proc. The First China-Japan-Korea Joint Conference on Numerical Mathematics, Vol. 36,745-757 (2007).

    T. Sakurai, H. Tadano, T. Ikegami, U. Nagashima, A parallel eigensolver using contour integration for generalized eigenvalue problems in molecular simulation, CS-TR-08-14, Tsukuba(2008)

    K. Senzaki, H. Tadano, T. Sakurai and Z. Bai, A method for profiling the distribution of eigenvalues using the AS method,RANMEP 2008 Taiwanese J. Math, Vol. 14,839-853 (2010).

    M. Suzuki, I. Suzuk, Bloch theorem and Energy band, Lecture Note on Solid State Physics(2006)

    P.-T. Tang, E. Polizzi,FEAST as a subspace iteration eigensolver Aacelerated by approximate spectral projection,SIMAX, (accepted)

    I. Yamazaki, T. Ikegami, H. Tadano, T. Sakurai, Performance comparison of parallel eigensolvers based on a contour integral method and a Lanczos method, Parallel Computing ,39,280–290.(2013)

    C. Yang, W.-G. Gao, Z.-J. Bai, X.-S. Li, L.-Q. Lee,
    P. Husbands. An algebraic substructuring method for large-scale eigenvalue calculation. SIAM J. SCI. COMPUT. Vol. 27, No. 3. 873–892(2005)

    K.-S. Yee, Numerical solution of initial boundary value problems involving Maxwell equations in isotropic media,IEEE Trans. Antennas Propagat.AP-14(3), 302-307. (1966)

    下載圖示
    QR CODE