Archive

Archive for the ‘Cell’ Category

薄型PS3が出るぞ! だがしかしっ!

8月 19th, 2009

薄くなった新 PS3 正式発表、2万9980円 (動画、ギャラリー追加)


ちょっと前からあちこちで話題になってはいたけど,どうやら公式に発表されたらしい.薄い!比較表は本家engadgetがわかりやすい.


PS3 Slim sized up: smaller, deeper, no Linux or PS2 compatibility


……あ……れ……PS2互換はともかく,Linux互換が削れてるぞ……?

待て待て,Linuxを入れればCellがいじれる,ってのはPS3の購入動機の一つとして市民権を得ているモノだと勝手に思っていたんだが,そんなことは無かったかチクショウ.

安価に手に入る,気軽に使えるってのは重要なはずなんだけどなあ?


薄型PS3が出たら(弄る時間があるかどうかはともかくとして)Cellを弄って遊ぶために私物として一台購入しようと思っていた私が涙目です.ファッキン.

tgbt *nix・サーバ関連, Cell, ゲーム

このエントリをはてなブックマークに追加このエントリをdel.icio.usに追加このエントリをLivedoor Clipに追加このエントリをYahoo!ブックマークに追加このエントリをFC2ブックマークに追加このエントリをNifty Clipに追加このエントリをPOOKMARK. Airlinesに追加このエントリをBuzzurl(バザール)に追加このエントリをChoixに追加このエントリをnewsingに追加


CellChallenge予選終了

2月 27th, 2009

もう少し時間があるけど,やめた.


昨夜,狙ったかのように講座のPS3Linuxが不調に.SPEが応答しない.

仕方がないので主催者側で用意しているCellマシンを利用したが,混雑していて自由にプログラムを実行できない.ところがCell上でのプログラム実行時間はばらつきがそれなりにあるので,繰り返し実行しないと最適化の優位性が確認できない.

というわけで,アンローリングとか細かい最適化をいろいろやってみたかったけど,成果がよくわからない環境状況(?)なので諦めた.

まぁ正しく動くものは出したつもりなので,決勝に残れたらもう少しやるか.


昨夜DMAダブルバッファリングにはまりまくったんだけど,原因は自分で作成した関数の引数順序を勘違いしていたのが原因だった.ばーかばーか.まぁ正しく動いても全然速くなってないんだけどね! orz

tgbt Cell

このエントリをはてなブックマークに追加このエントリをdel.icio.usに追加このエントリをLivedoor Clipに追加このエントリをYahoo!ブックマークに追加このエントリをFC2ブックマークに追加このエントリをNifty Clipに追加このエントリをPOOKMARK. Airlinesに追加このエントリをBuzzurl(バザール)に追加このエントリをChoixに追加このエントリをnewsingに追加


CellChallenge提出日

2月 27th, 2009

17時だったかな.

26日の夜に研究室のPS3がいかれたのはお約束。


とりあえずSIMD化はした.PS3上で倍くらいの性能になった.

あとはDMAダブルバッファリングでもやってみようと思って,調整しやすいように静的スケジューリングを作成した.でも肝心のDMAが,非同期にした時点でアウト.どの処理が正しくない結果を出しちゃうのかはわかるのに,どうして正しい結果が得られないのかがわからないという素晴らしい状態.もうどうにでもなあーれっ.


とある配列のDMA転送を同期か非同期か変更すると駄目になるんだけど,非同期化させる方法が間違ってるんだろうなぁ.問題が起きる場所がピンポイント過ぎて意味がわからない.


とりあえず提出状態にはしたので,あと朝起きたら夕方までループアンローリングを行う仕事でもします.

予選をまともに通過したらDMAダブルバッファリングに再挑戦するかもね.

tgbt Cell

このエントリをはてなブックマークに追加このエントリをdel.icio.usに追加このエントリをLivedoor Clipに追加このエントリをYahoo!ブックマークに追加このエントリをFC2ブックマークに追加このエントリをNifty Clipに追加このエントリをPOOKMARK. Airlinesに追加このエントリをBuzzurl(バザール)に追加このエントリをChoixに追加このエントリをnewsingに追加


ほぼ初めてみたいなSIMD(ry 2

2月 10th, 2009

ほぼ初めてみたいなSIMDプログラミング,略してはじしむ?


とりあえずCellのSPE向けにSIMDってみた.

サンプルとしてメディアンフィルタ……みたいなものを意識して,ある点の周囲の点の座標を合算するものを作ってみた.


まずはCで書いた元プログラム.

C:

  1. for(m=0; m<LOOP; m++){
  2.     margin = m*128*128;
  3.     for(i=1; i<128-1; i++){
  4.       for(k=1; k<128-1; k++){
  5.         data2[margin + i*128+k] =
  6.           (data1[margin + (i-1)*128+(k-1)]
  7.            + data1[margin + (i-1)*128+(k+1)]
  8.            + data1[margin + (i+1)*128+(k-1)]
  9.            + data1[margin + (i+1)*128+(k+1)]);
  10.       }
  11.     }
  12.   }


マージンが取ってあるのは,繰り返し計算して実行時間を引き延ばすため.めんどいから端点の処理はしない.ただ足しているだけなのは処理を最低限にしないと正しく実行できているかのチェックがツライため.チキンコード.

ちなみにデータ型はunsigned int.


続いてSPE向けのコード.

C:

  1. for(i=1; i<128-1; i++){
  2.       for(k=1; k<128-1; k++){
  3.         outdata[i*128+k] =
  4.           (
  5.            indata[(i-1)*128+(k-1)] +
  6.            indata[(i-1)*128+(k+1)] +
  7.            indata[(i+1)*128+(k-1)] +
  8.            indata[(i+1)*128+(k+1)]
  9.            );
  10.       }
  11.     }


データを128*128ごとにDMAってるのでこんなんになる.


続いてSIMD化コード.

C:

  1. for(i=1; i<128-1; i++){
  2.       int x, y;
  3.       vector unsigned int *v[6];
  4.       vector unsigned int *t;
  5.       vector unsigned int w[6];
  6.       vector unsigned char s0 = (vector unsigned char){
  7.         0x0c, 0x0d, 0x0e, 0x0f,
  8.         0x10, 0x11, 0x12, 0x13,
  9.         0x14, 0x15, 0x16, 0x17,
  10.         0x18, 0x19, 0x1a, 0x1b
  11.       };
  12.       vector unsigned char s1 = (vector unsigned char){
  13.         0x04, 0x05, 0x06, 0x07,
  14.         0x08, 0x09, 0x0a, 0x0b,
  15.         0x0c, 0x0d, 0x0e, 0x0f,
  16.         0x10, 0x11, 0x12, 0x13
  17.       };
  18.       for(x=2; x<30; x++){
  19.         v[0] = (vector unsigned int *)&indata[(i-1)*128+(x-1)*4];
  20.         v[1] = (vector unsigned int *)&indata[(i-1)*128+()*4];
  21.         v[2] = (vector unsigned int *)&indata[(i-1)*128+(x+1)*4];
  22.         v[3] = (vector unsigned int *)&indata[(i+1)*128+(x-1)*4];
  23.         v[4] = (vector unsigned int *)&indata[(i+1)*128+()*4];
  24.         v[5] = (vector unsigned int *)&indata[(i+1)*128+(x+1)*4];
  25.         t = (vector unsigned int *)&outdata[i*128+x*4];
  26.  
  27.         w[0] = spu_shuffle(*v[0], *v[1], s0);
  28.         w[1] = spu_shuffle(*v[1], *v[2], s1);
  29.         w[2] = spu_shuffle(*v[3], *v[4], s0);
  30.         w[3] = spu_shuffle(*v[4], *v[5], s1);
  31.         w[4] = spu_add(w[0], w[1]);
  32.         w[5] = spu_add(w[2], w[3]);
  33.         *t = spu_add(w[4], w[5]);
  34.       }
  35.       for(k=1; k<8; k++){
  36.         outdata[i*128+k] =
  37.           (
  38.            indata[(i-1)*128+(k-1)] +
  39.            indata[(i-1)*128+(k+1)] +
  40.            indata[(i+1)*128+(k-1)] +
  41.            indata[(i+1)*128+(k+1)]
  42.            );
  43.       }
  44.       for(k=120; k<128-1; k++){
  45.         outdata[i*128+k] =
  46.           (
  47.            indata[(i-1)*128+(k-1)] +
  48.            indata[(i-1)*128+(k+1)] +
  49.            indata[(i+1)*128+(k-1)] +
  50.            indata[(i+1)*128+(k+1)]
  51.            );
  52.       }
  53.  
  54.     }


やりやすいところだけSIMD化した.kループの真ん中あたりだけ.もっとはじっこまで行ける気がするのはとりあえず置いておく.


計算対象の配列のアドレスをvector型のポインタで参照しておいて,vector型で計算すると,計算対象にちゃんと反映されるよ!みたいな感じなのかね.まぁポインタだから当然だけど.なんとなくわかった気がする.

C:

  1. 正:t = (vector unsigned int *)&outdata[i*128+x*4];
  2.         誤:t = (vector unsigned int *)outdata[i*128+x*4];


みたいなミスを繰り返しまくったんだけど,vector型へのキャストが強力すぎるのかなんなのか,コンパイラが何も言わないのでしばらく頭を抱えた.これはツライ.


で,とりあえずSPE1つで実行時間を測定してみたんだけど……5回繰り返した最短時間を見ると

  • SIMD化する前: 13.36145091 sec
  • SIMD化した後: 11.95487118 sec

……うーむ,まぁ速くはなったっぽい?

というか,実行時間が結構ばらつきます.実行環境はPS3なんだけど,ハードウェアの特性かな?最長時間を見ると

  • SIMD化する前: 16.74284220 sec
  • SIMD化した後: 15.84491611 sec

……最短との差は3割.まぁそんなものかしら?


SIMDの効果ってこんなものかなあと疑問を持ったので同じものをSSEあたりでやってみようと思ったけど,何を参考にどうくめばいいのかわからず.1時間以上webを漁ったけどよくわからない.また明日以降に試す.

というか,おとなしくPPE向けのSIMD化でもやろう,うん.

tgbt Cell, プログラミング一般

このエントリをはてなブックマークに追加このエントリをdel.icio.usに追加このエントリをLivedoor Clipに追加このエントリをYahoo!ブックマークに追加このエントリをFC2ブックマークに追加このエントリをNifty Clipに追加このエントリをPOOKMARK. Airlinesに追加このエントリをBuzzurl(バザール)に追加このエントリをChoixに追加このエントリをnewsingに追加