float Vec::length() const
{
float len = x*x+y*y+z*z;
#ifdef _MSC_VER
__asm
{
fld len
fsqrt
fstp len
}
#else
len = sqrt(len);
#endif
return len;
}
速度重視のため、アセンブラ部分はかなりの部分をベタ書きのハードコーディングにしていたが、インライン関数に置き換えるなどしないとメンテしづらいことがわかった。
#include <iostream>
#include <vector>
#include <string>
typedef std::string String;
typedef std::vector<String> StringList;
/**
* 文字列を分割してコンテナに格納する
* @param dest 分割後の文字列格納先へのポインタ
* @param src 分割したい文字列への参照
* @param delim デリミタ文字列への参照
* @return 分割後の文字列格納先へのポインタ
*/
StringList *SpritString(
StringList *dest,
const String &src,
const String &delim
)
{
String::size_type start = 0;//デリミタを検索するインデクス番号
while(true){
//デリミタが現れる最初のインデクスを求める
String::size_type end = src.find(delim, start);
//デリミタが見つかった場合
if(end != String::npos){
dest->push_back(src.substr(start, end - start));
}
//デリミタが見つからなかった場合
else{
//文末までを格納して返す
dest->push_back(src.substr(start, src.length() - start));
break;
}
//次の開始地点へ移動
start = end + delim.length();
}
return dest;
}
/**
* 呼び出し側
*/
int main(){
using namespace std;
String src("hoge piyo fuga");//区切りたい文字列
String delim(" ");//デリミタ（区切り文字）
StringList out;//区切られた文字が配列となって帰ってくる
SpritString(&out, src, delim);//空白文字で区切る
//出力
for(String::size_type i=0; i<out.size(); ++i){
cout << out[i] << endl;
}
/* 出力結果---------------------------
> hoge
> piyo
> fuga
*/
return 0;
}
参考先のサイトでは関数内部で作成したリストオブジェクトをそのまま返していたので、軽くする目的でポインタ経由での受け渡しにしてみた。そもそも軽くしたいならstd::stringなんか使うなって話だが…。便利なんだもん。/**
* floatであらわされる色構造体
*/
struct FRGB{
float r,b,g;
static FRGB *Set(FRGB *out, unsigned char r, unsigned char g, unsigned char b);
};
/**
* 自作マテリアル構造体
*/
struct Material{
FRGB diffuse;//拡散反射光係数
FRGB ambient;//環境光係数
FRGB specular;//鏡面反射光係数
FRGB transmit;//透過係数
};
透過係数は、0~1で定義される透明度。0の時は光を一切通さず、1の時は完全に光を通す（＝透明となる）ものとする。
入射光をIpこの時、自然数nの値を大きくすればするほど、ハイライトの広がりは抑えられて、シャープな印象になるそうです。
鏡面反射係数をks
レイの正反射方向と視線のなす角をa
レイの正反射方向ベクトルをL'
視線へのベクトルをV
と置いたときの表面輝度Isは
Is = ks*Ip*cosna = ks*Ip*(dot(L', V))n
となる
入射ベクトルをL最終的な表面の輝度は【環境光＋拡散反射光＋鏡面反射光】になるので、ランバートシェーダの処理のあとに上記の処理を入れるだけでフォンのモデルは実装できちゃうわけですね。
表面の法線をNとすると
L' = 2*(dot(N,V))*N - L
Vector3 &Vector3::normalize()
{
float len = x*x+y*y+z*z;//ノルムの平方
//以下len = sqrt(len);と同じ
__asm{
fld len
fsqrt
fstp len//お行儀が悪いがlenを使いまわす
}
//0割防止
if(len < FLT_EPSILON){
x=y=z=0;
}
else{
float denom = 1.f / len;
x *= denom;
y *= denom;
z *= denom;
}
return *this;
}
上記のif文の評価式をlen >= EPSILONにして、ブロックを入れ替えた方がジャンプが減る分確率的に速くなる…気がする。てかここでごっそり速くなった。
if(len >= FLT_EPSILON){
float denom = 1.f / len;
x *= denom;
y *= denom;
z *= denom;
}
else{
x=y=z=0;
}