有没有可以解决以下特征值的Fortran解算器？

Question

A是N乘N矩阵。 I是（N-2）乘以（N-2）的单位矩阵。 B是另一个N乘N矩阵，其被定义为

B=[I 0 0; 
     0 0 0; 
     0 0 0] 

。 x是一个有N个元素的数组。如何解决以下形式的本征值 甲X = C 乙X，其中c是本征值， 通过使用特征值求解器？

Answer 1

您可以查看Lapack库，它提供了广义矩阵的特征值问题的解决方案。根据您的数据类型和矩阵类型，您将需要使用不同的子程序。 在这方面看看here。也可以看看here作为用作矩阵类型的术语。

最后，我之前写了this module来给出一些Lapack功能的实现示例，其中包括特征值问题。在那里你可以找到一个单值的广义矩阵（sgeev）。