 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 9 | 1/1 | 返回列表 |
| 查看: 1972 | 回復(fù): 8 | |||
[求助]
fortran到底怎么弄稀疏矩陣啊？
|
|
大家好，我對(duì)稀疏矩陣很模糊，在網(wǎng)上和論壇上搜了下一直沒(méi)找到該怎么去構(gòu)建稀疏矩陣。我這里有個(gè)圖形，是六角格子，沿著x方向具有平移對(duì)稱性，也就是說(shuō)x方向是一樣的，y方向我取8個(gè)點(diǎn)為一個(gè)原包，也就是y方向8個(gè)點(diǎn)為一個(gè)最小單元，我現(xiàn)在把這8個(gè)點(diǎn)的哈密頓矩陣寫出來(lái)： 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的近鄰兩點(diǎn)之間相互作用，從哈密頓形式上看，是稀疏矩陣，我有一個(gè)問(wèn)題：怎么用fortran把它寫成稀疏矩陣？比如我把上面的哈密頓寫出來(lái)之后，怎么把它轉(zhuǎn)成稀疏矩陣？ 這個(gè)小程序就是上面這個(gè)哈密頓矩陣的寫法。 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 我想把它轉(zhuǎn)成稀疏矩陣，該怎么辦呢？sparse命令到底是怎么用的呢？糾結(jié)�。�！ |
回復(fù)此樓» 猜你喜歡
材料與化工（0856）304求B區(qū)調(diào)劑
已經(jīng)有7人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有8人回復(fù)
-
285求調(diào)劑
已經(jīng)有7人回復(fù)
-
317求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿華中師范化學(xué)332分求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑
已經(jīng)有25人回復(fù)
-
一志愿南昌大學(xué)324求調(diào)劑
已經(jīng)有9人回復(fù)
-
材料科學(xué)與工程 317求調(diào)劑
已經(jīng)有7人回復(fù)
-
085701環(huán)境工程求調(diào)劑
已經(jīng)有11人回復(fù)
-
286分調(diào)劑
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- FORTRAN中比雙精度（real（kind=8））更高的精度的數(shù)據(jù)類型是什么？
已經(jīng)有10人回復(fù)
- 想用Fortran語(yǔ)言實(shí)現(xiàn)矩陣連乘，并且矩陣?yán)锩嬗幸粋�(gè)變量
已經(jīng)有6人回復(fù)
- fortran中如何創(chuàng)建由二維矩陣為元素的1維矩陣
已經(jīng)有4人回復(fù)
- 請(qǐng)問(wèn)在fortran里面怎么求復(fù)數(shù)矩陣的行列式�。�！
已經(jīng)有4人回復(fù)
- 求助：利用fortran寫一個(gè)矩陣求逆的程序
已經(jīng)有14人回復(fù)
- fortran 怎么按行讀取數(shù)據(jù)（有小例子）
已經(jīng)有21人回復(fù)
- 編一個(gè)fortran程序
已經(jīng)有17人回復(fù)
- 關(guān)于6*6剛度系數(shù)矩陣和9*9剛度系數(shù)矩陣的問(wèn)題
已經(jīng)有4人回復(fù)
- fortran real(kind=16)
已經(jīng)有10人回復(fù)
- 稀疏矩陣的處理
已經(jīng)有12人回復(fù)
- 怎么將此for循環(huán)改為矩陣形式？
已經(jīng)有5人回復(fù)
- fortran輸出求助
已經(jīng)有6人回復(fù)
- fortran如何實(shí)現(xiàn)張量相乘（并乘）
已經(jīng)有9人回復(fù)
- 用fortran存儲(chǔ)稀疏矩陣
已經(jīng)有4人回復(fù)
- fortran中到底該如何比較兩個(gè)實(shí)數(shù)是否相等
已經(jīng)有8人回復(fù)
- 【求助】fortran求解矩陣
已經(jīng)有7人回復(fù)
- 【求助】用fortran怎么描述下面這樣的情況，我是新手，請(qǐng)大家?guī)兔?/a>
已經(jīng)有8人回復(fù)
- 【求助】如何學(xué)習(xí)FORTRAN的LAPACK?
已經(jīng)有4人回復(fù)
- 【求助】Fortran 如何讀取文件中特定行的內(nèi)容
已經(jīng)有16人回復(fù)
virtualzx
木蟲(chóng) (著名寫手)
- 應(yīng)助: 263 (大學(xué)生)
- 金幣: 7161.3
- 紅花: 54
- 帖子: 1605
- 在線: 317.6小時(shí)
- 蟲(chóng)號(hào): 2069080
- 注冊(cè): 2012-10-18
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
【答案】應(yīng)助回帖
★
感謝參與，應(yīng)助指數(shù) +1
jjdg: 金幣+1, 感謝參與 2012-11-10 01:05:08
感謝參與，應(yīng)助指數(shù) +1
jjdg: 金幣+1, 感謝參與 2012-11-10 01:05:08
| 可以試試這個(gè) http://math.nist.gov/~KRemington/fspblas/ |
2樓2012-11-09 03:03:32
3樓2012-11-09 19:21:31
weiwei43
至尊木蟲(chóng) (正式寫手)
- 應(yīng)助: 16 (小學(xué)生)
- 金幣: 18316.1
- 紅花: 6
- 帖子: 737
- 在線: 94.7小時(shí)
- 蟲(chóng)號(hào): 250282
- 注冊(cè): 2006-05-13
- 專業(yè): 機(jī)械動(dòng)力學(xué)
4樓2012-11-10 08:22:12
5樓2012-11-10 08:31:44
weiwei43
至尊木蟲(chóng) (正式寫手)
- 應(yīng)助: 16 (小學(xué)生)
- 金幣: 18316.1
- 紅花: 6
- 帖子: 737
- 在線: 94.7小時(shí)
- 蟲(chóng)號(hào): 250282
- 注冊(cè): 2006-05-13
- 專業(yè): 機(jī)械動(dòng)力學(xué)
★
jjdg: 金幣+1, 感謝參與 2012-11-11 01:12:48
jjdg: 金幣+1, 感謝參與 2012-11-11 01:12:48
|
我這里主要是利用稀疏矩陣特點(diǎn)存儲(chǔ)，至于求解方法什么的適當(dāng)?shù)母鶕?jù)存儲(chǔ)特點(diǎn)修改了一下，主要是減少計(jì)算量。稀疏矩陣可以求特征值和特征向量我感覺(jué)是和一般的方陣求解是一樣的。 關(guān)于書(shū)，我現(xiàn)在手頭沒(méi)有，沒(méi)法給你提出具體數(shù)目。你可以看一下70-80年代的有限元書(shū)或結(jié)構(gòu)分析書(shū)，很多書(shū)后面都附有程序的。 |
6樓2012-11-10 08:49:45
7樓2012-11-10 20:52:35
mariohollis
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 5
- 帖子: 1
- 在線: 3小時(shí)
- 蟲(chóng)號(hào): 3605349
- 注冊(cè): 2014-12-21
8樓2015-06-15 10:19:10
9樓2015-06-15 10:41:36
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 田山東 的主題更新
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +7 | 十閑wx 2026-03-24 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分，前三科223分 +6 | dongfang59 2026-03-27 | 6/300 |
|
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +14 | 北道巷 2026-03-24 | 14/700 |
|
|
|
[考研] 一志愿鄭大材料工程290求調(diào)劑 +9 | Youth_ 2026-03-30 | 9/450 |
|
|
|
[考研] 287求調(diào)劑 +14 | land xuxu 2026-03-26 | 14/700 |
|
|
|
[考研] 277跪求調(diào)劑 +7 | 1915668 2026-03-27 | 11/550 |
|
|
|
[碩博家園] 求調(diào)劑 有機(jī)化學(xué)考研356分 +10 | Nadiums 2026-03-25 | 11/550 |
|
|
|
[考研] 085600，材料與化工321分求調(diào)劑 +10 | 大饞小子 2026-03-28 | 10/500 |
|
|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +7 | @taotao 2026-03-29 | 7/350 |
|
|
|
[考研] 一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
|
[考研] 347求調(diào)劑 +3 | 山頂見(jiàn)α 2026-03-25 | 3/150 |
|
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么！ 2026-03-26 | 7/350 |
|
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
|
[考博] 26申博 +3 | 加油沖��！ 2026-03-26 | 3/150 |
|
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
|
[考研]
材料調(diào)劑
 5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|