#!/bin/bash

nucleus=Ba132
pN=1
nN=1
jsp="example/sdg7h11/CQi50-82.jsp"
xpn="example/sdg7h11/CQi50-82.pn"
pjpair="output/best-proton.pair0"
njpair="output/best-neutron.pair0"

echo "======================================================================="
echo "             variation starts "
echo "-----------------------------------------------------------------------"
./pvpc.x << input
v-wo-Pandya               # 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 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-wo-Pandya                     // 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 "======================================================================="