24小时热门版块排行榜    

查看: 1443  |  回复: 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的回帖

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的回帖
查看全部 5 个回答

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 PbwVEMK5haaE 2026-08-25 4/200 2026-08-26 03:20 by cNXvBfCpiZOM
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 7K1CJE38xLG4 2026-08-25 4/200 2026-08-26 02:04 by cNXvBfCpiZOM
[基金申请] 今天务委会开完了,明天出结果吗 +18 angus9576 2026-08-25 22/1100 2026-08-26 00:41 by merchancy
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见