This commit is contained in:
kossoski 2022-02-18 18:24:42 +01:00
parent 367d6159ba
commit 7485e2b9b0
13 changed files with 515 additions and 67 deletions

View File

@ -156,14 +156,14 @@ paste det.dat a.dat D.dat > det_aD.dat
grep 'FCI' det_aD.dat | sed 's/FCI/1 /g' > det_aD_FCI.dat
grep 'FCI' det_aD.dat | sed 's/FCI/1.0E9 /g' >> det_aD_FCI.dat
grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_CI.dat
grep 'hfCIS ' det_aD.dat | sed 's/hfCIS /1 /g' >> det_aD_CI.dat
grep 'hfCISD ' det_aD.dat | sed 's/hfCISD /2 /g' >> det_aD_CI.dat
#grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_CI.dat
#grep 'hfCIS ' det_aD.dat | sed 's/hfCIS /1 /g' >> det_aD_CI.dat
grep 'hfCISD ' det_aD.dat | sed 's/hfCISD /2 /g' > det_aD_CI.dat
grep 'hfCISDT ' det_aD.dat | sed 's/hfCISDT /3 /g' >> det_aD_CI.dat
grep 'hfCISDTQ' det_aD.dat | sed 's/hfCISDTQ/4 /g' >> det_aD_CI.dat
grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_ooCI.dat
grep 'ooCIS ' det_aD.dat | sed 's/ooCIS /1 /g' >> det_aD_ooCI.dat
#grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_ooCI.dat
grep 'ooCIS ' det_aD.dat | sed 's/ooCIS /1 /g' > det_aD_ooCI.dat
grep 'ooCISD ' det_aD.dat | sed 's/ooCISD /2 /g' >> det_aD_ooCI.dat
grep 'ooCISDT ' det_aD.dat | sed 's/ooCISDT /3 /g' >> det_aD_ooCI.dat
@ -174,8 +174,8 @@ grep 'hfCIs4 ' det_aD.dat | sed 's/hfCIs4/4 /g' >> det_aD_CIs.dat
grep 'ooCIs0 ' det_aD.dat | sed 's/ooCIs0/0 /g' > det_aD_ooCIs.dat
grep 'ooCIs2 ' det_aD.dat | sed 's/ooCIs2/2 /g' >> det_aD_ooCIs.dat
grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_CIo.dat
grep 'hfCIo1 ' det_aD.dat | sed 's/hfCIo1/1 /g' >> det_aD_CIo.dat
#grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_CIo.dat
grep 'hfCIo1 ' det_aD.dat | sed 's/hfCIo1/1 /g' > det_aD_CIo.dat
grep 'hfCIo1.5' det_aD.dat | sed 's/hfCIo1.5/1.5 /g' >> det_aD_CIo.dat
grep 'hfCIo2 ' det_aD.dat | sed 's/hfCIo2/2 /g' >> det_aD_CIo.dat
grep 'hfCIo2.5' det_aD.dat | sed 's/hfCIo2.5/2.5 /g' >> det_aD_CIo.dat
@ -183,8 +183,8 @@ grep 'hfCIo3 ' det_aD.dat | sed 's/hfCIo3/3 /g' >> det_aD_CIo.dat
grep 'hfCIo3.5' det_aD.dat | sed 's/hfCIo3.5/3.5 /g' >> det_aD_CIo.dat
grep 'hfCIo4 ' det_aD.dat | sed 's/hfCIo4/4 /g' >> det_aD_CIo.dat
grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_ooCIo.dat
grep 'ooCIo1 ' det_aD.dat | sed 's/ooCIo1/1 /g' >> det_aD_ooCIo.dat
#grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_ooCIo.dat
grep 'ooCIo1 ' det_aD.dat | sed 's/ooCIo1/1 /g' > det_aD_ooCIo.dat
grep 'ooCIo1.5' det_aD.dat | sed 's/ooCIo1.5/1.5 /g' >> det_aD_ooCIo.dat
grep 'ooCIo2 ' det_aD.dat | sed 's/ooCIo2/2 /g' >> det_aD_ooCIo.dat
grep 'ooCIo2.5' det_aD.dat | sed 's/ooCIo2.5/2.5 /g' >> det_aD_ooCIo.dat
@ -255,14 +255,14 @@ paste det.dat xe.dat D.dat > det_xe.dat
grep 'FCI' det_xe.dat | sed 's/FCI/1 /g' > det_xe_FCI.dat
grep 'FCI' det_xe.dat | sed 's/FCI/1E9 /g' >> det_xe_FCI.dat
grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_CI.dat
grep 'hfCIS ' det_xe.dat | sed 's/hfCIS /1 /g' >> det_xe_CI.dat
grep 'hfCISD ' det_xe.dat | sed 's/hfCISD /2 /g' >> det_xe_CI.dat
#grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_CI.dat
#grep 'hfCIS ' det_xe.dat | sed 's/hfCIS /1 /g' >> det_xe_CI.dat
grep 'hfCISD ' det_xe.dat | sed 's/hfCISD /2 /g' > det_xe_CI.dat
grep 'hfCISDT ' det_xe.dat | sed 's/hfCISDT /3 /g' >> det_xe_CI.dat
grep 'hfCISDTQ' det_xe.dat | sed 's/hfCISDTQ/4 /g' >> det_xe_CI.dat
grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_ooCI.dat
grep 'ooCIS ' det_xe.dat | sed 's/ooCIS /1 /g' >> det_xe_ooCI.dat
#grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_ooCI.dat
grep 'ooCIS ' det_xe.dat | sed 's/ooCIS /1 /g' > det_xe_ooCI.dat
grep 'ooCISD ' det_xe.dat | sed 's/ooCISD /2 /g' >> det_xe_ooCI.dat
grep 'ooCISDT ' det_xe.dat | sed 's/ooCISDT /3 /g' >> det_xe_ooCI.dat
@ -273,8 +273,8 @@ grep 'hfCIs4 ' det_xe.dat | sed 's/hfCIs4/4 /g' >> det_xe_CIs.dat
grep 'ooCIs0 ' det_xe.dat | sed 's/ooCIs0/0 /g' > det_xe_ooCIs.dat
grep 'ooCIs2 ' det_xe.dat | sed 's/ooCIs2/2 /g' >> det_xe_ooCIs.dat
grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_CIo.dat
grep 'hfCIo1 ' det_xe.dat | sed 's/hfCIo1/1 /g' >> det_xe_CIo.dat
#grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_CIo.dat
grep 'hfCIo1 ' det_xe.dat | sed 's/hfCIo1/1 /g' > det_xe_CIo.dat
grep 'hfCIo1.5' det_xe.dat | sed 's/hfCIo1.5/1.5 /g' >> det_xe_CIo.dat
grep 'hfCIo2 ' det_xe.dat | sed 's/hfCIo2/2 /g' >> det_xe_CIo.dat
grep 'hfCIo2.5' det_xe.dat | sed 's/hfCIo2.5/2.5 /g' >> det_xe_CIo.dat
@ -282,8 +282,8 @@ grep 'hfCIo3 ' det_xe.dat | sed 's/hfCIo3/3 /g' >> det_xe_CIo.dat
grep 'hfCIo3.5' det_xe.dat | sed 's/hfCIo3.5/3.5 /g' >> det_xe_CIo.dat
grep 'hfCIo4 ' det_xe.dat | sed 's/hfCIo4/4 /g' >> det_xe_CIo.dat
grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_ooCIo.dat
grep 'ooCIo1 ' det_xe.dat | sed 's/ooCIo1/1 /g' >> det_xe_ooCIo.dat
#grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_ooCIo.dat
grep 'ooCIo1 ' det_xe.dat | sed 's/ooCIo1/1 /g' > det_xe_ooCIo.dat
grep 'ooCIo1.5' det_xe.dat | sed 's/ooCIo1.5/1.5 /g' >> det_xe_ooCIo.dat
grep 'ooCIo2 ' det_xe.dat | sed 's/ooCIo2/2 /g' >> det_xe_ooCIo.dat
grep 'ooCIo2.5' det_xe.dat | sed 's/ooCIo2.5/2.5 /g' >> det_xe_ooCIo.dat

292
Be2_cc-pvtz/fit_eq_0.sh Executable file
View File

@ -0,0 +1,292 @@
#!/bin/bash
#######################################################################
if [ -f "fit.log" ]; then
rm "fit.log"
fi
gnuplot fit_eq.gnu
#######################################################################
grep ^a fit.log | grep = | sed 's/=//g' | cut -c -30 |
sed 's/a0 /FCI /g' |
sed 's/a1 /RHF /g' |
sed 's/a2 /hfCISD /g' |
sed 's/a3 /hfCISDT /g' |
sed 's/a4 /hfCISDTQ /g' |
sed 's/a5 /hfCIs0 /g' |
sed 's/a6 /hfCIs2 /g' |
sed 's/a7 /hfCIs4 /g' |
sed 's/a8 /hfCIo1 /g' |
sed 's/a9 /hfCIo1.5 /g' |
sed 's/a10/hfCIo2 /g' |
sed 's/a11/hfCIo2.5 /g' |
sed 's/a12/hfCIo3 /g' |
sed 's/a13/hfCIo3.5 /g' |
sed 's/a14/hfCIo4 /g' |
sed 's/a15/ooCISD /g' |
sed 's/a16/ooCISDT /g' |
sed 's/a17/ooCIs0 /g' |
sed 's/a18/ooCIs2 /g' |
sed 's/a19/ooCIo1 /g' |
sed 's/a20/ooCIo1.5 /g' |
sed 's/a21/ooCIo2 /g' |
sed 's/a22/ooCIo2.5 /g' |
sed 's/a23/ooCIo3 /g' > a.dat
grep ^D fit.log | grep = | sed 's/=//g' | cut -c -30 |
sed 's/D0 /FCI /g' |
sed 's/D1 /RHF /g' |
sed 's/D2 /hfCISD /g' |
sed 's/D3 /hfCISDT /g' |
sed 's/D4 /hfCISDTQ /g' |
sed 's/D5 /hfCIs0 /g' |
sed 's/D6 /hfCIs2 /g' |
sed 's/D7 /hfCIs4 /g' |
sed 's/D8 /hfCIo1 /g' |
sed 's/D9 /hfCIo1.5 /g' |
sed 's/D10/hfCIo2 /g' |
sed 's/D11/hfCIo2.5 /g' |
sed 's/D12/hfCIo3 /g' |
sed 's/D13/hfCIo3.5 /g' |
sed 's/D14/hfCIo4 /g' |
sed 's/D15/ooCISD /g' |
sed 's/D16/ooCISDT /g' |
sed 's/D17/ooCIs0 /g' |
sed 's/D18/ooCIs2 /g' |
sed 's/D19/ooCIo1 /g' |
sed 's/D20/ooCIo1.5 /g' |
sed 's/D21/ooCIo2 /g' |
sed 's/D22/ooCIo2.5 /g' |
sed 's/D23/ooCIo3 /g' > D.dat
paste a.dat D.dat > aD.dat
grep 'FCI' aD.dat | sed 's/FCI/-0.1 /g' > aD_FCI.dat
grep 'FCI' aD.dat | sed 's/FCI/4.1 /g' >> aD_FCI.dat
grep ^a fit.log | grep = | sed 's/=//g' | cut -c -30 > a0.dat
grep 'a0 ' a0.dat | sed 's/a0 /FCI /g' > a.dat
grep 'a1 ' a0.dat | sed 's/a1 /RHF /g' >> a.dat
grep 'a1 ' a0.dat | sed 's/a1 /hfCIS /g' >> a.dat
grep 'a2 ' a0.dat | sed 's/a2 /hfCISD /g' >> a.dat
grep 'a3 ' a0.dat | sed 's/a3 /hfCISDT /g' >> a.dat
grep 'a4 ' a0.dat | sed 's/a4 /hfCISDTQ /g' >> a.dat
grep 'a5 ' a0.dat | sed 's/a5 /hfCIs0 /g' >> a.dat
grep 'a6 ' a0.dat | sed 's/a6 /hfCIs2 /g' >> a.dat
grep 'a7 ' a0.dat | sed 's/a7 /hfCIs4 /g' >> a.dat
grep 'a8 ' a0.dat | sed 's/a8 /hfCIo1 /g' >> a.dat
grep 'a9 ' a0.dat | sed 's/a9 /hfCIo1.5 /g' >> a.dat
grep 'a10 ' a0.dat | sed 's/a10/hfCIo2 /g' >> a.dat
grep 'a11 ' a0.dat | sed 's/a11/hfCIo2.5 /g' >> a.dat
grep 'a12 ' a0.dat | sed 's/a12/hfCIo3 /g' >> a.dat
grep 'a13 ' a0.dat | sed 's/a13/hfCIo3.5 /g' >> a.dat
grep 'a14 ' a0.dat | sed 's/a14/hfCIo4 /g' >> a.dat
grep 'a1 ' a0.dat | sed 's/a1 /ooCIS /g' >> a.dat
grep 'a15 ' a0.dat | sed 's/a15/ooCISD /g' >> a.dat
grep 'a16 ' a0.dat | sed 's/a16/ooCISDT /g' >> a.dat
grep 'a17 ' a0.dat | sed 's/a17/ooCIs0 /g' >> a.dat
grep 'a18 ' a0.dat | sed 's/a18/ooCIs2 /g' >> a.dat
grep 'a19 ' a0.dat | sed 's/a19/ooCIo1 /g' >> a.dat
grep 'a20 ' a0.dat | sed 's/a20/ooCIo1.5 /g' >> a.dat
grep 'a21 ' a0.dat | sed 's/a21/ooCIo2 /g' >> a.dat
grep 'a22 ' a0.dat | sed 's/a22/ooCIo2.5 /g' >> a.dat
grep 'a23 ' a0.dat | sed 's/a23/ooCIo3 /g' >> a.dat
grep ^D fit.log | grep = | sed 's/=//g' | cut -c -30 > D0.dat
grep 'D0 ' D0.dat | sed 's/D0 /FCI /g' > D.dat
grep 'D1 ' D0.dat | sed 's/D1 /RHF /g' >> D.dat
grep 'D1 ' D0.dat | sed 's/D1 /hfCIS /g' >> D.dat
grep 'D2 ' D0.dat | sed 's/D2 /hfCISD /g' >> D.dat
grep 'D3 ' D0.dat | sed 's/D3 /hfCISDT /g' >> D.dat
grep 'D4 ' D0.dat | sed 's/D4 /hfCISDTQ /g' >> D.dat
grep 'D5 ' D0.dat | sed 's/D5 /hfCIs0 /g' >> D.dat
grep 'D6 ' D0.dat | sed 's/D6 /hfCIs2 /g' >> D.dat
grep 'D7 ' D0.dat | sed 's/D7 /hfCIs4 /g' >> D.dat
grep 'D8 ' D0.dat | sed 's/D8 /hfCIo1 /g' >> D.dat
grep 'D9 ' D0.dat | sed 's/D9 /hfCIo1.5 /g' >> D.dat
grep 'D10 ' D0.dat | sed 's/D10/hfCIo2 /g' >> D.dat
grep 'D11 ' D0.dat | sed 's/D11/hfCIo2.5 /g' >> D.dat
grep 'D12 ' D0.dat | sed 's/D12/hfCIo3 /g' >> D.dat
grep 'D13 ' D0.dat | sed 's/D13/hfCIo3.5 /g' >> D.dat
grep 'D14 ' D0.dat | sed 's/D14/hfCIo4 /g' >> D.dat
grep 'D1 ' D0.dat | sed 's/D1 /ooCIS /g' >> D.dat
grep 'D15 ' D0.dat | sed 's/D15/ooCISD /g' >> D.dat
grep 'D16 ' D0.dat | sed 's/D16/ooCISDT /g' >> D.dat
grep 'D17 ' D0.dat | sed 's/D17/ooCIs0 /g' >> D.dat
grep 'D18 ' D0.dat | sed 's/D18/ooCIs2 /g' >> D.dat
grep 'D19 ' D0.dat | sed 's/D19/ooCIo1 /g' >> D.dat
grep 'D20 ' D0.dat | sed 's/D20/ooCIo1.5 /g' >> D.dat
grep 'D21 ' D0.dat | sed 's/D21/ooCIo2 /g' >> D.dat
grep 'D22 ' D0.dat | sed 's/D22/ooCIo2.5 /g' >> D.dat
grep 'D23 ' D0.dat | sed 's/D23/ooCIo3 /g' >> D.dat
if [[ $(grep 'a0 ' a0.dat) ]]; then grep '0 ' det_FCI.dat | cut -c 4- > det.dat; fi
if [[ $(grep 'a1 ' a0.dat) ]]; then grep '0 ' det_CI.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a1 ' a0.dat) ]]; then grep '1 ' det_CI.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a2 ' a0.dat) ]]; then grep '2 ' det_CI.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a3 ' a0.dat) ]]; then grep '3 ' det_CI.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a4 ' a0.dat) ]]; then grep '4 ' det_CI.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a5 ' a0.dat) ]]; then grep '0 ' det_CIs.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a6 ' a0.dat) ]]; then grep '2 ' det_CIs.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a7 ' a0.dat) ]]; then grep '4 ' det_CIs.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a8 ' a0.dat) ]]; then grep '1 ' det_CIo.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a9 ' a0.dat) ]]; then grep '1.5 ' det_CIo.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a10 ' a0.dat) ]]; then grep '2 ' det_CIo.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a11 ' a0.dat) ]]; then grep '2.5 ' det_CIo.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a12 ' a0.dat) ]]; then grep '3 ' det_CIo.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a13 ' a0.dat) ]]; then grep '3.5 ' det_CIo.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a14 ' a0.dat) ]]; then grep '4 ' det_CIo.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a1 ' a0.dat) ]]; then grep '1 ' det_CI.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a15 ' a0.dat) ]]; then grep '2 ' det_CI.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a16 ' a0.dat) ]]; then grep '3 ' det_CI.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a17 ' a0.dat) ]]; then grep '0 ' det_CIs.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a18 ' a0.dat) ]]; then grep '2 ' det_CIs.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a19 ' a0.dat) ]]; then grep '1 ' det_CIo.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a20 ' a0.dat) ]]; then grep '1.5 ' det_CIo.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a21 ' a0.dat) ]]; then grep '2 ' det_CIo.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a22 ' a0.dat) ]]; then grep '2.5 ' det_CIo.dat | cut -c 4- >> det.dat; fi
if [[ $(grep 'a23 ' a0.dat) ]]; then grep '3 ' det_CIo.dat | cut -c 4- >> det.dat; fi
paste det.dat a.dat D.dat > det_aD.dat
grep 'FCI' det_aD.dat | sed 's/FCI/1 /g' > det_aD_FCI.dat
grep 'FCI' det_aD.dat | sed 's/FCI/1.0E8 /g' >> det_aD_FCI.dat
#grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_CI.dat
#grep 'hfCIS ' det_aD.dat | sed 's/hfCIS /1 /g' >> det_aD_CI.dat
grep 'hfCISD ' det_aD.dat | sed 's/hfCISD /2 /g' > det_aD_CI.dat
grep 'hfCISDT ' det_aD.dat | sed 's/hfCISDT /3 /g' >> det_aD_CI.dat
grep 'hfCISDTQ' det_aD.dat | sed 's/hfCISDTQ/4 /g' >> det_aD_CI.dat
#grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_ooCI.dat
#grep 'ooCIS ' det_aD.dat | sed 's/ooCIS /1 /g' >> det_aD_ooCI.dat
grep 'ooCISD ' det_aD.dat | sed 's/ooCISD /2 /g' > det_aD_ooCI.dat
grep 'ooCISDT' det_aD.dat | sed 's/ooCISDT/3 /g' >> det_aD_ooCI.dat
grep 'hfCIs0 ' det_aD.dat | sed 's/hfCIs0/0 /g' > det_aD_CIs.dat
grep 'hfCIs2 ' det_aD.dat | sed 's/hfCIs2/2 /g' >> det_aD_CIs.dat
grep 'hfCIs4 ' det_aD.dat | sed 's/hfCIs4/4 /g' >> det_aD_CIs.dat
grep 'ooCIs0 ' det_aD.dat | sed 's/ooCIs0/0 /g' > det_aD_ooCIs.dat
grep 'ooCIs2 ' det_aD.dat | sed 's/ooCIs2/2 /g' >> det_aD_ooCIs.dat
#grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_CIo.dat
grep 'hfCIo1 ' det_aD.dat | sed 's/hfCIo1/1 /g' > det_aD_CIo.dat
grep 'hfCIo1.5' det_aD.dat | sed 's/hfCIo1.5/1.5 /g' >> det_aD_CIo.dat
grep 'hfCIo2 ' det_aD.dat | sed 's/hfCIo2/2 /g' >> det_aD_CIo.dat
grep 'hfCIo2.5' det_aD.dat | sed 's/hfCIo2.5/2.5 /g' >> det_aD_CIo.dat
grep 'hfCIo3 ' det_aD.dat | sed 's/hfCIo3/3 /g' >> det_aD_CIo.dat
grep 'hfCIo3.5' det_aD.dat | sed 's/hfCIo3.5/3.5 /g' >> det_aD_CIo.dat
grep 'hfCIo4 ' det_aD.dat | sed 's/hfCIo4/4 /g' >> det_aD_CIo.dat
#grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_ooCIo.dat
grep 'ooCIo1 ' det_aD.dat | sed 's/ooCIo1/1 /g' > det_aD_ooCIo.dat
grep 'ooCIo1.5' det_aD.dat | sed 's/ooCIo1.5/1.5 /g' >> det_aD_ooCIo.dat
grep 'ooCIo2 ' det_aD.dat | sed 's/ooCIo2/2 /g' >> det_aD_ooCIo.dat
grep 'ooCIo2.5' det_aD.dat | sed 's/ooCIo2.5/2.5 /g' >> det_aD_ooCIo.dat
grep 'ooCIo3 ' det_aD.dat | sed 's/ooCIo3/3 /g' >> det_aD_ooCIo.dat
#######################################################################
#######################################################################
grep ^xe fit.log | grep = | sed 's/=//g' | cut -c -30 > xe0.dat
grep ^xe fit.log | grep = | sed 's/=//g' | cut -c -30 |
sed 's/xe0 /FCI /g' |
sed 's/xe1 /RHF /g' |
sed 's/xe2 /hfCISD /g' |
sed 's/xe3 /hfCISDT /g' |
sed 's/xe4 /hfCISDTQ /g' |
sed 's/xe5 /hfCIs0 /g' |
sed 's/xe6 /hfCIs2 /g' |
sed 's/xe7 /hfCIs4 /g' |
sed 's/xe8 /hfCIo1 /g' |
sed 's/xe9 /hfCIo1.5 /g' |
sed 's/xe10/hfCIo2 /g' |
sed 's/xe11/hfCIo2.5 /g' |
sed 's/xe12/hfCIo3 /g' |
sed 's/xe13/hfCIo3.5 /g' |
sed 's/xe14/hfCIo4 /g' |
sed 's/xe15/ooCISD /g' |
sed 's/xe16/ooCISDT /g' |
sed 's/xe17/ooCIs0 /g' |
sed 's/xe18/ooCIs2 /g' |
sed 's/xe19/ooCIo1 /g' |
sed 's/xe20/ooCIo1.5 /g' |
sed 's/xe21/ooCIo2 /g' |
sed 's/xe22/ooCIo2.5 /g' |
sed 's/xe23/ooCIo3 /g' > xe.dat
grep ^xe fit.log | grep = | sed 's/=//g' | cut -c -30 > xe0.dat
grep 'xe0 ' xe0.dat | sed 's/xe0 /FCI /g' > xe.dat
grep 'xe1 ' xe0.dat | sed 's/xe1 /RHF /g' >> xe.dat
grep 'xe1 ' xe0.dat | sed 's/xe1 /hfCIS /g' >> xe.dat
grep 'xe2 ' xe0.dat | sed 's/xe2 /hfCISD /g' >> xe.dat
grep 'xe3 ' xe0.dat | sed 's/xe3 /hfCISDT /g' >> xe.dat
grep 'xe4 ' xe0.dat | sed 's/xe4 /hfCISDTQ /g' >> xe.dat
grep 'xe5 ' xe0.dat | sed 's/xe5 /hfCIs0 /g' >> xe.dat
grep 'xe6 ' xe0.dat | sed 's/xe6 /hfCIs2 /g' >> xe.dat
grep 'xe7 ' xe0.dat | sed 's/xe7 /hfCIs4 /g' >> xe.dat
grep 'xe8 ' xe0.dat | sed 's/xe8 /hfCIo1 /g' >> xe.dat
grep 'xe9 ' xe0.dat | sed 's/xe9 /hfCIo1.5 /g' >> xe.dat
grep 'xe10 ' xe0.dat | sed 's/xe10/hfCIo2 /g' >> xe.dat
grep 'xe11 ' xe0.dat | sed 's/xe11/hfCIo2.5 /g' >> xe.dat
grep 'xe12 ' xe0.dat | sed 's/xe12/hfCIo3 /g' >> xe.dat
grep 'xe13 ' xe0.dat | sed 's/xe13/hfCIo3.5 /g' >> xe.dat
grep 'xe14 ' xe0.dat | sed 's/xe14/hfCIo4 /g' >> xe.dat
grep 'xe1 ' xe0.dat | sed 's/xe1 /ooCIS /g' >> xe.dat
grep 'xe15 ' xe0.dat | sed 's/xe15/ooCISD /g' >> xe.dat
grep 'xe16 ' xe0.dat | sed 's/xe16/ooCISDT /g' >> xe.dat
grep 'xe17 ' xe0.dat | sed 's/xe17/ooCIs0 /g' >> xe.dat
grep 'xe18 ' xe0.dat | sed 's/xe18/ooCIs2 /g' >> xe.dat
grep 'xe19 ' xe0.dat | sed 's/xe19/ooCIo1 /g' >> xe.dat
grep 'xe20 ' xe0.dat | sed 's/xe20/ooCIo1.5 /g' >> xe.dat
grep 'xe21 ' xe0.dat | sed 's/xe21/ooCIo2 /g' >> xe.dat
grep 'xe22 ' xe0.dat | sed 's/xe22/ooCIo2.5 /g' >> xe.dat
grep 'xe23 ' xe0.dat | sed 's/xe23/ooCIo3 /g' >> xe.dat
paste det.dat xe.dat D.dat > det_xe.dat
grep 'FCI' det_xe.dat | sed 's/FCI/1 /g' > det_xe_FCI.dat
grep 'FCI' det_xe.dat | sed 's/FCI/1E8 /g' >> det_xe_FCI.dat
#grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_CI.dat
#grep 'hfCIS ' det_xe.dat | sed 's/hfCIS /1 /g' >> det_xe_CI.dat
grep 'hfCISD ' det_xe.dat | sed 's/hfCISD /2 /g' > det_xe_CI.dat
grep 'hfCISDT ' det_xe.dat | sed 's/hfCISDT /3 /g' >> det_xe_CI.dat
grep 'hfCISDTQ' det_xe.dat | sed 's/hfCISDTQ/4 /g' >> det_xe_CI.dat
#grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_ooCI.dat
#grep 'ooCIS ' det_xe.dat | sed 's/ooCIS /1 /g' >> det_xe_ooCI.dat
grep 'ooCISD ' det_xe.dat | sed 's/ooCISD /2 /g' > det_xe_ooCI.dat
grep 'ooCISDT' det_xe.dat | sed 's/ooCISDT/3 /g' >> det_xe_ooCI.dat
grep 'hfCIs0 ' det_xe.dat | sed 's/hfCIs0/0 /g' > det_xe_CIs.dat
grep 'hfCIs2 ' det_xe.dat | sed 's/hfCIs2/2 /g' >> det_xe_CIs.dat
grep 'hfCIs4 ' det_xe.dat | sed 's/hfCIs4/4 /g' >> det_xe_CIs.dat
grep 'ooCIs0 ' det_xe.dat | sed 's/ooCIs0/0 /g' > det_xe_ooCIs.dat
grep 'ooCIs2 ' det_xe.dat | sed 's/ooCIs2/2 /g' >> det_xe_ooCIs.dat
#grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_CIo.dat
grep 'hfCIo1 ' det_xe.dat | sed 's/hfCIo1/1 /g' > det_xe_CIo.dat
grep 'hfCIo1.5' det_xe.dat | sed 's/hfCIo1.5/1.5 /g' >> det_xe_CIo.dat
grep 'hfCIo2 ' det_xe.dat | sed 's/hfCIo2/2 /g' >> det_xe_CIo.dat
grep 'hfCIo2.5' det_xe.dat | sed 's/hfCIo2.5/2.5 /g' >> det_xe_CIo.dat
grep 'hfCIo3 ' det_xe.dat | sed 's/hfCIo3/3 /g' >> det_xe_CIo.dat
grep 'hfCIo3.5' det_xe.dat | sed 's/hfCIo3.5/3.5 /g' >> det_xe_CIo.dat
grep 'hfCIo4 ' det_xe.dat | sed 's/hfCIo4/4 /g' >> det_xe_CIo.dat
#grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_ooCIo.dat
grep 'ooCIo1 ' det_xe.dat | sed 's/ooCIo1/1 /g' > det_xe_ooCIo.dat
grep 'ooCIo1.5' det_xe.dat | sed 's/ooCIo1.5/1.5 /g' >> det_xe_ooCIo.dat
grep 'ooCIo2 ' det_xe.dat | sed 's/ooCIo2/2 /g' >> det_xe_ooCIo.dat
grep 'ooCIo2.5' det_xe.dat | sed 's/ooCIo2.5/2.5 /g' >> det_xe_ooCIo.dat
grep 'ooCIo3 ' det_xe.dat | sed 's/ooCIo3/3 /g' >> det_xe_ooCIo.dat
#######################################################################
#epspdf fit_eq.eps
#okular fit_eq.pdf

View File

@ -31,3 +31,43 @@
3.4 -2.06260088
3.5 -2.05451700
3.6 -2.04724897
3.7 -2.04075589
3.8 -2.03498844
3.9 -2.02989402
4.0 -2.02541500
4.1 -2.02149231
4.2 -2.01807018
4.3 -2.01509383
4.4 -2.01251351
4.5 -2.01028245
4.6 -2.00835935
4.7 -2.00670418
4.8 -2.00528242
4.9 -2.00406073
5.0 -2.00301776
5.1 -2.00213173
5.2 -2.00137685
5.3 -2.00073304
5.4 -2.00018601
5.5 -1.99972098
5.6 -1.99932659
5.7 -1.99899193
5.8 -1.99870814
5.9 -1.99846786
6.0 -1.99826413
6.1 -1.99809160
6.2 -1.99794551
6.3 -1.99782184
6.4 -1.99771717
6.5 -1.99762841
6.6 -1.99755320
6.7 -1.99748942
6.8 -1.99743528
6.9 -1.99738933
7.0 -1.99735029
7.5 -1.99722636
8.0 -1.99716905
8.5 -1.99714119
9.0 -1.99712732
9.5 -1.99712026
10.0 -1.99711687

73
H4_cc-pvdz/pes_ooCIs4.dat Normal file
View File

@ -0,0 +1,73 @@
1.0 -1.89416406
1.1 -2.03006332
1.2 -2.12225090
1.3 -2.18372568
1.35 -2.20580385
1.4 -2.22325188
1.45 -2.23675483
1.5 -2.24689119
1.55 -2.25415055
1.6 -2.25895296
1.65 -2.26165601
1.7 -2.26256902
1.75 -2.26195628
1.8 -2.26004722
1.85 -2.25704006
1.9 -2.25310597
2.0 -2.24303439
2.1 -2.23080806
2.2 -2.21716895
2.3 -2.20268609
2.4 -2.18779608
2.5 -2.17283680
2.6 -2.15806971
2.7 -2.14369739
2.8 -2.12987642
2.9 -2.11672522
3.0 -2.10433183
3.1 -2.09275639
3.2 -2.08203592
3.3 -2.07218609
3.4 -2.06320395
3.5 -2.05507098
3.6 -2.04775524
3.7 -2.04121517
3.8 -2.03540184
3.9 -2.03026136
4.0 -2.02573774
4.1 -2.02177422
4.2 -2.01831523
4.3 -2.01530714
4.4 -2.01269952
4.5 -2.01044552
4.6 -2.00850205
4.7 -2.00683018
4.8 -2.00539489
4.9 -2.00416498
5.0 -2.00311285
5.1 -2.00221416
5.2 -2.00144765
5.3 -2.00079468
5.4 -2.00023911
5.5 -1.99976688
5.6 -1.99936588
5.7 -1.99902562
5.8 -1.99873710
5.9 -1.99849257
6.0 -1.99828540
6.1 -1.99810994
6.2 -1.99796134
6.3 -1.99783549
6.4 -1.99772890
6.5 -1.99763859
6.6 -1.99756205
6.7 -1.99749714
6.8 -1.99744206
6.9 -1.99739529
7.0 -1.99735552
7.5 -1.99722213
8.0 -1.99716834
8.5 -1.99714166
9.0 -1.99712801
9.5 -1.99711922
10.0 -1.99711666

View File

@ -49,10 +49,10 @@
5.2 -3.93653352
5.3 -3.93724673
5.4 -3.93790329
#5.5 -3.87997693
5.5 -3.93851531
#5.6 -3.97676976
#5.7 -3.89997467
#5.8 -3.94015306
5.7 -3.93961604
5.8 -3.94012202
#5.9 -3.94038278
6.0 -3.94106391
#6.1 -3.93960957
@ -61,9 +61,9 @@
#6.4 -4.03895963
#6.5 -3.93781661
6.6 -3.94353415
#6.7 -3.88251916
6.7 -3.94390984
#6.8 -3.94486071
#6.9 -3.94131095
6.9 -3.94462960
#7.0 -3.94076026
7.5 -3.94652098
8.0 -3.94774171

View File

@ -44,3 +44,5 @@
4.7 -3.79187796
4.8 -3.77798402
4.9 -3.76474094
5.0 -4.00645858
5.1 -4.00462218

View File

@ -0,0 +1 @@
1.0 -3.44375043

View File

@ -0,0 +1,4 @@
1.0 -3.44264191
1.1 -3.80394822
1.2 -4.05321075
1.3 -4.22379281

View File

@ -1,12 +1,6 @@
0.5 -99.51299514
0.55 -99.76770787
0.6 -99.93901362
0.65 -100.05395002
0.7 -100.12857177
0.75 -100.17623091
0.8 -100.20523177
0.85 -100.22126897
0.9 -100.22645291
0.85 -100.22127021
0.9 -100.22731775
0.95 -100.22684695
1.0 -100.22189555
1.05 -100.21391078

View File

@ -40,31 +40,31 @@ set ylabel 'Energy (Hartree)'
# check if file exists
file_exists(file) = int(system("[ -f '".file."' ] && echo '1' || echo '0'"))
D0=0.0;a0=1.0;xe0=1.0;C0=0.0
D1=0.0;a1=1.0;xe1=1.0;C1=0.0
D2=0.0;a2=1.0;xe2=1.0;C2=0.0
D3=0.0;a3=1.0;xe3=1.0;C3=0.0
D4=0.0;a4=1.0;xe4=1.0;C4=0.0
D5=0.0;a5=1.0;xe5=1.0;C5=0.0
D6=0.0;a6=1.0;xe6=1.0;C6=0.0
D7=0.0;a7=1.0;xe7=1.0;C7=0.0
D8=0.0;a8=1.0;xe8=1.0;C8=0.0
D9=0.0;a9=1.0;xe9=1.0;C9=0.0
D10=0.0;a10=1.0;xe10=1.0;C10=0.0
D11=0.0;a11=1.0;xe11=1.0;C11=0.0
D12=0.0;a12=1.0;xe12=1.0;C12=0.0
D13=0.0;a13=1.0;xe13=1.0;C13=0.0
D14=0.0;a14=1.0;xe14=1.0;C14=0.0
D15=0.0;a15=1.0;xe15=1.0;C15=0.0
D16=0.0;a16=1.0;xe16=1.0;C16=0.0
D17=0.0;a17=1.0;xe17=1.0;C17=0.0
D18=0.0;a18=1.0;xe18=1.0;C18=0.0
D19=0.0;a19=1.0;xe19=1.0;C19=0.0
D20=0.0;a20=1.0;xe20=1.0;C20=0.0
D21=0.0;a21=1.0;xe21=1.0;C21=0.0
D22=0.0;a22=1.0;xe22=1.0;C22=0.0
D23=0.0;a23=1.0;xe23=1.0;C23=0.0
D24=0.0;a24=1.0;xe24=1.0;C24=0.0
D0=0.3;a0=1.0;xe0=1.0;C0=0.0
D1=0.3;a1=1.0;xe1=1.0;C1=0.0
D2=0.3;a2=1.0;xe2=1.0;C2=0.0
D3=0.3;a3=1.0;xe3=1.0;C3=0.0
D4=0.3;a4=1.0;xe4=1.0;C4=0.0
D5=0.3;a5=1.0;xe5=1.0;C5=0.0
D6=0.3;a6=1.0;xe6=1.0;C6=0.0
D7=0.3;a7=1.0;xe7=1.0;C7=0.0
D8=0.3;a8=1.0;xe8=1.0;C8=0.0
D9=0.3;a9=1.0;xe9=1.0;C9=0.0
D10=0.3;a10=1.0;xe10=1.0;C10=0.0
D11=0.3;a11=1.0;xe11=1.0;C11=0.0
D12=0.3;a12=1.0;xe12=1.0;C12=0.0
D13=0.3;a13=1.0;xe13=1.0;C13=0.0
D14=0.3;a14=1.0;xe14=1.0;C14=0.0
D15=0.3;a15=1.0;xe15=1.0;C15=0.0
D16=0.3;a16=1.0;xe16=1.0;C16=0.0
D17=0.3;a17=1.0;xe17=1.0;C17=0.0
D18=0.3;a18=1.0;xe18=1.0;C18=0.0
D19=0.3;a19=1.0;xe19=1.0;C19=0.0
D20=0.3;a20=1.0;xe20=1.0;C20=0.0
D21=0.3;a21=1.0;xe21=1.0;C21=0.0
D22=0.3;a22=1.0;xe22=1.0;C22=0.0
D23=0.3;a23=1.0;xe23=1.0;C23=0.0
D24=0.3;a24=1.0;xe24=1.0;C24=0.0
# function for fitting

View File

@ -154,7 +154,7 @@ if [[ $(grep 'a24 ' a0.dat) ]]; then grep '3 ' det_CIo.dat | cut -c 4- >> det.
paste det.dat a.dat D.dat > det_aD.dat
grep 'FCI' det_aD.dat | sed 's/FCI/1 /g' > det_aD_FCI.dat
grep 'FCI' det_aD.dat | sed 's/FCI/1.0E9 /g' >> det_aD_FCI.dat
grep 'FCI' det_aD.dat | sed 's/FCI/1.0E11/g' >> det_aD_FCI.dat
grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_CI.dat
grep 'hfCIS ' det_aD.dat | sed 's/hfCIS /1 /g' >> det_aD_CI.dat
@ -253,7 +253,7 @@ grep 'xe24 ' xe0.dat | sed 's/xe24/ooCIo3 /g' >> xe.dat
paste det.dat xe.dat D.dat > det_xe.dat
grep 'FCI' det_xe.dat | sed 's/FCI/1 /g' > det_xe_FCI.dat
grep 'FCI' det_xe.dat | sed 's/FCI/1E9 /g' >> det_xe_FCI.dat
grep 'FCI' det_xe.dat | sed 's/FCI/1E11/g' >> det_xe_FCI.dat
grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_CI.dat
grep 'hfCIS ' det_xe.dat | sed 's/hfCIS /1 /g' >> det_xe_CI.dat

View File

@ -2,17 +2,17 @@
1.6 -77.29332495
1.7 -77.53168130
1.8 -77.70802692
1.9 -77.83677134
2.0 -77.92890141
2.1 -77.99223605
2.2 -78.03465007
1.9 -77.84159820
2.0 -77.93610860
2.1 -78.00252424
2.2 -78.04732377
2.3 -78.07553049
2.4 -78.07339214
2.5 -78.07667227
2.4 -78.09104792
2.5 -78.09691316
2.52053 -78.09714826
2.6 -78.07260476
2.6 -78.09548975
2.7 -78.08861735
2.8 -78.04937855
2.8 -78.07772924
2.9 -78.06394457
3.0 -78.04814002
3.1 -78.03100749
@ -25,4 +25,22 @@
3.8 -77.90893039
3.9 -77.89353420
4.0 -77.87879484
4.2 -77.64245355
4.2 -77.85124327
4.4 -77.82613461
4.6 -77.80327734
4.8 -77.78246546
5.0 -77.76349711
5.2 -77.74618578
5.4 -77.73036272
5.6 -77.71587711
5.8 -77.70259501
6.0 -77.69039650
6.2 -77.76688989
6.4 -77.76929692
6.6 -77.77128725
6.8 -77.77294261
7.0 -77.77432459
7.2 -77.77548143
7.4 -77.77645064
7.6 -77.77726332
7.8 -77.77794502

View File

@ -1,3 +1,27 @@
1.5 -77.04903276
1.6 -77.37200507
1.7 -77.61389011
1.8 -77.79366129
1.9 -77.92573191
2.0 -78.02113169
2.1 -78.08832259
2.2 -78.13390033
2.3 -78.16289654
2.4 -78.17926547
2.5 -78.18608883
2.52053 -78.18653736
2.6 -78.18576775
2.7 -78.18017151
2.8 -78.17075380
2.9 -78.15864566
3.0 -78.14472842
3.1 -78.12968920
3.2 -78.11406141
3.3 -78.09825578
3.4 -78.08258197
3.5 -78.06727022
3.6 -78.05248331
3.7 -78.03834858
3.8 -78.02492837
3.9 -78.01226506
4.0 -78.00037586