/*****************  m_amatvec.c  (in u1.a) *****************************
*									*
*  void mult_adj_u1_mat_vec( u1_matrix *a, u1_vector *b,*c )		*
*  C  <-  A_adjoint * B							*
*/
#include "complex.h"
#include "u1.h"

void mult_adj_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;
}