需要从Basis Set Exchange制作自己的基组,以下简单记录基组制作过程。
# 1. 查找数据集冲所有的几何结构文件
find . -type f -iname "geometry.in"
# 2. 打印文件的最后一列
awk '{print $NF}' file
# 3. 结合以上两个命令进行循环
for i in `find . -type f -iname "geometry.in"`;do awk '{print $NF}' $i; done
# 4. 将结果保存到文件并进行去重处理
for i in `find . -type f -iname "geometry.in"`; \
do awk '{print $NF}' $i >> all_elements; done; \
sort all_elements | uniq > uniq_elements
统计多个数据集中唯一的元素
#!/bin/sh
# Author: LittleTan24
# Created Time : 2021年12月21日 星期二 00时28分15秒
# File Name: uniq_elements_subsets.sh
# Description:
echo "" > "$1".all_elements
echo "" > "$1".all_uniq_elements
while read -r line
do
if [ -d "$line" ]
then
echo "Entering $line"
cd "$line" || exit
cat uniq_elements >> ../"$1".all_elements
cd ../ || exit
echo "Leaving $line"
fi
done < "$1"
sort "$1".all_elements | uniq > "$1".all_uniq_elements