Minor mods to Fsim, Nicola fixed twin reconstruction in AssembleVol

git-svn-id: https://svn.code.sf.net/p/dct/code/trunk@767 4c865b51-4357-4376-afb4-474e03ccb993

5_reconstruction/GtAssembleVol3D.m

@@ -516,9 +516,16 @@ classdef GtAssembleVol3D < handle
             fprintf('Phase %02d: ', phaseID);
 
             % want to work in reverse size order
-            for ii = 1:length(list)
-                vol = obj.localPar.cache.get('grain', {phaseID, list(ii)}, 'seg');
-                grainsize(ii) = length(find(vol(:)));
+            grainsize = zeros(grainsNum, 1);
+            for ii = 1:grainsNum
+                try
+                    vol = obj.localPar.cache.get('grain', {phaseID, list(ii)}, 'seg');
+                    grainsize(ii) = length(find(vol(:)));
+                catch mexc
+                    % We ignore the error because it will come up again
+                    % later
+                    grainsize(ii) = 0;
+                end
             end
             %sort - to do largest grain first
             [~, index] = sort(grainsize, 'descend');

5_reconstruction/gtFsimCheckMissingGrains.m

@@ -26,6 +26,10 @@ load(fullfile( '4_grains',['phase_' sprintf('%02d',phase_num)],'index.mat' ),'gr
 % total number of grains
 nof_grains = length(grain);
 
+% predefine is_uptodate and isgrain
+is_uptodate=zeros(1,length(grain));
+is_grain=zeros(1, length(grain));
+
 list_grains=dir(fullfile( '4_grains',['phase_' sprintf('%02d',phase_num)],'grain_*.mat'));
 for ii=1:length(list_grains)
     if datenum(list_grains(ii).date,'dd-mmm-yyyy')>=datemax