我对求解稀疏复线性系统Ax=b感兴趣,其中A是复数的方阵,b是复数的向量。
如果可能的话,我希望对这样的库进行模板化(为了便于安装和使用)本着艾根精神
我检查了 Eigen,但我认为它看起来不支持解决具有复杂稀疏矩阵的线性方程(尽管可以在复杂矩阵上创建和执行基本运算。
有人向我建议的另一个技巧是,通过使用 (A1 + iA2)(x1 + ix2) = (b1 + ib2) 的事实来解决两倍维度的扩展系统来解决这个问题但我更喜欢一些简单的黑匣子来完成工作。
有什么建议吗?
将其转移到维度两倍的实值系统可能是最直接的方法。您可以编写一个适配器来打包转换逻辑。也可以试试这个:http://trilinos.sandia.gov/packages/docs/r4.0/packages/komplex/doc/html/