RF: use faster isfield
OPT: use faster generation of spaces
OPT: do not check dataset when stacking or slicing
NF: provide values in output