%************************************************
% Define elastic stiffness tensor (GPa)
%************************************************
%
% Reference Elastic constants
% Alexsandrov et al(1974) An57-60 XIIa YII(Z x X) ZIIc*
% Elastic constants measured as monoclinic called "pseudo-monoclinic"
% 
%
% Reference Crystal Structure
% LABRADORITE  An 66 Wenk and Kroll(1985)
%
% Define density (g/cm3)
 rho= 2.6800;
% N.B. Crystal structure is triclinic
% Define Cartesian Tensor crystal symmetry & frame
cs_Tensor = crystalSymmetry('-1',[  8.1747 12.8708 14.2030],...
 [  93.4610 116.0860  90.5140]*degree,'x||a','z||c*',...
 'mineral','An57-60 1974');
%
% elastic Cij stiffness tensor (GPa) as matrix M
 M =.... 
 [[   99.40   62.80   48.70    0.00   -0.25    0.00];...
 [    62.80  158.40   26.70    0.00  -10.70    0.00];...
 [    48.70   26.70  150.30    0.00  -12.40    0.00];...
 [     0.00    0.00    0.00   21.70    0.00   -5.40];...
 [    -0.25  -10.70  -12.40    0.00   34.50    0.00];...
 [     0.00    0.00    0.00   -5.40    0.00   37.10]];
%
% M as stiffness tensor C with MTEX tensor command
C = tensor(M,cs_Tensor);
%
