Fix bugs in relatedOperations for NW,GT and GTprime.

python/damask/orientation.py

@@ -875,7 +875,7 @@ class Lattice:
       [[ 17, 12,  5],[ 17,  7, 17]],
       [[  5, 17, 12],[ 17, 17,  7]],
       [[ 12, -5,-17],[  7,-17,-17]],
-      [[-17,-12,  5],[-17,  7, 17]]],dtype='float')}
+      [[-17,-12,  5],[-17,-7, 17]]],dtype='float')}
  
   # Greninger--Troiano' orientation relationship for fcc <-> bcc transformation
   # from Y. He et al., Journal of Applied Crystallography 39:72-81, 2006 
@@ -901,7 +901,7 @@ class Lattice:
       [[-17,-17,  7],[-17, -5, 12]],
       [[  7,-17,-17],[ 12,-17, -5]],
       [[ 17, -7,-17],[ 5, -12,-17]],
-      [[ 17,-17,  7],[ 17, -5,-12]],
+      [[ 17,-17, -7],[ 17, -5,-12]],
       [[ -7, 17,-17],[-12, 17, -5]],
       [[-17,  7,-17],[ -5, 12,-17]],
       [[-17, 17, -7],[-17,  5,-12]]],dtype='float'),
@@ -957,7 +957,7 @@ class Lattice:
       [[  2,  1, -1],[  0, -1,  1]],
       [[ -1, -2, -1],[  0, -1,  1]],
       [[ -1,  1,  2],[  0, -1,  1]],
-      [[ -1,  2,  1],[  0, -1,  1]],
+      [[  2, -1,  1],[  0, -1,  1]], #It is wrong in the paper, but matrix is correct
       [[ -1,  2,  1],[  0, -1,  1]],
       [[ -1, -1, -2],[  0, -1,  1]]],dtype='float')}
                

python/tests/reference/Rotation/PoleFigures_OR.m

@@ -0,0 +1,99 @@
+%% Import Script for EBSD Data
+%
+% Use MTEX
+clear ; clear all
+
+%% Specify Crystal and Specimen Symmetries
+
+% crystal symmetry
+CS_bcc = {... 
+  crystalSymmetry('m-3m', [2.8665 2.8665 2.8665], 'mineral', 'Iron-alpha', 'color', 'light blue'),...
+  crystalSymmetry('m-3m', [1 1 1], 'color', 'light blue')};
+
+CS_fcc = {... 
+  crystalSymmetry('m-3m', [3.662 3.662 3.662], 'mineral', 'Iron', 'color', 'light blue'),...
+  crystalSymmetry('m-3m', [1 1 1], 'color', 'light blue')};
+
+% plotting convention
+setMTEXpref('xAxisDirection','north');
+setMTEXpref('zAxisDirection','outOfPlane');
+
+%% path to files
+pname = 'L:\f.gallardo\DAMASK\python\tests\reference\Rotation'; % has to be changed
+
+% which files to be imported
+fname1 = [pname '\bcc_Bain.txt']; fname2 = [pname '\bcc_GT.txt']; fname3 = [pname '\bcc_GT_prime.txt'];
+fname4 = [pname '\bcc_KS.txt']; fname5 = [pname '\bcc_NW.txt']; fname6 = [pname '\bcc_Pitsch.txt'];
+fname7 = [pname '\fcc_Bain.txt']; fname8 = [pname '\fcc_GT.txt']; fname9 = [pname '\fcc_GT_prime.txt'];
+fname10 = [pname '\fcc_KS.txt']; fname11 = [pname '\fcc_NW.txt']; fname12 = [pname '\fcc_Pitsch.txt'];
+
+
+%% Import the Data
+
+% create an EBSD variable containing the data
+ebsd1 = loadEBSD(fname1,CS_bcc,'interface','generic',...
+  'ColumnNames', { 'phi1' 'Phi' 'phi2' 'x' 'y'}, 'Bunge', 'Passive Rotation');
+ebsd2 = loadEBSD(fname2,CS_bcc,'interface','generic',...
+  'ColumnNames', { 'phi1' 'Phi' 'phi2' 'x' 'y'}, 'Bunge', 'Passive Rotation');
+ebsd3 = loadEBSD(fname3,CS_bcc,'interface','generic',...
+  'ColumnNames', { 'phi1' 'Phi' 'phi2' 'x' 'y'}, 'Bunge', 'Passive Rotation');
+ebsd4 = loadEBSD(fname4,CS_bcc,'interface','generic',...
+  'ColumnNames', { 'phi1' 'Phi' 'phi2' 'x' 'y'}, 'Bunge', 'Passive Rotation');
+ebsd5 = loadEBSD(fname5,CS_bcc,'interface','generic',...
+  'ColumnNames', { 'phi1' 'Phi' 'phi2' 'x' 'y'}, 'Bunge', 'Passive Rotation');
+ebsd6 = loadEBSD(fname6,CS_bcc,'interface','generic',...
+  'ColumnNames', { 'phi1' 'Phi' 'phi2' 'x' 'y'}, 'Bunge', 'Passive Rotation');
+
+ebsd7 = loadEBSD(fname7,CS_fcc,'interface','generic',...
+  'ColumnNames', { 'phi1' 'Phi' 'phi2' 'x' 'y'}, 'Bunge', 'Passive Rotation');
+ebsd8 = loadEBSD(fname8,CS_fcc,'interface','generic',...
+  'ColumnNames', { 'phi1' 'Phi' 'phi2' 'x' 'y'}, 'Bunge', 'Passive Rotation');
+ebsd9 = loadEBSD(fname9,CS_fcc,'interface','generic',...
+  'ColumnNames', { 'phi1' 'Phi' 'phi2' 'x' 'y'}, 'Bunge', 'Passive Rotation');
+ebsd10 = loadEBSD(fname10,CS_fcc,'interface','generic',...
+  'ColumnNames', { 'phi1' 'Phi' 'phi2' 'x' 'y'}, 'Bunge', 'Passive Rotation');
+ebsd11 = loadEBSD(fname11,CS_fcc,'interface','generic',...
+  'ColumnNames', { 'phi1' 'Phi' 'phi2' 'x' 'y'}, 'Bunge', 'Passive Rotation');
+ebsd12 = loadEBSD(fname12,CS_fcc,'interface','generic',...
+  'ColumnNames', { 'phi1' 'Phi' 'phi2' 'x' 'y'}, 'Bunge', 'Passive Rotation');
+
+%% Plot Data 1stpart_bcc
+h1 = [Miller(1,0,0,ebsd1.CS),Miller(1,1,0,ebsd1.CS),Miller(1,1,1,ebsd1.CS)]; % 3 pole figures
+plotPDF(ebsd1.orientations,h1,'MarkerSize',5,'MarkerColor','r','MarkerEdgeColor','r','DisplayName','BCC-Bain')
+hold on
+plotPDF(ebsd2.orientations,h1,'MarkerSize',5,'MarkerColor','w','MarkerEdgeColor','b','DisplayName','BCC-GT')
+plotPDF(ebsd3.orientations,h1,'MarkerSize',5,'MarkerColor','w','MarkerEdgeColor','g','DisplayName','BCC-GT_Prime')
+legend('show','location','southoutside') 
+cd 'L:\f.gallardo\DAMASK\python\tests\reference\Rotation'; % has to be changed 
+orient('landscape')
+print('-bestfit','1_BCC.pdf','-dpdf')
+
+%% Plot Data 2nd part_bcc
+close
+plotPDF(ebsd4.orientations,h1,'MarkerSize',5,'MarkerColor','r','MarkerEdgeColor','w','DisplayName','BCC-KS')
+hold on
+plotPDF(ebsd5.orientations,h1,'MarkerSize',5,'MarkerColor','w','MarkerEdgeColor','m','DisplayName','BCC-NW')
+plotPDF(ebsd6.orientations,h1,'MarkerSize',5,'MarkerColor','y','MarkerEdgeColor','w','DisplayName','BCC-Pitsch')
+legend('show','location','southoutside')
+print('-bestfit','2_BCC.pdf','-dpdf')
+
+%% Plot Data 1stpart_fcc
+close
+h2 = [Miller(1,0,0,ebsd7.CS),Miller(1,1,0,ebsd7.CS),Miller(1,1,1,ebsd7.CS)]; % 3 pole figures
+plotPDF(ebsd7.orientations,h2,'MarkerSize',5,'MarkerColor','r','MarkerEdgeColor','r','DisplayName','FCC-Bain')
+hold on
+plotPDF(ebsd8.orientations,h2,'MarkerSize',5,'MarkerColor','w','MarkerEdgeColor','b','DisplayName','FCC-GT')
+plotPDF(ebsd9.orientations,h2,'MarkerSize',5,'MarkerColor','w','MarkerEdgeColor','g','DisplayName','FCC-GT_Prime')
+legend('show','location','southoutside')
+print('-bestfit','1_FCC.pdf','-dpdf')
+
+%% Plot Data 2nd part_bcc
+close
+plotPDF(ebsd10.orientations,h2,'MarkerSize',5,'MarkerColor','r','MarkerEdgeColor','w','DisplayName','FCC-KS')
+hold on
+plotPDF(ebsd11.orientations,h2,'MarkerSize',5,'MarkerColor','w','MarkerEdgeColor','m','DisplayName','FCC-NW')
+plotPDF(ebsd12.orientations,h2,'MarkerSize',5,'MarkerColor','y','MarkerEdgeColor','w','DisplayName','FCC-Pitsch')
+legend('show','location','southoutside')
+print('-bestfit','2_FCC.pdf','-dpdf')
+close
+

python/tests/reference/Rotation/bcc_GT.txt

@@ -23,4 +23,4 @@
 352.1156357053931 43.82007387041961 14.074783631236542 1 21
 77.82610341510008 43.397849654402556 273.4002228089796 1 22
 193.60440567265297 9.976439066337806 123.24637065555939 1 23
-153.65751914298576 65.6559553854118 185.90444335627936 1 24
+172.11563570539317 43.82007387041961 194.07478363123653 1 24

python/tests/reference/Rotation/bcc_GT_prime.txt

@@ -20,7 +20,7 @@
 183.40022280897963 43.397849654402556 167.8261034151001 1 18
 255.92521636876344 43.82007387041961 97.88436429460687 1 19
 33.24637065555936 9.976439066337804 283.60440567265294 1 20
-26.291675350407385 65.60048732963618 354.34378938496315 1 21
+356.59977719102034 43.39784965440254 12.173896584899929 1 21
 75.92521636876346 43.82007387041961 277.8843642946069 1 22
 213.24637065555936 9.976439066337804 103.604405672653 1 23
 176.59977719102034 43.397849654402556 192.17389658489986 1 24

python/tests/reference/Rotation/bcc_NW.txt

@@ -9,6 +9,6 @@
 134.58444405678858 83.13253115922213 96.91733794010702 1 7
 225.41555594321142 83.13253115922213 173.082662059893 1 8
 0.0 9.735610317245317 135.0 1 9
-260.40196970123213 45.81931182053556 283.6387072794765 1 10
+99.59803029876785 45.81931182053557 166.36129272052355 1 10
 260.40196970123213 45.81931182053556 283.6387072794765 1 11
 180.0 99.73561031724535 225.0 1 12

python/tests/reference/Rotation/fcc_GT.txt

@@ -23,4 +23,4 @@
 165.92521636876344 43.82007387041961 187.88436429460683 1 21
 266.59977719102034 43.39784965440254 102.17389658489992 1 22
 56.75362934444064 9.976439066337804 346.395594327347 1 23
-354.0955566437206 65.6559553854118 26.342480857014277 1 24
+345.9252163687635 43.82007387041961 7.884364294606862 1 24

python/tests/reference/Rotation/fcc_GT_prime.txt

@@ -20,7 +20,7 @@
 12.173896584899929 43.39784965440254 356.59977719102034 1 18
 82.11563570539313 43.82007387041961 284.0747836312365 1 19
 256.395594327347 9.976439066337804 146.75362934444064 1 20
-185.65621061503683 65.60048732963617 153.70832464959264 1 21
+167.8261034151001 43.397849654402556 183.40022280897963 1 21
 262.1156357053931 43.82007387041961 104.07478363123654 1 22
 76.39559432734703 9.976439066337806 326.75362934444064 1 23
 347.8261034151001 43.39784965440255 3.400222808979685 1 24

python/tests/reference/Rotation/fcc_NW.txt

@@ -9,6 +9,6 @@
 83.082662059893 83.13253115922213 45.415555943211444 1 7
 6.917337940106983 83.13253115922211 314.5844440567886 1 8
 45.0 9.73561031724532 180.0 1 9
-256.36129272052347 45.81931182053556 279.59803029876775 1 10
+13.638707279476469  45.81931182053557 80.40196970123216 1 10
 256.36129272052347 45.81931182053556 279.59803029876775 1 11
 315.0 99.73561031724536 0.0 1 12