在cc中为有限差分法创建五对角稀疏矩阵的最佳方法是什么？

概述在MATLAB中,使用如下命令创建五对角稀疏矩阵非常方便： I = eye(m); % create identity matrixe = ones(m,1); % create an array of all 1'sT = spdiags([e -4*e e],[-1 0 1],m,m); S = spdiags([e e],[-1 1],m,m);A = (kron(I,T) + k 在MATLAB中,使用如下命令创建五对角稀疏矩阵非常方便：

I = eye(m);  % create IDentity matrixe = ones(m,1); % create an array of all 1'sT = spdiags([e -4*e e],[-1 0 1],m,m); S = spdiags([e e],[-1 1],m);A = (kron(I,T) + kron(S,I))/hˆ2;

我想知道是否有任何巧妙的技巧在c / c中做同样的事情.

解决方法 C中没有稀疏的Matrix类型.但是网络上有很多开源代数库(或者你可以编写自己的代码库).

Boost uBLAS支持稀疏矩阵,如果你只想“试验”有限差分,它可能是最好的选择.
如果您需要更高级的求解器,您应该看一下GSL,或者考虑一下LAPACK的C版本.

至于你原来的问题,据我所知,这些库都没有实现kron函数,因为它只是一个“方便”的例程.

总结

以上是内存溢出为你收集整理的在c / c中为有限差分法创建五对角稀疏矩阵的最佳方法是什么？全部内容，希望文章能够帮你解决在c / c中为有限差分法创建五对角稀疏矩阵的最佳方法是什么？所遇到的程序开发问题。

如果觉得内存溢出网站内容还不错，欢迎将内存溢出网站推荐给程序员好友。