/*******************  m_mat_nn.c  (in u1.a) ****************************
*									*
* void mult_u1_nn(u1_matrix *a, u1_matrix *b, u1_matrix *c)		*
* matrix multiply, no adjoints 						*
* C  <-  A*B								*
*/

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

void mult_u1_nn(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];
}