From 72e84657a75eb67f64f4885ea3536f867fc8ceee Mon Sep 17 00:00:00 2001
From: Johannes Pfeifer <jpfeifer@gmx.de>
Date: Sat, 20 Aug 2016 20:43:20 +0200
Subject: [PATCH] Condition use of eig to compute left eigenvector on Matlab
 version

Functionality is only supported in Matlab from 2014a/8.3 onwards. Closes #1266
---
 matlab/identification_analysis.m | 9 ++++++++-
 1 file changed, 8 insertions(+), 1 deletion(-)

diff --git a/matlab/identification_analysis.m b/matlab/identification_analysis.m
index 46613c5c5..68329c744 100644
--- a/matlab/identification_analysis.m
+++ b/matlab/identification_analysis.m
@@ -194,7 +194,14 @@ if info(1)==0,
 %             [V,D,W]=eig(cmm);
             sd=sqrt(diag(cmm));
             cc=cmm./(sd*sd');
-            [V,D,W]=eig(cc);
+            if isoctave && octave_ver_less_than('8.3')                 
+                [V,D]=eig(cc);
+                %fix for older Matlab versions that do not support computing left eigenvalues, see http://mathworks.com/help/releases/R2012b/matlab/ref/eig.html
+                [W,junk] = eig(cc.'); 
+                W = conj(W);
+            else
+                [V,D,W]=eig(cc);
+            end
             id=find(diag(D)>1.e-8);
             siTMP=siJ./repmat(sd,[1 nparam]);
             MIM=(siTMP'*V(:,id))*(D(id,id)\(W(:,id)'*siTMP));