大对称正带矩阵

计算科学 线性代数 有限差分 带状矩阵 英特尔-mkl
2021-12-16 12:46:57

我使用Intel MKL的gpbsv命令来解决对称正带系统。但不幸的是,当系统很大时,我收到错误访问违反VisualStudio 中的写入位置。

有人可以给我一个建议,我该如何解决这个系统。也许我可以使用更优化的库?也许为大问题编写自己的求解器会更简单(我在哪里可以阅读有关实现的信息)?

PS:我只有一台Core i3笔记本。

最良好的祝愿。

1个回答

如果这是一个大波段矩阵,我认为它是稀疏的,您可以很好地使用 MKL 中的稀疏求解器(例如 Pardiso 或迭代,例如共轭梯度)

http://software.intel.com/sites/products/documentation/hpc/mkl/mklman/GUID-78889273-7E77-426A-9B5E-23A7C2378D78.htm