/******************  s_m_s_vec.c  (in u1.a) ******************************
*									  *
* void scalar_mult_sub_u1_vector(u1_vector *a, u1_vector *b,           *
                                                 float s, u1_vector *c); *
* C <- A - s*B,  A, B and C vectors 						  *
*/
#include "complex.h"
#include "u1.h"

void scalar_mult_sub_u1_vector(u1_vector *a, u1_vector *b,
                                float s, u1_vector *c) {
      c->c.real = a->c.real - s*b->c.real;
      c->c.imag = a->c.imag - s*b->c.imag;
}