The order of operator may not be consistent between eigSolver and flda.

[nennigb]

In eigSolver.py, the matrices are passed to eigenvalue solvers depending on their position in the K list. This assume that operators are well ordered in K, dK and in flda. A more robust approach could be to use flda to sort and to pass the matrices to solvers.

The modification will concern all solve methods of EigSolver concrete class. For instance

        # FIXME need to modify matrix order based on flda
        if self.pb_type=='std':
            self.Lda,Vec = sp.linalg.eig(self.K[0],b=None) 
        elif self.pb_type=='gen':            
            self.Lda,Vec = sp.linalg.eig(self.K[0],b=-self.K[1])
        elif self.pb_type=='PEP':
            self.Lda,Vec = self._pep(self.K)
        else:
            raise NotImplementedError('The pb_type {} is not yet implemented'.format(self.pb_type))