/*******************  cs_m_a_mat.c (in u1.a)  **************************
*									*
*  c_scalar_mult_add_u1mat(u1_matrix *a, u1_matrix *b,               * 
*				  complex phase, u1_matrix *c):        *
*  multiply an u1 matrix by a complex scalar and add it to another     *
*  matrix:   c <- a + phase*b                                           *
*/

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

void c_scalar_mult_add_u1_matrix(u1_matrix *a, u1_matrix *b, 
				              complex phase, u1_matrix *c) {
  c->e[0] = a->e[0] + phase.real*b->e[0] - phase.imag*b->e[1];
  c->e[1] = a->e[1] + phase.imag*b->e[0] + phase.real*b->e[1];
}