/****************  m_amatvec_ns.c  (in u1.a) *******************************
*				  					   *
* void mult_adj_u1_mat_vec_nsum(u1_matrix *a, u1_vector *b, u1_vector *c)
* adj u1_matrix times u1_vector multiply and subtract from another       *
*  C  <- C - A_adjoint*B						   *
*/
#include "complex.h"
#include "u1.h"

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