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

/* a <- a + s*b, vectors */
#ifdef PROTO
void scalar_mult_sum_su3_vector(su3_vector *a, su3_vector *b, float s){
#else
void scalar_mult_sum_su3_vector(a,b,s) su3_vector *a,*b; float s;{
#endif
register int i;
    for(i=0;i<3;i++){
	a->c[i].real += s*b->c[i].real;
	a->c[i].imag += s*b->c[i].imag;
    }
}