我对求解稀疏线性系统Ax = b感兴趣,其中A是复数的方阵,b是复数的 vector 。

如果可能,我希望将此类库作为模板(为了易于安装和使用)
本着Eigen的精神

我检查了Eigen,但我认为它似乎不支持使用复杂的稀疏矩阵求解线性方程组(尽管可以在复杂的矩阵上创建和执行基本运算)。

有人向我建议的另一个技巧是,可以解决这一问题,方法是使用(A1 + iA2)(x1 + ix2)=(b1 + ib2)来求解两倍尺寸的扩展系统
但是我更喜欢一些简单的黑盒子来完成工作。

有什么建议么?

最佳答案

将其转移到二维的实值系统可能是最直接的方法。您可以编写一个适配器来打包转换逻辑。也可以尝试以下一种:http://trilinos.sandia.gov/packages/docs/r4.0/packages/komplex/doc/html/

关于c++ - 解决复杂线性系统Ax = b的C++库,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/20672637/

10-14 19:49
查看更多