PVPC kb3g pn/upn Ti44 LAP 脚本

随便记在这吧
LAPkb3g.sh

#!/bin/bash
# ./LAPkb3g.sh 1 1 example/kb3g.jsp example/kb3g.pn y 10 0.001 0.001 projection_pf0202_p.dat projection_pf0202_n.dat

nucleus=test
pN=$1
nN=$2
jsp=$3
pn=$4
Scale=$5
Jmax=$6
fJtol=$7
lineartol=$8
pjpair=$9
njpair=${10}

echo delete output/$nucleus.eig*, output/$nucleus.wave*
rm output/$nucleus.eig*; rm output/$nucleus.wave*

./pvpc.x << Input
LAP-wo-Pandya                     // mode
$nucleus                // nucleus
$pN                     // # of proton pairs
$nN                     // # of neutron pairs
$jsp                    // file for j-scheme single particle space
pn                     // format of interactions
$pn                    // file for interactions
$Scale                       // scale or not
1                       // number of pair condensates
$pjpair   $njpair
$Jmax                      // Jmax to be projected, not doubled
$fJtol   
$lineartol     
Input

LAPupn.sh

#!/bin/bash
# ./LAPupn.sh 1 1 example/kb3g.jsp example/kb3g.upn y 10 0.001 0.001 projection_pf0202_p.dat projection_pf0202_n.dat

nucleus=test
pN=$1
nN=$2
jsp=$3
upn=$4
Scale=$5
Jmax=$6
fJtol=$7
lineartol=$8
pjpair=$9
njpair=${10}

echo delete output/$nucleus.eig*, output/$nucleus.wave*
rm output/$nucleus.eig*; rm output/$nucleus.wave*

./pvpc.x << Input
LAP-wo-Pandya                     // mode
$nucleus                // nucleus
$pN                     // # of proton pairs
$nN                     // # of neutron pairs
$jsp                    // file for j-scheme single particle space
upn                     // format of interactions
$upn                    // file for interactions
$Scale                       // scale or not
1                       // number of pair condensates
$pjpair   $njpair
$Jmax                      // Jmax to be projected, not doubled
$fJtol   
$lineartol     
Input

又又又又又复习了一遍 PandaWarrior 里面 iso -> pn 格式转换，真的难用。真的应该找个时间把那个代码调一下，变成一个好用的可执行程序，然后放进工具箱。