/****************  s_m_a_mat.c  (in su3.a) ******************************
*									*
* void scalar_mult_add_su3_matrix(a,b,s,c) su3_matrix *a,*b,*c; float s	*
* C <- A + s*B								*
*/
#include "complex.h"
#include "su3.h"

/* c <- a + s*b, matrices */
#ifdef PROTO
void scalar_mult_add_su3_matrix(su3_matrix *a,su3_matrix *b,float s,
	su3_matrix *c){
#else
void scalar_mult_add_su3_matrix(a,b,s,c) su3_matrix *a,*b,*c; float s;{
#endif

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