void bench_fp()
{
// const char * ch = "../TMP/f3Reduc1.tif";
const char * ch = "/home/data/mpd/Andalousie/andaReduc1.tif";
Tiff_Im tiff(ch);
ELISE_fp fp(ch, ELISE_fp::READ);
tiff.show();
INT NbTx = tiff.sz().x/ tiff.SzTile()[0];
INT NbTy = tiff.sz().y/ tiff.SzTile()[1];
REAL ttot =0;
INT k=0;
INT f = 1;
char * c = new char [tiff.SzTile()[0] * tiff.SzTile()[1]*f*3];
for (INT x=0; x< NbTx -1; x++)
for (INT y=0; y< NbTy ; y++)
{
fp.seek_begin(tiff.offset_tile(x,y,0));
ElTimer tim;
cout << tiff.offset_tile(x,y,0) << " " << tiff.byte_count_tile(x,y,0) << "\n";
fp.read(c,1,tiff.byte_count_tile(x,y,0)/1);
REAL t = tim.sval();
ttot += t;
k++;
cout << "Read Time moy " << ttot/k << " " << x << " " << y << "\n";
}
REAL m1 = ttot/k ;
ttot = 0;
k=0;
{
for (INT x=0; x< NbTx ; x++)
for (INT y=0; y< NbTy ; y++)
{
fp.seek_begin(0);
ElTimer tim;
fp.seek_begin(tiff.offset_tile(x,y,0));
REAL t = tim.sval();
ttot += t;
k++;
cout << "Time moy " << ttot/k << " " << x << " " << y << "\n";
}
}
REAL m2 = ttot/k ;
cout << "Moy read = " << m1 << " Moy Seek = " << m2 << "\n";
}
void Optim_L1FormLin::One_bench_craig()
{
/*
ElMatrix<REAL> Smpd = MpdSolve();
ElMatrix<REAL> SBar = BarrodaleSolve();
for (INT v=0; v<_NbVar ; v++)
{
cout << SBar(0,v) << " "
<< Smpd(0,v) << " "
<< (SBar(0,v)-Smpd(0,v)) *1e10 <<"\n";
}
*/
ElTimer t;
BarrodaleSolve();
cout << "Time = " << t.sval() << " " << _NbVar << " " << _NbForm << "\n";
}