/******************  cs_m_s_vec.c  (in u1.a) ***************************
*									*
*  c_scalar_mult_sub_u1vec():						*
*  multiply an u1 vector by a complex scalar and subtract from another	*
*  vector:   v1 <- v1 - number*v2 					*
*  NOTE: non-standard order of target and source v1 and v2!             *
*/
#include "complex.h"
#include "u1.h"

void c_scalar_mult_sub_u1vec(u1_vector *v1, complex *phase, u1_vector *v2) {
  complex t;

    CMUL(v2->c, *phase, t);
    CSUB(v1->c, t, v1->c);
}