/*******************  m_mat_aa.c  (in u1.a) ****************************
*									*
* void mult_u1_aa(u1_matrix *a, u1_matrix *b, u1_matrix *c)		*
* matrix multiply, with both adjoints 						*
* C  <-  A^dagger * B^dagger						*
*/

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

void mult_u1_aa(u1_matrix *a, u1_matrix *b, u1_matrix *c) {
  c->e[0] =  a->e[0]*b->e[0] - a->e[1]*b->e[1];
  c->e[1] = -a->e[0]*b->e[1] - a->e[1]*b->e[0];
}