/****************  m_mv_s_4dir.c  (in su2.a) *****************************
*									 *
* void mult_su2_mat_vec_sum_4dir(su2_matrix *a, su2_vector *b0,          *
*        su2_vector *b1, su2_vector *b2, su2_vector *b3, su2_vector *c); *
* Multiply the elements of an array of four su2_matrices by the          *
* four su2_vectors, and add the results to produce a single              *
* su2_vector.                                                            *
*                                                                        *
* C <- A[0]*B[0]+A[1]*B[1]+A[2]*B[2]+A[3]*B[3]                           * 
*/
#include "complex.h"
#include "globaldefs.h"
#include "su2.h"

void mult_su2_mat_vec_sum_4dir(su2_matrix *a, su2_vector *b0,
       su2_vector *b1, su2_vector *b2, su2_vector *b3, su2_vector *c) {

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