1template<class T> inline T min(T a, T b) { return a < b ? a : b; }
2double cornerbound(double *P, double (*m)(double, double))
3{
4  double b=m(P[0],P[3]);
5  return m(b,P[12]);
6}
7void bound(double *P, double (*m)(double, double), double b)
8{
9  m(b,cornerbound(P,m));
10}
11void bounds(double fuzz, unsigned maxdepth)
12{
13  double Px[]={};
14  double bx=Px[0];
15  bound(Px,min,bx);
16}
17