/********************  s_m_s_mat.c (in u1.a)  **************************
*									*
* void scalar_mult_sub_u1_matrix(a,b,s,c) u1_matrix *a,*b,*c; float s *
* C <- A + s*B                                                          *
*/

#include "complex.h"
#include "u1.h"

void scalar_mult_sub_u1_matrix(u1_matrix *a, u1_matrix *b, float s, 
		    u1_matrix *c) {
  register int i;
  for(i=0;i<2;i++) c->e[i] = a->e[i] - s*b->e[i];
}