| 查看: 1833 | 回复: 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命令到底是怎么用的呢？纠结啊！！ |
回复此楼» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有17人回复
-
不自信的我
已经有5人回复
-
磺酰氟产物，毕不了业了！
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
3楼2012-11-09 19:21:31
5楼2012-11-10 08:31:44
7楼2012-11-10 20:52:35
9楼2015-06-15 10:41:36