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

/* c <- a - s*b, vectors */
#ifdef PROTO
void scalar_mult_sub_su3_vector(su3_vector *a,su3_vector *b,float s,
	su3_vector *c){
#else
void scalar_mult_sub_su3_vector(a,b,s,c) su3_vector *a,*b,*c; float s;{
#endif

register int i;
    for(i=0;i<3;i++){
	c->c[i].real = a->c[i].real - s*b->c[i].real;
	c->c[i].imag = a->c[i].imag - s*b->c[i].imag;
    }
}