#!/usr/bin/env python
import sys
block='Population'
eigen='eigenvalues'
occst='occ'
virst='vir'
sepst='--'
for filename in sys.argv[1:] :
  f=open(filename,'r')
  txt=f.read().splitlines()
  f.close()
  for line in txt :
    if line.find(block) > -1 or line.find(eigen) > -1 :
      if line.find(block) > -1 :
        eocc=[]
        evir=[]
      elif line.find(occst) > -1 :
        data=line.split(sepst)[1]
        eocc=eocc+[float(i) for i in data.split()]
      elif line.find(virst) > -1 :
        data=line.split(sepst)[1]
        evir=evir+[float(i) for i in data.split()]
  print "File :",filename
  if len(eocc) > 0 : print "  HOMO :",max(eocc)
  if len(evir) > 0 : print "  LUMO :",min(evir)

9楼: Originally posted by asaka at 2011-12-09 04:13:39:
说一下我的想法啊：

首先用grep把有用的输出信息过滤一下：
一个是分割各个计算结果，免得把好几个计算的信息混在一起，
另外也是楼主要求的，只输出最后的结果，
另一个就是各个轨道能的数据。

然后把占 ...

12楼: Originally posted by 秋雨子 at 2011-12-09 21:38:52:
非常感谢您的不吝赐教！思路很清晰 ，不过我还有一点小问题，麻烦您再解答一下，我只会用一点shell，python不知道怎么用，我把您的命令粘帖到text里面，运行了一下 没有输出文件，我把后面的python写的那两行删了 ...

11楼: Originally posted by asaka at 2011-12-09 17:40:15:
shell版实在是太慢，既然要用到python，那就纯python吧，速度快多了。
[code]
#!/usr/bin/env python
import sys
block='Population'
eigen='eigenvalues'
occst='occ'
virst='vir'
sepst='--'
for file ...

#!/bin/bash
# print homo lumo energy of gaussian output
# usage: print_homo_lumo filename(s)
for file in "$@"
do
echo $file
tac $file | tr "\n" "@" | egrep -o "@ Alpha virt\. eigenvalues.{0,75}@ Alpha  occ\. eigenvalues.{0,75}@ Alpha  occ\. eigenvalues" tio2.rev| tr "@" "\n" | head -n3 | awk '{if ($2 == "virt.") print $1" LUMO "$5; if($2 == "occ."){i=NF; print $1" HOMO "$i} }'
done

18楼: Originally posted by asaka at 2011-12-09 20:33:17:
以上作者只考虑闭壳层计算的情况，你有开壳层的计算吗？那时候你需要总的HOMO,LUMO还是分别alpha以及beta的HOMO和LUMO?
他们的程序都只有一点小问题，相信LZ可以搞定。
mchen10 的程序少了一个 |，多了一个奇怪 ...

#!/bin/bash
# print homo lumo energy of gaussian output
# usage: print_homo_lumo filename(s)
for file in "$@"
do
echo $file
tac $file | tr "\n" "@" | egrep -o "@ Alpha virt\. eigenvalues.{0,75}@ Alpha  occ\. eigenvalues.{0,75}@ Alpha  occ\. eigenvalues" | tr "@" "\n" | head -n3 | awk '{if ($2 == "virt.") print $1" LUMO "$5; if($2 == "occ."){i=NF; print $1" HOMO "$i} }'
done

16楼: Originally posted by ZDBWHZ at 2011-12-10 13:31:18:
直接shell:
#!/bin/sh
# ZDB, HNU, 2011/11/15.
# ZDB, HNU, 2011/11/23 modified to be more pratical.
`ls -l *.log | awk '{print $9}' | sed 's/.log//g' > file_list`
cat file_list | while read  ...