Code:

int main( void ){BEGIN_PROFILE( "NON_SSE_CROSS" );v = cross( v1, v2 );END_PROFILE();BEGIN_PROFILE( "SSE_CROSS" ); v = cross_sse( v1, v2 );END_PROFILE();if ( time_non_sse < time_sse )  printf( "I did a pretty good job!" );return 0;}

Code:

#pragma oncestatic const float SZero = 0; __declspec(align(16)) struct Vector3D{#pragma region Constructorpublic: Vector3D(){X = 0;Y = 0;Z = 0;W = 0;}public: Vector3D(float x, float y, float z){X = x;Y = y;Z = z;W = 0;}#pragma endregion#pragma region Destructorpublic: ~Vector3D(void){}#pragma endregion#pragma region Fieldspublic: float X;public: float Y;public: float Z;private: float W;#pragma endregion#pragma region Properties#pragma endregion#pragma region Operator Overloads#pragma endregion#pragma region Methodspublic: void NormalizePrecise(){_asm{//get lengthmoveax, thismovdxmm0, dword ptr ds:[eax]movdxmm1, dword ptr ds:[eax + 4]movdxmm2, dword ptr ds:[eax + 8]mulssxmm0, xmm0mulssxmm1, xmm1mulssxmm2, xmm2addssxmm0, xmm1addssxmm0, xmm2sqrtssxmm0, xmm0ucomissxmm0, SZerojeDone//duplicate length across registerpshufdxmm0, xmm0, 0//divide by lengthmovdqaxmm1, xmmword ptr ds:[eax]divpsxmm1, xmm0//store resultmovdqaxmmword ptr ds:[eax], xmm1Done:}}public: void Normalize(){_asm{//get recipricol lengthmoveax, thismovdxmm0, dword ptr ds:[eax]movdxmm1, dword ptr ds:[eax + 4]movdxmm2, dword ptr ds:[eax + 8]mulssxmm0, xmm0mulssxmm1, xmm1mulssxmm2, xmm2addssxmm0, xmm1addssxmm0, xmm2ucomissxmm0, SZerojeDonersqrtssxmm0, xmm0//duplicate length across registerpshufdxmm0, xmm0, 0//multiply by recipricol length (divide by length)movdqaxmm1, xmmword ptr ds:[eax]mulpsxmm1, xmm0//store resultmovdqaxmmword ptr ds:[eax], xmm1Done:}}public: void Absolute(){_asm{moveax, thisanddword ptr ds:[eax], 07FFFFFFFhanddword ptr ds:[eax + 4], 07FFFFFFFhanddword ptr ds:[eax + 8], 07FFFFFFFh}}public: void Maximize(const Vector3D &v1){_asm{//get parameter informationmoveax, v1movdqaxmm0, xmmword ptr ds:[eax]moveax, this//compute maximummaxpsxmm0, xmmword ptr ds:[eax]//store resultmovdqaxmmword ptr ds:[eax], xmm0}}public: void Minimize(const Vector3D &v1){_asm{//get parameter informationmoveax, v1movdqaxmm0, xmmword ptr ds:[eax]moveax, this//compute minimumminpsxmm0, xmmword ptr ds:[eax]//store resultmovdqaxmmword ptr ds:[eax], xmm0}}public: void Cross(const Vector3D &v1){_asm{//get parameter informationmoveax, v1movdqaxmm0, xmmword ptr ds:[eax]moveax, thismovdqaxmm1, xmmword ptr ds:[eax]//align vectors to be multipliedpshufdxmm2, xmm0, 11001001b//(v1.Y, v1.Z, v1.X)pshufdxmm4, xmm1,11001001b//(Y, Z, X)pshufdxmm3, xmm0, 11010010b//(v1.Z, v1.X, v1.Y)pshufdxmm5, xmm1,11010010b//(Z, X, Y)//perform cross-productmulpsxmm4, xmm3mulpsxmm5, xmm2subpsxmm4, xmm5//store resultmovdqaxmmword ptr ds:[eax], xmm4}}public: void Lerp(const Vector3D &v1, float interpolator){_asm{//get parameter informationmoveax, v1movdqaxmm0, xmmword ptr ds:[eax]moveax, thismovdqaxmm1, xmmword ptr ds:[eax]//duplicate interpolator across registerpshufdxmm2, interpolator, 0//interpolatesubpsxmm0, xmm1mulpsxmm0, xmm2addpsxmm0, xmm1//store resultmovdqaxmmword ptr ds:[eax], xmm0}}#pragma endregion};