%************************************************
% Define elastic stiffness tensor (GPa)
%************************************************
%
% Elastic constants reference
% Low Albite M.Brown et al.2006 PCM An0 An0 x||*a, y||(Z x X) , z||c
%
% Crystal structure reference
% Low Albite An0 Harlow and Brown (1980) Am.Min.65,986-995.
% s.g. C-1
% Define density (g/cm3)
rho= 2.6223;
%
% Define Cartesian Tensor crystal symmetry & frame
cs_Tensor = crystalSymmetry('-1',[  8.1420 12.7850  7.1590],...
 [  94.1900 116.6100  87.6800]*degree,'x||*a','z||c',...
 'mineral','Low Albite An0 2007');
%
% elastic stiffness tensor (GPa) as matrix M
 M =.... 
 [[   69.90   33.24   31.56    5.28   -2.46   -0.72];...
 [    33.24  183.28    7.53    5.31   -7.60   -4.23];...
 [    31.56    7.53  175.65  -17.48    5.86  -11.29];...
 [     5.28    5.31  -17.48   26.93   -3.94   -6.56];...
 [    -2.46   -7.60    5.86   -3.94   26.91    0.98];...
 [    -0.72   -4.23  -11.29   -6.56    0.98   33.39]];
%
% M as stiffness tensor C with MTEX tensor command
C = tensor(M,cs_Tensor);
%
