$B>pJs2J3X9V5f(B - MPI

$B0p3@(B $BNI0l(B, $BNkLZ(B $BE/=((B, $BCfB<(B $B9%0l(B

• $B3FC
• $B9TNs1i;;$N9bB.2=(B
• $B%W%m%0%i%`(B

$B3FC

$B"%La$k(B

$B9TNs7W;;$N9bB.2=(B

C=A*B$B$N9TNs7W;;$r(BMPI$B$rMQ$$$FJBNs2=!#4D6-$d(BPE$B?t$G7W;;B.EY$r7WB,!&Hf3S$7$^$7$?!#(B

$B

• $B9TNs7W;;(B
  • $B9TNs$N%V%m%C%/2=(B
  • $B9TNs$O0l
• $B$=$NJBNs2=(B
  • $B;HMQ(BPE$B?t(B(n$B$H$9$k(B)$B$K$"$o$;$F9TNs(BA$B$r(Ba₁, a₂,...a_n$B$H9TJ}8~$KJ,3d!#(B
  • a_i$B$H(BB$B$r3F(BPE$B$KG[$j!"(Bc_i=a_i*B$B$r7W;;!#(B
  • $B3F(BPE$B$O7W;;7k2L(Bc₁, c₂,...c_n$B$r?F(BPE$B$KJV$7$F=8Ls$9$k!#(B

$B7WB,J}K!(B

1$B$D$N>r7o$K$D$$$F!"(B20$B2s!JBh(B4$BC$BIC(B$B$G$9!#(B

$B;HMQ$7$?%3%s%Q%$%k%*%W%7%g%s$O(B-O3 -funroll-loops$B$NFs$D!#(B

$B7k2L(B

edam

PE$B?t(B	500*500	1000*1000	1500*1500
1	7.52146180	57.830102	195.8328042
2	3.98120615	29.303396	99.6881416
3	2.70881360	21.979148	66.8359671
4	2.03025785	15.926105	56.6987160
5	1.81121250	13.191404	41.5771440
6	1.48576565	12.009133	37.1108106
7	1.34372030	10.846369	32.1309268
8	1.25148895	9.478561	28.5594552
9	1.19568575	8.233586	26.1275945
10	1.17231610	7.728422	24.5318186

$B#1%W%m%;%9;~$N.$5$/$J$C$F$*$j!"(B $BJ,C4=hM}$N8z2L$,$"$i$o$l$F$$$k!#(B

PE$B?t(B	500*500	1000*1000	1500*1500
1	1	1	1
2	0.529312819	0.506715274	0.509047205
3	0.360144567	0.38006414	0.341290967
4	0.269928626	0.275394724	0.289526141
5	0.240805916	0.228106186	0.212309394
6	0.197536821	0.207662317	0.189502524
7	0.178651482	0.187555765	0.164073261
8	0.166389059	0.163903584	0.14583591
9	0.158969863	0.142375436	0.133417864
10	0.1558628	0.133640124	0.1252692

Folon$B-7(B

$B"((BMyrinet2k$B$r;H$C$F$$$^$9!#6u=j$O%G!<%?L$:N=8(B

PE$B?t(B	500	1000	1500	2000	3000	4000
1	1.13526235	9.02041140	30.26844630	71.59615130	241.1294860	566.5069245
2	0.59692580	4.63417445	15.41836675	36.33357125	121.7571952	286.4863531
3	0.42209055	3.18271645	11.05881840	24.63382525	82.1141123	192.2656200
4	0.34422660	2.52148800	8.14342180	18.98036210	62.8521370	146.2858617
5	0.30704455
6	0.27810975
7	0.26648115
8	0.26197680	1.61582020	5.07055825	10.98437880	34.6697289	79.4249104

Folon$B-7$G$N.$5$$$H(Bedam$B$h$j8zN($,0-$/$J$C$F$7$^$C$F$$$k!#(B $B5U$K(Bedam$B$O%a%b%j$r6&M-$7$F$$$k$+$i8zN($,$h$$$N$+$b$7$l$J$$!#(B $B9TNs%5%$%:$,Bg$-$/$J$k$H!"$[$\(Bedam$B$HF1$8$/$i$$$N%Q%U%)!<%^%s%9$rH/4x$7$F$$$k$H8+$k$3$H$,$G$-$k!#(B

PE$B?t(B	500	1000	1500	2000	3000	4000
1	1	1	1	1	1	1
2	0.525804278	0.513743137	0.509387452	0.507479391	0.504945277	0.50570671
3	0.371800007	0.352834955	0.36535798	0.344066333	0.340539491	0.339387943
4	0.303213262	0.279531375	0.269039967	0.265103106	0.260657201	0.258224313
8	0.230763224	0.179129324	0.167519608	0.153421359	0.143780545	0.140201129

$BBh#4C

$B"(;29M$G$9!#Bf?t$,A}$($k$H%Q%U%)!<%^%s%9Mn$A$k$7!&!&!&(B

PE$B?t(B	500	1000
1	4.0723943	29.613245
2	5.3833844	21.223144
3	12.39142405	44.710156
4	14.55225079	44.831185
5	16.79908344	43.363055
6	22.24323065
7	34.10413025
8	41.6756541
9	46.93410175
10	61.63975955

$B"%La$k(B

$B%W%m%0%i%`(B

$B;HMQ$7$?%W%m%0%i%`$r(B$BCV$$$F$*$-$^$9(B$B!#:#$N$H$3$m9TNs1i;;$N%W%m%0%i%`0l<0$N$_!#(B $B%=!<%9$O1x$$$N$G!"2K$,$"$l$P@6=q$9$k$+$b$7$l$^$;$s!#(B

matrix.c$B$H(Bmatrix.h$B$O9TNs1i;;$N$?$a$N%W%m%0%i%`$G$9!#HFMQ@-$OHyL/$K$"$k$+$b!#(B $B2.ED$5$s$N(B$B9TNs>h;;9bB.2=%W%m%8%'%/%H(B$B$r;29M$K$7$^$7$?!#(B

$B"%La$k(B

<< $B%H%C%W$KLa$k(B

---

Copyright© 2002,$B0p3@(B $BNI0l(B<inagaki@ueda.info.waseda.ac.jp>
$B:G=*99?7F|(B 2002/2/11