/****************  m_mv_s_4dir.c  (in u1.a) *****************************
*									 *
* void mult_u1_mat_vec_sum_4dir(u1_matrix *a, u1_vector *b0,          *
*        u1_vector *b1, u1_vector *b2, u1_vector *b3, u1_vector *c); *
* Multiply the elements of an array of four u1_matrices by the          *
* four u1_vectors, and add the results to produce a single              *
* u1_vector.                                                            *
*                                                                        *
* C <- A[0]*B[0]+A[1]*B[1]+A[2]*B[2]+A[3]*B[3]                           * 
*/
#include "complex.h"
#include "u1.h"

void mult_u1_mat_vec_sum_4dir(u1_matrix *a, u1_vector *b0,
       u1_vector *b1, u1_vector *b2, u1_vector *b3, u1_vector *c) {

   mult_u1_mat_vec(a+0,b0,c);
   mult_u1_mat_vec_sum( a+1,b1,c);
   mult_u1_mat_vec_sum( a+2,b2,c);
   mult_u1_mat_vec_sum( a+3,b3,c);
}