/****************  cs_m_mat.c  (in su3.a) *******************************
*									*
* void c_scalar_mult_su3mat(b,s,c) su3_matrix *b,*c; complex *s		*
* C <- s*B,   B and C matrices 						*
*/
#include "complex.h"
#include "su3.h"

/* c <- s*b, matrices */
#ifdef PROTO
void c_scalar_mult_su3mat(su3_matrix *b,complex *s, su3_matrix *c){
#else
void c_scalar_mult_su3mat(b,s,c) su3_matrix *b,*c; complex *s;{
#endif

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