#include <TGraph.h>
#include <TRandom3.h>
#include <TStopwatch.h>
#include <algorithm>

int profile_eval()
{
  double x[] = {0,1,2,3,4,5,6,7,8,9};
  double y[] = {4,6,3,7,8,34,2,6,8,3};

  TGraph g(10,x,y);
  const int seed = 12345;
  const int N = 100000000;

  TStopwatch ts1;

  ts1.Start();
  TRandom3 tr1(seed);
  volatile double xval0;
  for(int i = 0; i < N;i++)
  {
    xval0 = g.Eval(tr1.Rndm());
  }
  ts1.Stop();
  ts1.Print();

  TStopwatch ts2;
  ts2.Start();
  TRandom3 tr2(seed);
  for(int i = 0;i<N;i++)
  {
    const double * const xarr = g.GetX();
    const int idx = std::upper_bound(xarr,xarr+g.GetN(),tr2.Rndm())-1-xarr;
    double xval,yval;
    g.GetPoint(idx,xval,yval);
  }
  ts2.Stop();
  ts2.Print();
  return 0;
}