/****************  m_matvec.c  (in u1.a) *******************************
*									*
* void mult_u1_mat_vec(u1_matrix *A, u1_vector *b, u1_vector *c)	*
* matrix times vector multiply, no adjoints 				*
*  c  <-  A*b								*
*/
#include "complex.h"
#include "u1.h"

void mult_u1_mat_vec(u1_matrix *a, u1_vector *b, u1_vector *c) {

   c->c.real = a->e[0]*b->c.real - a->e[1]*b->c.imag;
   c->c.imag = a->e[1]*b->c.real + a->e[0]*b->c.imag;
}