/******************  u1_rdot.c  (in u1.a) ******************************
*                                                                       *
* float u1_rdot(u1_vector *a, u1_vector *b);                        *
* return dot product of two u1_vectors 
     = a^* b

     same as u1_dot
*/
#include "complex.h"
#include "u1.h"

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