Small changes

func/OR/getVariants.m

@@ -1,4 +1,4 @@
-function [v, ind] = getVariants(o, ORmat, CS)
+function [v, ind] = getVariants(o, ORmat, CS, varargin)
 % Get all crystallographic variant of transformation products.
 %   Array 'v' contain all variants. Variants groups by symmetry, first
 %   variants for one rotation.
@@ -21,9 +21,18 @@
 % Set misorientation
 ORr = rotation('matrix', ORmat);
 
-% Some variants
-v0 = o*CS*ORr;
-v = orientation(v0,get(o,'CS'),get(o,'SS'));
-ind = repmat(1:length(o),1,size(CS,1))'; % ???
+if check_option(varargin, 'someVariants')
+    i = get_option(varargin, 'someVariants');
+    % Some variants
+    v00 = ORr*o;
+    v0 = inverse(rotation(CS(i)))*v00(:);
+    v = orientation(v0,get(o,'CS'),get(o,'SS'));
+    ind = ones(length(v),1);
+else
+    % Some variants
+    v0 = o*CS*ORr;
+    v = orientation(v0,get(o,'CS'),get(o,'SS'));
+    ind = repmat(1:length(o),1,size(CS,1))'; % ???
+end
 
 end

func/misc/exportData.m

@@ -6,6 +6,9 @@
         str = exportEuler(data, style);
     case 'matrix'
         str = exportMatrix(data, style);
+    case 'ori'
+        eul = Euler(data);
+        str = exportEuler(eul/degree, style);
     otherwise
         error('Unknown type.')
 end

func/recon/findUniqueParent.m

@@ -98,19 +98,14 @@
 [op, opi] = getVariants(ori, ORmat, CS);
 np = length(op);
 
-% Prevent out of memory
-if length(op) > getpref('ebsdam', 'maxProbParents')
-    fprintf('Too much probable parent\n'); return;
-end
-
-
 %% Calculation of probabilities of potential parents
 
 % Angle between potential parents
 %  aaa = repmat(op, 1, np);
 %  bbb = repmat(reshape(op,1,np), np,1);
 %  ma = angle(aaa,bbb);
-ma = angle(op\op);
+
+ma = smartAngle(op);
 
 % Check potential parents
 %   Pp  - probability of parent
@@ -183,6 +178,35 @@
 end
 
 
+function ma = smartAngle(op)
+% Prevent out of memory
+if length(op) > getpref('ebsdam', 'maxProbParents')
+    fprintf('Too much probable parent\n');
+
+    ma = [];
+
+    n0 = getpref('ebsdam', 'maxProbParents');
+    n1 = length(op);
+    n = ceil(n1*n1/(n0*n0));
+    ni = ceil(n1/n);
+    is = 1;
+    ie = ni;
+    for i = 1:n
+        ind = is:ie;
+        mai = angle(op(ind)\op);
+        ma = [ma; mai]; %#ok<AGROW>
+        is = is + ni;
+        ie = ie + ni;
+        if (ie > n1)
+            ie = n1;
+        end
+    end
+else
+    ma = angle(op\op);
+end
+end
+
+
 function [bool, PR] = checkRatio(PS, PRmin, Pmin, vn, VNmin)
 % Calculate the ratio of parent probabilities
 %