#coding=utf-8import re
# 请在此处填写代码#********** Begin **********#
f =open('src/step1/constants.txt','r')
constants ={
}
count =0for i in f:
count=count+1if count>2:
j=re.split('\s{2,}',i)
constants[j[0]]=float(j[1])print(constants)#********** End **********#
字符串操作 - 物质名称缩写
#coding=utf-8import re
defread_densities1(filename):# 请在此处填写代码# ********** Begin **********#
f=open(filename,'r')
dirt1={
}for i in f:
list1=re.split('\s{2,}',i.rstrip('\n'),2,re.U)
dirt1[list1[0]]=eval(list1[1])return dirt1
# ********** End **********#defread_densities2(filename):# 请在此处填写代码# ********** Begin **********#
f=open(filename,'r')
dirt1={
}for i in f:
list1=i.rstrip('\n').split()
str1=' '.join(list1[0:len(list1)-1])
dirt1[str1]=eval(list1[-1])return dirt1
# ********** End **********#deftest_densities(filename):# 请在此处填写代码# ********** Begin **********#if read_densities2(filename)==read_densities1(filename):returnTrueelse:returnFalse# ********** End **********#
结构文件解析与绘图 - 解释程序的输出
#coding=utf-8import warnings
warnings.filterwarnings("ignore")import matplotlib
matplotlib.use('Agg')import matplotlib.pyplot as plt
# 请在此处填写代码#********** Begin **********#
f =open('src/step3/insum.txt','r')
epsilon=[]
n=[]for i in f:
i=i.rstrip('\n')
epsilon.append(eval(i[9:15]))
n.append(eval(i[i.find('=')+1:len(i)]))
plt.semilogy(n, epsilon)#********** End **********#
plt.savefig('src/step3/output/result.jpg')