#include <stdio.h>
#include <sys/time.h>
#define _SHOW_SSE
//#define _SHOW_BASIC_ASM


double getcurrenttime()
{
    struct timeval timestruct;
    gettimeofday(&timestruct,NULL);
    return (double)( 1.0 *timestruct.tv_sec + (timestruct.tv_usec/1000000.0) );		
}

struct float4{
	float w;
	float x;
	float y;
	float z;
};

float A[4] = {10, 20, 30, 40};
float B[4] = {100, 200, 300, 400};
float C[4];



inline void mulfvec4(float *a, float *b, float *c)
{
	int i;

	for (i=0;i<4;i++)
	{
		c[i] = a[i] * b[i];
	}
}

#define mulfvec4_SSE(a, b, c) \
{ \
__asm__ __volatile__ ("movups %1, %%xmm0 \n\t" \
                      "movups %2, %%xmm1 \n\t" \
                      "mulps %%xmm0, %%xmm1 \n\t" \
                      "movups %%xmm1, %0 \n\t" \
                      :"=m" (c) \
                      :"m" (a), "m" (b)); \
}

#define MUL_TIMES 10000000

#ifdef _SHOW_SSE
int main(void)
{
    double non_SSE_start_time = getcurrenttime();

	for (long i=0;i<MUL_TIMES;i++)
	{
		for (int j=0;j<4;j++)
		{
			A[j] += 0.001;
			B[j] += 0.001;
		}
		mulfvec4(A, B, C);
	}
    double non_SSE_done_time = getcurrenttime();

	printf("done with mul in %f sec : %f, %f, %f, %f \n", non_SSE_done_time - non_SSE_start_time, C[0], C[1], C[2], C[3]);

	struct float4 AA, BB, CC;

	AA.w =  10;
	AA.x =  20;
	AA.y =  30;
	AA.z =  40;

	BB.w =  100;
	BB.x =  200;
	BB.y =  300;
	BB.z =  400;

    double SSE_start_time = getcurrenttime();

	for (long i=0;i<MUL_TIMES;i++)
	{
		AA.w +=  0.001;
		AA.x +=  0.001;
		AA.y +=  0.001;
		AA.z +=  0.001;

		BB.w +=  0.001;
		BB.x +=  0.001;
		BB.y +=  0.001;
		BB.z +=  0.001;
	
		mulfvec4_SSE(AA, BB, CC);
	}

	double SSE_done_time = getcurrenttime();

	printf("done with mul SSE in %f sec: %f, %f, %f, %f\n", SSE_done_time - SSE_start_time, CC.w, CC.x, CC.y, CC.z);

}
#endif


#ifdef _SHOW_BASIC_ASM
int main(void)
{
	int foo = 10, bar = 15;
	int clear = 0;
	__asm__ __volatile__("addl %%ebx, %%eax"
			     :"=b"(foo)
			     :"a"(foo), "b"(bar)
			    );

	printf("foo+bar = %d\n", foo);
	return 0;
}
#endif // _SHOW_BASIC_ASM