Специализированные методы для комплексных симметричных трехдиагональных обобщенных задач на собственные значения

Я должен решить обобщенные задачи на собственные значения где A и B являются трехдиагональными, B симметрично положительно определенным и вещественным, но A является только комплексно симметричным (не определенным или эрмитовым). Кроме того, мне нужно полное собственное разложение. В настоящее время я просто называю обобщенный eigensolver Лапака, но мне интересно, есть ли лучшие методы для этой конкретной, высоко структурированной проблемы. В частности, лучше иметь свободно доступный код (C ++).$Ax = \lambda Bx$$A$$B$$B$$A$ZGGEV

Метод расширения полюсов и выбранной инверсии ( PEXSI ) может быть ответом. Я не использовал этот метод, но он предлагает процедуру инверсии для сложных симметричных матриц. Он не специфичен для трехдиагональных матриц, но использует разреженность.