应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » PWSCF在哪里设置lapack？再问一句声子谱对角化出错是lapack的问题吗？
51/1返回列表
查看: 1242  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

gemucai

木虫 (正式写手)

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

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

   Band Structure Calculation
      Davidson diagonalization with overlap

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

还有，这个问题是lapack的问题吗？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-07-11 14:47:48
已阅   回复此楼   关注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 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[考研] 石河子大学（211、双一流）硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 08工科 320总分 求调剂 +11 梨花珞晚风 2026-03-17 11/550 2026-03-22 17:42 by luoyongfeng
[考研] 一志愿华中科技大学071000，求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 材料工程专硕 348分求调剂 +3 冬辞. 2026-03-17 5/250 2026-03-21 18:47 by 学员8dgXkO
[考研] 278求调剂 +9 烟火先于春 2026-03-17 9/450 2026-03-21 17:47 by 学员8dgXkO
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 306求0703调剂一志愿华中师范 +5 纸鱼ly 2026-03-21 5/250 2026-03-21 17:11 by 学员8dgXkO
[考研] 346求调剂[0856] +4 WayneLim327 2026-03-16 7/350 2026-03-21 04:02 by JourneyLucky
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 一志愿中国海洋大学，生物学，301分，求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
信息提示
关闭
请填处理意见
关闭