/*****************  u1_anti_proj.c  (in u1.a) *************************
*									*
* void u1_antiherm_projector(u1_vector *u, u1_vector *v, u1_matrix *m);
* M  <- outer product of U and V, projected onto traceless antihermitian M
*  M_ij = U_i * V_adjoint_j |_projected					*
*/
#include "complex.h"
#include "u1.h"

void u1_antiherm_projector(u1_vector *u, u1_vector *v, u1_matrix *m) {
  
  u1_projector(u,v,m);
  m->e[0] = 0;
}