关于pardiso处理复数稀疏矩阵的几点疑惑

hello_lxl 2019-04-24 03:48:25

恳请各位大佬懂得能给予点帮助，谢谢~ 首先我处理的复数矩阵是三万阶左右的复数方阵，我在调用pardiso时已经把前面复数矩阵中非零元素转化成CSR存储格式，设置的mtype=13，感觉各个方面都已经准备好了，但最后在运行的时候总是说访问冲突？我也试着把iparm（60）设置成2，但是没作用，弄了一天多还是没结果，所以来此请各位大佬看看错在哪里？我其他iparm设置的是默认

...全文

445 2 打赏收藏转发到动态举报

写回复

2 条回复

切换为时间正序

请发表友善的回复…

发表回复

qq_25424105 2021-05-11

打赏
举报

回复

非对称的可以试试双共轭梯度或者稳定的双共轭梯度 gmres也是可以的

hello_lxl 2019-04-25

打赏
举报

回复

或者说对于处理复数非对称的稀疏矩阵应该用什么方法？我之前用的是共轭梯度法，但是我想知道有没有什么其他方法？我尝试过Zgetrf+Zgetri，其他预处理的共轭梯度，但是不是不符合要求就是水平不够，改不出程序，恳请大家帮帮忙。。。

1、稀疏矩阵定义及存储 稀疏矩阵即含0较多的矩阵。 MATLAB为了节省内存空间，对与实数稀疏矩阵采用3个内部数组来保存非零元素及其对应的位置。对与复数稀疏矩阵，采用第4个内部数组保存非0复数的虚部。说明:M_full和M_sparse其实是同一个矩阵，但是内存却相差了一万倍，可怕。 2、稀疏矩阵的创建 MATLAB不会自己创建稀疏矩阵。一般由用户根据商来确定是否要采用稀疏矩阵。商也叫密度...

1、稀疏矩阵的储存方式对于稀疏矩阵，MATLAB仅储存矩阵所有非零元素的值及其位置（行号和列号）。2、稀疏矩阵的生成1）利用sparse函数从满矩阵转换得到稀疏矩阵函数名称表示意义sparse(A)由非零元素和下标建立稀疏矩阵A。如果A已是一个稀疏矩阵，则返回A本身。sparse(u,v,a)生成大小为max(u) x max(v)的稀疏矩阵。其中u和v是整数向量，a是实数或者复数向量sparse...

Copyright 2008说明：这一段时间用Matlab做了LDPC码的性能仿真，过程中涉及了大量的矩阵运算，本文记录了Matlab中矩阵的相关知识，特别的说明了稀疏矩阵和有限域中的矩阵。Matlab的运算是在矩阵意义下进行的，这里所提到的是狭义上的矩阵，即通常意义上的矩阵。目录内容第一部分：矩阵基本知识(只作基本介绍，详细说明请参考Matlab帮助文档)矩阵是进行数据处理和运算的基本元素。在M...

本文对稀疏矩阵SparseMatrix的主要操作进行了总结。首先，建议先阅读《Eigen教程2 - 稀疏矩阵操作》。关于稀疏矩阵，最重要的一点是：稀疏矩阵的存储方式，是按列优先储存，还是按行优先存储。绝大多数的稀疏矩阵的算术操作都会断言（判断）操作数的存储方式是否相同。 稀疏矩阵初始化构造函数默认列优先。 SparseMatrix<double> sm1(1000,1000); //创建一个

关于稀疏矩阵CSR存储方式的一点理解

69,371

社区成员

243,082

社区内容

发帖

与我相关

我的任务

社区管理员

加入社区

近7日
近30日
至今

加载中

查看更多榜单

社区公告

暂无公告

试试用AI创作助手写篇文章吧

+ 用AI写文章