/****************  m_amatvec_s.c  (in u1.a) ********************************
*									    *
* void mult_adj_u1_mat_vec_sum(u1_matrix *a, u1_vector *b, u1_vector *c)*
* adj u1_matrix times u1_vector multiply and add to another u1_vector    *
*  C  <- C + A_adjoint*B						    *
*/
#include "complex.h"
#include "u1.h"

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