もっと高速な方法はあるでしょうか? 島村1243 2007/12/13 14:21:55 (修正1回) └続き 島村1243 2007/12/13 14:23:05 (修正1回) ├(1)について(処理手順の改良) 山中和義 2007/12/14 08:55:07
Re: 続き |
返事を書く ノートメニュー |
山中和義 <drdlxujciw> 2007/12/14 08:55:07 | |
(1)について(処理手順の改良)
等高線の数だけ複素平面を走査して毎回wを算出すると時間を費やします。 FOR文のネストをかえると多少速くなります。 「1回の走査で該当する等高線を描く」ように改修する場合(FOR文の部分抜粋) SET POINT STYLE 1 FOR k=0 TO c_Div*b-1 !Y座標 ※ドットをもとに LET y=Ya+dy*k ![Ya,Yb] FOR i=0 TO c_Div*a-1 !X座標 LET x=Xa+dx*i ![Xa,Xb] LET z=x+y*j !空間上の点を計算する IF z<>a1 AND z<>a2 THEN LET w=q1/ABS(z-a1)+q2/ABS(z-a2) !点zにおける合成電位 FOR wb=0 TO 0.5 STEP 0.02 !等電位幅値の設定 IF wb<=0 THEN SET POINT COLOR "black" ELSE SET POINT COLOR "red" END IF IF ABS(w-wb) < cEps THEN PLOT POINTS: x,y !点画 NEXT wb END IF NEXT i NEXT k |
│├!もったいない、整理すれば、4倍速。 SECOND 2007/12/14 21:34:04 ││└確かに速いです。 島村1243 2007/12/14 22:33:52 ││ └どっちでしょう? SECOND 2007/12/14 22:50:11 (修正1回) │└20秒速くなりました。 島村1243 2007/12/14 22:06:46 └(3)について各点での流線 山中和義 2007/12/15 09:43:05 ├!等ポテンシャルを色の濃淡で表現する 山中和義 2007/12/15 11:03:55 └これを利用すると流れの方向がわかるので、... 山中和義 2007/12/15 20:11:13 └とても数理的で良いですね。 島村1243 2007/12/15 23:15:50 (修正1回) └接地導体球の作画に問題が有る様 島村1243 2007/12/16 07:37:08 └対策案できました。 島村1243 2007/12/16 10:44:04 (修正1回) └ポテンシャル関数f(x,y)なら、電位分布は等... 山中和義 2007/12/16 13:56:21 └完璧に目的達成です! 島村1243 2007/12/16 19:04:03 (修正1回) └プログラムの続き 島村1243 2007/12/16 19:05:45