use output methods (defined once) instead of coding it each time
error remaining: a voxel being in a subset should not appear in a bigset