///////////////////////////////////////////////////////////////////////////////
//
// 	FILE: dotprod.c
//
///////////////////////////////////////////////////////////////////////////////

const int N = 20;

///////////////////////////////////////////////////////////////////////////////
// double a_dot_b( double*, double* )
///////////////////////////////////////////////////////////////////////////////

double a_dot_b( double *a, double *b )
{
    int i;
    double output = 0.0;

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

    return( output );
}

///////////////////////////////////////////////////////////////////////////////
// double a_dot_c( double*, double* )
///////////////////////////////////////////////////////////////////////////////

double a_dot_c( double *a, double *c )
{
    int i;
    double output = 0.0;

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

    return( output );
}

///////////////////////////////////////////////////////////////////////////////
// double a_dot_d( double*, double* )
///////////////////////////////////////////////////////////////////////////////

double a_dot_d( double *a, double *d )
{
    int i;
    double output = 0.0;

    for (i=0; i<N/2; i++) 
    {
		output += ( a[2*i] * d[i] );
	}

    return( output );
}