sbool closeToOne(double d, double tolerance) {
  ret d <= 1+tolerance && d >= 1/(1+tolerance);
}