| 查看: 2298 | 回复: 8 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
fortran到底怎么弄稀疏矩阵啊？
|
||
 大家好，我对稀疏矩阵很模糊，在网上和论坛上搜了下一直没找到该怎么去构建稀疏矩阵。我这里有个图形，是六角格子，沿着x方向具有平移对称性，也就是说x方向是一样的，y方向我取8个点为一个原包，也就是y方向8个点为一个最小单元，我现在把这8个点的哈密顿矩阵写出来： 0 t 0 0 0 0 0 t t 0 t 0 0 0 0 0 0 t 0 t 0 0 0 0 0 0 t 0 t 0 0 0 0 0 0 t 0 t 0 0 0 0 0 0 t 0 t 0 0 0 0 0 0 t 0 t t 0 0 0 0 0 t 0 t的近邻两点之间相互作用，从哈密顿形式上看，是稀疏矩阵，我有一个问题：怎么用fortran把它写成稀疏矩阵？比如我把上面的哈密顿写出来之后，怎么把它转成稀疏矩阵？ 这个小程序就是上面这个哈密顿矩阵的写法。 program xishu implicit none real(8) :: t=1.0 integer :: i1,ny real(8),allocatable :: h(:,: ) ny=8 allocate(h(ny,ny)) h=0.0 do i1=1,ny-1 h(i1,i1+1)=t h(i1+1,i1)=t enddo h(1,ny)=t h(ny,1)=t end 我想把它转成稀疏矩阵，该怎么办呢？sparse命令到底是怎么用的呢？纠结啊！！ |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- FORTRAN中比双精度（real（kind=8））更高的精度的数据类型是什么？
已经有10人回复
- 想用Fortran语言实现矩阵连乘，并且矩阵里面有一个变量
已经有6人回复
- fortran中如何创建由二维矩阵为元素的1维矩阵
已经有4人回复
- 请问在fortran里面怎么求复数矩阵的行列式！！！
已经有4人回复
- 求助：利用fortran写一个矩阵求逆的程序
已经有14人回复
- fortran 怎么按行读取数据（有小例子）
已经有21人回复
- 编一个fortran程序
已经有17人回复
- 关于6*6刚度系数矩阵和9*9刚度系数矩阵的问题
已经有4人回复
- fortran real(kind=16)
已经有10人回复
- 稀疏矩阵的处理
已经有12人回复
- 怎么将此for循环改为矩阵形式？
已经有5人回复
- fortran输出求助
已经有6人回复
- fortran如何实现张量相乘（并乘）
已经有9人回复
- 用fortran存储稀疏矩阵
已经有4人回复
- fortran中到底该如何比较两个实数是否相等
已经有8人回复
- 【求助】fortran求解矩阵
已经有7人回复
- 【求助】用fortran怎么描述下面这样的情况，我是新手，请大家帮忙
已经有8人回复
- 【求助】如何学习FORTRAN的LAPACK?
已经有4人回复
- 【求助】Fortran 如何读取文件中特定行的内容
已经有16人回复
5楼2012-11-10 08:31:44
virtualzx
木虫 (著名写手)
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2012-11-10 01:05:08
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2012-11-10 01:05:08
2楼2012-11-09 03:03:32
weiwei43
至尊木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 18316.1
- 红花: 6
- 帖子: 737
- 在线: 94.7小时
- 虫号: 250282
- 注册: 2006-05-13
- 专业: 机械动力学
4楼2012-11-10 08:22:12
weiwei43
至尊木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 18316.1
- 红花: 6
- 帖子: 737
- 在线: 94.7小时
- 虫号: 250282
- 注册: 2006-05-13
- 专业: 机械动力学
6楼2012-11-10 08:49:45
40