| 612 | } |
| 613 | |
| 614 | volume empty_volume(ndim dim) { |
| 615 | volume out(dim); |
| 616 | LOOP_OVER_DIRECTIONS(dim, d) { |
| 617 | out.set_direction_max(d, 0.0); |
| 618 | out.set_direction_min(d, 0.0); |
| 619 | } |
| 620 | return out; |
| 621 | } |
| 622 | |
| 623 | volume grid_volume::dV(const ivec &here, double diameter) const { |
| 624 | const double hinva = 0.5 * inva * diameter; |