#include <stdio.h>    
#include <stdlib.h>
#include <omp.h>

/* This is a program that multiplies a square matrix A times a
   vector x to get a vector b.   
*/

int main(int argc, char* argv[])
{
    float  A[100][100];
    float  x[100];
    float  b[100];
    int mymax = 30000;
    int i, j, n;

/* INITIALIZATION: */

    n = 2;

/* Use rand and RAND_MAX from stdlib.h to generate random numbers for
   A and x.
   The function rand() returns a random integer between 0 and RAND_MAX.
   Divide rand() by RAND_MAX to get a number between 0 and 1,
   but tell C to convert from integers to floating point result by using
   the "casts" "(float)" in the statements below.
*/

    for (i=0; i<n; i++)
    {  x[i] = (float) rand() / (float) RAND_MAX;
       for (j=0; j<n; j++)
       {   A[i][j] = (float) rand() / (float) RAND_MAX;  
       }
    }

/* FORM THE PRODUCT b = A * x.  */

    for (i=0; i<n; i++)
    {   for (j=0; j<n; j++)
        {   b[i] += A[i][j]*x[j];
        }
    }

/* DISPLAY THE RESULTS.        */

    printf("\n");    /* Here is an example of how to display the results. */
    printf("A:\n");
    for (i=0; i<n; i++)
    {   for (j=0; j<n; j++)
        {   printf(" %f ",A[i][j]);
        }
        printf("\n");  /* end of line */
    }

    printf("\n");
    printf("x:\n");
    for (i=0; i<n; i++)
    {   printf(" %f \n",x[i]);
    }

    printf("\n");
    printf("b = A*x:\n");
    for (i=0; i<n; i++)
    {   printf(" %f \n",b[i]);
    }

}