/******************  u1_adjoint.c  (in u1.a) **************************
*                                                                       *
* void u1_adjoint(a,b) u1_matrix *a,*b;                               *
* B  <- A_adjoint,  adjoint of an U1 matrix                            *
*/

#include "complex.h"
#include "u1.h"

/* adjoint of an U1 matrix */
void u1_adjoint(u1_matrix *a, u1_matrix *b) {
   b->e[0] = a->e[0]; 
   b->e[1] = -a->e[1];
}