/******************  cs_m_a_mat.c  (in su3.a) ***************************
*									*
*  c_scalar_mult_add_su3mat():					*
*  multiply an su3 matrix by a complex scalar and add it to another	*
*  matrix:   m3 <- m1 + number*m2 					*
*/
#include "complex.h"
#include "su3.h"

void c_scalar_mult_add_su3mat(m1,m2,phase,m3)
su3_matrix *m1,*m2,*m3; complex *phase; {
register int i,j;
complex t;
    for(i=0;i<3;i++)for(j=0;j<3;j++){
	t = cmul(&m2->e[i][j],phase);
	m3->e[i][j] = cadd(&m1->e[i][j],&t);
    }
}