วิธีเฉพาะสำหรับปัญหาค่าลักษณะเฉพาะทั่วไปที่ซับซ้อนแบบสมมาตร

ฉันต้องแก้ปัญหาค่าลักษณะเฉพาะทั่วไปโดยที่AและBมีทั้งสามด้าน, Bเป็นสมมาตรบวกแน่นอนและจริง แต่Aเป็นสมมาตรเชิงซ้อนเท่านั้น (ไม่ใช่ชัดเจนหรือ Hermitian) นอกจากนี้ฉันต้องการ eigendecomposition เต็ม ขณะนี้ฉันเพิ่งเรียกอีเจนซอลเวอร์ทั่วไปของLapack แต่ฉันสงสัยว่ามีวิธีที่ดีกว่าสำหรับปัญหาที่มีโครงสร้างสูงนี้โดยเฉพาะหรือไม่ โดยเฉพาะอย่างยิ่งการมีโค้ดที่ใช้ได้อย่างอิสระ (C ++) จะดีที่สุด$Ax = \lambda Bx$$A$$B$$B$$A$ZGGEV

วิธีการขยายเสาและวิธีผกผันที่เลือก ( PEXSI ) อาจเป็นคำตอบ ฉันไม่ได้ใช้วิธีนี้ แต่มีวิธีการผกผันสำหรับเมทริกซ์สมมาตรที่ซับซ้อน มันไม่เฉพาะเจาะจงกับเมทริกซ์แบบไตรภาคี แต่ใช้ประโยชน์จาก sparsity