/*****************  m_amv_4dir.c  (in u1.a) *****************************
*									*
*  void mult_adj_u1_mat_vec_4dir(u1_matrix *mat, u1_vector *src,     *
*                                         u1_vector *dest);            *
*  Multiply an su3_vector by an array of four adjoint u1_matrices,	*
*  result in an array of four u1_vectors.				*
*  dest[i]  <-  A_adjoint[i] * src					*
*/
#include "complex.h"
#include "u1.h"

void mult_adj_u1_mat_vec_4dir(u1_matrix *mat, u1_vector *src,
                                                     u1_vector *dest) {
    mult_adj_u1_mat_vec( mat+0, src, dest+0 );
    mult_adj_u1_mat_vec( mat+1, src, dest+1 );
    mult_adj_u1_mat_vec( mat+2, src, dest+2 );
    mult_adj_u1_mat_vec( mat+3, src, dest+3 );
}