#!/bin/bash # compute Egap = E(n+1) - E(n) # use: egap.sh n file.r declare -i nocc=$1 shift declare -i nl=nocc/5+1 #echo "nl=" $nl declare -i nfrac=nocc-5*\(nocc/5\) #echo "nfrac=" $nfrac grep -h -A$nl '<eigenvalues ' ${*}| \ awk -v nl=$nl -v nfrac=$nfrac \ ' NR%(nl+2)==nl {e[0] = $5;} \ NR%(nl+2)==(nl+1) \ { e[1]=$1; e[2]=$2; e[3]=$3; e[4]=$4; e[5]=$5; \ e_homo = e[nfrac]; e_lumo = e[nfrac+1]; \ print "E_HOMO=",e_homo, "E_LUMO=",e_lumo, "Eg=",e_lumo-e_homo}' -
