Bug fix in the IVP functions - sorry...

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

6_rendering/gtIVPcubCmap.m

@@ -15,10 +15,10 @@ map = zeros(length(r_vectors+1),3);
 
 for i=1:length(r_vectors)
     
-    g=Rod2g(r_vectors(i, :));
+    crys2sam=Rod2g(r_vectors(i, :));
     
-    % z axis in cryst axes
-    Zc = g*Zlab';
+    % z axis in cryst axes \ to go sample->crystal
+    Zc = crys2sam\Zlab';
     
     % all equivilents by sym
     Zc_sym=[];
@@ -44,8 +44,8 @@ for i=1:length(r_vectors)
     
     % these shouldn't come out negative!   
     red=(45-chi)/45;
-    green=phi/45;
-    blue=(45-phi)/45;
+    green=(chi/45)*phi/45;
+    blue=(chi/45)*(45-phi)/45;
  
 % saturate colours    
 rgb=[red green blue];

6_rendering/gtIVPhexCmap.m

@@ -29,13 +29,13 @@ map = zeros(length(r_vectors+1),3);
 
 for i=1:length(r_vectors)
     
-    g=Rod2g(r_vectors(i, :));
+    crys2sam=Rod2g(r_vectors(i, :));
     
     %     disp([num2str(bodge) ' degrees bodge here!'])
     %     g=[cosd(bodge) sind(bodge) 0; -sind(bodge) cosd(bodge) 0; 0 0 1]*g;
     %
-    % RD in cryst axes
-    LDc = g*LD';
+    % RD in cryst axes - \ to do inverse, sample->crystal
+    LDc = crys2sam\LD';
     % all equivilents by sym
     LDc_sym=[];
     for jj=1:12