/******************  magsq_u1vec.c  (in u1.a) *************************
*                                                                       *
* float magsq_u1vec(u1_vector *a);                                    *
* return squared magnitude of an U1 vector
   same as u1_dot(a,a)
*/
#include "complex.h"
#include "u1.h"

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