using System;

using System.Diagnostics;

​

            int MFLOP = 50;

            int SIZEA = 25;

            int SIZEB = 40;

            Stopwatch  _timer = Stopwatch.StartNew();

            _timer.Stop();

            Random autoRand = new Random( );

​

​

            double[] vector_out = new double[SIZEA];

            double[] vector = new double[SIZEB];

            double[,] matrix = new double[SIZEA, SIZEB];

​

            for (int n = 1; n <= MFLOP*1000; n++)

            {

                for (int i = 0; i < SIZEA; i++)

                {

                    vector_out[i] = autoRand.NextDouble();

                    vector[i] = autoRand.NextDouble();

                    for (int j = 0; j < SIZEB; j++)

                        matrix[i, j] = autoRand.NextDouble();

                }

​

                _timer.Start();

                for (int i = 0; i < SIZEA; i++)

                    for (int j = 0; j < SIZEB; j++)

                        vector_out[i] = vector_out[i] + vector[j] * matrix[i, j];

                _timer.Stop();

            }

            Console.WriteLine("Msecs: " + _timer.ElapsedMilliseconds);

            var secs = Convert.ToDouble(_timer.ElapsedMilliseconds)/1000;

            Console.WriteLine("MFLOPS: " + MFLOP/secs);