'Error: Index list out of range' following computation of eigenvalues during continuation and/or bifurcation detection

[Maurizio De Pitta']

I keep on hitting this error, no matter what tolerances I choose. It seems that the bifurcation is detected but the code exits because of a bug (?) in the detection/computation of the unit vectors. Why is this occurring and how could I avoid it?

Test function going crazy:  class Hopf_Bor(BorderMethod, BiAltMethod)

Test function going crazy:  class Hopf_Bor(BorderMethod, BiAltMethod)

---------------------------------------------------------------------------
IndexError                                Traceback (most recent call last)
<ipython-input-77-324faa53a51c> in <module>()
     23 print "Computing EQ1"
     24 # PC1['EQ1'].backward()
---> 25 PC1['EQ1'].forward()
     26 print "done!"
     27

/usr/lib/python2.7/site-packages/PyDSTool-0.90.2-py2.7.egg/PyDSTool/PyCont/Continuation.pyc in forward(self)
   1244         self.CurveInfo = PointInfo()
   1245         if self.sol is None:
-> 1246             self._compute(v0=self.initdirec)
   1247             self.sol = self._curveToPointset()
   1248

/usr/lib/python2.7/site-packages/PyDSTool-0.90.2-py2.7.egg/PyDSTool/PyCont/Continuation.pyc in _compute(self, x0, v0, direc)
   1193                 # update self.loc for Corrector's reference
   1194                 self.loc = loc
-> 1195                 if self.BifPoints != {} and self._checkForBifPoints():
   1196                     stop = True
   1197

/usr/lib/python2.7/site-packages/PyDSTool-0.90.2-py2.7.egg/PyDSTool/PyCont/Continuation.pyc in _checkForBifPoints(self)
    533                 Xval, Vval = bfinfo.locate((curve[loc-1], V[loc-1]),
    534                                            (curve_loc, V_loc), self)
--> 535                 found = bfinfo.process(Xval, Vval, self)
    536 
    537                 if found:

/usr/lib/python2.7/site-packages/PyDSTool-0.90.2-py2.7.egg/PyDSTool/PyCont/BifPoint.pyc in process(self, X, V, C)
    217         W, VR = linalg.eig(A)
    218         W0 = [ind for ind, eig in enumerate(W) if abs(eig) < 5e-5] #5e-5
--> 219         V1 = real(VR[:,W0[0]])
    220 
    221         H = C.CorrFunc.hess(X, C.coords+C.params, C.coords+C.params)

IndexError: list index out of range

[Maurizio De Pitta']

Uhm I actually posted already on this topic a while ago. Sorry, about it sometimes short-term memory becomes an issue...