// mat_vec.c  compute matrix times vector multiply
#include <stdio.h>

void mat_vec(long int n, double A[][n], double x[], double y[])
{
  long int i, j;

  printf("mat_vec n=%ld \n",n);
  for(i=0; i<n; i++)
  {
    y[i] = 0.0;
    for(j=0; j<n; j++)
    {
      y[i] = y[i] + A[i][j]*x[j];
    } // end j
  } // end i
} // end mat_vec

int main()
{
  long int n3 = 3;
  double A1[3][3] = {{1.0,2.0,3.0},{4.0,5.0,6.0},{7.0,8.0,9.0}}; // A matrix
  double x1[3] = {1.0,2.0,3.0};          // x vector
  double y1[3];

  printf("mat_vec_64.asm running \n"); 
  mat_vec(n3, A1, x1, y1);
  printf("mat_vec returned =%f %f %f \n",y1[0], y1[1], y1[2]); 
  return 0;
} // end mat_vec.c