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

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

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