有没有能够解决以下特征值的fortran求解器?

时间:2017-09-22 01:19:15

标签: fortran eigenvalue

A 是N×N矩阵。 I 是(N-2)×(N-2)的单位矩阵。 B是另一个N乘N矩阵,定义为

auto& i00 = array[0][0];

x 是一个包含N个元素的数组。如何解决以下形式的特征值 A x = c B x ,其中c是特征值, 通过使用特征值求解器?

1 个答案:

答案 0 :(得分:1)

您可以查看Lapack库,它提供了广义矩阵特征值问题的解决方案。根据您的数据类型和矩阵类型,您将需要使用不同的子例程。 在这方面看看here。对于用作矩阵类型的命名法,还要看here

最后,前段时间我写了this module来给出一些Lapack功能的实现示例,包括特征值问题。您可以找到的那个是单值通用矩阵(sgeev)。