24小时热门版块排行榜    

查看: 1239  |  回复: 4

gemucai

木虫 (正式写手)

[求助] PWSCF在哪里设置lapack?再问一句声子谱对角化出错是lapack的问题吗?

大家好,最近我算声子谱,算了好几天,在最后两步出错,报错如下:

   Band Structure Calculation
      Davidson diagonalization with overlap

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
      from cdiaghg : error #        66
      diagonalization (ZHEGV*) failed
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
有人告诉我说是lapack本身的bug,让我下载最新版本,我下载之后应该如何让程序指向我新加的lapack呢?

还有,这个问题是lapack的问题吗?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

★ ★
xueht987: 金币+2, 专家考核, 鼓励交流! 2012-07-11 17:30:46
在make.sys里修改Lapack库,然后重新编译。
使用david方法进行对角化,的确有容易失败的问题~具体算是lapack的bug还是什么就不清楚了,但是通过调节一些初始的输入参数,这个是可以避免的~
实在不成,可以尝试修改一下代码,把ph.x中david对角化方法修改成cg法……
思想重于技巧,内涵重于表象
2楼2012-07-11 16:00:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

引用回帖:
2楼: Originally posted by souledge at 2012-07-11 16:00:40
在make.sys里修改Lapack库,然后重新编译。
使用david方法进行对角化,的确有容易失败的问题~具体算是lapack的bug还是什么就不清楚了,但是通过调节一些初始的输入参数,这个是可以避免的~
实在不成,可以尝试修 ...

改成cg,确实是个好想法。我碰到的情况很诡异,算三维和一维的体系不会报错,只要算二维的石墨烯的声子就挂了,算二维的别的体系又没问题。
3楼2012-08-20 09:31:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
uuv2010: 金币+3, 专家考核, 多谢详细指导 2012-08-20 11:32:39
gemucai: 金币+5, ★★★很有帮助, 非常感谢。另外,你指的是增加电子计算时候的nbnd吧?可我这是声子计算。 2012-08-20 20:36:35
引用回帖:
3楼: Originally posted by gemucai at 2012-08-20 09:31:24
改成cg,确实是个好想法。我碰到的情况很诡异,算三维和一维的体系不会报错,只要算二维的石墨烯的声子就挂了,算二维的别的体系又没问题。...

不会啊,我算石墨烯从来没挂过啊~
并且,cg虽然基本不出对角化问题,但是计算速度比david慢很多,用goldenfish的话说,就是cg的并行很差~
而且,修改ph.x中的对角化部分并不是太容易,牵扯到不少代码的修改。
个人建议,更换赝势是首选的,其次,适当增加空带数很多时候能够避免对角化出错。
思想重于技巧,内涵重于表象
4楼2012-08-20 11:06:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)


uuv2010: 金币+1, 专家考核, 多谢补充 2012-08-20 11:33:01
顺便多说几句,虽然个人没搞过LAPACK的编程,但是从RSPt的manual看来,使用LAPACK对角化时,有很多种选择:ZHPEV,ZHEEV,ZHEEVD……个人猜测,QE使用ZHEGV……至于有什么区别有什么联系有什么优缺点,偶就不得而知了~等高手来解答。
思想重于技巧,内涵重于表象
5楼2012-08-20 11:11:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gemucai 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 08工科 320总分 求调剂 +11 梨花珞晚风 2026-03-17 11/550 2026-03-22 17:42 by luoyongfeng
[考研] 324求调剂 +6 lucky呀呀呀鸭 2026-03-20 6/300 2026-03-22 16:01 by ColorlessPI
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
[考研] 求调剂 +7 Auroracx 2026-03-22 7/350 2026-03-22 12:38 by 素颜倾城1988
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 材料工程(专)一志愿985 初试335求调剂 +3 hiloiy 2026-03-17 4/200 2026-03-21 03:04 by JourneyLucky
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 一志愿华中科技大学,080502,354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 广西大学家禽遗传育种课题组2026年硕士招生(接收计算机专业调剂) +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
信息提示
请填处理意见