6D-reconstruction: merged code paths for num_interp and shape-functions reconstructions

Num_interp with shape-functions is not working yet, and multi-ospace-region reconstructions, like twins might be broken.

Signed-off-by: Nicola Vigano <vigano@yoda.esrf.fr>