Krylov A04 October 8, 2010 T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

Size: px
Start display at page:

Download "Krylov A04 October 8, 2010 T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48"

Transcription

1 Krylov A04 October 8, 2010 T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

2 Krylov QCD, RSDFT, Shell model Block Krylov MATLAB Scilab T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

3 Krylov Krylov n A n v K k (A, v) span(v, Av, A 2 v,..., A k 1 v). Krylov Krylov Alexei Krylov (Wikipedia ) Krylov T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

4 Krylov Krylov 1 Ax = b x 0 2 r 0 = b Ax 0 3 k = 1, 2,... x k x 0 K k (A, r 0 ) x k ρ (k) j x k = x 0 + k 1 j=0 ρ (k) j A j r 0 T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

5 Krylov x k r k = b Ax k r k = I k 1 ρ (k) j A j+1 j=0 r 0 = R k (A)r 0 K k+1 (A, r 0 ) R k (A) A k r k Krylov r k x k T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

6 Krylov Krylov 4 1 Ritz-Galerkin b Ax k K k (A, r 0 ) x k 2 b Ax k 2 K k (A, r 0 ) x k 3 Petrov-Galerkin b Ax k k x k 4 x k x 0 2 x k A T K k (A T, r 0 ) T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

7 (CG) Krylov 1952 Hestenes, Stiefel x e = x x f(x) f(x) = 1 2 (e, Ae) = 1 2 (x, Ax) (x, b) (x, b) T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

8 (CG) p k k + 1 x k+1 x k+1 = x k + α k p k α k f(x k+1 ) x k+1 r k+1 r k+1 = r k α k Ap k. p k+1 p k+1 = r k+1 + β k p k p k+1 (p i, Ap j ) = 0, i j β k+1 T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

9 Krylov A A + σi Krylov K k (A + σi, v) = K k (A, v). A Krylov (A + σi)x σ = b x σ T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

10 CG r 0 = b Ax 0 r k = R k (A)r 0 K k+1 (A, r 0 ) R k (z) Lanczos Lanczos R k+1 (z) R 0 (z) = 1, P 0 (z) = 1 R k+1 (z) = R k (z) α k zp k (z) P k+1 (z) = R k+1 (z) + β k P k (z) P k (z) T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

11 CG Ax = b x 0 = 0 r 0 = b r k = R k (A)b K k+1 (A, b). (A + σi)x σ = b x σ 0 = 0 rσ 0 = b r σ k = Rσ k (A + σi)b K k+1(a + σi, b). Krylov K k (A, b) = K k (A + σi, b) R σ k (z + σ) = ξσ k R k(z). T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

12 CG Lanczos R σ k (z + σ) = ξσ k R k(z) ξ σ k+1 = ξ σ k ξσ k 1 α k 1 (ξ σ k 1 ξσ k )α kβ k 1 + ξ σ k 1 α k 1(1 + σα k ) α σ k = ξσ k+1 α ξk σ k β σ k = ( ξ σ k+1 ξ σ k ) 2 β k r σ k = ξσ k r k T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

13 SS(Sakurai-Sugiura) N ( ) zj γ k+1 (z j I A) 1 ˆF k = 1 N j=1 ρ span{ŝ 0, ŝ 0,, ŝ m 1 } Sakurai, et.al. (2003) T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

14 ŝ k = 1 N N j=1 ( ) zj γ k w j (z j I A) 1 v ρ Krylov N (z 1 I A)y 1 = v (z 2 I A)y 2 = v. (z N I A)y N = v 1/N T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

15 A R n n A λ 1,..., λ n x 1,..., x n A n A = λ j P j, j=1 P j = x j x T j. P j λ j T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

16 (zi A) 1 re(a) (zi A) 1, z re(a) A A n A = λ j P j j=1 n (zi A) 1 P j = z λ j j=1 T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

17 Γ α Jordan f(z) Γ f(α) = 1 f(z) 2πi Γ z α dz. r(z) Γ α 1, α 2,..., α n Γ r(z)dz = 2πi n Res z=αj r(z) j=1 Res z=αj r(z) z = α j r(z) T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

18 α 1,..., α n R ν 1,..., ν n R 0 1 2πi Γ r(z)dz = r(z) = n j=1 n j=1 1 2πi ν j z α j Γ ν j z α j dz = n j=1 ν j 1 2πi Γ z k r(z)dz = n j=1 1 2πi Γ ν j z k z α j dz = n ν j α k j. j=1 T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

19 Γ A Jordan 1 (zi A) 1 dz = 2πi Γ n P j. j=1 f(z) f(a) = 1 2πi Γ f(z)(zi A) 1 dz. v R n A k v = 1 z k (zi A) 1 dz v. 2πi Γ T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

20 Γ Jordan λ 1,..., λ m Γ (zi A) 1 P Γ (A) = 1 (zi A) 1 dz 2πi Γ m P Γ (A) = j=1 P j P Γ (A) x 1,..., x m T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

21 P Γ (A) v s = P Γ (A)v s = 1 (zi A) 1 dz v = 2πi Γ m P j v. j=1 v F (x) = 1 1 2πi z x dz Γ T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

22 Γ γ ρ z j = γ + ρ e 2πi N (j+1/2), j = 0, 1,..., N 1 ŝ = 1 N N 1 j=0 e 2πi N (j+1/2) (z j I A) 1 v ˆF (x) = 1 N N 1 j=0 e 2πi N (j+1/2) z j x T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

23 QCD QCD all-to-all propagator A = γ 5 M Low mode part: (λ j, v j ) High mode part: A 1 High = 1 m A 1 Low = n j=1 1 m φ j 1 j=1 λ j v j v j n v j v j η j j=1 A 1 = A 1 Low + A 1 High T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

24 QCD O(a) κ = (497, , 380, 864) 20 shifted CG (Ohno, et.al. accepted) T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

25 QCD T2K-Tsukuba 1 node SS ( ) SS ( ) PARPACK [sec] [sec] ( ) ( ) Arnoldi SS QCD Implicitly restarted Arnoldi T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

26 48 Cr in the model space consisting of single-particle orbits f 7/2, p 3/2, f 5/2, p 1/2. M-scheme dimension for M = 0 : : double-lanczos : SS with Shifted COCG (Mizusaki, et.al. 2010) T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

27 Arnoldi SS Jordan Γ µ Γ = 1 2πi Γ tr((zb A) 1 B)dz. trace tr((zb A) 1 ) 1 s s v T j (z jb A) 1 v j, j=1 v j 1 1 n v j T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

28 (RSDFT) Si ,248, 2000 (Futamura, et.al. (submitted)) T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

29 T2K cores time[sec] speedup Linear Speedup NVIDIA Tesla C1060( ) GPUs time[sec] speedup GPU GPU 8 GPU Speedup T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

30 Matrix Market: olm100 ( et.al (2010)) T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

31 Block Krylov B n L AX = B Block BiCG O Leary(1980) Block GMRES Vital(1990) Block QMR Freund(1997) Block BiCGSTAB Guennouni(2003) 1 1 Block Krylov K k (A, R 0 ) = (R 0, AR 0, A 2 R 0,..., A k 1 R 0 ), R 0 = B AX 0 X 0 T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

32 Block Krylov 1 R k X k B AX k R k T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

33 BiCGGR Tadano, et.al (2009) Block Krylov T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

34 Block Krylov QCD , 572, , 216, 064 Block BiCGSTAB Block BiCGGR T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

35 (NEPs) A 0, A 1,..., A p C n n, A p O p T (z) = z i A i, z C p i=0 T (z) T (λ)x = 0 λ x 0 p = 1 A 1 = I T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

36 (NEPs) T (z) Γ T (z) T (λ)x = 0 λ x 0 T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

37 : T (z) = zi + A 0 + A 1 e τ z E. Jarlebring, 2008 : T (z) = K z 2 M + i l j=1 z 2 σ 2 j W j Acoustic field computations: Rich Lee, et. al., 2008 T (z) = K(z) + zd + z 2 M, K(z) J. Mehrmann, 2008 T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

38 : f(x, y, z) = 0, g(x, y, z) = 0, h(x, y, z) = 0. Dixon δ 1 f(x, y, z) g(x, y, z) h(x, y, z) δ = f(α, y, z) g(α, y, z) h(α, y, z) (x α)(y β) f(α, β, z) g(α, β, z) h(α, β, z) 1 x = [1, α,, α n 1 β n 1 ] T (z). x n 1 y n 1 T (z)u = 0 Sakurai, et. al, 2008 T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

39 m n n m 3 T (z) = A 0 + za 1 + z 2 A 2 + z 3 A 3 C T = O O A 0 I O A 1 O I A 2 T (z) T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

40 : T (z) = K z 2 M + i tx j=1 q z 2 σ 2 j W j, t = 1, σ 1 = 0. : 1,100,242 1,100,242 : T2K-Tsukuba System Quad-core AMD Opteron 2.3GHz (Barcelona) 4 sockets, 32GB of memory, 10,368 cores (95 TFlops), Infiniband Fat-tree Network Compiler : Intel C Library : Intel MKL PARDISO in MKL (shared memory) SuperLU DIST(distributed memory) Rich Lee, Stanford Linear Accelerator Center T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

41 Bottom Level PARDISO PARDISO 1 sec. Solution time in seconds sec Number of cores SS Top Level, Second Level MPI T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

42 4 (Γ 1 Γ 4 ) PARDISO 8 ( = 1024 ) T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

43 (Linear Solver) z shift-invert T (z) (SuperLU DIST) LU LU factorization (sec) T2K Franklin #procs T2K 1 (16 ) Franklin T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

44 (GEP) (GEP) 40 SS Lanczos (TRLan) SS sec Franklin(LBNL) TRLan 100 SS #procs 128 (TRLan) (SS method) T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

45 (NEP) (NEP) SS 40 sec 5000 T2K Franklin(LBNL) T2K 100 Franklin #procs Franklin 1 2 SS T2K Franklin T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

46 I Krylov,,,, Shifted Linear Systems Krylov,, Vol. 14, No. 3, pp (2004)., (2010), ISBN T. Mizusaki, K. Kaneko, M. Honma, and T. Sakurai, Filter diagonalization of shell-model calculations, Phys. Rev. C 82, (2010). [10 pages] H. Ohno, Y. Kuramashi, H. Tadano, T. Sakurai, A quadrature-based eigensolver with a Krylov subspace method for shifted linear systems for Hermitian eigenproblems in lattice QCD, JSIAM Letters (accepted). Y. Futamura, H. Tadano, T. Sakurai, Parallel stochastic estimation method of eigenvalue distribution, JSIAM Letters, submitted. Block Krylov H. Tadano, T. Sakurai, On single precision preconditioners for Krylov subspace iterative methods, Lecture Notes in Computer Science, No. 4818, pp (2007). H. Tadano, T. Sakurai and Y. Kuramashi, Block BiCGGR: A new block Krylov subspace method for computing high accuracy solutions, JSIAM Letters, Vol.1, pp (2009). T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

47 II,, Block Krylov,, Vol. 2 No. 2, pp (2009). T. Sakurai, H. Tadano and Y. Kuramashi, Application of block Krylov subspace algorithms to the Wilson-Dirac equation with multiple right-hand sides in lattice QCD, Comput. Phys. Commun., Vol. 181, pp (2010). H. Tadano, Y. Kuramashi and T. Sakurai, Application of preconditioned block BiCGGR to the Wilson-Dirac equation with multiple right-hand sides in lattice QCD, Comput. Phys. Commun., Vol. 181, pp (2010). SS T. Sakurai and H. Sugiura, A projection method for generalized eigenvalue problems, J. Comput. Appl. Math. Vol. 159, pp (2003). 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, pp (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, pp (2009). T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

48 III 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, pp (2009). T. Ikegami, T. Sakurai and U. Nagashima, A filter diagonalization for generalized eigenvalue problems based on the Sakurai-Sugiura projection method, J. Comput, Appl. Math., Vol. 233, (2010). J. Asakura, T. Sakurai, H. Tadano, T. Ikegami and K. Kimura, A numerical method for polynomial eigenvalue problems using contour integral, Japan J. Indust. Appl. Math., Vol. 27, pp (2010). I. Yamazaki, M. Okada, H. Tadano, T. Sakurai and K. Teranishi, A block sparse approximate inverse with cutoff preconditioner for semi-sparse linear systems derived from Molecular Orbital calculations, JSIAM Letters, Vol. 2, pp (2010).,,,,, Cutoff 2,, Vol. 18, No. 4, pp (2008).,,,,, Vol. 17, No. 3, pp (2007). T. Sakurai (Univ. Tsukuba) Krylov October 8, / 48

untitled

untitled 1 1 Ax = b A R m m A b R m x R m A shift-and invert Lanczos - LU CG A = LU LU Ly = b Ux = y A LU A A = LL T 1 LU b,, Vol. 11, No. 4, pp. 14 18 (2006). x * x (0), x (1), x (2), A Ap A # x (n+1) = Cx (n)

More information

2 (March 13, 2010) N Λ a = i,j=1 x i ( d (a) i,j x j ), Λ h = N i,j=1 x i ( d (h) i,j x j ) B a B h B a = N i,j=1 ν i d (a) i,j, B h = x j N i,j=1 ν i

2 (March 13, 2010) N Λ a = i,j=1 x i ( d (a) i,j x j ), Λ h = N i,j=1 x i ( d (h) i,j x j ) B a B h B a = N i,j=1 ν i d (a) i,j, B h = x j N i,j=1 ν i 1. A. M. Turing [18] 60 Turing A. Gierer H. Meinhardt [1] : (GM) ) a t = D a a xx µa + ρ (c a2 h + ρ 0 (0 < x < l, t > 0) h t = D h h xx νh + c ρ a 2 (0 < x < l, t > 0) a x = h x = 0 (x = 0, l) a = a(x,

More information

医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 第 2 版 1 刷発行時のものです.

医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます.   このサンプルページの内容は, 第 2 版 1 刷発行時のものです. 医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/009192 このサンプルページの内容は, 第 2 版 1 刷発行時のものです. i 2 t 1. 2. 3 2 3. 6 4. 7 5. n 2 ν 6. 2 7. 2003 ii 2 2013 10 iii 1987

More information

橡固有値セミナー2_棚橋改.PDF

橡固有値セミナー2_棚橋改.PDF 1 II. 2003 5 14 2... Arnoldi. Lanczos. Jacobi-Davidson . 3 4 Ax = x A A Ax = Mx M: M 5 Householder ln ln-1 0 l3 0 l2 l1 6 Lanczos Lanczos, 1950 Arnoldi Arnoldi, 1951 Hessenberg Jacobi-Davidson Sleijpen

More information

2012年度HPCサマーセミナー_多田野.pptx

2012年度HPCサマーセミナー_多田野.pptx ! CCS HPC! I " tadano@cs.tsukuba.ac.jp" " 1 " " " " " " " 2 3 " " Ax = b" " " 4 Ax = b" A = a 11 a 12... a 1n a 21 a 22... a 2n...... a n1 a n2... a nn, x = x 1 x 2. x n, b = b 1 b 2. b n " " 5 Gauss LU

More information

untitled

untitled A = QΛQ T A n n Λ Q A = XΛX 1 A n n Λ X GPGPU A 3 T Q T AQ = T (Q: ) T u i = λ i u i T {λ i } {u i } QR MR 3 v i = Q u i A {v i } A n = 9000 Quad Core Xeon 2 LAPACK (4/3) n 3 O(n 2 ) O(n 3 ) A {v i }

More information

untitled

untitled A = QΛQ T A n n Λ Q A = XΛX 1 A n n Λ X GPGPU A 3 T Q T AQ = T (Q: ) T u i = λ i u i T {λ i } {u i } QR MR 3 v i = Q u i A {v i } A n = 9000 Quad Core Xeon 2 LAPACK (4/3) n 3 O(n 2 ) O(n 3 ) A {v i }

More information

微分積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.

微分積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます.   このサンプルページの内容は, 初版 1 刷発行時のものです. 微分積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. ttp://www.morikita.co.jp/books/mid/00571 このサンプルページの内容は, 初版 1 刷発行時のものです. i ii 014 10 iii [note] 1 3 iv 4 5 3 6 4 x 0 sin x x 1 5 6 z = f(x, y) 1 y = f(x)

More information

IDRstab(s, L) GBiCGSTAB(s, L) 2. AC-GBiCGSTAB(s, L) Ax = b (1) A R n n x R n b R n 2.1 IDR s L r k+1 r k+1 = b Ax k+1 IDR(s) r k+1 = (I ω k A)(r k dr

IDRstab(s, L) GBiCGSTAB(s, L) 2. AC-GBiCGSTAB(s, L) Ax = b (1) A R n n x R n b R n 2.1 IDR s L r k+1 r k+1 = b Ax k+1 IDR(s) r k+1 = (I ω k A)(r k dr 1 2 IDR(s) GBiCGSTAB(s, L) IDR(s) IDRstab(s, L) GBiCGSTAB(s, L) Verification of effectiveness of Auto-Correction technique applied to preconditioned iterative methods Keiichi Murakami 1 Seiji Fujino 2

More information

I A A441 : April 15, 2013 Version : 1.1 I Kawahira, Tomoki TA (Shigehiro, Yoshida )

I A A441 : April 15, 2013 Version : 1.1 I   Kawahira, Tomoki TA (Shigehiro, Yoshida ) I013 00-1 : April 15, 013 Version : 1.1 I Kawahira, Tomoki TA (Shigehiro, Yoshida) http://www.math.nagoya-u.ac.jp/~kawahira/courses/13s-tenbou.html pdf * 4 15 4 5 13 e πi = 1 5 0 5 7 3 4 6 3 6 10 6 17

More information

ペタスケール計算環境に向けたFFTライブラリ

ペタスケール計算環境に向けたFFTライブラリ A01 高橋班 大規模並列環境における 数値計算アルゴリズム 研究代表者 : 高橋大介 筑波大学大学院システム情報工学研究科 研究組織 研究代表者 高橋大介 ( 筑波大学 ): 研究統括および高速アルゴリズム 研究分担者 今村俊幸 ( 電気通信大学 ): 性能チューニング 多田野寛人 ( 筑波大学 ): 大規模線形計算 連携研究者 佐藤三久 ( 筑波大学 ): 並列システムの性能評価 朴泰祐 ( 筑波大学

More information

²ÄÀÑʬΥ»¶ÈóÀþ·¿¥·¥å¥ì¡¼¥Ç¥£¥ó¥¬¡¼ÊýÄø¼°¤ÎÁ²¶á²òÀÏ Asymptotic analysis for the integrable discrete nonlinear Schrödinger equation

²ÄÀÑʬΥ»¶ÈóÀþ·¿¥·¥å¥ì¡¼¥Ç¥£¥ó¥¬¡¼ÊýÄø¼°¤ÎÁ²¶á²òÀÏ  Asymptotic analysis for the integrable discrete nonlinear Schrödinger equation Asymptotic analysis for the integrable discrete nonlinear Schrödinger equation ( ) ( ) 2016 12 17 1. Schrödinger focusing NLS iu t + u xx +2 u 2 u = 0 u(x, t) =2ηe 2iξx 4i(ξ2 η 2 )t+i(ψ 0 +π/2) sech(2ηx

More information

201711grade1ouyou.pdf

201711grade1ouyou.pdf 2017 11 26 1 2 52 3 12 13 22 23 32 33 42 3 5 3 4 90 5 6 A 1 2 Web Web 3 4 1 2... 5 6 7 7 44 8 9 1 2 3 1 p p >2 2 A 1 2 0.6 0.4 0.52... (a) 0.6 0.4...... B 1 2 0.8-0.2 0.52..... (b) 0.6 0.52.... 1 A B 2

More information

EGunGPU

EGunGPU Super Computing in Accelerator simulations - Electron Gun simulation using GPGPU - K. Ohmi, KEK-Accel Accelerator Physics seminar 2009.11.19 Super computers in KEK HITACHI SR11000 POWER5 16 24GB 16 134GFlops,

More information

量子力学 問題

量子力学 問題 3 : 203 : 0. H = 0 0 2 6 0 () = 6, 2 = 2, 3 = 3 3 H 6 2 3 ϵ,2,3 (2) ψ = (, 2, 3 ) ψ Hψ H (3) P i = i i P P 2 = P 2 P 3 = P 3 P = O, P 2 i = P i (4) P + P 2 + P 3 = E 3 (5) i ϵ ip i H 0 0 (6) R = 0 0 [H,

More information

Dirac 38 5 Dirac 4 4 γ µ p µ p µ + m 2 = ( p µ γ µ + m)(p ν γ ν + m) (5.1) γ = p µ p ν γ µ γ ν p µ γ µ m + mp ν γ ν + m 2 = 1 2 p µp ν {γ µ, γ ν } + m

Dirac 38 5 Dirac 4 4 γ µ p µ p µ + m 2 = ( p µ γ µ + m)(p ν γ ν + m) (5.1) γ = p µ p ν γ µ γ ν p µ γ µ m + mp ν γ ν + m 2 = 1 2 p µp ν {γ µ, γ ν } + m Dirac 38 5 Dirac 4 4 γ µ p µ p µ + m 2 p µ γ µ + mp ν γ ν + m 5.1 γ p µ p ν γ µ γ ν p µ γ µ m + mp ν γ ν + m 2 1 2 p µp ν {γ µ, γ ν } + m 2 5.2 p m p p µ γ µ {, } 10 γ {γ µ, γ ν } 2η µν 5.3 p µ γ µ + mp

More information

I

I I 6 4 10 1 1 1.1............... 1 1................ 1 1.3.................... 1.4............... 1.4.1.............. 1.4................. 1.4.3........... 3 1.4.4.. 3 1.5.......... 3 1.5.1..............

More information

HPC (pay-as-you-go) HPC Web 2

HPC (pay-as-you-go) HPC Web 2 ,, 1 HPC (pay-as-you-go) HPC Web 2 HPC Amazon EC2 OpenFOAM GPU EC2 3 HPC MPI MPI Courant 1 GPGPU MPI 4 AMAZON EC2 GPU CLUSTER COMPUTE INSTANCE EC2 GPU (cg1.4xlarge) ( N. Virgina ) Quadcore Intel Xeon 5570

More information

x (x, ) x y (, y) iy x y z = x + iy (x, y) (r, θ) r = x + y, θ = tan ( y ), π < θ π x r = z, θ = arg z z = x + iy = r cos θ + ir sin θ = r(cos θ + i s

x (x, ) x y (, y) iy x y z = x + iy (x, y) (r, θ) r = x + y, θ = tan ( y ), π < θ π x r = z, θ = arg z z = x + iy = r cos θ + ir sin θ = r(cos θ + i s ... x, y z = x + iy x z y z x = Rez, y = Imz z = x + iy x iy z z () z + z = (z + z )() z z = (z z )(3) z z = ( z z )(4)z z = z z = x + y z = x + iy ()Rez = (z + z), Imz = (z z) i () z z z + z z + z.. z

More information

数学の基礎訓練I

数学の基礎訓練I I 9 6 13 1 1 1.1............... 1 1................ 1 1.3.................... 1.4............... 1.4.1.............. 1.4................. 3 1.4.3........... 3 1.4.4.. 3 1.5.......... 3 1.5.1..............

More information

t = h x z z = h z = t (x, z) (v x (x, z, t), v z (x, z, t)) ρ v x x + v z z = 0 (1) 2-2. (v x, v z ) φ(x, z, t) v x = φ x, v z

t = h x z z = h z = t (x, z) (v x (x, z, t), v z (x, z, t)) ρ v x x + v z z = 0 (1) 2-2. (v x, v z ) φ(x, z, t) v x = φ x, v z I 1 m 2 l k 2 x = 0 x 1 x 1 2 x 2 g x x 2 x 1 m k m 1-1. L x 1, x 2, ẋ 1, ẋ 2 ẋ 1 x = 0 1-2. 2 Q = x 1 + x 2 2 q = x 2 x 1 l L Q, q, Q, q M = 2m µ = m 2 1-3. Q q 1-4. 2 x 2 = h 1 x 1 t = 0 2 1 t x 1 (t)

More information

nsg02-13/ky045059301600033210

nsg02-13/ky045059301600033210 φ φ φ φ κ κ α α μ μ α α μ χ et al Neurosci. Res. Trpv J Physiol μ μ α α α β in vivo β β β β β β β β in vitro β γ μ δ μδ δ δ α θ α θ α In Biomechanics at Micro- and Nanoscale Levels, Volume I W W v W

More information

Design of highly accurate formulas for numerical integration in weighted Hardy spaces with the aid of potential theory 1 Ken ichiro Tanaka 1 Ω R m F I = F (t) dt (1.1) Ω m m 1 m = 1 1 Newton-Cotes Gauss

More information

163 KdV KP Lax pair L, B L L L 1/2 W 1 LW = ( / x W t 1, t 2, t 3, ψ t n ψ/ t n = B nψ (KdV B n = L n/2 KP B n = L n KdV KP Lax W Lax τ KP L ψ τ τ Cha

163 KdV KP Lax pair L, B L L L 1/2 W 1 LW = ( / x W t 1, t 2, t 3, ψ t n ψ/ t n = B nψ (KdV B n = L n/2 KP B n = L n KdV KP Lax W Lax τ KP L ψ τ τ Cha 63 KdV KP Lax pair L, B L L L / W LW / x W t, t, t 3, ψ t n / B nψ KdV B n L n/ KP B n L n KdV KP Lax W Lax τ KP L ψ τ τ Chapter 7 An Introduction to the Sato Theory Masayui OIKAWA, Faculty of Engneering,

More information

( ) 5 Reduction ( ) A M n (C) Av = λv (v 0) (11.1) λ C A (eigenvalue) v C n A λ (eigenvector) M n (R) A λ(a) A M n (R) n A λ

( ) 5 Reduction ( ) A M n (C) Av = λv (v 0) (11.1) λ C A (eigenvalue) v C n A λ (eigenvector) M n (R) A λ(a) A M n (R) n A λ 125 11 ( ) 5 Reduction 11.1 11.1.1 ( ) A M n (C) Av = λv (v 0) (11.1) λ C A (eigenvalue) v C n A λ (eigenvector) M n (R) A λ(a) 11.1.2 A M n (R) n A λi = 0 A C n 5 126 11 A n λ 1 (A) λ 2 (A) λ n (A) A

More information

Note.tex 2008/09/19( )

Note.tex 2008/09/19( ) 1 20 9 19 2 1 5 1.1........................ 5 1.2............................. 8 2 9 2.1............................. 9 2.2.............................. 10 3 13 3.1.............................. 13 3.2..................................

More information

2 1 1 α = a + bi(a, b R) α (conjugate) α = a bi α (absolute value) α = a 2 + b 2 α (norm) N(α) = a 2 + b 2 = αα = α 2 α (spure) (trace) 1 1. a R aα =

2 1 1 α = a + bi(a, b R) α (conjugate) α = a bi α (absolute value) α = a 2 + b 2 α (norm) N(α) = a 2 + b 2 = αα = α 2 α (spure) (trace) 1 1. a R aα = 1 1 α = a + bi(a, b R) α (conjugate) α = a bi α (absolute value) α = a + b α (norm) N(α) = a + b = αα = α α (spure) (trace) 1 1. a R aα = aα. α = α 3. α + β = α + β 4. αβ = αβ 5. β 0 6. α = α ( ) α = α

More information

1 4 1 ( ) ( ) ( ) ( ) () 1 4 2

1 4 1 ( ) ( ) ( ) ( ) () 1 4 2 7 1995, 2017 7 21 1 2 2 3 3 4 4 6 (1).................................... 6 (2)..................................... 6 (3) t................. 9 5 11 (1)......................................... 11 (2)

More information

°ÌÁê¿ô³ØII

°ÌÁê¿ô³ØII July 14, 2007 Brouwer f f(x) = x x f(z) = 0 2 f : S 2 R 2 f(x) = f( x) x S 2 3 3 2 - - - 1. X x X U(x) U(x) x U = {U(x) x X} X 1. U(x) A U(x) x 2. A U(x), A B B U(x) 3. A, B U(x) A B U(x) 4. A U(x),

More information

A Feasibility Study of Direct-Mapping-Type Parallel Processing Method to Solve Linear Equations in Load Flow Calculations Hiroaki Inayoshi, Non-member

A Feasibility Study of Direct-Mapping-Type Parallel Processing Method to Solve Linear Equations in Load Flow Calculations Hiroaki Inayoshi, Non-member A Feasibility Study of Direct-Mapping-Type Parallel Processing Method to Solve Linear Equations in Load Flow Calculations Hiroaki Inayoshi, Non-member (University of Tsukuba), Yasuharu Ohsawa, Member (Kobe

More information

最新耐震構造解析 ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 第 3 版 1 刷発行時のものです.

最新耐震構造解析 ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます.   このサンプルページの内容は, 第 3 版 1 刷発行時のものです. 最新耐震構造解析 ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/052093 このサンプルページの内容は, 第 3 版 1 刷発行時のものです. i 3 10 3 2000 2007 26 8 2 SI SI 20 1996 2000 SI 15 3 ii 1 56 6

More information

t χ 2 F Q t χ 2 F 1 2 µ, σ 2 N(µ, σ 2 ) f(x µ, σ 2 ) = 1 ( exp (x ) µ)2 2πσ 2 2σ 2 0, N(0, 1) (100 α) z(α) t χ 2 *1 2.1 t (i)x N(µ, σ 2 ) x µ σ N(0, 1

t χ 2 F Q t χ 2 F 1 2 µ, σ 2 N(µ, σ 2 ) f(x µ, σ 2 ) = 1 ( exp (x ) µ)2 2πσ 2 2σ 2 0, N(0, 1) (100 α) z(α) t χ 2 *1 2.1 t (i)x N(µ, σ 2 ) x µ σ N(0, 1 t χ F Q t χ F µ, σ N(µ, σ ) f(x µ, σ ) = ( exp (x ) µ) πσ σ 0, N(0, ) (00 α) z(α) t χ *. t (i)x N(µ, σ ) x µ σ N(0, ) (ii)x,, x N(µ, σ ) x = x+ +x N(µ, σ ) (iii) (i),(ii) z = x µ N(0, ) σ N(0, ) ( 9 97.

More information

ohpmain.dvi

ohpmain.dvi fujisawa@ism.ac.jp 1 Contents 1. 2. 3. 4. γ- 2 1. 3 10 5.6, 5.7, 5.4, 5.5, 5.8, 5.5, 5.3, 5.6, 5.4, 5.2. 5.5 5.6 +5.7 +5.4 +5.5 +5.8 +5.5 +5.3 +5.6 +5.4 +5.2 =5.5. 10 outlier 5 5.6, 5.7, 5.4, 5.5, 5.8,

More information

I , : ~/math/functional-analysis/functional-analysis-1.tex

I , : ~/math/functional-analysis/functional-analysis-1.tex I 1 2004 8 16, 2017 4 30 1 : ~/math/functional-analysis/functional-analysis-1.tex 1 3 1.1................................... 3 1.2................................... 3 1.3.....................................

More information

,, Andrej Gendiar (Density Matrix Renormalization Group, DMRG) 1 10 S.R. White [1, 2] 2 DMRG ( ) [3, 2] DMRG Baxter [4, 5] 2 Ising 2 1 Ising 1 1 Ising

,, Andrej Gendiar (Density Matrix Renormalization Group, DMRG) 1 10 S.R. White [1, 2] 2 DMRG ( ) [3, 2] DMRG Baxter [4, 5] 2 Ising 2 1 Ising 1 1 Ising ,, Andrej Gendiar (Density Matrix Renormalization Group, DMRG) 1 10 S.R. White [1, 2] 2 DMRG ( ) [3, 2] DMRG Baxter [4, 5] 2 Ising 2 1 Ising 1 1 Ising Model 1 Ising 1 Ising Model N Ising (σ i = ±1) (Free

More information

tomocci ,. :,,,, Lie,,,, Einstein, Newton. 1 M n C. s, M p. M f, p d ds f = dxµ p ds µ f p, X p = X µ µ p = dxµ ds µ p. µ, X µ.,. p,. T M p.

tomocci ,. :,,,, Lie,,,, Einstein, Newton. 1 M n C. s, M p. M f, p d ds f = dxµ p ds µ f p, X p = X µ µ p = dxµ ds µ p. µ, X µ.,. p,. T M p. tomocci 18 7 5...,. :,,,, Lie,,,, Einstein, Newton. 1 M n C. s, M p. M f, p d ds f = dxµ p ds µ f p, X p = X µ µ p = dxµ ds µ p. µ, X µ.,. p,. T M p. M F (M), X(F (M)).. T M p e i = e µ i µ. a a = a i

More information

TOP URL 1

TOP URL   1 TOP URL http://amonphys.web.fc.com/ 3.............................. 3.............................. 4.3 4................... 5.4........................ 6.5........................ 8.6...........................7

More information

Trapezoidal Rule θ = 1/ x n x n 1 t = 1 [f(t n 1, x n 1 ) + f(t n, x n )] (6) 1. dx dt = f(t, x), x(t 0) = x 0 (7) t [t 0, t 1 ] f t [t 0, t 1 ], x x

Trapezoidal Rule θ = 1/ x n x n 1 t = 1 [f(t n 1, x n 1 ) + f(t n, x n )] (6) 1. dx dt = f(t, x), x(t 0) = x 0 (7) t [t 0, t 1 ] f t [t 0, t 1 ], x x University of Hyogo 8 8 1 d x(t) =f(t, x(t)), dt (1) x(t 0 ) =x 0 () t n = t 0 + n t x x n n x n x 0 x i i = 0,..., n 1 x n x(t) 1 1.1 1 1 1 0 θ 1 θ x n x n 1 t = θf(t n 1, x n 1 ) + (1 θ)f(t n, x n )

More information

: , 2.0, 3.0, 2.0, (%) ( 2.

: , 2.0, 3.0, 2.0, (%) ( 2. 2017 1 2 1.1...................................... 2 1.2......................................... 4 1.3........................................... 10 1.4................................. 14 1.5..........................................

More information

第5章 偏微分方程式の境界値問題

第5章 偏微分方程式の境界値問題 October 5, 2018 1 / 113 4 ( ) 2 / 113 Poisson 5.1 Poisson ( A.7.1) Poisson Poisson 1 (A.6 ) Γ p p N u D Γ D b 5.1.1: = Γ D Γ N 3 / 113 Poisson 5.1.1 d {2, 3} Lipschitz (A.5 ) Γ D Γ N = \ Γ D Γ p Γ N Γ

More information

Macdonald, ,,, Macdonald. Macdonald,,,,,.,, Gauss,,.,, Lauricella A, B, C, D, Gelfand, A,., Heckman Opdam.,,,.,,., intersection,. Macdona

Macdonald, ,,, Macdonald. Macdonald,,,,,.,, Gauss,,.,, Lauricella A, B, C, D, Gelfand, A,., Heckman Opdam.,,,.,,., intersection,. Macdona Macdonald, 2015.9.1 9.2.,,, Macdonald. Macdonald,,,,,.,, Gauss,,.,, Lauricella A, B, C, D, Gelfand, A,., Heckman Opdam.,,,.,,., intersection,. Macdonald,, q., Heckman Opdam q,, Macdonald., 1 ,,. Macdonald,

More information

I, II 1, A = A 4 : 6 = max{ A, } A A 10 10%

I, II 1, A = A 4 : 6 = max{ A, } A A 10 10% 1 2006.4.17. A 3-312 tel: 092-726-4774, e-mail: hara@math.kyushu-u.ac.jp, http://www.math.kyushu-u.ac.jp/ hara/lectures/lectures-j.html Office hours: B A I ɛ-δ ɛ-δ 1. 2. A 1. 1. 2. 3. 4. 5. 2. ɛ-δ 1. ɛ-n

More information

(τ τ ) τ, σ ( ) w = τ iσ, w = τ + iσ (w ) w, w ( ) τ, σ τ = (w + w), σ = i (w w) w, w w = τ w τ + σ w σ = τ + i σ w = τ w τ + σ w σ = τ i σ g ab w, w

(τ τ ) τ, σ ( ) w = τ iσ, w = τ + iσ (w ) w, w ( ) τ, σ τ = (w + w), σ = i (w w) w, w w = τ w τ + σ w σ = τ + i σ w = τ w τ + σ w σ = τ i σ g ab w, w S = 4π dτ dσ gg ij i X µ j X ν η µν η µν g ij g ij = g ij = ( 0 0 ) τ, σ (+, +) τ τ = iτ ds ds = dτ + dσ ds = dτ + dσ δ ij ( ) a =, a = τ b = σ g ij δ ab g g ( +, +,... ) S = 4π S = 4π ( i) = i 4π dτ dσ

More information

2 1 κ c(t) = (x(t), y(t)) ( ) det(c (t), c x (t)) = det (t) x (t) y (t) y = x (t)y (t) x (t)y (t), (t) c (t) = (x (t)) 2 + (y (t)) 2. c (t) =

2 1 κ c(t) = (x(t), y(t)) ( ) det(c (t), c x (t)) = det (t) x (t) y (t) y = x (t)y (t) x (t)y (t), (t) c (t) = (x (t)) 2 + (y (t)) 2. c (t) = 1 1 1.1 I R 1.1.1 c : I R 2 (i) c C (ii) t I c (t) (0, 0) c (t) c(i) c c(t) 1.1.2 (1) (2) (3) (1) r > 0 c : R R 2 : t (r cos t, r sin t) (2) C f : I R c : I R 2 : t (t, f(t)) (3) y = x c : R R 2 : t (t,

More information

K 2 X = 4 MWG(f), X P 2 F, υ 0 : X P 2 2,, {f λ : X λ P 1 } λ Λ NS(X λ ), (υ 0 ) λ : X λ P 2 ( 1) X 6, f λ K X + F, f ( 1), n, n 1 (cf [10]) X, f : X

K 2 X = 4 MWG(f), X P 2 F, υ 0 : X P 2 2,, {f λ : X λ P 1 } λ Λ NS(X λ ), (υ 0 ) λ : X λ P 2 ( 1) X 6, f λ K X + F, f ( 1), n, n 1 (cf [10]) X, f : X 2 E 8 1, E 8, [6], II II, E 8, 2, E 8,,, 2 [14],, X/C, f : X P 1 2 3, f, (O), f X NS(X), (O) T ( 1), NS(X), T [15] : MWG(f) NS(X)/T, MWL(f) 0 (T ) NS(X), MWL(f) MWL(f) 0, : {f λ : X λ P 1 } λ Λ NS(X λ

More information

x () g(x) = f(t) dt f(x), F (x) 3x () g(x) g (x) f(x), F (x) (3) h(x) = x 3x tf(t) dt.9 = {(x, y) ; x, y, x + y } f(x, y) = xy( x y). h (x) f(x), F (x

x () g(x) = f(t) dt f(x), F (x) 3x () g(x) g (x) f(x), F (x) (3) h(x) = x 3x tf(t) dt.9 = {(x, y) ; x, y, x + y } f(x, y) = xy( x y). h (x) f(x), F (x [ ] IC. f(x) = e x () f(x) f (x) () lim f(x) lim f(x) x + x (3) lim f(x) lim f(x) x + x (4) y = f(x) ( ) ( s46). < a < () a () lim a log xdx a log xdx ( ) n (3) lim log k log n n n k=.3 z = log(x + y ),

More information

II A A441 : October 02, 2014 Version : Kawahira, Tomoki TA (Kondo, Hirotaka )

II A A441 : October 02, 2014 Version : Kawahira, Tomoki TA (Kondo, Hirotaka ) II 214-1 : October 2, 214 Version : 1.1 Kawahira, Tomoki TA (Kondo, Hirotaka ) http://www.math.nagoya-u.ac.jp/~kawahira/courses/14w-biseki.html pdf 1 2 1 9 1 16 1 23 1 3 11 6 11 13 11 2 11 27 12 4 12 11

More information

x V x x V x, x V x = x + = x +(x+x )=(x +x)+x = +x = x x = x x = x =x =(+)x =x +x = x +x x = x ( )x = x =x =(+( ))x =x +( )x = x +( )x ( )x = x x x R

x V x x V x, x V x = x + = x +(x+x )=(x +x)+x = +x = x x = x x = x =x =(+)x =x +x = x +x x = x ( )x = x =x =(+( ))x =x +( )x = x +( )x ( )x = x x x R V (I) () (4) (II) () (4) V K vector space V vector K scalor K C K R (I) x, y V x + y V () (x + y)+z = x +(y + z) (2) x + y = y + x (3) V x V x + = x (4) x V x + x = x V x x (II) x V, α K αx V () (α + β)x

More information

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63>

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63> 通信方式第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/072662 このサンプルページの内容は, 第 2 版発行当時のものです. i 2 2 2 2012 5 ii,.,,,,,,.,.,,,,,.,,.,,..,,,,.,,.,.,,.,,.. 1990 5 iii 1 1

More information

Microsoft Word - 11問題表紙(選択).docx

Microsoft Word - 11問題表紙(選択).docx A B A.70g/cm 3 B.74g/cm 3 B C 70at% %A C B at% 80at% %B 350 C γ δ y=00 x-y ρ l S ρ C p k C p ρ C p T ρ l t l S S ξ S t = ( k T ) ξ ( ) S = ( k T) ( ) t y ξ S ξ / t S v T T / t = v T / y 00 x v S dy dx

More information

磁性物理学 - 遷移金属化合物磁性のスピンゆらぎ理論

磁性物理学 - 遷移金属化合物磁性のスピンゆらぎ理論 email: takahash@sci.u-hyogo.ac.jp May 14, 2009 Outline 1. 2. 3. 4. 5. 6. 2 / 262 Today s Lecture: Mode-mode Coupling Theory 100 / 262 Part I Effects of Non-linear Mode-Mode Coupling Effects of Non-linear

More information

n ξ n,i, i = 1,, n S n ξ n,i n 0 R 1,.. σ 1 σ i .10.14.15 0 1 0 1 1 3.14 3.18 3.19 3.14 3.14,. ii 1 1 1.1..................................... 1 1............................... 3 1.3.........................

More information

Z: Q: R: C: sin 6 5 ζ a, b

Z: Q: R: C: sin 6 5 ζ a, b Z: Q: R: C: 3 3 7 4 sin 6 5 ζ 9 6 6............................... 6............................... 6.3......................... 4 7 6 8 8 9 3 33 a, b a bc c b a a b 5 3 5 3 5 5 3 a a a a p > p p p, 3,

More information

S I. dy fx x fx y fx + C 3 C vt dy fx 4 x, y dy yt gt + Ct + C dt v e kt xt v e kt + C k x v k + C C xt v k 3 r r + dr e kt S Sr πr dt d v } dt k e kt

S I. dy fx x fx y fx + C 3 C vt dy fx 4 x, y dy yt gt + Ct + C dt v e kt xt v e kt + C k x v k + C C xt v k 3 r r + dr e kt S Sr πr dt d v } dt k e kt S I. x yx y y, y,. F x, y, y, y,, y n http://ayapin.film.s.dendai.ac.jp/~matuda n /TeX/lecture.html PDF PS yx.................................... 3.3.................... 9.4................5..............

More information

D = [a, b] [c, d] D ij P ij (ξ ij, η ij ) f S(f,, {P ij }) S(f,, {P ij }) = = k m i=1 j=1 m n f(ξ ij, η ij )(x i x i 1 )(y j y j 1 ) = i=1 j

D = [a, b] [c, d] D ij P ij (ξ ij, η ij ) f S(f,, {P ij }) S(f,, {P ij }) = = k m i=1 j=1 m n f(ξ ij, η ij )(x i x i 1 )(y j y j 1 ) = i=1 j 6 6.. [, b] [, d] ij P ij ξ ij, η ij f Sf,, {P ij } Sf,, {P ij } k m i j m fξ ij, η ij i i j j i j i m i j k i i j j m i i j j k i i j j kb d {P ij } lim Sf,, {P ij} kb d f, k [, b] [, d] f, d kb d 6..

More information

7 π L int = gψ(x)ψ(x)φ(x) + (7.4) [ ] p ψ N = n (7.5) π (π +,π 0,π ) ψ (σ, σ, σ )ψ ( A) σ τ ( L int = gψψφ g N τ ) N π * ) (7.6) π π = (π, π, π ) π ±

7 π L int = gψ(x)ψ(x)φ(x) + (7.4) [ ] p ψ N = n (7.5) π (π +,π 0,π ) ψ (σ, σ, σ )ψ ( A) σ τ ( L int = gψψφ g N τ ) N π * ) (7.6) π π = (π, π, π ) π ± 7 7. ( ) SU() SU() 9 ( MeV) p 98.8 π + π 0 n 99.57 9.57 97.4 497.70 δm m 0.4%.% 0.% 0.8% π 9.57 4.96 Σ + Σ 0 Σ 89.6 9.46 K + K 0 49.67 (7.) p p = αp + βn, n n = γp + δn (7.a) [ ] p ψ ψ = Uψ, U = n [ α

More information

index calculus

index calculus index calculus 2008 3 8 1 generalized Weil descent p :, E/F p 3 : Y 2 = f(x), where f(x) = X 3 + AX + B, A F p, B F p 3 E(F p 3) 3 : Generalized Weil descent E(F p 4) 2 Index calculus Plain version Double-large-prime

More information

4 Mindlin -Reissner 4 δ T T T εσdω= δ ubdω+ δ utd Γ Ω Ω Γ T εσ (1.1) ε σ u b t 3 σ ε. u T T T = = = { σx σ y σ z τxy τ yz τzx} { εx εy εz γ xy γ yz γ

4 Mindlin -Reissner 4 δ T T T εσdω= δ ubdω+ δ utd Γ Ω Ω Γ T εσ (1.1) ε σ u b t 3 σ ε. u T T T = = = { σx σ y σ z τxy τ yz τzx} { εx εy εz γ xy γ yz γ Mindlin -Rissnr δ εσd δ ubd+ δ utd Γ Γ εσ (.) ε σ u b t σ ε. u { σ σ σ z τ τ z τz} { ε ε εz γ γ z γ z} { u u uz} { b b bz} b t { t t tz}. ε u u u u z u u u z u u z ε + + + (.) z z z (.) u u NU (.) N U

More information

20 4 20 i 1 1 1.1............................ 1 1.2............................ 4 2 11 2.1................... 11 2.2......................... 11 2.3....................... 19 3 25 3.1.............................

More information

[1] convention Minkovski i Polchinski [2] 1 Clifford Spin 1 2 Euclid Clifford 2 3 Euclid Spin 6 4 Euclid Pin Clifford Spin 10 A 12 B 17 1 Cliffo

[1] convention Minkovski i Polchinski [2] 1 Clifford Spin 1 2 Euclid Clifford 2 3 Euclid Spin 6 4 Euclid Pin Clifford Spin 10 A 12 B 17 1 Cliffo [1] convention Minkovski i Polchinski [2] 1 Clifford Spin 1 2 Euclid Clifford 2 3 Euclid Spin 6 4 Euclid Pin + 8 5 Clifford Spin 10 A 12 B 17 1 Clifford Spin D Euclid Clifford Γ µ, µ = 1,, D {Γ µ, Γ ν

More information

並列計算の数理とアルゴリズム サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.

並列計算の数理とアルゴリズム サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます.  このサンプルページの内容は, 初版 1 刷発行時のものです. 並列計算の数理とアルゴリズム サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/080711 このサンプルページの内容は, 初版 1 刷発行時のものです. Calcul scientifique parallèle by Frédéric Magoulès and François-Xavier

More information

JKR Point loading of an elastic half-space 2 3 Pressure applied to a circular region Boussinesq, n =

JKR Point loading of an elastic half-space 2 3 Pressure applied to a circular region Boussinesq, n = JKR 17 9 15 1 Point loading of an elastic half-space Pressure applied to a circular region 4.1 Boussinesq, n = 1.............................. 4. Hertz, n = 1.................................. 6 4 Hertz

More information

it-ken_open.key

it-ken_open.key 深層学習技術の進展 ImageNet Classification 画像認識 音声認識 自然言語処理 機械翻訳 深層学習技術は これらの分野において 特に圧倒的な強みを見せている Figure (Left) Eight ILSVRC-2010 test Deep images and the cited4: from: ``ImageNet Classification with Networks et

More information

ohpr.dvi

ohpr.dvi 2003/12/04 TASK PAF A. Fukuyama et al., Comp. Phys. Rep. 4(1986) 137 A. Fukuyama et al., Nucl. Fusion 26(1986) 151 TASK/WM MHD ψ θ ϕ ψ θ e 1 = ψ, e 2 = θ, e 3 = ϕ ϕ E = E 1 e 1 + E 2 e 2 + E 3 e 3 J :

More information

.2 ρ dv dt = ρk grad p + 3 η grad (divv) + η 2 v.3 divh = 0, rote + c H t = 0 dive = ρ, H = 0, E = ρ, roth c E t = c ρv E + H c t = 0 H c E t = c ρv T

.2 ρ dv dt = ρk grad p + 3 η grad (divv) + η 2 v.3 divh = 0, rote + c H t = 0 dive = ρ, H = 0, E = ρ, roth c E t = c ρv E + H c t = 0 H c E t = c ρv T NHK 204 2 0 203 2 24 ( ) 7 00 7 50 203 2 25 ( ) 7 00 7 50 203 2 26 ( ) 7 00 7 50 203 2 27 ( ) 7 00 7 50 I. ( ν R n 2 ) m 2 n m, R = e 2 8πε 0 hca B =.09737 0 7 m ( ν = ) λ a B = 4πε 0ħ 2 m e e 2 = 5.2977

More information

ʪ¼Á¤Î¥È¥Ý¥í¥¸¥«¥ë¸½¾Ý (2016ǯ¥Î¡¼¥Ù¥ë¾Þ¤Ë´ØÏ¢¤·¤Æ)

ʪ¼Á¤Î¥È¥Ý¥í¥¸¥«¥ë¸½¾Ý  (2016ǯ¥Î¡¼¥Ù¥ë¾Þ¤Ë´ØÏ¢¤·¤Æ) (2016 ) Dept. of Phys., Kyushu Univ. 2017 8 10 1 / 59 2016 Figure: D.J.Thouless F D.M.Haldane J.M.Kosterlitz TOPOLOGICAL PHASE TRANSITIONS AND TOPOLOGICAL PHASES OF MATTER 2 / 59 ( ) ( ) (Dirac, t Hooft-Polyakov)

More information

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63>

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63> 単純適応制御 SAC サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/091961 このサンプルページの内容は, 初版 1 刷発行当時のものです. 1 2 3 4 5 9 10 12 14 15 A B F 6 8 11 13 E 7 C D URL http://www.morikita.co.jp/support

More information

X G P G (X) G BG [X, BG] S 2 2 2 S 2 2 S 2 = { (x 1, x 2, x 3 ) R 3 x 2 1 + x 2 2 + x 2 3 = 1 } R 3 S 2 S 2 v x S 2 x x v(x) T x S 2 T x S 2 S 2 x T x S 2 = { ξ R 3 x ξ } R 3 T x S 2 S 2 x x T x S 2

More information

IPSJ SIG Technical Report Vol.2011-HPC-131 No /10/6 1 1 Parareal-in-Time Applicability of Time-domain Parallelism to Iterative Linear Calculus T

IPSJ SIG Technical Report Vol.2011-HPC-131 No /10/6 1 1 Parareal-in-Time Applicability of Time-domain Parallelism to Iterative Linear Calculus T 1 1 Parareal-in-Time Applicability of Time-domain Parallelism to Iterative Linear Calculus Toshiya Taami 1 and Aira Nishida 1 The time-domain parallelism, nown as Parareal-in-Time algorithm, has been applied

More information

IA hara@math.kyushu-u.ac.jp Last updated: January,......................................................................................................................................................................................

More information

x = a 1 f (a r, a + r) f(a) r a f f(a) 2 2. (a, b) 2 f (a, b) r f(a, b) r (a, b) f f(a, b)

x = a 1 f (a r, a + r) f(a) r a f f(a) 2 2. (a, b) 2 f (a, b) r f(a, b) r (a, b) f f(a, b) 2011 I 2 II III 17, 18, 19 7 7 1 2 2 2 1 2 1 1 1.1.............................. 2 1.2 : 1.................... 4 1.2.1 2............................... 5 1.3 : 2.................... 5 1.3.1 2.....................................

More information

: u i = (2) x i Smagorinsky τ ij τ [3] ij u i u j u i u j = 2ν SGS S ij, (3) ν SGS = (C s ) 2 S (4) x i a u i ρ p P T u ν τ ij S c ν SGS S csgs

: u i = (2) x i Smagorinsky τ ij τ [3] ij u i u j u i u j = 2ν SGS S ij, (3) ν SGS = (C s ) 2 S (4) x i a u i ρ p P T u ν τ ij S c ν SGS S csgs 15 C11-4 Numerical analysis of flame propagation in a combustor of an aircraft gas turbine, 4-6-1 E-mail: tominaga@icebeer.iis.u-tokyo.ac.jp, 2-11-16 E-mail: ntani@iis.u-tokyo.ac.jp, 4-6-1 E-mail: itoh@icebeer.iis.u-tokyo.ac.jp,

More information

untitled

untitled 18 1 2,000,000 2,000,000 2007 2 2 2008 3 31 (1) 6 JCOSSAR 2007pp.57-642007.6. LCC (1) (2) 2 10mm 1020 14 12 10 8 6 4 40,50,60 2 0 1998 27.5 1995 1960 40 1) 2) 3) LCC LCC LCC 1 1) Vol.42No.5pp.29-322004.5.

More information

01.Œk’ì/“²fi¡*

01.Œk’ì/“²fi¡* AIC AIC y n r n = logy n = logy n logy n ARCHEngle r n = σ n w n logσ n 2 = α + β w n 2 () r n = σ n w n logσ n 2 = α + β logσ n 2 + v n (2) w n r n logr n 2 = logσ n 2 + logw n 2 logσ n 2 = α +β logσ

More information

S I. dy fx x fx y fx + C 3 C dy fx 4 x, y dy v C xt y C v e kt k > xt yt gt [ v dt dt v e kt xt v e kt + C k x v + C C k xt v k 3 r r + dr e kt S dt d

S I. dy fx x fx y fx + C 3 C dy fx 4 x, y dy v C xt y C v e kt k > xt yt gt [ v dt dt v e kt xt v e kt + C k x v + C C k xt v k 3 r r + dr e kt S dt d S I.. http://ayapin.film.s.dendai.ac.jp/~matuda /TeX/lecture.html PDF PS.................................... 3.3.................... 9.4................5.............. 3 5. Laplace................. 5....

More information

211 kotaro@math.titech.ac.jp 1 R *1 n n R n *2 R n = {(x 1,..., x n ) x 1,..., x n R}. R R 2 R 3 R n R n R n D D R n *3 ) (x 1,..., x n ) f(x 1,..., x n ) f D *4 n 2 n = 1 ( ) 1 f D R n f : D R 1.1. (x,

More information

: 1g99p038-8

: 1g99p038-8 16 17 : 1g99p038-8 1 3 1.1....................................... 4 1................................... 5 1.3.................................. 5 6.1..................................... 7....................................

More information

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63>

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63> MATLAB/Simulink による現代制御入門 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/9241 このサンプルページの内容は, 初版 1 刷発行当時のものです. i MATLAB/Simulink MATLAB/Simulink 1. 1 2. 3. MATLAB/Simulink

More information

1 (Contents) (1) Beginning of the Universe, Dark Energy and Dark Matter Noboru NAKANISHI 2 2. Problem of Heat Exchanger (1) Kenji

1 (Contents) (1) Beginning of the Universe, Dark Energy and Dark Matter Noboru NAKANISHI 2 2. Problem of Heat Exchanger (1) Kenji 8 4 2018 6 2018 6 7 1 (Contents) 1. 2 2. (1) 22 3. 31 1. Beginning of the Universe, Dark Energy and Dark Matter Noboru NAKANISHI 2 2. Problem of Heat Exchanger (1) Kenji SETO 22 3. Editorial Comments Tadashi

More information

2005 2006.2.22-1 - 1 Fig. 1 2005 2006.2.22-2 - Element-Free Galerkin Method (EFGM) Meshless Local Petrov-Galerkin Method (MLPGM) 2005 2006.2.22-3 - 2 MLS u h (x) 1 p T (x) = [1, x, y]. (1) φ(x) 0.5 φ(x)

More information

φ 4 Minimal subtraction scheme 2-loop ε 2008 (University of Tokyo) (Atsuo Kuniba) version 21/Apr/ Formulas Γ( n + ɛ) = ( 1)n (1 n! ɛ + ψ(n + 1)

φ 4 Minimal subtraction scheme 2-loop ε 2008 (University of Tokyo) (Atsuo Kuniba) version 21/Apr/ Formulas Γ( n + ɛ) = ( 1)n (1 n! ɛ + ψ(n + 1) φ 4 Minimal subtraction scheme 2-loop ε 28 University of Tokyo Atsuo Kuniba version 2/Apr/28 Formulas Γ n + ɛ = n n! ɛ + ψn + + Oɛ n =,, 2, ψn + = + 2 + + γ, 2 n ψ = γ =.5772... Euler const, log + ax x

More information

* 1 1 (i) (ii) Brückner-Hartree-Fock (iii) (HF, BCS, HFB) (iv) (TDHF,TDHFB) (RPA) (QRPA) (v) (vi) *

* 1 1 (i) (ii) Brückner-Hartree-Fock (iii) (HF, BCS, HFB) (iv) (TDHF,TDHFB) (RPA) (QRPA) (v) (vi) * * 1 1 (i) (ii) Brückner-Hartree-Fock (iii) (HF, BCS, HFB) (iv) (TDHF,TDHFB) (RPA) (QRPA) (v) (vi) *1 2004 1 1 ( ) ( ) 1.1 140 MeV 1.2 ( ) ( ) 1.3 2.6 10 8 s 7.6 10 17 s? Λ 2.5 10 10 s 6 10 24 s 1.4 ( m

More information

xia2.dvi

xia2.dvi Journal of Differential Equations 96 (992), 70-84 Melnikov method and transversal homoclinic points in the restricted three-body problem Zhihong Xia Department of Mathematics, Harvard University Cambridge,

More information

TOP URL 1

TOP URL   1 TOP URL http://amonphys.web.fc.com/ 1 19 3 19.1................... 3 19.............................. 4 19.3............................... 6 19.4.............................. 8 19.5.............................

More information

/02/18

/02/18 3 09/0/8 i III,,,, III,?,,,,,,,,,,,,,,,,,,,,?,?,,,,,,,,,,,,,,!!!,? 3,,,, ii,,,!,,,, OK! :!,,,, :!,,,,,, 3:!,, 4:!,,,, 5:!,,! 7:!,,,,, 8:!,! 9:!,,,,,,,,, ( ),, :, ( ), ( ), 6:!,,, :... : 3 ( )... iii,,

More information

January 27, 2015

January 27, 2015 e-mail : kigami@i.kyoto-u.ac.jp January 27, 205 Contents 2........................ 2.2....................... 3.3....................... 6.4......................... 2 6 2........................... 6

More information

基礎数学I

基礎数学I I & II ii ii........... 22................. 25 12............... 28.................. 28.................... 31............. 32.................. 34 3 1 9.................... 1....................... 1............

More information

( ; ) C. H. Scholz, The Mechanics of Earthquakes and Faulting : - ( ) σ = σ t sin 2π(r a) λ dσ d(r a) =

( ; ) C. H. Scholz, The Mechanics of Earthquakes and Faulting : - ( ) σ = σ t sin 2π(r a) λ dσ d(r a) = 1 9 8 1 1 1 ; 1 11 16 C. H. Scholz, The Mechanics of Earthquakes and Faulting 1. 1.1 1.1.1 : - σ = σ t sin πr a λ dσ dr a = E a = π λ σ πr a t cos λ 1 r a/λ 1 cos 1 E: σ t = Eλ πa a λ E/π γ : λ/ 3 γ =

More information

Part () () Γ Part ,

Part () () Γ Part , Contents a 6 6 6 6 6 6 6 7 7. 8.. 8.. 8.3. 8 Part. 9. 9.. 9.. 3. 3.. 3.. 3 4. 5 4.. 5 4.. 9 4.3. 3 Part. 6 5. () 6 5.. () 7 5.. 9 5.3. Γ 3 6. 3 6.. 3 6.. 3 6.3. 33 Part 3. 34 7. 34 7.. 34 7.. 34 8. 35

More information

i 18 2H 2 + O 2 2H 2 + ( ) 3K

i 18 2H 2 + O 2 2H 2 + ( ) 3K i 18 2H 2 + O 2 2H 2 + ( ) 3K ii 1 1 1.1.................................. 1 1.2........................................ 3 1.3......................................... 3 1.4....................................

More information

COE-RES Discussion Paper Series Center of Excellence Project The Normative Evaluation and Social Choice of Contemporary Economic Systems Graduate Scho

COE-RES Discussion Paper Series Center of Excellence Project The Normative Evaluation and Social Choice of Contemporary Economic Systems Graduate Scho COE-RES Discussion Paper Series Center of Excellence Project The Normative Evaluation and Social Choice of Contemporary Economic Systems Graduate School of Economics and Institute of Economic Research

More information

2016 ǯ¥Î¡¼¥Ù¥ëʪÍý³Ø¾Þ²òÀ⥻¥ß¥Ê¡¼ Kosterlitz-Thouless ž°Ü¤È Haldane ͽÁÛ

2016 ǯ¥Î¡¼¥Ù¥ëʪÍý³Ø¾Þ²òÀ⥻¥ß¥Ê¡¼  Kosterlitz-Thouless ž°Ü¤È Haldane ͽÁÛ 2016 Kosterlitz-Thouless Haldane Dept. of Phys., Kyushu Univ. 2016 11 29 2016 Figure: D.J.Thouless F D.M.Haldane J.M.Kosterlitz TOPOLOGICAL PHASE TRANSITIONS AND TOPOLOGICAL PHASES OF MATTER ( ) ( ) (Dirac,

More information

スライド 1

スライド 1 GPU クラスタによる格子 QCD 計算 広大理尾崎裕介 石川健一 1.1 Introduction Graphic Processing Units 1 チップに数百個の演算器 多数の演算器による並列計算 ~TFLOPS ( 単精度 ) CPU 数十 GFLOPS バンド幅 ~100GB/s コストパフォーマンス ~$400 GPU の開発環境 NVIDIA CUDA http://www.nvidia.co.jp/object/cuda_home_new_jp.html

More information

1/2 ( ) 1 * 1 2/3 *2 up charm top -1/3 down strange bottom 6 (ν e, ν µ, ν τ ) -1 (e) (µ) (τ) 6 ( 2 ) 6 6 I II III u d ν e e c s ν µ µ t b ν τ τ (2a) (

1/2 ( ) 1 * 1 2/3 *2 up charm top -1/3 down strange bottom 6 (ν e, ν µ, ν τ ) -1 (e) (µ) (τ) 6 ( 2 ) 6 6 I II III u d ν e e c s ν µ µ t b ν τ τ (2a) ( August 26, 2005 1 1 1.1...................................... 1 1.2......................... 4 1.3....................... 5 1.4.............. 7 1.5.................... 8 1.6 GIM..........................

More information

三石貴志.indd

三石貴志.indd 流通科学大学論集 - 経済 情報 政策編 - 第 21 巻第 1 号,23-33(2012) SIRMs SIRMs Fuzzy fuzzyapproximate approximatereasoning reasoningusing using Lukasiewicz Łukasiewicz logical Logical operations Operations Takashi Mitsuishi

More information

keisoku01.dvi

keisoku01.dvi 2.,, Mon, 2006, 401, SAGA, JAPAN Dept. of Mechanical Engineering, Saga Univ., JAPAN 4 Mon, 2006, 401, SAGA, JAPAN Dept. of Mechanical Engineering, Saga Univ., JAPAN 5 Mon, 2006, 401, SAGA, JAPAN Dept.

More information