/********************  submat.c (in u1.a)  *****************************
*									*
*  Subtract two U1 matrices 						*
*/
#include "complex.h"
#include "u1.h"

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