Fortran求特征根问题

pasuka · 发表于 2016-9-29 08:03:04

首先，这个不是程序的问题：
1、LAPACK难道不能求解广义特征值问题吗？
2、有限元质量阵就一定是SPD吗？

pasuka · 发表于 2016-9-29 10:38:52

zzgjtx 发表于 2016-9-29 10:20
1. 在我的问题里，质量矩阵是正定对称矩阵，我验证过了。
2. 不管ivf还是NAG函数库，都是调用的LAPACK。 ...

一步到位改成DSYGVD不是挺好嘛，NAG网站都有现成例子
http://www.nag.com/lapack-ex/node104.html

pasuka · 发表于 2016-9-29 21:06:08

zzgjtx 发表于 2016-9-29 11:39
谢谢，我又仔细去看了这个例子。你说的是对的，最近被特征根问题搞的头大了。

问题好多，发现我求解的A* ...

计算方法和矩阵论补齐的话，再配上戈卢布的矩阵计算，就没有搞不明白的LAPACK程序调用问题

[线性代数] Fortran求特征根问题