import std.getopt;
import std.stdio;
import std.conv : to;
import std.datetime;

immutable int n = 1000, l=1000, m=1000;

auto A = new double [1000][m];
auto B = new double [1000][n];
auto C = new double [1000][n];
void main(string[] args)
{
  matrix!(typeof(A)).ini(A);
  matrix!(typeof(A)).ini(B);
  run_test();

  auto r = benchmark!(run_test)(500);

  auto res = to!Duration(r[0]);
  writefln(res.toString());

}

void run_test(){
    matrix!(typeof(A)).MatMul(A,B,C);
}



template matrix(T){
  void ini(T mtx){
    for(int i = 0; i <mtx.length; i++) {
      for(int j = 0; j< mtx[i].length; j++) {
    mtx[i][j] = i*j + 0.6*j +3.4;
      }
    }
  }

  void print(T mtx){
    for(int i = 0; i <mtx.length; i++) {
      writefln("%(%8.1f %)", mtx[i]);

    }
  }


  void MatMul(T A, T B, T C){
    for(int k = 0; k < C[0].length; k++) {
      for(int i = 0; i < C.length; i++) {
    C[i][k] = 0;
    for(int j = 0; j < A[0].length; j++) {
      C[i][k] += A[i][j] * B[j][k];
    }
      }
    }
  }

}