void print_3d_cell(FILE *wr, interval_t cell[])
  { fprintf(wr, 
      "[%+8.5f _ %+8.5f]×[%+8.5f _ %+8.5f]×[%+8.5f _ %+8.5f]",
      LO(cell[0]), HI(cell[0]), 
      LO(cell[1]), HI(cell[1]), 
      LO(cell[2]), HI(cell[2])
    );
  }