24小时热门版块排行榜    

查看: 1980  |  回复: 2

zjucy

金虫 (小有名气)

[求助] chgdiff.pl的CHGCAR文件顺序问题 已有1人参与

chgdiff.pl CHGCAR1 CHGCAR2意思究竟是 CHGCAR2-CHGCAR1?  还是CHGCAR1-CHGCAR2?还请大神们给个准确的说法。
回复此楼

» 猜你喜欢

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

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

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
fzx2008: 金币+3, 谢谢指导 2014-10-15 15:30:38
CHGCAR2 - CHGCAR1

proof 1
程式碼內寫到
usage: chgdiff.pl <reference CHGCAR> <CHGCAR2>
所以第一項是 reference

proof 2
程式碼內相減的部份只有
    $line1 = <IN1>;
    $line2 = <IN2>;
    @line1 = split(/\s+/,$line1);
    @line2 = split(/\s+/,$line2);
    for ($j=1; $j<@line1; $j++) {
        $line1[$j] = $line2[$j]-$line1[$j];
    }
這一小塊
可以看到它是用 line2 - line1 也就是 CHGCAR2 - CHGCAR1
2楼2014-10-14 15:25:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magicmonk

至尊木虫 (著名写手)

引用回帖:
2楼: Originally posted by jpchou at 2014-10-14 15:25:46
CHGCAR2 - CHGCAR1

proof 1
程式碼內寫到
usage: chgdiff.pl <reference CHGCAR> <CHGCAR2>
所以第一項是 reference

proof 2
程式碼內相減的部份只有
    $line1 = <IN1>;
    $lin ...

请问是 自洽-非自洽 吗?
我们最深的恐惧不是我们能力不够,我们最深的恐惧是能力超越了极限。
3楼2014-12-19 16:50:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zjucy 的主题更新
信息提示
请填处理意见