///////////////////////////////////////////////////////////////////////////////
//
// 	FILE: dotprod_main.c
//
///////////////////////////////////////////////////////////////////////////////

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

///////////////////////////////////////////////////////////////////////////////
// global variables
///////////////////////////////////////////////////////////////////////////////

double pm a[] = {
	
   9.9619469809174555e-01,
   9.8480775301220802e-01,
   9.6592582628906831e-01,
   9.3969262078590843e-01,
   9.0630778703664994e-01,
   8.6602540378443871e-01,
   8.1915204428899180e-01,
   7.6604444311897801e-01,
   7.0710678118654757e-01,
   6.4278760968653936e-01,
   8.7155742747658166e-02,
   1.7364817766693033e-01,
   2.5881904510252074e-01,
   3.4202014332566871e-01,
   4.2261826174069944e-01,
   4.9999999999999994e-01,
   5.7357643635104605e-01,
   6.4278760968653925e-01,
   7.0710678118654746e-01,
   7.6604444311897801e-01 
};

double b[] = {
	
	-8.7155742747658232e-03,
	-1.7364817766693030e-02,
	-2.5881904510252084e-02,
	-3.4202014332566873e-02,
	-4.2261826174069934e-02, 
	-4.9999999999999975e-02,
	-5.7357643635104580e-02,
	-6.4278760968653939e-02,  
	-7.0710678118654752e-02,  
	-7.6604444311897793e-02,  
	9.9619469809174549e-02,  
	9.8480775301220805e-02,  
	9.6592582628906826e-02,  
	9.3969262078590840e-02,  
	9.0630778703665010e-02,  
	8.6602540378443865e-02,  
	8.1915204428899199e-02,  
	7.6604444311897807e-02,  
	7.0710678118654752e-02,
	6.4278760968653953e-02 
};

double c[] = {
	
	6.4278760968653939e-02,  
	5.7357643635104608e-02,   
	5.0000000000000010e-02,   
	4.2261826174069941e-02,   
	3.4202014332566880e-02,   
	2.5881904510252074e-02,   
	1.7364817766693040e-02,   
	8.7155742747658145e-03,   
	6.1232339957367663e-18,  
	-8.7155742747658232e-03,   
	7.6604444311897807e-02,   
	8.1915204428899185e-02,   
	8.6602540378443865e-02,   
	9.0630778703664996e-02,   
	9.3969262078590826e-02,   
	9.6592582628906826e-02,   
	9.8480775301220805e-02,   
	9.9619469809174549e-02,   
	1.0000000000000001e-01,   
	9.9619469809174549e-02
};

double d[] = {
	
	-9.9619469809174549e-02,
	-9.6592582628906840e-02,
	-9.0630778703665010e-02,
	-8.1915204428899199e-02,
	-7.0710678118654766e-02,
	-8.7155742747657937e-03,
	-2.5881904510252036e-02,
	-4.2261826174069927e-02,
	-5.7357643635104580e-02,
	-7.0710678118654752e-02 
};

///////////////////////////////////////////////////////////////////////////////
// function prototypes
///////////////////////////////////////////////////////////////////////////////

extern double a_dot_b( double pm *, double * );
/*extern double a_dot_c( double pm *, double * );*/
extern double a_dot_c_asm( double pm *, double * ); 
extern double a_dot_d( double pm *, double * ); 

///////////////////////////////////////////////////////////////////////////////
// void main()
///////////////////////////////////////////////////////////////////////////////

void main()
{
	int i;
	double result[3] = {0};
    
    result[0] = a_dot_b( a, b );
/*     result[1] = a_dot_c( a, c );*/ 
    result[1] = a_dot_c_asm( a, c );
    result[2] = a_dot_d( a, d );

    for( i=0; i<3; i++ ) 
    {
		printf( "Dot product [%d] = %f\n", i, result[i] );
	}
}