double expert(const double *array, int len, double tolerance,
              double *scratch)
{
  double val;
  val = 0;
  /* assume length of scratch is >= len */
  /* do operation with input array, using scratch */

  return val;
}