SAD部屋 : Recent Changes
http://jyurai.ddo.jp/~amorita/SAD/?c=recent
SAD部屋 Recent Changes
ja
Copyright (C) Akio Morita
2024-11-18T13:52:14+00:00
-
SADベンチ(bench2.sad)
http://jyurai.ddo.jp/~amorita/SAD/?Benchmarks%2FSAD%2FBench2
2024-11-18T13:52:14+00:00
@@ -10,9 +10,9 @@
||Ryzen TR 2950X||1.166768||1.830007||0.547440||0.121408||0.053810||
!! amorita forked branch r7688
-
-||!CPU ||!Compiler ||!Function||!Optics||!Tracking||!Matching||!Overall||
-||Ryzen 7950X||GNU Fortran 12.2||0.608823||0.964649||0.444471||0.063279||0.029833||
-||Ryzen 5950X||GNU Fortran 12.2||0.943251||1.299135||0.588711||0.084467||0.043003||
-||Ryzen 3700X||GNU Fortran 12.2||1.073548||1.529004||0.716498||0.101149||0.050005||
-||Ryzen TR 2950X||GNU Fortran 12.2||1.204148||1.826535||0.541115||0.121333||0.054597||
+* GNU Fortran 13.2 + FreeBSD clang
+||!CPU ||!Function||!Optics||!Tracking||!Matching||!Overall||
+||Ryzen 7950X||0.608823||0.964649||0.444471||0.063279||0.029833||
+||Ryzen 5950X||0.943251||1.299135||0.588711||0.084467||0.043003||
+||Ryzen 3700X||1.073548||1.529004||0.716498||0.101149||0.050005||
+||Ryzen TR 2950X||1.204148||1.826535||0.541115||0.121333||0.054597||]]> -
SideMenu
http://jyurai.ddo.jp/~amorita/SAD/?SideMenu
2023-04-10T02:52:27+00:00
@@ -26,7 +26,7 @@
!!SAD
*[[本家|http://acc-physics.kek.jp/SAD/]]
*[[公式Wiki|SADWiki:FrontPage]]
-*[[掲示板|https://acc-physics.kek.jp/SAD/wforum/sad-sad/]]
+*[[掲示板|https://acc-physics.kek.jp/SAD/wforum/sad-sad/]]([[2020～|https://acc-physics.kek.jp/cgi-bin/SAD2/wforum.cgi]]/[[2014～|https://acc-physics.kek.jp/cgi-bin/SAD/SAD-BB/wforum.cgi?page=0&list=tree]]/[[2003～|https://acc-physics.kek.jp/SAD_legacy/ArchivedBBS/bbs.html]])
!!サイト内
*[[Orz日記|http://jyurai.ddo.jp/~amorita/diary/]]]]> -
開発メモ
http://jyurai.ddo.jp/~amorita/SAD/?%E9%96%8B%E7%99%BA%E3%83%A1%E3%83%A2
2023-03-24T14:55:01+00:00
@@ -20,9 +20,10 @@
* Besselの実装改良
** BesselYの精度不良
** 張り合わせ付近の精度不足
+*** $1/z \sim 0$近傍の漸近展開は、比較的安定に精度が出ている模様
** BesselKなので、一部の NaNが出る不良領域あり
** 定義域によっては、積分公式を二重指数変換型の数値積分を行うのも現実的
-*** BesselIの右半面などは、比較的良好な結果が得られる模様
+*** BesselKの右半面などは、比較的良好な結果が得られる模様
* GUI系のカスタマイズホックの増強
** KBFCursorEntryへの '''TextWidgetOption'''追加
** Canvasへのグラフ描画時のFont・Ticksカスタマイズ]]> -
FrontPage
http://jyurai.ddo.jp/~amorita/SAD/?FrontPage
2023-03-24T12:37:37+00:00
@@ -1,4 +1,4 @@
!このサイトは何？
-[[SAD|http://acc-physics.kek.jp/SAD/]]の個人的な開発forkに関わる
+[[SAD]]の個人的な開発forkに関わる
雑多な文書を整理するために立ち上げた[[Hiki|http://hikiwiki.org/ja/]]サイトです]]> -
InterWikiName
http://jyurai.ddo.jp/~amorita/SAD/?InterWikiName
2023-03-24T12:37:17+00:00
@@ -1,9 +1,10 @@
*[[GoogleJ|http://www.google.com/search?num=50&lr=lang_ja&q=]] sjis
*[[Hiki|http://hikiwiki.org/ja/?]] euc
*[[Image|http://images.google.co.jp/images?hl=ja;btnG=Google?hl=ja;btnG=Google+%8C%9F%8D%F5;q=]] sjis
*[[WikiWikiWeb|http://www.c2.com/cgi/wiki?]]
*[[Ruby|http://www.ruby-lang.org/]] alias
*[[SADWiki|http://www-kekb.kek.jp/Documentation/SAD/index.php?]]
+*[[SAD|https://acc-physics.kek.jp/SAD/]] alias
*[[self|http://jyurai.ddo.jp/~amorita/]] alias
// Alias for TextFormattingRules]]> -
姫野ベンチ
http://jyurai.ddo.jp/~amorita/SAD/?Benchmarks%2FHimeno
2023-02-03T16:02:20+00:00
@@ -23,13 +23,6 @@
||!Processor ||!Compiler ||!XS ||!S ||!M ||!L ||!XL ||!Options ||
'''2023/mm/dd測定'''
-||Processor ||Ways|| MHz||Architecture ||Compiler ||XS ||S ||M ||L ||XL ||Options
-||Xeon X5550 || 2 ||2667||AMD64-FreeBSD8 ||GNU Fortran 4.5.1 20100506||3874.3845||7833.4702||7518.5513||6616.7612||5510.0771||-O3
-||Xeon X5482 || 2 ||3200||i386-Darwin9 ||GNU Fortran 4.5.0 20100114||3460.2549||3869.0618||3075.3809||2776.5073||no memory||-O3
-||Xeon E5462 || 2 ||2800||i386-Darwin9 ||GNU Fortran 4.5.0 20100114||2821.7212||2857.2898||2437.5750||2208.3286||no memory||-O3
-||Xeon 5160 || 2 ||3000||i386-Darwin9 ||GNU Fortran 4.5.0 20100114||2164.6001||1748.8965||1538.7083||swaping ||no memory||-O3
-||PowerPC G5 || 2 ||2500||PowerPC-Darwin8||GNU Fortran 4.5.0 20100114||616.30756||621.87097||632.66132||466.56357||no memory||-O3
-
!!雑感
* C dynamic allocがC staticに対して性能が出てないのは、添字計算マクロの効率もしくは配列次元が定数畳み込み出来ていないためか？
** 原理的には、配列次元の定数畳み込みを除けば同じところまでは行けるはず…]]> -
SAD乱数生成ベンチ
http://jyurai.ddo.jp/~amorita/SAD/?Benchmarks%2FSAD%2FPRNG
2023-02-03T05:27:14+00:00
@@ -3,4 +3,4 @@
||!CPU ||!SAD||!MT19936/53bit||!MT19936/32bit||!SFMT19936/53bit||!SFMT19936/32bit
||Ryzen 7950X|| 0.378|| 5.966|| 3.178|| 1.354|| 0.669||
||Ryzen 5950X|| 0.452|| 6.156|| 3.312|| 1.641|| 0.825||
-||Ryzen TR 2950X||
+||Ryzen TR 2950X|| 0.668|| 7.073|| 3.792|| 2.066|| 1.177||]]> -
開発メモ/64bit環境
http://jyurai.ddo.jp/~amorita/SAD/?%E9%96%8B%E7%99%BA%E3%83%A1%E3%83%A2%2F64bit%E7%92%B0%E5%A2%83
2023-02-01T15:52:03+00:00
@@ -7,8 +7,8 @@
* rlist(*)をFortran pointerとして定義し、最初に割り当てたヒープブロックを基準に先頭アドレスを調整する
* mmap(2)の拡張機能を使い'''16GiB'''の仮想アドレス空間を予約し、rlist(*)の予約領域内に割り当てる
-!64bit環境のサポート法
-!! 1. データサイズを一回り大きくする(コードには手を付けず解釈だけ変えて誤魔化す)
+!64bit化
+!! データサイズを一回り大きくする(コードには手を付けず解釈だけ変えて誤魔化す)
* '''REAL'''を4倍精度実数(''128bit'')
* '''INTEGER'''や'''INTEGER*4'''を''64bit''整数
* '''BYTE'''、'''CHARACTER'''や'''INTEGER*1'''を''16bit''整数
@@ -17,12 +17,14 @@
** 最近の'''M17N'''のトレンドを考えると、8倍精度実数/128bit整数/UCS-4な文字表現の方が良いかも(IPv6アドレスが整数型に直接収容できるし)
Fortran I/O(READ/WRITE文)を対応させる必要があるので、コンパイラ側の対応、若しくは改造が必須
-!! 2. 真面目に、内部アドレス表現を 64bit化する(SADの Fortranメインコードを書き直す)
-* 内部データ構造が、コード中で隠蔽されていないのでとても大変
-* データ構造のまともな資料が現存しない
-* 改修後のコードの動作を検証するテストセットが存在しないので、検証が大変
+!! 真面目に、内部アドレス表現を 64bit化する(SADの Fortranメインコードを書き直す)
+# アドレス表現を含む内部APIをポインタベースのAPIに置き換える
+# アドレス表現を含む内部データ構造へのアクセスに抽象化層を挟み、すべて抽象化層経由に書き換える
+# 抽象化層で隠蔽されたアドレス表現を含む部データ構造を再設計・再実装する
+* データ構造のまともな資料が現存しない・実装コードからの操作が一貫していない(旧い・バグったコードが残っている)
+* 改修後のコードの動作を検証するテストセットが存在しない(変更後の動作検証が大変)
-!! 3. 古くさい Fortranコードを捨てて、もっと保守性の良いコードを書き起こす
+!! 古くさい Fortranコードを捨てて、もっと保守性の良いコードをfrom scratchで書き起こす
* 車輪の再発明
** J-PARCでSADScriptインタプリタ部の Javaによる再実装をやっているらしいので、注意しないと'''車輪の再再発明'''になる''罠''
* 仕様書やまともなテストセットが現存しないので、互換性保持が困難]]> -
開発メモ/Process
http://jyurai.ddo.jp/~amorita/SAD/?%E9%96%8B%E7%99%BA%E3%83%A1%E3%83%A2%2FProcess
2023-02-01T14:06:25+00:00
@@ -0,0 +1,14 @@
+!Process制御
+* プロセス間通信関連機能
+** Object Serializerを用意する(Network Byte Orderサポートも?)
+** Unibyte stream Encoder/Decoderを用意する(BASE64辺り？)
+* MPIサポート
+** Rankの問い合わせ
+** ノード間のオブジェクト転送 (Object Serializerが必要)
+** ノード間同期
+* User Level Multi Thread
+** Software Threadの場合
+*** インタプリタ状態の保存・復元・Context Dispatch・Thread Scheduler
+** Native Threadの場合
+*** インタプリタ本体および関連HelperコードのReentrant化/Multi Thread Safe化
+*** SigActionのシグナルハンドラ上からのスクリプト実行に Signal Safe化が必要]]> -
開発メモ/IO
http://jyurai.ddo.jp/~amorita/SAD/?%E9%96%8B%E7%99%BA%E3%83%A1%E3%83%A2%2FIO
2023-02-01T14:01:44+00:00
@@ -0,0 +1,7 @@
+!Fortran I/Oを捨てる
+* READ statementの置き換えは、Oide氏が実装ずみ
+** read(2)による itfgetbuf_実装は非常に遅い
+** stdioを使うには、lun -> FILE型のマップを自前で維持する必要が...
+** 自前で、read bufferを管理する?(FSEEKと干渉する/FSEEKをこちらで実装すると、WRITE側のFSEEKで問題が)
+* WRITE(lfno,)系が残っている
+* FSEEK(lfno,,)が残っている]]>