stopped using inv to get matrix inverse and use \ operator instead. Possible performance gains.

Author: dpmcsuss

classifiers/discriminant_classifiers/discriminant_classifiers.m

@@ -22,7 +22,7 @@
     if LDA
         e0=(preds(:,i)-params.mu0);
         e1=(preds(:,i)-params.mu1);
-        yhat.LDA(i)=e0'*params.InvSig*e0 - params.lnprior0 > e1'*params.InvSig*e1 - params.lnprior1;
+        yhat.LDA(i)=e0'*(params.Sig\e0) - params.lnprior0 > e1'*(params.Sig\e1) - params.lnprior1;
     end
 
     if dLDA
@@ -34,7 +34,7 @@
     if QDA
         e0=(preds(:,i)-params.mu0);
         e1=(preds(:,i)-params.mu1);
-        yhat.QDA(i)=e0'*params.InvSig0*e0 - params.lnprior0 > e1'*params.InvSig1*e1 - params.lnprior1;
+        yhat.QDA(i)=e0'*(params.Sig0\e0) - params.lnprior0 > e1'*(params.Sig1\e1) - params.lnprior1;
     end
 
     if dQDA

classifiers/discriminant_classifiers/get_discriminant_params.m

@@ -13,7 +13,8 @@
 end
 
 if isfield(discrim,'LDA')
-    params.InvSig=inv(Sig);
+    %params.InvSig=inv(Sig);
+    params.Sig=Sig;
 end
 
 if isfield(discrim,'dLDA')
@@ -37,8 +38,10 @@
     Sig0=cov(x(:,inds.y0trn)');
     Sig1=cov(x(:,inds.y1trn)');
     if isfield(discrim,'QDA') 
-        params.InvSig0=inv(Sig0);
-        params.InvSig1=inv(Sig1);
+        params.Sig0=Sig0;
+        params.Sig1=Sig1;
+        %params.InvSig0=inv(Sig0);
+        %params.InvSig1=inv(Sig1);
     end
     if isfield( discrim,'dQDA' )
         params.InvdSig0 = diag(diag(Sig0).^-1);