/******************  su3_dot.c  (in su3.a) ******************************
*									*
* complex su3_dot(a,b) su3_vector *a,*b;				*
* return dot product of two su3_vectors					*
*/
#include "complex.h"
#include "su3.h"

complex su3_dot(a,b) su3_vector *a,*b; {
complex temp1,temp2;
    CMULJ_(a->c[0],b->c[0],temp1) 
    CMULJ_(a->c[1],b->c[1],temp2)
    CSUM(temp1,temp2);
    CMULJ_(a->c[2],b->c[2],temp2)
    CSUM(temp1,temp2);
    return(temp1);
}