标签:脚本 PVPC xpn 代码 gitee pf example best
我们写的 PVPC 的代码在 gitee 上可以下载:https://gitee.com/luyi07/pvpc.git
这里摘录一个脚本,用来自动做变分+投影,在最新的版本里会放在 example/run.sh 中,运行方法:
bash run.sh
#!/bin/bash
nucleus=Cr48
pN=2
nN=2
jsp="example/pf/pf.jsp"
xpn="example/pf/pf+pairing.xpn"
pjpair="best-proton.pair01"
njpair="best-neutron.pair01"
:<<!
echo "======================================================================="
echo " variation starts "
echo "-----------------------------------------------------------------------"
./pvpc.x << input
v # mode: optimize
$nucleus # nucleus
$pN # number of proton pairs
$nN # number of neutron pairs
$jsp # file for j-scheme single particle orbits
xpn # format of interactions: J-scheme P+Q interactions
$xpn # file for j-scheme P+Q interactions
y # scaling
1 # number of optimized pair condensates ordered, though they may turn out to be linearly dependent.
y # start from random pairs
input
echo "-----------------------------------------------------------------------"
echo " variation ends "
echo "======================================================================="
!
echo "======================================================================="
echo " LAP starts "
echo "-----------------------------------------------------------------------"
./pvpc.x << Input
LAP // mode
$nucleus // nucleus
$pN // # of proton pairs
$nN // # of neutron pairs
$jsp // file for j-scheme single particle space
xpn // format of interactions
$xpn // file for interactions
y // scale or not
1 // number of pair condensates
$pjpair $njpair
16 // Jmax to be projected, not doubled
Input
echo "-----------------------------------------------------------------------"
echo " LAP ends "
echo "======================================================================="
标签:脚本,PVPC,xpn,代码,gitee,pf,example,best 来源: https://www.cnblogs.com/luyi07/p/14958726.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。