6D-Reconstruction/shape-functions/twins: fixes

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

zUtil_Deformation/Gt6DReconstructionAlgorithmFactory.m

@@ -73,7 +73,7 @@ classdef Gt6DReconstructionAlgorithmFactory < handle
             for ii_o_reg = 1:num_regions
                 gr = sampler(ii_o_reg).get_reference_grain();
                 fprintf('%d) Grainid %d:\n', ii_o_reg, gr.id)
-                shape_funcs = self.get_shape_functions(sampler);
+                shape_funcs = self.get_shape_functions(sampler(ii_o_reg));
                 algo_params(ii_o_reg) = self.get_algorithm_params( ...
                     sampler(ii_o_reg), num_interp, shape_funcs); %#ok<AGROW>
                 fprintf('\n')

zUtil_Deformation/gtDefShapeFunctionsFwdProj.m

@@ -84,7 +84,9 @@ function shape_funcs = gtDefShapeFunctionsFwdProj(sampler, varargin)
 
     gv.used_ind = find(gv.pow > 0);
 
-    sel_bls = find(sampler.selected);
+    % Adapting to the recent change in GtOrientationSampling, where we only
+    % compute the scattering info for the selected reflections
+    sel_bls = 1:numel(find(sampler.selected));
 
     % Only one lattice allowed for the moment!
     num_ors = numel(sampler.lattice(1).gr);