%************************************************
% Define elastic stiffness tensor (GPa)
%************************************************
%
% Reference Elastic constants
% Alexandrov et al(1974) Crystal No.1 X||a* Y||b Z||c
% Hornblende No.1
%
% Reference Crystal Structure
% Phillips M W, Draheim J E, Popp R K, Clowe C A, Pinkerton A A
% American Mineralogist 74 (1989) 764-773.
%
% Define density (g/cm3)
 rho= 3.1240;
%
% Define Cartesian Tensor crystal symmetry & frame
cs_Tensor = crystalSymmetry('2/m',[ 9.803 18.046 5.313],...
 [ 90, 105.05, 90 ]*degree,'x||a*','z||c',...
 'mineral','Hornblende No.1 1974');
%
% elastic Cij stiffness tensor (GPa) as matrix M
 M =.... 
 [[  115.80   49.30   63.00    0.00   -5.50    0.00];...
 [    49.30  159.00   65.30    0.00  -18.70    0.00];...
 [    63.00   65.30  191.30    0.00   -8.70    0.00];...
 [     0.00    0.00    0.00   58.80    0.00   -6.60];...
 [    -5.50  -18.70   -8.70    0.00   31.70    0.00];...
 [     0.00    0.00    0.00   -6.60    0.00   38.50]];
%
% M as stiffness tensor C with MTEX tensor command
C = tensor(M,cs_Tensor);
%
