#include <iostream>

#include "fflas-ffpack/ffpack.h"
#include "fflas-ffpack/modular-balanced.h"
#include "timer.h"
#include "Matio.h"
 
using namespace std;

int main(int argc, char** argv) {

  // parameter: p, n, iteration, file

  float    p    = atof(argv[1]);
  int n    = atoi(argv[2]);
  size_t iter = atoi(argv[3]);
  

  typedef ModularBalanced<double> Field;
  //  typedef ModularBalanced<float> Field;
  typedef Field::Element Element;
      
  Field F(p);

  Timer chrono;
  double time=0.0;
  int singular;
  
  Element *A;

  for (size_t i=0;i<iter;++i){
    
    A = new Element[n*n];
    Field::RandIter G(F);
    for (size_t i=0; i< n*n; ++i)
      G.random(*(A+i));      
    
    size_t * P = new size_t[n];
    size_t * Q = new size_t[n];
    
    chrono.clear();
    chrono.start();
    FFPACK::LUdivine (F, FFLAS::FflasNonUnit, FFLAS::FflasNoTrans, n, n, A, n,
		      P, Q);
    chrono.stop();
 
    time+=chrono.realtime();
    delete[] P;
    delete[] Q;
    delete[] A;
    
  }

  cerr<<"n: "<<n<<" p: "<<p<<std::endl
      <<" time:  "<<time/iter<<std::endl
      <<" speed: "<<2/3.0*n/1000.0*n/1000.0*n/1000.0/time*iter<<" Gffops"<<std::endl;


  return 0;
}