home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c072 / 1.ddi / PRG8_2.C < prev    next >
Encoding:
C/C++ Source or Header  |  1987-09-19  |  1.3 KB  |  48 lines
  1. /*Program 8_2 - Floating point benchmark
  2.   by Stephen R. Davis
  3.  
  4.   The matrix multiply is something of a floating point
  5.   performance standard.  The size may be varied at will, but
  6.   the default is 50 x 50.  To execute this program without
  7.   using the 8087 in a machine equipped with one, issue the
  8.   DOS command SET 87=N before starting.  SET 87=Y will reenable
  9.   the 8087.
  10. */
  11.  
  12. #include <stdio.h>
  13. #include <stdlib.h>
  14.  
  15. #define msize 50
  16.  
  17. /*define our data requirements*/
  18. unsigned i, j, k;
  19. float x [msize][msize], y [msize][msize], z[msize][msize];
  20. double accum;
  21.  
  22. /*Main - perform the benchmark*/
  23. main ()
  24. {
  25.      /*begin by initializing the matrices*/
  26.      printf ("Initializing matrices\n");
  27.      for (i = 0; i < msize; i++)
  28.           for (j = 0; j < msize; j++) {
  29.                x [i][j] = rand () / 10000.;
  30.                y [i][j] = rand () / 10000.;
  31.           }
  32.  
  33.      /*now begin the benchmark*/
  34.      printf ("Enter return and start stop watch at same time\n");
  35.      getchar ();
  36.      printf ("Start...");
  37.  
  38.      for (i = 0; i < msize; i++)
  39.           for (j = 0; j < msize; j++) {
  40.                accum = 0;
  41.                for (k = 0; k < msize; k++)
  42.                    accum += x [i][k] * y [k][j];
  43.                z [i][j] = accum;
  44.           }
  45.  
  46.       printf ("stop!\n");
  47. }
  48.