lapack矩阵存储

shanweikun 2014-04-08 11:18:17
我要求解一个特征值问题A*x=aB*x, a是特征值,在我形成两个矩阵A,B的时候我是直接存储的两个矩阵的非零元 按照i,j,k格式,即是说,数组i 存的是行号,j是列号,k是对应位置的矩阵元素值, 因为矩阵很大很稀疏 所以为了内存考虑 我就这么存储的 但是当我调用lapack函数求解特征值的时候 比如用函数dsygv_或者别的,因为lapack支持的矩阵存储就是带状形式 压缩形式 三角形式,所以我就想问问 怎么利用我的这种i,j,k格式 转化为函数可以计算的格式呢?lapack没有这样的转换函数吧? 我想的是能不能利用我最初的i,j,k对矩阵进行重新组合 使得带宽尽可能小 这样进而利用带状的形式带入函数进行运算~~这只是我想达到的效果 还请各位大神有什么思路和建议 请多多指教啊~~感谢感谢
...全文
145 回复 打赏 收藏 转发到动态 举报
写回复
用AI写文章
回复
切换为时间正序
请发表友善的回复…
发表回复

65,209

社区成员

发帖
与我相关
我的任务
社区描述
C++ 语言相关问题讨论,技术干货分享,前沿动态等
c++ 技术论坛(原bbs)
社区管理员
  • C++ 语言社区
  • encoderlee
  • paschen
加入社区
  • 近7日
  • 近30日
  • 至今
社区公告
  1. 请不要发布与C++技术无关的贴子
  2. 请不要发布与技术无关的招聘、广告的帖子
  3. 请尽可能的描述清楚你的问题,如果涉及到代码请尽可能的格式化一下

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