#include <iostream.h>


extern void MatrixVectorMult(double out[], double Val[], int Row[], int Col[],
				double in[], int m, int nnz);


int main(int argc, char ** argv)
{
   double Val[] = {11, 12, 14, 22, 23, 25, 31, 33, 34, 42, 45,
		46, 55, 65, 66, 67, 75, 77, 78, 87, 88};
   int Row[] = {0, 0, 0, 1, 1, 1, 2, 2, 2, 3, 3,
		3, 4, 5, 5, 5, 6, 6, 6, 7, 7};
   int Col[] = {0, 1, 3, 1, 2, 4, 0, 2, 3, 1, 4,
		5, 4, 4, 5, 6, 4, 6, 7, 6, 7};
   double u[8], v[8];
   int i, j;

   u[0] = 3; u[1] = 5; u[2] = 2; u[3] = 1;
   u[4] = 0; u[5] = 1; u[6] = 2; u[7] = 4;


   MatrixVectorMult(v, Val, Row, Col,  u, 8, 21);

   for (i = 0; i < 8; i = i + 1)
   {
       printf("      %4.1lf\n", v[i]);
   }
}