版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

田山東

捐助貴賓 (著名寫手)

應助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時
蟲號: 800307
注冊: 2009-06-28
性別: GG
專業(yè): 基礎(chǔ)物理學

引用回帖:

17樓: Originally posted by lurencyj at 2012-10-12 15:39:51
順便問一個問題，樓主為什么不用matlab來求本征值問題呢，如果內(nèi)存不夸張的話？

老大，內(nèi)存很夸張，體系大�。簄=3*24*160，復數(shù)矩陣H(n,n)，光這樣的矩陣就三個

贊一下

回復此樓

everythinghasitsseason.enjoyyourlife.

21樓2012-10-12 16:25:46

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

田山東

捐助貴賓 (著名寫手)

應助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時
蟲號: 800307
注冊: 2009-06-28
性別: GG
專業(yè): 基礎(chǔ)物理學

引用回帖:

20樓: Originally posted by snoopyzhao at 2012-10-12 16:10:08
其實沒有必要用 ln -s 形成軟鏈接的

在 gfortran 編譯的時候，記得把 libarpack.a 的路徑寫全就可以了……...

老大，怎么找到這個庫文件的路徑啊?超級菜鳥問題了。我裝好arpack之后在文件夾ARPACK里有個libarpack_SUN4.a,是它嗎？

贊一下

回復此樓

everythinghasitsseason.enjoyyourlife.

22樓2012-10-12 16:27:48

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

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學

【答案】應助回帖

引用回帖:

22樓: Originally posted by 田山東 at 2012-10-12 16:27:48
老大，怎么找到這個庫文件的路徑啊?超級菜鳥問題了。我裝好arpack之后在文件夾ARPACK里有個libarpack_SUN4.a,是它嗎？...

我認為就是它，當然，我不太清楚你是怎么安裝的，呵呵……

贊一下

回復此樓

23樓2012-10-12 17:24:40

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

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學

【答案】應助回帖

★ ★ ★ ★ ★
田山東: 金幣+5, ★★★★★最佳答案, 那個文件看到了，說，如果改變libarpack的路徑要干嘛干嘛的 2012-10-12 17:56:38

引用回帖:

你研究一下那個 ARPACK/example/simple/makefile 就知道了……

贊一下(1人)

回復此樓

24樓2012-10-12 17:26:46

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

田山東

捐助貴賓 (著名寫手)

應助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時
蟲號: 800307
注冊: 2009-06-28
性別: GG
專業(yè): 基礎(chǔ)物理學

引用回帖:

24樓: Originally posted by snoopyzhao at 2012-10-12 17:26:46
你研究一下那個 ARPACK/example/simple/makefile 就知道了……...

我又重新編譯了一下ARmake.inc,它給出的libarpack_linux.a的路徑是：
/home/tianhongyu/arpack/ARPACK/libarpack_LINUX.a，而我剛開始用的路徑是：/arpack/ARPACK/libarpack_LINUX.a，所以會出錯。
現(xiàn)在路徑找到了，發(fā)現(xiàn)還是有問題：
/home/tianhongyu/arpack/ARPACK/libarpack_SUN4.a(second.o): In function `second_':
second.f

.text+0x11): undefined reference to `etime_'
collect2: ld returned 1 exit status
問題正在一個一個浮出水面。阿門，上帝來幫下忙吧。

贊一下

回復此樓

everythinghasitsseason.enjoyyourlife.

25樓2012-10-12 18:16:22

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

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
田山東: 金幣+10, ★★★★★最佳答案 2012-10-12 19:31:37

引用回帖:

25樓: Originally posted by 田山東 at 2012-10-12 18:16:22
我又重新編譯了一下ARmake.inc,它給出的libarpack_linux.a的路徑是：
/home/tianhongyu/arpack/ARPACK/libarpack_LINUX.a，而我剛開始用的路徑是：/arpack/ARPACK/libarpack_LINUX.a，所以會出錯。
現(xiàn)在路徑找到 ...

你當時沒有修改 second.f 嗎？
實在不行，把那個 etime 改成 Fortran 95 內(nèi)置的函數(shù) cpu_time 算了……

贊一下

回復此樓

26樓2012-10-12 19:05:27

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

田山東

捐助貴賓 (著名寫手)

應助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時
蟲號: 800307
注冊: 2009-06-28
性別: GG
專業(yè): 基礎(chǔ)物理學

引用回帖:

26樓: Originally posted by snoopyzhao at 2012-10-12 19:05:27
你當時沒有修改 second.f 嗎？
實在不行，把那個 etime 改成 Fortran 95 內(nèi)置的函數(shù) cpu_time 算了……...

second.f這個程序在哪里我都不知道啊。
剛才又看了下ARmake.inc，發(fā)現(xiàn)沒有etime這項，這個選項在哪里修改呢？

贊一下

回復此樓

everythinghasitsseason.enjoyyourlife.

27樓2012-10-12 19:37:10

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

田山東

捐助貴賓 (著名寫手)

應助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時
蟲號: 800307
注冊: 2009-06-28
性別: GG
專業(yè): 基礎(chǔ)物理學

引用回帖:

second.f找到了。把etime改成cpu_time還是不行。    SUBROUTINE SECOND( T )
   REAL T
   CALL CPU_TIME(T)
   RETURN
   END
再一次編譯EXAMPLE下某個程序會出現(xiàn)這種錯誤：
/tmp/ccdi4uxo.o: In function `MAIN__':
dndrv2.f

.text+0x31b): undefined reference to `dgttrf_'
dndrv2.f

.text+0x693): undefined reference to `dgttrs_'
/home/tianhongyu/arpack/ARPACK/libarpack_LINUX.a(second.o): In function `second_':
second.f

.text+0x11): undefined reference to `etime_'
collect2: ld returned 1 exit status
windows下有人安裝了，這種方法，不知道是否可以用在linux下，因為有很多人在linux下裝過arpack都沒這么搞過，所以不確定啊：
先到ARPACK的官方網(wǎng)站去下載程序包，我選擇下載了 arpack96.tar.gz 和補丁 patch.tar.gz。
1. 把兩個壓縮文件分別解壓以后，第一步是先打上補丁。復制patch下的ARPACK 粘貼到 arpack96 里面去覆蓋同名文件，在Win7上會提示合并文件夾，選擇是，然后會提示復制和替換文件，選擇全部替換就行了。
2. 把 arpack96\ARPACK\BLAS，arpack96\ARPACK\LAPACK 和 arpack96\ARPACK\UTIL 中的 *.f 文件都復制到 arpack96\ARPACK\SRC 里面去。然后有三個源文件需要改動。首先，找到second.f 打開，把里面的內(nèi)容改成：
   SUBROUTINE SECOND( T )
   REAL T
   CALL CPU_TIME(T)
   RETURN
   END
保存。因為原來的程序要調(diào)用ETIME來返回時間，但是連接的時候會出錯提示找不到ETIME。我看了一下gfortran 里好像是用 CPU_TIME 來計時的，所以就改過來。還有就是 dnrm2.f 和 snrm2.f 使用了ASSIGN statement，這個語句在Fortran 95 的標準里面已經(jīng)被刪除了，所以gfortran 不能編譯這兩個文件。因此，需要到netlib上下載它們的更新版本： dnrm2.f 和 snrm2.f 并替換舊版文件。
3. 接下來編譯ARPACK和LAPACK基本原理是類似的。用DOS進入 arpack96\ARPACK\SRC 目錄，用gfortran 開始編譯：
E:\arpack96\ARPACK\SRC>gfortran -fimplicit-none -g -c *.f -O2
4. 用ar 連接成一個庫
E:\arpack96\ARPACK\SRC>ar rcs libarpack.a *.o
這樣就在SRC文件夾里得到ARPACK庫文件 libarpack.a。

贊一下

回復此樓

everythinghasitsseason.enjoyyourlife.

28樓2012-10-12 20:14:56

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

田山東

捐助貴賓 (著名寫手)

應助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時
蟲號: 800307
注冊: 2009-06-28
性別: GG
專業(yè): 基礎(chǔ)物理學

網(wǎng)上也有人碰到類似問題，他們是這么總結(jié)的，幫忙看下用在我的系統(tǒng)里該怎么處理？

dseupd.f

.text+0x1445): undefined reference to `dorm2r_'
dseupd.f

.text+0x1580): undefined reference to `dorm2r_'
dseupd.f

.text+0x167d): undefined reference to `dscal_'

從錯誤來看，是沒有找到 BLAS 和 LAPACK 的庫函數(shù)。
我重新修改了 ARmake.inc 處 BLAS 和 LAPACK 的路徑，不使用 ARPACK
自帶的BLAS 和把相對
路徑改成絕對路徑
BLASdir = /usr/local/share/BLAS
LAPACKdir = /usr/local/share/lapack-3.1.1
給出了庫的路徑
LAPACKLIB = /usr/local/lib/liblapack.a
BLASLIB = /usr/local/lib/libblas.a
再次編譯了一次 make lib

這一次我編譯了 bsbdar1
make bsbdr1
tlu@linux-1482:~/proj/arpack/EXAMPLES/BAND> make dsbdr1
Making dsbdr1.o from dsbdr1.f
f77 -O dsbdr1.o dsband.o /usr/local/share/ARPACK/libarpack_LINUX.a /usr/local/lib/liblapack.a /usr/local/lib/libblas.a -o dsbdr1

成功了，原來 dsband.f 是不能單獨成為可執(zhí)行程序的, 他是提供一個骨架，
具體的實現(xiàn)要有類似于 dsbdr1.f 的文件實現(xiàn)。

贊一下

回復此樓

everythinghasitsseason.enjoyyourlife.

29樓2012-10-12 20:17:26

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

田山東

捐助貴賓 (著名寫手)

應助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時
蟲號: 800307
注冊: 2009-06-28
性別: GG
專業(yè): 基礎(chǔ)物理學

引用回帖:

10樓: Originally posted by lurencyj at 2012-10-12 12:51:32
1. 關(guān)于安裝。建議你不要自己編譯源代碼，如果你沒有頭緒，沒有相應的經(jīng)驗。最好通過suse自帶的軟件管理來安裝arpack. 這里有個鏈接是介紹opensuse下如何安裝軟件的：http://zhidao.baidu.com/question/240291320. ...

您好，可以幫我解決個機械性的問題嗎？現(xiàn)在arpack基本上看懂了，只是有個問題一直沒跨過去。http://www.gaoyang168.com/bbs/viewthread.php?tid=5189740
怎么能從一系列本征值中找到一段本征值，中間一段，比如系統(tǒng)本征值從-9到9之間，而我取-3到3之間的一小段，有沒有相關(guān)程序�。恐x謝�。�！
太麻煩你們了

贊一下

回復此樓

everythinghasitsseason.enjoyyourlife.

30樓2012-11-16 16:54:52

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

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

返回列表

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

[考研] 280分求調(diào)劑一志愿085802 +4	PUMPT 2026-03-22	7/350	2026-03-22 22:13 by 星空星月
[考研] 一志愿西安交通大學材料工程專業(yè) 282分求調(diào)劑 +11	楓橋ZL 2026-03-18	13/650	2026-03-22 20:26 by edmund7
[考研] 一志愿武理材料工程348求調(diào)劑 +5	￣^￣゜汗 2026-03-19	7/350	2026-03-22 19:44 by 公瑾逍遙
[考研] 260求調(diào)劑 +3	朱芷琳 2026-03-20	4/200	2026-03-22 15:12 by 朱芷琳
[考研] 生物學調(diào)劑 +5	Surekei 2026-03-21	5/250	2026-03-22 14:39 by tcx007
[考研] 生物學一志愿985，分數(shù)349求調(diào)劑 +4	zxts12 2026-03-21	7/350	2026-03-22 09:57 by zxts12
[基金申請] 山東省面上項目限額評審 +4	石瑞0426 2026-03-19	4/200	2026-03-22 08:50 by Wei_ren
[考研] 297求調(diào)劑 +3	喜歡還是不甘心 2026-03-20	3/150	2026-03-21 18:33 by 學員8dgXkO
[考研] 296求調(diào)劑 +4	www_q 2026-03-20	4/200	2026-03-21 17:26 by 學員8dgXkO
[考研] 279分求調(diào)劑一志愿211 +14	chaojifeixia 2026-03-19	15/750	2026-03-21 13:24 by zhukairuo
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3	hiloiy 2026-03-17	4/200	2026-03-21 03:04 by JourneyLucky
[考研] 化學求調(diào)劑 +4	臨澤境llllll 2026-03-17	5/250	2026-03-21 02:23 by JourneyLucky
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 288求調(diào)劑 +16	于海海海海 2026-03-19	16/800	2026-03-20 22:28 by JourneyLucky
[考研] 350求調(diào)劑 +5	weudhdk 2026-03-19	5/250	2026-03-20 22:04 by luoyongfeng
[考研] 求調(diào)劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學校友
[考研] 材料與化工求調(diào)劑 +7	為學666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽w
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan