// a.cpp
#include <stdexcept>
char array[128 * 1024]; // 128 KB
main()
{
    std::runtime_error e("e");
}

$ spu-g++ -Os -D NDEBUG a.cpp
/usr/lib/gcc/spu/4.1.1/../../../../spu/bin/ld: .bss exceeds local store range

みたいな何の変哲もないコードで SPE のサイズをオーバーしちゃってハマってたのね。

　で、原因を調べてみたら、どうやら原因は std::runtime_error が依存している std::string にあるらしい事が判明。

// a2.cpp
main()
{
}

$ spu-g++ -Os -D NDEBUG a2.cpp
$ ll a.out
-rwxr-xr-x 1 yoffy yoffy 97K 2008-08-30 11:01 a.out

// a3.cpp
#include <string>
main()
{
    std::string s;
}

$ spu-g++ -Os -D NDEBUG a3.cpp
$ ll a.out
-rwxr-xr-x 1 yoffy yoffy 167K 2008-08-30 11:04 a.out

　std::allocator を少し疑ったんだけど、std::vector は問題ないので std::string 側に原因がありそうな気がする。面倒くさいので細かくは追ってない。
　別に stdexcept や string が無くても困らないけど、SPE 用にもう少し考えてくれてもいいんじゃないだろうか。

// a4.c
main()
{
}

$ spu-gcc -Os -D NDEBUG a4.c
$ ll a.out
-rwxr-xr-x 1 yoffy yoffy 4.9K 2008-08-30 11:34 a.out

　…。
　何もしてないのに 97 KB も使うっておかしいなあとは思ってたんだよ。これ、libstdc++ が怠慢じゃない？コンパイラオプションで例外とか止めたら小さくできるのかな。

NSColor * c = [NSColor whiteColor];
float r = [c recComponent];

　こんな単純なコードが動かない。
　今日重い腰を上げてドキュメントを読んでみたら、「RGBカラーでなければ例外を投げる」と書いてあった。whiteColor は RGB カラーではなくグレースケールカラーとして定義されているらしい。
　過去のアプリケーションとの互換性低すぎないか？これ。

　ということで、理由は分かったけれどもシステム標準の色からいつでも安全に RGB を得たい場合はどうすればいいんだろう。
　わざわざ変換しなくてはならないとか？
　それともシステムカラーは使えないものとして [NSColor colorWithCalibratedRed:green:blue:] で自力で解決するとか？

#include <altivec.h>

vector unsigned int test(vector unsigned int a, vector unsigned int b)
{
	vector unsigned int c;
	c = vec_cmpeq(a, b);
	return c;
}

$ g++ -maltivec test.cpp
test.cpp: In function 'unsigned int __vector__ test(unsigned int __vector__, unsigned int __vector__)':
test.cpp:6: error: cannot convert 'unsigned int __vector__' to 'unsigned int __vector__' in assignment

　日本語でおk

　そのまま自分の環境に当てはめると、これでいけた。

　これは畳み込みの一種であり、例えばよく使われている移動平均なんかは次のような感じで畳み込み

元画像と移動平均

1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0

ガウスぼかしは次のようなので畳み込む。

ガウスぼかし

0.37 0.6 0.37 0.6 1.0 0.6 0.37 0.6 0.37

　Non Local Means は畳み込む画素の周囲でマッチングを行って、その一致率でウエイトを作る。似通った領域では広く平均が取られ、似通っていない領域では平均されない。
　一致率の出し方が工夫のしどころではあるけど、NL-Means filter では差の二乗和を d、係数を h とすると

exp(-d*h)

というように指数に入れている。差が無ければウエイトが 1 になり、差が少しでも大きくなるとすぐ 0 に落ち込む。
　また、h が 0 であれば移動平均になり、h が大きくなるにつれ似通った領域のみ使われるようになる。最終的に h が無限大になれば元の画像に一致する。

h が小さい場合

h が大きい場合

　この手法の特徴は、パターンが一致している部分のみを使うので、エッジが完全に保存される。メディアンフィルタでは落ちてしまうような一本の線でさえも保存される。
　しかし、逆に一点のノイズであっても保存されてしまうので、結局強いノイズは予め他の手法で取り除いておかないといけない。

メディアンフィルタ 移動平均に比べてエッジが保存されている

　Non Local Means のみで強いノイズを除去するには h の値を小さくしないといけないが、それだと移動平均に近づくのでボケる。
　h を小さくしつつボケにくくする為には、畳み込みのサイズを小さくするしか無い。それでは広く平均を取れるこの手法の利点が薄くなる。

　そこで折衷案として、h を 0 にした時にガウスぼかしになるようにしてみた。
　こうすると領域は円形になるし、サイズを小数単位で微調整出来る。何より遠くにいく程弱くなるので、広過ぎる範囲で平均されない。

ガウスぼかしに近づく方法で h を小さくした場合

　うーん。SAMSUNG の文字は多少読みやすくなってるけど、あまり差が無いかな。