/*****************  sub4vecs.c  (in u1.a) ******************************
*									*
*  Subtract four u1_vectors from an u1_vector				*
* void sub_four_u1_vecs(u1_vector *a, u1_vector *b1, u1_vector *b2, *
*                                  u1_vector *b3, u1_vector *b4)	*
* A  <-  A - B1 - B2 - B3 - B4						*
*/
#include "complex.h"
#include "u1.h"

void sub_four_u1_vecs(u1_vector *a, u1_vector *b1, u1_vector *b2, 
		       u1_vector *b3, u1_vector *b4) {
     CSUB( a->c, b1->c, a->c );
     CSUB( a->c, b2->c, a->c );
     CSUB( a->c, b3->c, a->c );
     CSUB( a->c, b4->c, a->c );
}