/*****************  make_ahmat.c  (in u1.a) ****************************
*									*
* void make_anti_hermitian( u1_matrix *m2, anti_hermitmat *ah2)	*
* take the traceless and anti_hermitian part of an u1 matrix 		*
*/
#include "complex.h"
#include "u1.h"

void make_anti_hermitian(u1_matrix *m2, u1_anti_hermitmat *ah2) {
  *ah2 = *m2;
  ah2->e[0] = 0;  
}