Fortran90/95 2. (p 74) f g h x y z f x h x = f x + g x h y = f y + g y h z = f z + g z f x f y f y f h = f + g Fortran 1 3 a b c c(1) = a(1) + b(1) c(

Size: px
Start display at page:

Download "Fortran90/95 2. (p 74) f g h x y z f x h x = f x + g x h y = f y + g y h z = f z + g z f x f y f y f h = f + g Fortran 1 3 a b c c(1) = a(1) + b(1) c("

Transcription

1 Fortran90/ n n = 5 x1,x2,x3,,x4,x5 5 average = ( x1 + x2 + x3 + x4 + x5 )/5.0 n n x (subscript) x 1 x 2 average = 1 n n x i i=1 Fortran ( ) x(1) x(2) x(n) Fortran ( ) average = sum(x(1:n))/real(n) x x(1) x(2) (p 61) x 2 y 5 z 3 z(3) = real(x(2)) + y(5) (extent) 3 x(1) x(2) x(3) x 3 (p 68) Fortran a i j a(i,j) (dimension) Fortran 7 (size)

2 Fortran90/95 2. (p 74) f g h x y z f x h x = f x + g x h y = f y + g y h z = f z + g z f x f y f y f h = f + g Fortran 1 3 a b c c(1) = a(1) + b(1) c(2) = a(2) + b(2) c(3) = a(3) + b(3) c = a + b (p 76) 3 a b b = abs(a) Fortran b(1)=abs(a(1)) b(2)=abs(a(2)) b(3)=abs(a(3)) sum ( )

3 Fortran90/95 Fortran c c + 1 c 1 3. FORTRAN FORTRAN (array) (dimension) (tensor) (order) 0 = (scalar) 1 1 = (vector) 2 2 = (matrix) (extent) (size) (dimension) FORTRAN (i) 1 {x} x FORTRAN 1 1 y z z = y n 1 n FORTRAN (ii) (p 71) Fortran(90/95)

4 Fortran90/ n x n 10 Fortran90/95 :: ( ) integer real Fortran program example4_01 3 real :: x(10) 4 real :: average, variance 5 integer :: n, i 6 7 average = read(*,*) n 9 do i = 1, n 10 read(*,*) x(i) 11 average = average + x(i) 12 enddo 13 average = average/real(n) variance = do i = 1, n 17 variance = variance + (x(i)-average)**2 18 enddo 19 variance = variance/real(n) write(*,*) "average=", average, " variance=", variance stop 24 end program example4_

5 Fortran90/ n 2,1,0,-1,-2 n 1 1 n select case -2 2 ( ) ( : ) : program example4_02 3 integer :: k(-2:2) 4 integer :: i, j, n 5 real :: average 6 7 do i = -2, 2 8 k(i) = 0 9 enddo 10 average = read(*,*) n 13 do i = 1, n 14 read(*,*) j 15 k(j) = k(j) average= average + real(j) 17 enddo do i = -2, 2 20 write(*,*) i, " ->",k(i) 21 enddo 22 write(*,*) "average=", average/real(n) stop 25 end program example4_

6 Fortran90/ (, ) read(*,*) write(*,*) x 1 5 read(*,*) x read(*,*) x(1), x(2), x(3), x(4), x(5) program example4_03 3 integer :: ip(10), max_point, min_point 4 real :: average 5 integer :: i 6 7 read(*,*) ip 8 9 max_point = 0 10 min_point = average = do i = 1, if(ip(i) > max_point) max_point = ip(i) 14 if(ip(i) < min_point) min_point = ip(i) 15 average = average + real(ip(i)) 16 enddo write(*,*) "max point=",max_point 19 write(*,*) "min point=",min_point 20 write(*,*) "average= ", average/ stop 23 end program example4_

7 Fortran90/95 do n 1 1 do read(*,*) ( ( ), =,, ) 1, do (i) (ii) + (iii) + 2. do program example4_04 3 real :: x(10) 4 real :: average, variance 5 integer :: n, i 6 7 average = read(*,*) n, (x(i), i = 1, n) 9 10 do i = 1, n 11 average = average + x(i) 12 enddo 13 average = average/real(n) variance = do i = 1, n 17 variance = variance + (x(i)-average)**2 18 enddo 19 variance = variance/real(n) 20 write(*,*) "average=", average, " variance=", variance stop 23 end program example4_04 8 n do

8 Fortran90/ n, (/ /) (/ 1, 2,, n /) (/ /) program example4_05 3 integer :: month_day(12)=(/31,28,31,30,31,30,31,31,30,31,30,31/) 4 integer :: month, day, day_year, i 5 6 read(*,*) month, day 7 8 day_year = 0 9 do i = 1, month-1 10 day_year = day_year + month_day(i) 11 enddo 12 day_year = day_year + day write(*,*) "day of year=", day_year stop 17 end program example4_05 3 (/ 1, 2, /)

9 Fortran90/95 (i) ( ) (/ 2.2, x, 1.2 /) x (ii) ( ) do (/ (i, i=1,9,2) /) 5 i (iii) (/ a, b /) a b do (i ) ( ) real, parameter :: x = 1.7 real :: a(3) = (/ 2.2, x, 1.2 /) (ii ) ( ) do integer :: i integer :: k(5) = (/ (i, i=1,9,2) /) (iii ) integer, parameter :: ic(2) = (/ 1, 2 /) integer :: id(4) = (/ ic, (/ 3, 4 /) /) program example4_05a implicit none integer :: month_day(12)=(/31,28,31,30,31,30,31,31,30,31,30,31/) integer :: month, day read(*,*) day do month = 1, 12 if( day <= month_day(month) ) exit day = day - month_day(month) enddo write(*,*) "month=", month, " day=", day stop end program example4_05a

10 Fortran90/95 4. data 4-6 C H O C= 12 H= 1 O= 16 (i) read (ii) 3 (iii) data *1 data data 1 / 1 /, 2 / 2 /,, (a(i),i=5,10) data 1,2,0,0,0,0,3 * 1,2,4*0,3 data end program example4_06 3 real :: w(3), w_mole 4 integer :: k(3), i 5 data w/12.0, 1.0, 16.0/ 6 7 read(*,*) k 8 w_mole = do i = 1, 3 10 w_mole = w_mole + w(i)*real(k(i)) 11 enddo 12 write(*,*) "molecular weight=", w_mole stop 15 end program example4_06 *1-67 -

11 Fortran90/ A B C D Excel 1 :: ( 1, 2 ) ,75,78 ip(1,1) ip(1,2) ip(1,3) 65,90,70 -> ip(2,1) ip(2,2) ip(2,3) 70,73,85 ip(3,1) ip(3,2) ip(3,3) 66,72,71 ip(4,1) ip(4,2) ip(4,3) program example4_07 3 integer :: ip(4,3) 4 real :: ave_subject(3) 5 integer :: i, j 6 7 do i = 1, 4 8 read(*,*) (ip(i,j), j= 1, 3) 9 enddo do j = 1, 3 12 ave_subject(j) = do i = 1, 4 14 ave_subject(j) = ave_subject(j) + ip(i,j) 15 enddo 16 ave_subject(j) = ave_subject(j)/ enddo write(*,*) (ave_subject(j), j = 1, 3) stop 22 end program example4_

12 Fortran90/95 Fortran 1 1 ( ) 2 * ip * 3 ip(1,1), ip(2,1), ip(3,1), ip(4,1), ip(1,2),, ip(3,3), ip(4,3) read(*,*) ip write(*,*) ip Fortran (1) data data data integer :: ip(4,3) data ip/ 82, 65, 70, 66, 75, 90, 73, 72, 78, 70, 85, 71 / do data ((ip(i,j),j=1,3),i=1,4) / 82, 75, 78, 65, 90, 70, 70, 73, 85, 66, 72, 71 / (2) reshape reshape reshape( source, shape ) source shape (/ /) ( ) ( ) ip(4,3) = ip = reshape((/ 82, 65, 70, 66, 75, 90, 73, 72, 78, 70, 85, 71 /), (/ 4, 3 /)) *2 C *3-69 -

13 Fortran90/ ,2,4,7,8,10, n 3 1 1,4,7, 2 2,5,8, 3 3,6,9, do a a(n1:n2[:n3)] n1 a n2 a n3 n3=1 n3 1 do a n2 n1 n2 n3 (i) n1=1 1 (ii) n2 a n2 (iii) n3 1 :1 : 1 a(:) program example4_08 3 integer, parameter :: m = 14 4 integer :: n(m) 5 6 read(*,*) n 7 8 write(*,*) n(1:m:3)! n(::3) 9 write(*,*) n(2:m:3) 10 write(*,*) n(3:m:3) stop 13 end program example4_08 do

14 Fortran90/95 (n1,n2[,n3]) (i) (ii) (iii) a(4,3) a(1,1) a(1,2) a(1,3) a(2,1) a(2,2) a(2,3) a(3,1) a(3,2) a(3,3) a(4,1) a(4,2) a(4,3) a(2:4:2,:) a(2,1) a(2,2) a(2,3) a(4,1) a(4,2) a(4,3) a(3,:) a(3,1) a(3,2) a(3,3) a(n1) a(n1:n1) 1 a(n1,n2) 2 a(n1:n1,n2:n2) 1 a(2,2) = a(2:2,2:2) a(n1,n2) n1>n2 (size) n a m num a num a(num) a(num(1)),a(num(2))...a(num(m)) m a write(*,*) (a(num(i)), i=1, m) do write(*,*) a(num(:)) num num a read

15 Fortran90/ n x n 5 81, 72, 70, 64, n ( x(10000) ) (i) (ii) (iii) n n step1: allocatable, allocatable :: (:) (:,:) step2: allocate allocate ( ( )) n allocate deallocate step3: deallocate deallocate ( ) deallocate deallocate allocate

16 Fortran90/ program example4_09 3 real, allocatable :: x(:) 4 real :: sum 5 integer :: n, i 6 7 read(*,*) n 8 allocate (x(n)) 9 read(*,*) x(:) 10 sum = do i = 1, n 12 sum = sum + x(i) 13 enddo 14 write(*,*) "average =", sum/real(n) deallocate (x) 17 stop 18 end program example4_ do 9 read(*,*) x x n n program example4_09a implicit none integer, allocatable :: ip(:,:) real :: average integer :: i, j, n read(*,*) n; allocate (ip(n,3)) read(*,*) ((ip(i,j), j =1, 3), i =1, n) do j = 1, 3 average = 0.0 do i = 1, n average = average + ip(i,j) enddo write(*,*) average/real(n) enddo deallocate (ip) stop end program example4_09a

17 Fortran90/ a = (1) a + b (2) 2 a 3 b (3) a b ( ) 1.2 b = n Fortran n a b data read a b c ( ) c = a + b ( ) c = a - b ( ) c = a * b ( ) c = a / b *4 c(i) = a(i) + b(i) a b c x a y y = a * x y = x * a x y *4 c(:) = a(:) + b(:) c(:,:) = a(:,:) + b(:,:) (2 )

18 Fortran90/95 d e f f = d ** e program example4_10 3 real :: a(3)=(/2.0,1.5,3.2/), b(3)=(/1.2,-1.1,0.6/) 4 real :: c(3), sum 5 integer :: i 6 7 write(*,*) a(:) + b(:) 8 write(*,*) 2.0*a(:) - 3.0*b(:) 9 c(:) = a(:) * b(:) 10 sum = do i = 1, 3 12 sum = sum + c(i) 13 enddo 14 write(*,*) sum stop 17 end program example4_ program example4_10a implicit none real :: a(2,2), b(2,2), c(2,2) data a/1.0,3.0,2.0,4.0/! data b/2.5,0.5,1.0,-1.5/! c(:,:) = a(:,:) + b(:,:) write(*,*) c(1,1), c(1,2) write(*,*) c(2,1), c(2,2) stop end program example4_10a

19 Fortran90/ a = do real :: a(3)=(/1.5,0.5,-0.9/) integer :: b(3), i do i=1,3 b(i) = nint((a(i)) enddo b 1 ( nint) 2 = ( 1) 1 2 ( ) program example4_11 3 real :: a(3)=(/1.5,0.5,-0.9/) 4 integer :: b(3) 5 6 b(:) = nint(a(:)) 7 write(*,*) b(:) 8 9 stop 10 end program example4_

20 Fortran90/ where where where where Fortran90 if 4-12 n a 0.0 n a = (/2.5,1.2,-1.2,1.1,-0.5/) where where if where( mask ) where ().true.( ).false.( ) (mask ) ( a(:) < 1.0 ) where ( where ) where if * 5 where( ) 1 2 end where program example4_12a 3 real :: a(5)=(/2.5,1.2,-1.2,1.1,-0.5/) 4 5 where( a(:)< 0.0 ) a(:) = 0.0 *5 elesewhere( ) elsewhere if elseif else

21 Fortran90/95 6 write(*,*) a(:) 7 8 stop 9 end program example4_12a 2. forall forall (95) forall Fortran95 where do where forall forall forall forall(do 1 [ do 2 ]..[ ]) do = : : forall(i=1:n, a(i)>0) a(i)=1.0/a(i) forall forall(do 1 [ do 2 ]..[ ]) 1 2 end forall do loop do loop forall CPU forall program example4_12b 3 real :: a(5)=(/2.5,1.2,-1.2,1.1,-0.5/) 4 integer :: i 5 6 forall(i=1:5, a(i)<0.0) a(i) = write(*,*) a(:) 8 9 stop 10 end program example4_12b

22 Fortran90/ (90/95) fortran90/ Fortran(90/95) (1) 2 2 a = 1 1 b = 1 a b 3 2 (2) C = 1 2 C 2 C t C (C t C ) 0 1 (dot_product) u v (u,v) dot_product(u,v) u v 1 u v u v u v u v u v u u (conjugate) u(i),v(i) u v u(i).and. v(i) program example4_13a 2! 3 implicit none 4 integer :: a(3) = (/2, 1, 3/) 5 integer :: b(3) = (/1,-1, 2/) 6 integer :: p 7 8 p = dot_product(a,b) 9 write(*,*) p stop 12 end program example4_13a i

23 Fortran90/95 (matmul ) A B C=AB c = matmul(a,b) a k m 2 b m n 2 c k n 2 a b a b a b 2 b m 1 u v=au v = matmul(a,u) v k 1 2 a k 1 v u = t va v = matmul(v,a) u m 1 ( transpose ) A *6 C = t A c = transpose(a) a k m 2 m k program example4_13b 2! 3 implicit none 4 integer :: c(2,2), d(2,2), e(2,2) 5 integer :: i, j 6 data c/1,0,2,1/ 7 8 d = matmul(c,c) 9 write(*,*) d(1,:) ; write(*,*) d(2,:) ; write(*,*) e = matmul(transpose(c),c) 12 write(*,*) e(1,:) ; write(*,*) e(2,:) stop 15 end program example4_13b *6 C = t A A i j = C ji

24 Fortran90/ A n n m real :: a(2,-1:3,0:2) (shape) A shape(a) 1 A i A i (size) A size(a [,dim]) dim A dim dim A (ubound) A ubound(a [,dim]) dim A dim dim A A (lbound) A ubound lbound(a [,dim]) program example4_14 3 real :: a(2,-1:3,0:2) 4 integer :: n, m(7)! 7 5 n = size(shape(a)) 6 m(1:n) = shape(a) 7 write(*,*) n 8 write(*,*) m(1:n) 9 stop 10 end program example4_

25 Fortran90/ A (1) n (2) ave real :: a(7) = (/1.2, 2.1, 3.0, 999.0, 2.0, 1.1, 999.0/) sum if FORTRAN77 Fortran90/95 mask mask A sum(a,mask=(a>0)) mask mask= mask mask ( A ) 1 1 (sum) sum(a[,dim][,mask]) A A ( ) dim dim A -1 mask mask (product) product(a[,dim][,mask]) sum (maxval) maxval(a[,dim][,mask]) A A ( ) dim dim A -1 mask mask

26 Fortran90/95 (minval) minval(a[,dim][,mask]) maxval mask mask= (count) count(mask[,dim]) mask ( ) dim dim dim mask -1 (any) any(mask[,dim]) mask 1 ( ) dim dim 1 dim mask -1 (all) all(mask[,dim]) mask ( ) dim dim dim mask program example4_15 3 real :: a(7) = (/1.2, 2.1, 3.0, 999.0, 2.0, 1.1, 999.0/), ave 4 integer :: n 5 6 n = count(a<999.0) 7 write(*,*) n 8 ave = sum(a, mask=(a<999.0) ) 9 write(*,*) ave/real(n) stop 12 end program example4_

27 Fortran90/ A real :: a(-3:3) = (/1.2, 2.1, 3.0, 999.0, 2.0, 1.1, 999.0/) 1 (maxloc) maxloc(a[,mask]) A 1 mask mask (minloc) minloc(a[,mask]) maxloc program example4_16 3 real :: a(-3:3) = (/1.2, 2.1, 3.0, 999.0, 2.0, 1.1, 999.0/) 4 integer :: m(1) 5 6 m = maxloc(a, mask=(a<999.0)) 7 write(*,*) lbound(a,1) + m(1) stop 10 end program example4_1-84 -

28 Fortran90/ ( ) 1000 ( ) (1) 2 (2) 2 (3) 3 2 (4) ( ) 6 1,2,3, (10000 ) 4-3 ( ) 4-4 ( ) 5 25, 12, 50, 35, 29 ( ) A = A 2 (= AA) t AA a = 2.0 b = 0.0 a b a = 2.0 b = 1.0 (a b) = t ab a t b matmul

1F90/kouhou_hf90.dvi

1F90/kouhou_hf90.dvi Fortran90 3 33 1 2 Fortran90 FORTRAN 1956 IBM IBM704 FORTRAN(FORmula TRANslation ) 1965 FORTRAN66 1978 FORTRAN77 1991 Fortran90 Fortran90 Fortran Fortran90 6 Fortran90 77 90 90 Fortran90 [ ] Fortran90

More information

untitled

untitled Fortran90 ( ) 17 12 29 1 Fortran90 Fortran90 FORTRAN77 Fortran90 1 Fortran90 module 1.1 Windows Windows UNIX Cygwin (http://www.cygwin.com) C\: Install Cygwin f77 emacs latex ps2eps dvips Fortran90 Intel

More information

11042 計算機言語7回目 サポートページ:

11042 計算機言語7回目  サポートページ: 11042 7 :https://goo.gl/678wgm November 27, 2017 10/2 1(print, ) 10/16 2(2, ) 10/23 (3 ) 10/31( ),11/6 (4 ) 11/13,, 1 (5 6 ) 11/20,, 2 (5 6 ) 11/27 (7 12/4 (9 ) 12/11 1 (10 ) 12/18 2 (10 ) 12/25 3 (11

More information

all.dvi

all.dvi fortran 1996 4 18 2007 6 11 2012 11 12 1 3 1.1..................................... 3 1.2.............................. 3 2 fortran I 5 2.1 write................................ 5 2.2.................................

More information

(Basic Theory of Information Processing) Fortran Fortan Fortan Fortan 1

(Basic Theory of Information Processing) Fortran Fortan Fortan Fortan 1 (Basic Theory of Information Processing) Fortran Fortan Fortan Fortan 1 17 Fortran Formular Tranlator Lapack Fortran FORTRAN, FORTRAN66, FORTRAN77, FORTRAN90, FORTRAN95 17.1 A Z ( ) 0 9, _, =, +, -, *,

More information

情報処理概論(第二日目)

情報処理概論(第二日目) 情報処理概論 工学部物質科学工学科応用化学コース機能物質化学クラス 第 8 回 2005 年 6 月 9 日 前回の演習の解答例 多項式の計算 ( 前半 ): program poly implicit none integer, parameter :: number = 5 real(8), dimension(0:number) :: a real(8) :: x, total integer

More information

インテル(R) Visual Fortran Composer XE 2013 Windows版 入門ガイド

インテル(R) Visual Fortran Composer XE 2013 Windows版 入門ガイド Visual Fortran Composer XE 2013 Windows* エクセルソフト株式会社 www.xlsoft.com Rev. 1.1 (2012/12/10) Copyright 1998-2013 XLsoft Corporation. All Rights Reserved. 1 / 53 ... 3... 4... 4... 5 Visual Studio... 9...

More information

( ) 1 Windows HTML ( ) ( ) ( ) WWW 10 ( )

( ) 1 Windows HTML ( ) ( ) ( ) WWW 10 ( ) ( ) 1 Windows HTML ( ) ( ) ( ) 1. 2. 3. 4. WWW 10 ( ) 2 1. 2. 1 3. ( ) 4. 5. 3 Windows 2 7 8 MS Word MS Excel 1. MS Word 600 2. MS Excel 1 34 2 83 3 23 4 70 5 100 6 45 7 58 8 29 9 68 10 72 11 37 12 12

More information

01_OpenMP_osx.indd

01_OpenMP_osx.indd OpenMP* / 1 1... 2 2... 3 3... 5 4... 7 5... 9 5.1... 9 5.2 OpenMP* API... 13 6... 17 7... 19 / 4 1 2 C/C++ OpenMP* 3 Fortran OpenMP* 4 PC 1 1 9.0 Linux* Windows* Xeon Itanium OS 1 2 2 WEB OS OS OS 1 OS

More information

情報科学概論 第1回資料

情報科学概論 第1回資料 1. Excel (C) Hiroshi Pen Fujimori 1 2. (Excel) 2.1 Excel : 2.2Excel Excel (C) Hiroshi Pen Fujimori 2 256 (IV) :C (C 65536 B4 :2 (2 A3 Excel (C) Hiroshi Pen Fujimori 3 Tips: (1) B3 (2) (*1) (3) (4)Word

More information

2 Excel =sum( ) =average( ) B15:D20 : $E$26 E26 $ =A26*$E$26 $ $E26 E$26 E$26 $G34 $ E26 F4

2 Excel =sum( ) =average( ) B15:D20 : $E$26 E26 $ =A26*$E$26 $ $E26 E$26 E$26 $G34 $ E26 F4 1234567 0.1234567 = 2 3 =2+3 =2-3 =2*3 =2/3 =2^3 1:^, 2:*/, 3:+- () =2+3*4 =(2+3)*4 =3*2^2 =(3*2)^2 =(3+6)^0.5 A12 =A12+B12 ( ) ( )0.4 ( 100)0.9 % 1 2 Excel =sum( ) =average( ) B15:D20 : $E$26 E26 $ =A26*$E$26

More information

. (.8.). t + t m ü(t + t) + c u(t + t) + k u(t + t) = f(t + t) () m ü f. () c u k u t + t u Taylor t 3 u(t + t) = u(t) + t! u(t) + ( t)! = u(t) + t u(

. (.8.). t + t m ü(t + t) + c u(t + t) + k u(t + t) = f(t + t) () m ü f. () c u k u t + t u Taylor t 3 u(t + t) = u(t) + t! u(t) + ( t)! = u(t) + t u( 3 8. (.8.)............................................................................................3.............................................4 Nermark β..........................................

More information

プラズマ核融合学会誌5月号【81-5】/内外情報_ソフト【注:欧フォント特殊!】

プラズマ核融合学会誌5月号【81-5】/内外情報_ソフト【注:欧フォント特殊!】 PROGRAM PLOTDATA USE NUM_KINDS, ONLY : wp=>dp, i4b USE MYLIB, ONLY : GET_SIZE, GET_DATA INTEGER(i4b) :: ntime, nx REAL(wp), ALLOCATABLE :: time(:), x(:), Temp(:,:) Fortran Temp, temp, TEMP temporal REAL(wp)

More information

... 3... 3... 3... 3... 4... 7... 10... 10... 11... 12... 12... 13... 14... 15... 18... 19... 20... 22... 22... 23 2

... 3... 3... 3... 3... 4... 7... 10... 10... 11... 12... 12... 13... 14... 15... 18... 19... 20... 22... 22... 23 2 1 ... 3... 3... 3... 3... 4... 7... 10... 10... 11... 12... 12... 13... 14... 15... 18... 19... 20... 22... 22... 23 2 3 4 5 6 7 8 9 Excel2007 10 Excel2007 11 12 13 - 14 15 16 17 18 19 20 21 22 Excel2007

More information

コンピュータ概論

コンピュータ概論 4.1 For Check Point 1. For 2. 4.1.1 For (For) For = To Step (Next) 4.1.1 Next 4.1.1 4.1.2 1 i 10 For Next Cells(i,1) Cells(1, 1) Cells(2, 1) Cells(10, 1) 4.1.2 50 1. 2 1 10 3. 0 360 10 sin() 4.1.2 For

More information

<4D F736F F F696E74202D D F95C097F D834F E F93FC96E5284D F96E291E85F8DE391E52E >

<4D F736F F F696E74202D D F95C097F D834F E F93FC96E5284D F96E291E85F8DE391E52E > SX-ACE 並列プログラミング入門 (MPI) ( 演習補足資料 ) 大阪大学サイバーメディアセンター日本電気株式会社 演習問題の構成 ディレクトリ構成 MPI/ -- practice_1 演習問題 1 -- practice_2 演習問題 2 -- practice_3 演習問題 3 -- practice_4 演習問題 4 -- practice_5 演習問題 5 -- practice_6

More information

Fortran90/95 [9]! (1 ) " " 5 "Hello!"! 3. (line) Fortran Fortran 1 2 * (1 ) 132 ( ) * 2 ( Fortran ) Fortran ,6 (continuation line) 1

Fortran90/95 [9]! (1 )   5 Hello!! 3. (line) Fortran Fortran 1 2 * (1 ) 132 ( ) * 2 ( Fortran ) Fortran ,6 (continuation line) 1 Fortran90/95 2.1 Fortran 2-1 Hello! 1 program example2_01! end program 2! first test program ( ) 3 implicit none! 4 5 write(*,*) "Hello!"! write Hello! 6 7 stop! 8 end program example2_01 1 program 1!

More information

untitled

untitled I 9 MPI (II) 2012 6 14 .. MPI. 1-3 sum100.f90 4 istart=myrank*25+1 iend=(myrank+1)*25 0 1 2 3 mpi_recv 3 isum1 1 isum /tmp/120614/sum100_4.f90 program sum100_4 use mpi implicit none integer :: i,istart,iend,isum,isum1,ip

More information

(2-1) x, m, 2 N(m, 2 ) x REAL*8 FUNCTION NRMDST (X, M, V) X,M,V REAL*8 x, m, 2 X X N(0,1) f(x) standard-norm.txt normdist1.f x=0, 0.31, 0.5

(2-1) x, m, 2 N(m, 2 ) x REAL*8 FUNCTION NRMDST (X, M, V) X,M,V REAL*8 x, m, 2 X X N(0,1) f(x) standard-norm.txt normdist1.f x=0, 0.31, 0.5 2007/5/14 II II agata@k.u-tokyo.a.jp 0. 1. x i x i 1 x i x i x i x x+dx f(x)dx f(x) f(x) + 0 f ( x) dx = 1 (Probability Density Funtion 2 ) (normal distribution) 3 1 2 2 ( x m) / 2σ f ( x) = e 2πσ x m

More information

Microsoft Word - 資料 docx

Microsoft Word - 資料 docx y = Asin 2πt T t t = t i i 1 n+1 i i+1 Δt t t i = Δt i 1 ( ) y i = Asin 2πt i T 29 (x, y) t ( ) x = Asin 2πmt y = Asin( 2πnt + δ ) m, n δ (x, y) m, n 30 L A x y A L x 31 plot sine curve program sine implicit

More information

OpenMP¤òÍѤ¤¤¿ÊÂÎó·×»»¡Ê£±¡Ë

OpenMP¤òÍѤ¤¤¿ÊÂÎó·×»»¡Ê£±¡Ë 2012 5 24 scalar Open MP Hello World Do (omp do) (omp workshare) (shared, private) π (reduction) PU PU PU 2 16 OpenMP FORTRAN/C/C++ MPI OpenMP 1997 FORTRAN Ver. 1.0 API 1998 C/C++ Ver. 1.0 API 2000 FORTRAN

More information

OpenMP¤òÍѤ¤¤¿ÊÂÎó·×»»¡Ê£±¡Ë

OpenMP¤òÍѤ¤¤¿ÊÂÎó·×»»¡Ê£±¡Ë 2011 5 26 scalar Open MP Hello World Do (omp do) (omp workshare) (shared, private) π (reduction) scalar magny-cours, 48 scalar scalar 1 % scp. ssh / authorized keys 133. 30. 112. 246 2 48 % ssh 133.30.112.246

More information

●70974_100_AC009160_KAPヘ<3099>ーシス自動車約款(11.10).indb

●70974_100_AC009160_KAPヘ<3099>ーシス自動車約款(11.10).indb " # $ % & ' ( ) * +, -. / 0 1 2 3 4 5 6 7 8 9 : ; < = >? @ A B C D E F G H I J K L M N O P Q R S T U V W X Y " # $ % & ' ( ) * + , -. / 0 1 2 3 4 5 6 7 8 9 : ; < = > ? @ A B

More information

Microsoft Word - 03-数値計算の基礎.docx

Microsoft Word - 03-数値計算の基礎.docx δx f x 0 + δ x n=0 a n = f ( n) ( x 0 ) n δx n f x x=0 sin x = x x3 3 + x5 5 x7 7 +... x ( ) = a n δ x n ( ) = sin x ak = (-mod(k,2))**(k/2) / fact_k 10 11 I = f x dx a ΔS = f ( x)h I = f a h I = h b (

More information

D0020.PDF

D0020.PDF n 3 X n Y n = Z n 17 1995 300 n n 2 3 2 a b c c 2 a 2 b 2 600 2000 322 3 15 2 3 580 3 1 5 4 3 2 1 300 2 1 2 1 1 ExcelVBA 2 VBA 1 VBA 2 API Sleep ExcelVBA 2 100 60 80 50 ExcelVBA API Sleep 3 100 60 (80

More information

Microsoft Word - 資料 (テイラー級数と数値積分).docx

Microsoft Word - 資料 (テイラー級数と数値積分).docx δx δx n x=0 sin x = x x3 3 + x5 5 x7 7 +... x ak = (-mod(k,2))**(k/2) / fact_k ( ) = a n δ x n f x 0 + δ x a n = f ( n) ( x 0 ) n f ( x) = sin x n=0 58 I = b a ( ) f x dx ΔS = f ( x)h I = f a h h I = h

More information

情報活用資料

情報活用資料 y = Asin 2πt T t t = t i i 1 n+1 i i+1 Δt t t i = Δt i 1 ( ) y i = Asin 2πt i T 21 (x, y) t ( ) x = Asin 2πmt y = Asin( 2πnt + δ ) m, n δ (x, y) m, n 22 L A x y A L x 23 ls -l gnuplot gnuplot> plot "sine.dat"

More information

コンピュータ概論

コンピュータ概論 5.1 VBA VBA Check Point 1. 2. 5.1.1 ( bug : ) (debug) On Error On Error On Error GoTo line < line > 5.1.1 < line > Cells(i, j) i, j 5.1.1 MsgBox Err.Description Err1: GoTo 0 74 Visual Basic VBA VBA Project

More information

120802_MPI.ppt

120802_MPI.ppt CPU CPU CPU CPU CPU SMP Symmetric MultiProcessing CPU CPU CPU CPU CPU CPU CPU CPU CPU CPU CPU CPU CP OpenMP MPI MPI CPU CPU CPU CPU CPU CPU CPU CPU CPU CPU MPI MPI+OpenMP CPU CPU CPU CPU CPU CPU CPU CP

More information

i

i 14 i ii iii iv v vi 14 13 86 13 12 28 14 16 14 15 31 (1) 13 12 28 20 (2) (3) 2 (4) (5) 14 14 50 48 3 11 11 22 14 15 10 14 20 21 20 (1) 14 (2) 14 4 (3) (4) (5) 12 12 (6) 14 15 5 6 7 8 9 10 7

More information

Excel ではじめる数値解析 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.

Excel ではじめる数値解析 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます.   このサンプルページの内容は, 初版 1 刷発行時のものです. Excel ではじめる数値解析 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/009631 このサンプルページの内容は, 初版 1 刷発行時のものです. Excel URL http://www.morikita.co.jp/books/mid/009631 i Microsoft Windows

More information

koji07-02.dvi

koji07-02.dvi 007 I II III 1,, 3, 4, 5, 6, 7 5 4 1 ε-n 1 ε-n ε-n ε-n. {a } =1 a ε N N a a N= a a

More information

リファレンス

リファレンス ii iii iv v vi NEC Corporation 1998 vii C O N T E N T S PART 1 PART 2 viii ix C O N T E N T S PART 3 PART 4 x xi C O N T E N T S PART 5 xii xiii xiv P A R T 1 2 1 3 4 5 1 6 7 1 8 1 9 10 11 1 12 13 1 14

More information

橡統計担当者のためのエクセル表紙.PDF

橡統計担当者のためのエクセル表紙.PDF Num Lock(Pad Lock) ( Num Lock) [ ][ ] [ ][ ] Alt Alt 4 + ( ) 3+3 - ( ) 3-1 -1 * ( ) 3*3 / ( ) 3/3 % ( ) 20% ^ ( ) 3^2 (3*3 ) 2 TRUE FALSE = ( ) A1=B1 > ( ) A1>B1 < ( ) A1= ( ) A1>=B1

More information

Microsoft PowerPoint _MPI-03.pptx

Microsoft PowerPoint _MPI-03.pptx 計算科学演習 Ⅰ ( 第 11 回 ) MPI を いた並列計算 (III) 神戸大学大学院システム情報学研究科横川三津夫 yokokawa@port.kobe-u.ac.jp 2014/07/03 計算科学演習 Ⅰ:MPI を用いた並列計算 (III) 1 2014/07/03 計算科学演習 Ⅰ:MPI を用いた並列計算 (III) 2 今週の講義の概要 1. 前回課題の解説 2. 部分配列とローカルインデックス

More information

OHP.dvi

OHP.dvi 0 7 4 0000 5.. 3. 4. 5. 0 0 00 Gauss PC 0 Gauss 3 Gauss Gauss 3 4 4 4 4 3 4 4 4 4 3 4 4 4 4 3 4 4 4 4 u [] u [3] u [4] u [4] P 0 = P 0 (),3,4 (,), (3,), (4,) 0,,,3,4 3 3 3 3 4 4 4 4 0 3 6 6 0 6 3 6 0 6

More information

XMPによる並列化実装2

XMPによる並列化実装2 2 3 C Fortran Exercise 1 Exercise 2 Serial init.c init.f90 XMP xmp_init.c xmp_init.f90 Serial laplace.c laplace.f90 XMP xmp_laplace.c xmp_laplace.f90 #include int a[10]; program init integer

More information

FORTRAN文法の基礎

FORTRAN文法の基礎 FORTRAN 文法の基礎 ( 初級編 ) 2009-04-16 泉聡志 1 はじめに FORTRAN は数あるプログラム言語の中で最も数値計算に適した言語であり かつ最もかんたんである 加えて FORTRAN を使って数値計算プログラムを作成する工学者は 最小限のことを知っていれば良く 高度な知識は要求されない また 多くのプログラミングは scratch から作らず ベースとなるものを真似て改造して使う場合が多い

More information

広域防災拠点検討会報告書030723表紙_0829_.PDF

広域防災拠点検討会報告書030723表紙_0829_.PDF 15 3 i 15 3 ii iii iv ( ) ( ) ( ) ... i...iii... 1.... 1.... 1..... 1..... 2.... 3... 5.... 5..... 5..... 5.... 6..... 6..... 6.... 7..... 7..... 8... 12... 13.... 13..... 13..... 16..... 16.... 17....

More information

I? 3 1 3 1.1?................................. 3 1.2?............................... 3 1.3!................................... 3 2 4 2.1........................................ 4 2.2.......................................

More information

20 15 14.6 15.3 14.9 15.7 16.0 15.7 13.4 14.5 13.7 14.2 10 10 13 16 19 22 1 70,000 60,000 50,000 40,000 30,000 20,000 10,000 0 2,500 59,862 56,384 2,000 42,662 44,211 40,639 37,323 1,500 33,408 34,472

More information

- 2 -

- 2 - - 2 - - 3 - (1) (2) (3) (1) - 4 - ~ - 5 - (2) - 6 - (1) (1) - 7 - - 8 - (i) (ii) (iii) (ii) (iii) (ii) 10 - 9 - (3) - 10 - (3) - 11 - - 12 - (1) - 13 - - 14 - (2) - 15 - - 16 - (3) - 17 - - 18 - (4) -

More information

2 1980 8 4 4 4 4 4 3 4 2 4 4 2 4 6 0 0 6 4 2 4 1 2 2 1 4 4 4 2 3 3 3 4 3 4 4 4 4 2 5 5 2 4 4 4 0 3 3 0 9 10 10 9 1 1

2 1980 8 4 4 4 4 4 3 4 2 4 4 2 4 6 0 0 6 4 2 4 1 2 2 1 4 4 4 2 3 3 3 4 3 4 4 4 4 2 5 5 2 4 4 4 0 3 3 0 9 10 10 9 1 1 1 1979 6 24 3 4 4 4 4 3 4 4 2 3 4 4 6 0 0 6 2 4 4 4 3 0 0 3 3 3 4 3 2 4 3? 4 3 4 3 4 4 4 4 3 3 4 4 4 4 2 1 1 2 15 4 4 15 0 1 2 1980 8 4 4 4 4 4 3 4 2 4 4 2 4 6 0 0 6 4 2 4 1 2 2 1 4 4 4 2 3 3 3 4 3 4 4

More information

1 (1) (2)

1 (1) (2) 1 2 (1) (2) (3) 3-78 - 1 (1) (2) - 79 - i) ii) iii) (3) (4) (5) (6) - 80 - (7) (8) (9) (10) 2 (1) (2) (3) (4) i) - 81 - ii) (a) (b) 3 (1) (2) - 82 - - 83 - - 84 - - 85 - - 86 - (1) (2) (3) (4) (5) (6)

More information

a n a n ( ) (1) a m a n = a m+n (2) (a m ) n = a mn (3) (ab) n = a n b n (4) a m a n = a m n ( m > n ) m n 4 ( ) 552

a n a n ( ) (1) a m a n = a m+n (2) (a m ) n = a mn (3) (ab) n = a n b n (4) a m a n = a m n ( m > n ) m n 4 ( ) 552 3 3.0 a n a n ( ) () a m a n = a m+n () (a m ) n = a mn (3) (ab) n = a n b n (4) a m a n = a m n ( m > n ) m n 4 ( ) 55 3. (n ) a n n a n a n 3 4 = 8 8 3 ( 3) 4 = 8 3 8 ( ) ( ) 3 = 8 8 ( ) 3 n n 4 n n

More information

=

= 2. 2.1 2.2 kuri@ice.uec.ac.jp ( 2007/10/30/16:46) 1 . 1. 1 + 2 = 5. 2. 180. 3. 3 3. 4.. 5.. 2 2.1 1.,,,,. 2., ( ) ( ).,,,, 3.,. 4.,,,. 3 1.,. 1. 1 + 2 = 5. (, ) 2. 180. (, ) 3. 3, 3. (, ) 4.. (, ) 5..

More information

都道府県別経済財政モデル(平成27年度版)_02

都道府県別経済財政モデル(平成27年度版)_02 -1 (--- 10-2 ---- 4.- 5-3 () 10 13 3 5-4 () 13 16 14-5 () 11 30-1 10 1. 1() Cw j C SNA 47 47 Chi LikL i k1 47 Chi k1 ij Cw j Ch i C SNA L ij j i SNA i j - 2 - -2 5-5 19-3 4 3 4-5 - 3 - 茨 - 4 - -1 (---

More information

インターネット入門

インターネット入門 PART 7 PART 8 ii PART 1 PART 2 PART 3 PART 4 PART 7 PART 4 PART 5 PART 6 iii iv v vi NEC Corporation 1998 vii C O N T E N T S PART 1 PART 2 viii PART 3 ix C O N T E N T S PART 4 PART 5 x PART 6 xi C O

More information

D0050.PDF

D0050.PDF Excel VBA 6 3 3 1 Excel BLOCKGAME.xls Excel 1 OK 2 StepA D B1 B4 C1 C2 StepA StepA Excel Workbook Open StepD BLOCKGAME.xls VBEditor ThisWorkbook 3 1 1 2 2 3 5 UserForm1 4 6 UsorForm2 StepB 3 StepC StepD

More information

untitled

untitled RIKEN AICS Summer School 3 4 MPI 2012 8 8 1 6 MPI MPI 2 allocatable 2 Fox mpi_sendrecv 3 3 FFT mpi_alltoall MPI_PROC_NULL 4 FX10 /home/guest/guest07/school/ 5 1 A (i, j) i+j x i i y = Ax A x y y 1 y i

More information

Excel Excel Excel 20132 20 = 1048576 Excel 201316 100 III 7 (2014 11 18 ) 1

Excel Excel Excel 20132 20 = 1048576 Excel 201316 100 III 7 (2014 11 18 ) 1 III 7 VBA / III 7 (2014 11 18 ) Excel Excel Excel 20132 20 = 1048576 Excel 201316 100 III 7 (2014 11 18 ) 1 Excel VBA Excel Excel 2 20 Excel QR Excel R QR QR BLASLAPACK III 7 (2014 11 18 ) 2 VBA VBA (Visual

More information

OpenMP¤òÍѤ¤¤¿ÊÂÎó·×»»¡Ê£²¡Ë

OpenMP¤òÍѤ¤¤¿ÊÂÎó·×»»¡Ê£²¡Ë 2013 5 30 (schedule) (omp sections) (omp single, omp master) (barrier, critical, atomic) program pi i m p l i c i t none integer, parameter : : SP = kind ( 1. 0 ) integer, parameter : : DP = selected real

More information

My関数の作成演習問題集

My関数の作成演習問題集 Excel Sum,Average,Max 330 BMI Excel My Excel VBA Visual BASIC Editor AltF11 Visual BASIC Editor My Function Function -1- Function ( As Single, As Single) As Double Function Funciton Funciton As Single

More information

2X Y Y X θ 1, θ 2,... Y = f(x, θ 1, θ 2,...) θ k III 8 (2013 05 28 ) 1 / 39

2X Y Y X θ 1, θ 2,... Y = f(x, θ 1, θ 2,...) θ k III 8 (2013 05 28 ) 1 / 39 III 8 (3) VBA, R / III 8 (2013 05 28 ) / 39 2X Y Y X θ 1, θ 2,... Y = f(x, θ 1, θ 2,...) θ k III 8 (2013 05 28 ) 1 / 39 Y X 1, X 2,..., X n Y = f(x 1, X 2,..., X n, θ 1, θ 2,...) (y k, x 1,k, x 2,k,...)

More information

CUDA 連携とライブラリの活用 2

CUDA 連携とライブラリの活用 2 1 09:30-10:00 受付 10:00-12:00 Reedbush-H ログイン GPU 入門 13:30-15:00 OpenACC 入門 15:15-16:45 OpenACC 最適化入門と演習 17:00-18:00 OpenACC の活用 (CUDA 連携とライブラリの活用 ) CUDA 連携とライブラリの活用 2 3 OpenACC 簡単にGPUプログラムが作成できる それなりの性能が得られる

More information

2 2.1 Mac OS CPU Mac OS tar zxf zpares_0.9.6.tar.gz cd zpares_0.9.6 Mac Makefile Mekefile.inc cp Makefile.inc/make.inc.gfortran.seq.macosx make

2 2.1 Mac OS CPU Mac OS tar zxf zpares_0.9.6.tar.gz cd zpares_0.9.6 Mac Makefile Mekefile.inc cp Makefile.inc/make.inc.gfortran.seq.macosx make Sakurai-Sugiura z-pares 26 9 5 1 1 2 2 2.1 Mac OS CPU......................................... 2 2.2 Linux MPI............................................ 2 3 3 4 6 4.1 MUMPS....................................

More information

Gauss

Gauss 15 1 LU LDL T 6 : 1g00p013-5 1 6 1.1....................................... 7 1.2.................................. 8 1.3.................................. 8 2 Gauss 9 2.1.....................................

More information

2/ 土 :30 11:20 似通った科目名がありますので注意してください. 受験許可されていない科目を解答した場合は無効 整理番号と科目コードは受験許可証とよく照合し正確に記入

2/ 土 :30 11:20 似通った科目名がありますので注意してください. 受験許可されていない科目を解答した場合は無効 整理番号と科目コードは受験許可証とよく照合し正確に記入 2/ 土 28 9 10 10:30 11:20 似通った科目名がありますので注意してください. 受験許可されていない科目を解答した場合は無効 整理番号と科目コードは受験許可証とよく照合し正確に記入 30 10 11 12 00011 00016 01101 02607 02703 (1) AB AB 100 cm 2 3.00 cm 2 9.80 m/s 2 AB A B A 10.0 kg A

More information

n 第1章 章立ての部分は、書式(PC入門大見出し)を使います

n 第1章 章立ての部分は、書式(PC入門大見出し)を使います FORTRAN FORTRAN FORTRAN ) DO DO IF IF FORTRAN FORTRAN(FORmula TRANslator)1956 IBM FORTRAN IV FORTRAN77 Fortran90 FORTRAN77 FORTRAN FORTARN IF, DO C UNIX FORTRAN PASCAL COBOL PL/I BASIC Lisp PROLOG Lisp

More information

エクセルカバー入稿用.indd

エクセルカバー入稿用.indd i 1 1 2 3 5 5 6 7 7 8 9 9 10 11 11 11 12 2 13 13 14 15 15 16 17 17 ii CONTENTS 18 18 21 22 22 24 25 26 27 27 28 29 30 31 32 36 37 40 40 42 43 44 44 46 47 48 iii 48 50 51 52 54 55 59 61 62 64 65 66 67 68

More information

3. :, c, ν. 4. Burgers : u t + c u x = ν 2 u x 2, (3), ν. 5. : u t + u u x = ν 2 u x 2, (4), c. 2 u t 2 = c2 2 u x 2, (5) (1) (4), (1 Navier Stokes,.,

3. :, c, ν. 4. Burgers : u t + c u x = ν 2 u x 2, (3), ν. 5. : u t + u u x = ν 2 u x 2, (4), c. 2 u t 2 = c2 2 u x 2, (5) (1) (4), (1 Navier Stokes,., B:,, 2017 12 1, 8, 15, 22 1,.,,,,.,.,,,., 1,. 1. :, ν. 2. : u t = ν 2 u x 2, (1), c. u t + c u x = 0, (2), ( ). 1 3. :, c, ν. 4. Burgers : u t + c u x = ν 2 u x 2, (3), ν. 5. : u t + u u x = ν 2 u x 2,

More information

3. :, c, ν. 4. Burgers : t + c x = ν 2 u x 2, (3), ν. 5. : t + u x = ν 2 u x 2, (4), c. 2 u t 2 = c2 2 u x 2, (5) (1) (4), (1 Navier Stokes,., ν. t +

3. :, c, ν. 4. Burgers : t + c x = ν 2 u x 2, (3), ν. 5. : t + u x = ν 2 u x 2, (4), c. 2 u t 2 = c2 2 u x 2, (5) (1) (4), (1 Navier Stokes,., ν. t + B: 2016 12 2, 9, 16, 2017 1 6 1,.,,,,.,.,,,., 1,. 1. :, ν. 2. : t = ν 2 u x 2, (1), c. t + c x = 0, (2). e-mail: iwayama@kobe-u.ac.jp,. 1 3. :, c, ν. 4. Burgers : t + c x = ν 2 u x 2, (3), ν. 5. : t +

More information

働く女性の母性健康管理、母性保護に関する法律のあらまし

働く女性の母性健康管理、母性保護に関する法律のあらまし 17 1 3 3 12 3 13 10 19 21 22 22 23 26 28 33 33 35 36 38 39 1 I 23 2435 36 4/2 4/3 4/30 12 13 14 15 16 (1) 1 2 3 (2) 1 (1) (2)(1) 13 3060 32 3060 38 10 17 20 12 22 22 500 20 2430m 12 100 11 300m2n 2n

More information

1. A0 A B A0 A : A1,...,A5 B : B1,...,B

1. A0 A B A0 A : A1,...,A5 B : B1,...,B 1. A0 A B A0 A : A1,...,A5 B : B1,...,B12 2. 3. 4. 5. A0 A B f : A B 4 (i) f (ii) f (iii) C 2 g, h: C A f g = f h g = h (iv) C 2 g, h: B C g f = h f g = h 4 (1) (i) (iii) (2) (iii) (i) (3) (ii) (iv) (4)

More information

Pascal Pascal Free Pascal CPad for Pascal Microsoft Windows OS Pascal

Pascal Pascal Free Pascal CPad for Pascal Microsoft Windows OS Pascal Pascal Pascal Pascal Free Pascal CPad for Pascal Microsoft Windows OS 2010 10 1 Pascal 2 1.1.......................... 2 1.2.................. 2 1.3........................ 3 2 4 2.1................................

More information

01_.g.r..

01_.g.r.. I II III IV V VI VII VIII IX X XI I II III IV V I I I II II II I I YS-1 I YS-2 I YS-3 I YS-4 I YS-5 I YS-6 I YS-7 II II YS-1 II YS-2 II YS-3 II YS-4 II YS-5 II YS-6 II YS-7 III III YS-1 III YS-2

More information

橡Taro9-生徒の活動.PDF

橡Taro9-生徒の活動.PDF 3 1 4 1 20 30 2 2 3-1- 1 2-2- -3- 18 1200 1 4-4- -5- 15 5 25 5-6- 1 4 2 1 10 20 2 3-7- 1 2 3 150 431 338-8- 2 3 100 4 5 6 7 1-9- 1291-10 - -11 - 10 1 35 2 3 1866 68 4 1871 1873 5 6-12 - 1 2 3 4 1 4-13

More information

i ii iii iv v vi vii ( ー ー ) ( ) ( ) ( ) ( ) ー ( ) ( ) ー ー ( ) ( ) ( ) ( ) ( ) 13 202 24122783 3622316 (1) (2) (3) (4) 2483 (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 11 11 2483 13

More information

II 2 3.,, A(B + C) = AB + AC, (A + B)C = AC + BC. 4. m m A, m m B,, m m B, AB = BA, A,, I. 5. m m A, m n B, AB = B, A I E, 4 4 I, J, K

II 2 3.,, A(B + C) = AB + AC, (A + B)C = AC + BC. 4. m m A, m m B,, m m B, AB = BA, A,, I. 5. m m A, m n B, AB = B, A I E, 4 4 I, J, K II. () 7 F 7 = { 0,, 2, 3, 4, 5, 6 }., F 7 a, b F 7, a b, F 7,. (a) a, b,,. (b) 7., 4 5 = 20 = 2 7 + 6, 4 5 = 6 F 7., F 7,., 0 a F 7, ab = F 7 b F 7. (2) 7, 6 F 6 = { 0,, 2, 3, 4, 5 },,., F 6., 0 0 a F

More information

B 5 (2) VBA R / B 5 ( ) / 34

B 5 (2) VBA R / B 5 ( ) / 34 B 5 (2) VBAR / B 5 (2014 11 17 ) / 34 VBA VBA (Visual Basic for Applications) Visual Basic VBAVisual Basic Visual BasicC B 5 (2014 11 17 ) 1 / 34 VBA 2 Excel.xlsm 01 Sub test() 02 Dim tmp As Double 03

More information

sim98-8.dvi

sim98-8.dvi 8 12 12.1 12.2 @u @t = @2 u (1) @x 2 u(x; 0) = (x) u(0;t)=u(1;t)=0fort 0 1x, 1t N1x =1 x j = j1x, t n = n1t u(x j ;t n ) Uj n U n+1 j 1t 0 U n j =1t=(1x) 2 = U n j+1 0 2U n j + U n j01 (1x) 2 (2) U n+1

More information

2 X Y Y X θ 1,θ 2,... Y = f (X,θ 1,θ 2,...) θ k III 8 ( ) 1 / 39

2 X Y Y X θ 1,θ 2,... Y = f (X,θ 1,θ 2,...) θ k III 8 ( ) 1 / 39 III 8 (3) VBA, R / III 8 (2013 11 26 ) / 39 2 X Y Y X θ 1,θ 2,... Y = f (X,θ 1,θ 2,...) θ k III 8 (2013 11 26 ) 1 / 39 Y X 1, X 2,..., X n Y = f (X 1, X 2,..., X n,θ 1,θ 2,...) (y k, x k,1, x k,2,...)

More information

2 36 41 41 42 44 44 1 2 16 17 18 19 20 25 26 27 28 29 4 4.12 32 4.2 4.2.1 36 4.2.2 41 4.2.3 41 4.2.4 42 4.3 4.3.1 44 4.3.2 44 31 1 32 33 < 2 x 1 x x 2 < x 1 x1x 2 x1x 2 34 36 4.2 (1) (4) (1)

More information

1

1 005 11 http://www.hyuki.com/girl/ http://www.hyuki.com/story/tetora.html http://www.hyuki.com/ Hiroshi Yuki c 005, All rights reserved. 1 1 3 (a + b)(a b) = a b (x + y)(x y) = x y a b x y a b x y 4 5 6

More information

3 5 18 3 5000 1 2 7 8 120 1 9 1954 29 18 12 30 700 4km 1.5 100 50 6 13 5 99 93 34 17 2 2002 04 14 16 6000 12 57 60 1986 55 3 3 3 500 350 4 5 250 18 19 1590 1591 250 100 500 20 800 20 55 3 3 3 18 19 1590

More information

Sae x Sae x 1: 1. {x (i) 0 0 }N i=1 (x (i) 0 0 p(x 0) ) 2. = 1,, T a d (a) i (i = 1,, N) I, II I. v (i) II. x (i) 1 = f (x (i) 1 1, v(i) (b) i (i = 1,

Sae x Sae x 1: 1. {x (i) 0 0 }N i=1 (x (i) 0 0 p(x 0) ) 2. = 1,, T a d (a) i (i = 1,, N) I, II I. v (i) II. x (i) 1 = f (x (i) 1 1, v(i) (b) i (i = 1, ( ) 1 : ( ) Sampling/Imporance resampling (SIR) Kiagawa (1993, 1996), Gordon(1993) EnKF EnKF EnKF 1CPU 1core 2 x = f (x 1, v ) y = h (x, w ) (1a) (1b) PF p(x y 1 ) {x (i) 1 }N i=1, p(x y ) {x (i) }N i=1

More information

困ったときのQ&A

困ったときのQ&A ii iii iv NEC Corporation 1997 v P A R T 1 vi vii P A R T 2 viii P A R T 3 ix x xi 1P A R T 2 1 3 4 1 5 6 1 7 8 1 9 1 2 3 4 10 1 11 12 1 13 14 1 1 2 15 16 1 2 1 1 2 3 4 5 17 18 1 2 3 1 19 20 1 21 22 1

More information

178 5 I 1 ( ) ( ) 10 3 13 3 1 8891 8 3023 6317 ( 10 1914 7152 ) 16 5 1 ( ) 6 13 3 13 3 8575 3896 8 1715 779 6 (1) 2 7 4 ( 2 ) 13 11 26 12 21 14 11 21

178 5 I 1 ( ) ( ) 10 3 13 3 1 8891 8 3023 6317 ( 10 1914 7152 ) 16 5 1 ( ) 6 13 3 13 3 8575 3896 8 1715 779 6 (1) 2 7 4 ( 2 ) 13 11 26 12 21 14 11 21 I 178 II 180 III ( ) 181 IV 183 V 185 VI 186 178 5 I 1 ( ) ( ) 10 3 13 3 1 8891 8 3023 6317 ( 10 1914 7152 ) 16 5 1 ( ) 6 13 3 13 3 8575 3896 8 1715 779 6 (1) 2 7 4 ( 2 ) 13 11 26 12 21 14 11 21 4 10 (

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション ループ ループとは? ある条件を満たすまで 指定の命令を繰り返す Do... Loop For Next For Each Next While WEnd ループの種類 Do Loop Do While 条件 ステートメント Loop Do ステートメント Loop While 条件 Do Until 条件 ステートメント Loop Do ステートメント Until Loop 条件 Do Loop

More information

02_C-C++_osx.indd

02_C-C++_osx.indd C/C++ OpenMP* / 2 C/C++ OpenMP* OpenMP* 9.0 1... 2 2... 3 3OpenMP*... 5 3.1... 5 3.2 OpenMP*... 6 3.3 OpenMP*... 8 4OpenMP*... 9 4.1... 9 4.2 OpenMP*... 9 4.3 OpenMP*... 10 4.4... 10 5OpenMP*... 11 5.1

More information

リファレンス

リファレンス ii iii iv v vi NEC Corporation 1998 vii C O N T E N T S PART 1 viii ix C O N T E N T S x PART 2 xi C O N T E N T S PART 3 PART 4 xii PART 5 xiii C O N T E N T S xiv PART 6 xv xvi 2 3 4 5 6 7 P A R T

More information

<4D6963726F736F667420576F7264202D208DEC90AC837D836A83858341838B81698F4390B394C5816A2E646F63>

<4D6963726F736F667420576F7264202D208DEC90AC837D836A83858341838B81698F4390B394C5816A2E646F63> + = 付 録 Ⅰ マクロのフローチャート 開 始 キーワード 数 と 文 章 数 をカウントする Yes 全 ての 文 章 番 号 (1 0) を 照 合 したか No 全 ての 文 章 番 号 (data) を 照 合 したか Yes No 文 章 番 号 (1 0) が 文 章 番 号 (data) と 等 しいか Yes No 値 を 一 行 全 て 0にする 値 を 一

More information

MPI usage

MPI usage MPI (Version 0.99 2006 11 8 ) 1 1 MPI ( Message Passing Interface ) 1 1.1 MPI................................. 1 1.2............................... 2 1.2.1 MPI GATHER.......................... 2 1.2.2

More information

演習2

演習2 神戸市立工業高等専門学校電気工学科 / 電子工学科専門科目 数値解析 2017.6.2 演習 2 山浦剛 (tyamaura@riken.jp) 講義資料ページ h t t p://clim ate.aic s. riken. jp/m embers/yamaura/num erical_analysis. html 曲線の推定 N 次多項式ラグランジュ補間 y = p N x = σ N x x

More information

1. 入力画面

1. 入力画面 指定した時刻に指定したマクロ (VBA) を実行するプログラム (VBA) 益永八尋 様々な業務を行っている場合には 指定した時刻に指定したマクロ (Macro VBA) を実行したくなる場合がある たとえば 9:00 17: 00 や 1 時間 6 時間間隔に指定したマクロ (Macro VBA) を実行する この様な場合に対応できるように汎用性の高いプログラムを作成した この場合に注意する必要があるのは

More information