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

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