/******************  reunit_u1.c  (in u1.a) ***************************
*									*
* void reunit_u1( u1_matrix *a )					*
* reunitarize u1_matrix so that  it is a phase
*/

#include "complex.h"
#include "u1.h"
double sqrt(double);

void reunit_u1(u1_matrix *a) {
  register int i;
  float asq=0.;

  asq = (double) 1./sqrt((a->e[0])*(a->e[0]) + (a->e[1])*(a->e[1]));
  for(i=0; i<2; ++i) a->e[i] *= asq;
}