Use GaussLegendre quadrature for computation of Taylor factor

GaussLegendre quadrature needs less evaluation points but is slightly more expansive in NSOFT.
In the case of the taylor model the evaluation process is very expansive.

SO3Fun/@SO3FunHarmonic/quadrature.m

@@ -60,8 +60,12 @@
 % Curtis quadrature grid in fundamental region. 
 % Therefore adjust the bandwidth to crystal and specimen symmetry.
 bw = adjustBandwidth(N,SRight,SLeft);
-SO3G = quadratureSO3Grid(bw,'ClenshawCurtis',SRight,SLeft,'ABG');
-
+if check_option(varargin,'GaussLegendre')
+  SO3G = quadratureSO3Grid(bw,'GaussLegendre',SRight,SLeft,'ABG');
+else
+  SO3G = quadratureSO3Grid(bw,'ClenshawCurtis',SRight,SLeft,'ABG');
+end
+
 % Only evaluate unique orientations
 values = f.eval(SO3G);
 

TensorAnalysis/@strainTensor/calcTaylor.m

@@ -46,7 +46,8 @@
   bw = get_option(varargin,'bandwidth',32);
   numOut = nargout;
   F = SO3FunHandle(@(rot) calcTaylorFun(rot,eps,sS,numOut,varargin{:}),sS.CS,eps.CS);
-  SO3F = SO3FunHarmonic(F,'bandwidth',bw);
+  % Use Gauss-Legendre quadrature, since the evaluation process is very expansive
+  SO3F = SO3FunHarmonic(F,'bandwidth',bw,'GaussLegendre');
   M = SO3F(1);
   if nargout>1
     b = [];