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

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

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