 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 468 | 回復(fù): 2 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
haiqing1987銅蟲 (初入文壇)
|
[交流]
【求助】IMSL調(diào)用出錯，求解！(已解決) 已有2人參與
|
||
|
先上子程序： SUBROUTINE VIBSOL(NQ1max,NQ1,RQ1,TMS,eigenvaluor,Vpot,eigenvalue) use imsl IMPLICIT REAL*8 (A-H,O-Z) DIMENSION RQ1(NQ1),HAR(NQ1,NQ1),Vpot(NQ1),eigenvalue(NQ1), $WK1(1),WK2(1),eigenvaluor(NQ1,NQ1) RL=RQ1(NQ1)-RQ1(1)+2.D0*(RQ1(2)-RQ1(1)) NDVR=NQ1+1 DO I=1,NQ1 HAR(I,I)=Vpot(I)+DVRKE1(I,NDVR,RL,TMS) DO J=1,I-1 HAR(J,I)=DVRKE2(I,J,NDVR,RL,TMS) HAR(I,J)=HAR(J,I) enddo enddo eigenvalue=eig(HAR,v=eigenvaluor) RETURN END 運行時提示堆棧出錯，程序崩潰。調(diào)試時發(fā)現(xiàn)崩在了 EIG函數(shù)上，提示如截圖。第一次用IMSL庫，懂的不多。請教高手，代碼哪里有問題？ 謝謝！ 附截圖：  [ Last edited by haiqing1987 on 2010-12-28 at 19:47 ] |
回復(fù)此樓» 猜你喜歡
283分材料與化工求調(diào)劑
已經(jīng)有12人回復(fù)
-
349求調(diào)劑
已經(jīng)有7人回復(fù)
-
296材料專碩求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿上海海洋大學(xué)083200食品學(xué)碩，求調(diào)劑，接受其他專業(yè)
已經(jīng)有7人回復(fù)
-
材料340分調(diào)劑
已經(jīng)有8人回復(fù)
-
316求調(diào)劑
已經(jīng)有13人回復(fù)
-
321求調(diào)劑
已經(jīng)有11人回復(fù)
-
315求調(diào)劑
已經(jīng)有6人回復(fù)
-
330分求調(diào)劑
已經(jīng)有13人回復(fù)
-
化工調(diào)劑303分，過四級
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 高斯09中的L502出錯的問題
已經(jīng)有5人回復(fù)
- 請熱心人能提供一個帶imsl庫的intel visual fortran編譯器下載鏈接
已經(jīng)有6人回復(fù)
- LC-MS 求解
已經(jīng)有23人回復(fù)
- 求解實驗室發(fā)酵酸奶分層了怎么回事？
已經(jīng)有10人回復(fù)
- 出錯 warning#2 #180
已經(jīng)有7人回復(fù)
- 質(zhì)量監(jiān)督系統(tǒng)考試 色弱可以報嗎求解
已經(jīng)有4人回復(fù)
- 請教：orca 中用CCSD(T)計算單點能出錯--提示內(nèi)存不足 該怎么辦？
已經(jīng)有10人回復(fù)
- MS軟件老出錯，求指導(dǎo)
已經(jīng)有3人回復(fù)
- 【求助】matlab 求解微分方程中的未知參數(shù)
已經(jīng)有20人回復(fù)
- 【求助】在進行fluent模擬求解過程中，突然提醒出錯，然后進入死循環(huán)，怎么回事？
已經(jīng)有5人回復(fù)
- 【求助】fortran：如何求解 復(fù)數(shù)矩陣的秩【已完結(jié)】
已經(jīng)有5人回復(fù)
- 【求助】用fortran求解大型線性方程組時出現(xiàn)的錯誤【已解決】
已經(jīng)有11人回復(fù)
- 【轉(zhuǎn)帖】IMSL在IVF中的安裝和應(yīng)用 (Windows系統(tǒng))【搜索無重復(fù)】
已經(jīng)有4人回復(fù)
haiqing1987
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 293.9
- 帖子: 36
- 在線: 39.8小時
- 蟲號: 1155042
- 注冊: 2010-11-24
- 性別: GG
- 專業(yè): 原子和分子物理
3樓2010-12-18 22:21:37
maomao1210
金蟲 (正式寫手)
- 程序強帖: 5
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1432.8
- 散金: 242
- 紅花: 16
- 沙發(fā): 1
- 帖子: 991
- 在線: 442.1小時
- 蟲號: 253215
- 注冊: 2006-05-20
- 性別: MM
- 專業(yè): 考古理論
★ ★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
余澤成(金幣+3):謝謝參與應(yīng)助！ 2010-12-19 19:29:09
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
余澤成(金幣+3):謝謝參與應(yīng)助！ 2010-12-19 19:29:09
|
你這沒有問題呀，你搞錯了吧？ 你的意圖不就是構(gòu)建一個矩陣，調(diào)用庫函數(shù)將之對角化嗎？ 求出本征值與本征矢嘛？ 這是我簡單的測試結(jié)果，用IMSL庫還有另外一個庫的測試結(jié)果，用你的程序：  下面是簡單的測試代碼 可以參考，從你的程序意圖，好使構(gòu)建hermian算符的矩陣，應(yīng)該是實對稱的吧。   [ Last edited by maomao1210 on 2010-12-18 at 14:10 ] |
2樓2010-12-17 19:56:49
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿山東大學(xué)，085600，344 +3 | 魏子per 2026-04-02 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 336求調(diào)劑 +6 | kiyy 2026-04-01 | 6/300 |
|
|
|
[考研] 求調(diào)劑求調(diào)劑 +6 | 121. 2026-04-02 | 6/300 |
|
|
|
[考博] 材料工程專業(yè)碩士申博 +3 | 麟正宇 2026-03-30 | 3/150 |
|
|
|
[考研] 材料工程322分 +8 | 哈哈哈吼吼吼哈 2026-04-01 | 8/400 |
|
|
|
[考研] 292求調(diào)劑 +14 | 是妍子也是研子 2026-03-30 | 14/700 |
|
|
|
[考研] 348環(huán)境工程調(diào)劑 +3 | 吳彥祖24k 2026-04-01 | 3/150 |
|
|
|
[考研] 302求調(diào)劑一志愿北航070300，本科鄭大化學(xué) +8 | 圣日耳曼條 2026-04-01 | 11/550 |
|
|
|
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +14 | 一定上岸aaa 2026-03-27 | 18/900 |
|
|
|
[考研] 085600 一志愿9 總分351 求調(diào)劑學(xué)校 +7 | czhcz 2026-03-31 | 9/450 |
|
|
|
[考研] 349求調(diào)劑 +6 | 吃的不少 2026-04-01 | 6/300 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-30 | 7/350 |
|
|
|
[考研] 326求調(diào)劑 +4 | 崽崽仔 2026-03-31 | 4/200 |
|
|
|
[考研] 375求調(diào)劑 +7 | 雨夏整夜 2026-03-29 | 7/350 |
|
|
|
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +6 | 大饞小子 2026-03-29 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +10 | 家佳佳佳佳佳 2026-03-29 | 10/500 |
|
|
|
[考研] 2026年華南師范大學(xué)歡迎化學(xué)，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入！ +3 | llss0711 2026-03-28 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +7 | 爭取九點睡 2026-03-28 | 8/400 |
|
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|