PRACE NPT

#include <stdio.h>

#include <stdlib.h>

double cclock( void );

int    check( int m, int l, int n, double *c );

void   gendat( int, int, int, int, double *a, double *b );

void   mvm( int, int, int, int, double *a, double *b, double *c );

void   prthead( void );

void   prtspeed( int, int, int, double, int, int );

void   state( char * );

int main()

{

   int      lda, m, l, n;

   int      ok, nops, nrep;

   int      i;

   double   *a, *b, *c;

   double   time;

   FILE     *inl;

// ------------------------------------------------------------------------

   state( "mod2a" );

   prthead();

   inl = fopen( "mod2am.in", "r" );

   while( ( fscanf( inl, "%d%d%d%d\n", &m, &l, &n, &nrep ) != EOF ) ){

      lda = l + 1;

      a = calloc( m*lda, sizeof( double ) );

      b = calloc( l*n, sizeof( double ) );

      c = calloc( m*n, sizeof( double ) );

      gendat( lda, m, l, n, a, b );

      time = cclock();

      for( i = 0; i < nrep; i++ ){

         mxm( lda, m, l, n, a, b, c );

      }

      time = cclock() - time;

      ok   = check( m, l, n, c );

      time = time/nrep;

      nops  = 2*m*l*n;

      prtspeed( m, l, n, time, ok, nops );

      free( a ); free( b ); free( c);

   }

   printf( "-------------------------------------------------\n" );

}