/******************  u1_dot.c  (in u1.a) ******************************
*                                                                       *
* complex su3_dot(u1_vector *a, u1_vector *b);                        *
* return dot product of two u1_vectors 
     = a^* b
*/
#include "complex.h"
#include "u1.h"

float u1_dot(u1_vector *a, u1_vector *b) {
   return(a->c.real * b->c.real + a->c.imag * b->c.imag);
}