/*******************  trace_u1.c  (in u1.a) ***************************
*									*
* float trace_u1(a) u1_matrix *a;					*
* return trace of a U1 matrix
*/

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

/* trace of an U1 matrix */
complex trace_u1(u1_matrix *a) {
    return(cmplx(a->e[0], a->e[1]));
}