/******************  realtr.c  (in u1.a) *******************************
*									*
* float realtrace_u1(u1_matrix *a, u1_matrix *b)			*
* return Re( Tr( A_adjoint * B )  					*
*/

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

float realtrace_u1(u1_matrix *a, u1_matrix *b) {
  return( a->e[0]*b->e[0] + a->e[1]*b->e[1] );
}