add the input type of rod vectors

Signed-off-by: Zheheng Liu <zheheng.liu@esrf.fr>

zUtil_Twins/gtFindAnnealingTwinsAmongNeighboringGrains_v2.m

@@ -77,11 +77,21 @@ function [twin_clusters, neighboring_twin_clusters, twin_cands, neighbors] = gtF
         'merge_dist_deg', 1, ...
         'neighbors', [], ...
         'neighbor_strategy', 'dil_vol', ...
+        'input', 'gr_id', ...
         'save', false, ...
         'pad', 1);
     conf = parse_pv_pairs(conf, varargin);
-    %% Load the quaternions of grains
-    Quaternions = gtMathsRod2Quat(sample.phases{phaseID}.R_vector.');
+
+    switch conf.input
+        case 'gr_id'
+            %% Load the quaternions of grains
+            Quaternions = gtMathsRod2Quat(sample.phases{phaseID}.R_vector.');
+        case 'rod'
+            Quaternions = gtMathsRod2Quat(grain_ids);
+            grain_ids = 1:size(grain_ids, 2);
+            conf.save = false;
+            conf.neighbors = 1:10;
+    end
 
     %% Find neighboring grains
     if (size(conf.neighbors, 1) == tot_grains && size(conf.neighbors, 2) == tot_grains)