static double roundUpToPowerOfTen_double(double d) {
  regt iround(exp10(ceil(log10(d))));
}