版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

徐誼

木蟲 (著名寫手)

應(yīng)助: 25 (小學(xué)生)
金幣: 3693.4
散金: 2438
紅花: 9
帖子: 1184
在線: 603.9小時
蟲號: 1214236
注冊: 2011-02-26
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] fortran求矩陣本征值出錯已有2人參與

在用Fortran求矩陣本征值的時候，調(diào)用了MKL里面的DSYEVD函數(shù)，程序可以運(yùn)行，但是好像并沒有執(zhí)行這一語句。
程序如下：

CODE:

program test_mkl

      IMPLICIT NONE

      INTEGER          N

      PARAMETER        ( N = 2 )

      INTEGER          LDA

      PARAMETER        ( LDA = N )

      INTEGER          LWMAX

      PARAMETER        ( LWMAX = 1000 )

      INTEGER          INFO, LWORK, LIWORK

      INTEGER          IWORK( LWMAX )

      DOUBLE PRECISION A( LDA, N ), W( N ), WORK( LWMAX )

      DATA A/1,3,3,4/

      WRITE(*,*)'DSYEVD Example Program Results'

      LWORK = -1

      LIWORK = -1

      CALL DSYEVD ( 'V', 'U', N, A, LDA, W, WORK, LWORK,IWORK, LIWORK, INFO )

     ! Print eigenvalues.

      PRINT *, 'Eigenvalues', W

     ! Print eigenvectors.

      PRINT *, 'Eigenvectors (stored columnwise)', A

      STOP

    end

程序運(yùn)行沒問題！但是CALL DSYEVD ( 'V', 'U', N, A, LDA, W, WORK, LWORK,IWORK, LIWORK, INFO ) 這一句貌似沒有執(zhí)行，因?yàn)榉祷氐慕Y(jié)果還是原來的矩陣！
請問這是什么問題呢？謝謝！

1.jpg

回復(fù)此樓

» 猜你喜歡

一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑已經(jīng)有3人回復(fù)
南昌大學(xué)材料專碩311分求調(diào)劑已經(jīng)有6人回復(fù)
316求調(diào)劑已經(jīng)有6人回復(fù)
346求調(diào)劑[0856] 已經(jīng)有7人回復(fù)
一志愿山大07化學(xué) 332分四六級已過本科山東雙非求調(diào)劑！已經(jīng)有3人回復(fù)
310求調(diào)劑已經(jīng)有3人回復(fù)
機(jī)械專碩299求調(diào)劑至材料已經(jīng)有4人回復(fù)
070300化學(xué)319求調(diào)劑已經(jīng)有7人回復(fù)
08工科 320總分求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有12人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

（急急急）求助matlab大型矩陣求本征值 eigs命令使用后，提示錯誤已經(jīng)有7人回復(fù)
fortran中如何調(diào)用matlab中求解特征值函數(shù)？已經(jīng)有9人回復(fù)
求助高斯輸出文件中坐標(biāo)系怎么定義的？已經(jīng)有11人回復(fù)
有誰會求解復(fù)矩陣（且本征值是虛數(shù)）的幺正矩陣��？求大神幫助！！已經(jīng)有3人回復(fù)
求助fortran特征值算法！有大禮！已經(jīng)有4人回復(fù)
計算一特殊對稱矩陣的行列式或者本征值已經(jīng)有3人回復(fù)
fortran能求帶字母矩陣的本征值嗎？（矩陣元素帶有字母，能得到形式解嗎）已經(jīng)有5人回復(fù)
【求助】fortran求解矩陣已經(jīng)有7人回復(fù)
【求助】fortran：如何求解復(fù)數(shù)矩陣的秩【已完結(jié)】已經(jīng)有5人回復(fù)
【求助】C語言如何求解復(fù)數(shù)矩陣的本征值和本征向量已經(jīng)有7人回復(fù)

1樓 2014-11-05 19:49:28

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wlscj

木蟲之王 (文壇精英)

應(yīng)助: 170 (高中生)
金幣: 211474
散金: 28056
紅花: 320
帖子: 32581
在線: 1886.2小時
蟲號: 601878
注冊: 2008-09-13
專業(yè): 應(yīng)用地球物理學(xué)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

調(diào)用了MKL里面的DSYEVD函數(shù)?

DSYEVD 是個啥?

贊一下

回復(fù)此樓

2樓2014-11-06 19:17:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

徐誼

木蟲 (著名寫手)

應(yīng)助: 25 (小學(xué)生)
金幣: 3693.4
散金: 2438
紅花: 9
帖子: 1184
在線: 603.9小時
蟲號: 1214236
注冊: 2011-02-26
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

2樓: Originally posted by wlscj at 2014-11-06 19:17:22
調(diào)用了MKL里面的DSYEVD函數(shù)?

DSYEVD 是個啥?

嗯，是的。求本征值的
?syevd
Computes all eigenvalues and (optionally) all
eigenvectors of a real symmetric matrix using divide
and conquer algorithm.

贊一下

回復(fù)此樓

3樓2014-11-06 19:23:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wlscj

木蟲之王 (文壇精英)

應(yīng)助: 170 (高中生)
金幣: 211474
散金: 28056
紅花: 320
帖子: 32581
在線: 1886.2小時
蟲號: 601878
注冊: 2008-09-13
專業(yè): 應(yīng)用地球物理學(xué)

引用回帖:

3樓: Originally posted by 徐誼 at 2014-11-06 19:23:50
嗯，是的。求本征值的
?syevd
Computes all eigenvalues and (optionally) all
eigenvectors of a real symmetric matrix using divide
and conquer algorithm....

沒看到這個子程序啊

subroutine DSYEVD ( 'V', 'U', N, A, LDA, W, WORK, LWORK,IWORK, LIWORK, INFO )

贊一下

回復(fù)此樓

4樓2014-11-07 06:03:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wlscj

木蟲之王 (文壇精英)

應(yīng)助: 170 (高中生)
金幣: 211474
散金: 28056
紅花: 320
帖子: 32581
在線: 1886.2小時
蟲號: 601878
注冊: 2008-09-13
專業(yè): 應(yīng)用地球物理學(xué)

【答案】應(yīng)助回帖

也不明白MKL是個啥？是你用fortran軟件平臺上自帶的程序庫模塊？
MKL里面的DSYEVD

是不是要

USE MKL?

贊一下

回復(fù)此樓

5樓2014-11-07 06:09:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

徐誼

木蟲 (著名寫手)

應(yīng)助: 25 (小學(xué)生)
金幣: 3693.4
散金: 2438
紅花: 9
帖子: 1184
在線: 603.9小時
蟲號: 1214236
注冊: 2011-02-26
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

5樓: Originally posted by wlscj at 2014-11-07 06:09:51
也不明白MKL是個啥？是你用fortran軟件平臺上自帶的程序庫模塊？
MKL里面的DSYEVD

是不是要

USE MKL?

對的，是vs自帶的函數(shù)庫！應(yīng)該不用加use mkl吧，因?yàn)榱硗獾暮瘮?shù)運(yùn)行就沒問題，不知道為啥

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

6樓2014-11-07 09:37:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fish.yfyh

銅蟲 (小有名氣)

應(yīng)助: 42 (小學(xué)生)
金幣: 607.9
紅花: 4
帖子: 108
在線: 35.9小時
蟲號: 2867094
注冊: 2013-12-12
專業(yè): 化學(xué)生物學(xué)與生物有機(jī)化學(xué)

【答案】應(yīng)助回帖

調(diào)用mkl庫函數(shù)，只需要在編譯的時候，添加調(diào)用靜態(tài)庫或動態(tài)庫即可。例如:

$ ifort test.f90 -L/path/to/mkl/lib.a

贊一下

回復(fù)此樓

7樓2014-11-19 14:04:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主徐誼的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 265求調(diào)劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 307求調(diào)劑 +10	冷笙123 2026-03-17	10/500	2026-03-21 01:54 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 296求調(diào)劑 +6	www_q 2026-03-18	10/500	2026-03-20 23:56 by JourneyLucky
[考研] 295求調(diào)劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] 材料與化工 322求調(diào)劑 +4	然11 2026-03-19	4/200	2026-03-20 22:12 by luoyongfeng
[考研] 265求調(diào)劑 +12	梁梁校校 2026-03-19	13/650	2026-03-20 21:01 by 無際的草原
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +4	一瓶苯甲酸 2026-03-14	4/200	2026-03-20 20:36 by fen_rao
[考研] 319求調(diào)劑 +3	小力氣珂珂 2026-03-20	3/150	2026-03-20 19:47 by JourneyLucky
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 0856調(diào)劑，是學(xué)校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無懈可擊111
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	3/150	2026-03-20 11:58 by 呱呱呱呱叫
[考研] 0703化學(xué)調(diào)劑 +10	妮妮ninicgb 2026-03-15	14/700	2026-03-19 22:59 by 學(xué)員8dgXkO
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 085600材料與化工求調(diào)劑 +6	緒幸與子 2026-03-17	6/300	2026-03-19 13:27 by houyaoxu
[考研] 283求調(diào)劑 +3	聽風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629

24小時熱門版塊排行榜

[求助] fortran求矩陣本征值出錯 已有2人參與

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] fortran求矩陣本征值出錯已有2人參與

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: