24小时热门版块排行榜    

CyRhmU.jpeg
查看: 289  |  回复: 4
当前主题已经存档。

xiaojie7783

金虫 (小有名气)

[交流] 【讨论】测试计算

对于VASP的计算参数,大家都是怎么做测试的呢?是用个脚本测试吗?我见到这个脚本:#!/bin/sh
rm WAVECAR
for i in 5 7 9 11 13 15
do
cat > KPOINTS < Automatic generation
0
Monhkorst-Pack
$i $i $i
0.0 0.0 0.0
!
echo "k mesh = $i x $i x $i" ; time ./vasp
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n",$5 }'`
KP=`grep "irreducible" OUTCAR | tail -1 | awk '{printf "%5i \n",$2 }'`
echo $i $KP $E >>comment
done

但是测试出来comment文件里面的数据不全:
5 18 -107.791627
5 0 0.000000
7 40 -107.793812
7 75 -107.793818
9 75 -107.793736
5 18 -107.791626
7 40 -107.793812
9 75 -107.793736
9 126 -107.793746
11 126 -107.793743
13 196 -107.793743
15 288 -107.793746
11 288 -107.793752
13 196 -107.793750
15 288 -107.793754
这是怎么回事呢?
再测试一遍之后更加不全了

希望大家给些建议啊
回复此楼
成功其实很简单,就是当你坚持不住的时候,再坚持一下!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaojie7783

金虫 (小有名气)

是K点测试
成功其实很简单,就是当你坚持不住的时候,再坚持一下!
2楼2009-06-19 10:34:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

condensed

木虫 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
直接for循环计算完,grep一下就完了呗
3楼2009-06-19 11:06:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaojie7783

金虫 (小有名气)

引用回帖:
Originally posted by condensed at 2009-6-19 11:06:
直接for循环计算完,grep一下就完了呗

那为什么我计算出来的结果不全呢?是循环的不对吗?
成功其实很简单,就是当你坚持不住的时候,再坚持一下!
4楼2009-06-19 11:16:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

condensed

木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zdhlover(金币+1,VIP+0):谢谢,鼓励新虫,欢迎常来 6-18 23:23
你可以再计算(./vasp)后面,加上
cp OUTCAT OUTCAR_$i
这样,计算的结果文件就保留下来了
5楼2009-06-19 12:21:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaojie7783 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见