Python SW.fhsの例

コード例 #1

def SW_Fexc(T, n):
    #return SW.fhs(T,n)
    #return SW.fdisp(T,n)
    return SW.fdisp(T, n) + SW.fhs(T, n)

コード例 #2

def SW_Finfinity(n):
    T = 1.0e8
    return SW.fhs(T, n) / T

コード例 #3

colors = ['ro', 'go', 'bo', 'rx', 'gx', 'bx']
plt.figure()
size = 15.0
for i in range(0, len(allFinfinity)):
    plt.plot(allFinfinity[i][0],
             allFinfinity[i][1],
             colors[i % len(colors)],
             label="L=%f" % allCellLength[i])
plt.xlabel("filling fraction")
plt.ylabel("dataF-saftF")
plt.legend(loc='best')
plt.savefig("deltaF.pdf")
plt.close()

colors = ['ro', 'go', 'bo', 'rx', 'gx', 'bx']
plt.figure()
size = 15.0
for i in range(0, len(allFinfinity)):  #ms=size*(0.7)**i
    plt.plot(allFinfinityPercent[i][0],
             allFinfinityPercent[i][1],
             colors[i % len(colors)],
             label="L=%f" % allCellLength[i])
plt.xlabel("filling fraction")
plt.ylabel("(dataF-saftF)/dataF*100.0")
plt.legend(loc='best')
plt.savefig("deltaFpercent.pdf")
plt.close()

print SW.fhs(0, 0.1)
print SW.fhs(0, 0.2)

コード例 #4

ファイル: comparisonWithSaftV4.py プロジェクト: droundy/deft

def SW_Finfinity(n):
	T=1.0e8
	return SW.fhs(T,n)/T

コード例 #5

ファイル: comparisonWithSaftV4.py プロジェクト: droundy/deft

def SW_Fexc(T,n):
	#return SW.fhs(T,n)
	#return SW.fdisp(T,n)
	return SW.fdisp(T,n)+SW.fhs(T,n)

コード例 #6

ファイル: comparisonWithSaft.py プロジェクト: droundy/deft

			allFinfinityPercent.append([[],[]])
			for i in range(0,len(dataFreeEnergy)):
				allFinfinityPercent[-1][0].append(ffData[i])
				allFinfinityPercent[-1][1].append(dataFreeEnergy[i])
			
###############--------------all ww=1.5  T=infinity  dataF-saftF------##
colors=['ro','go','bo','rx','gx','bx']
plt.figure()
size=15.0
for i in range(0,len(allFinfinity)):
	plt.plot(allFinfinity[i][0],allFinfinity[i][1],colors[i%len(colors)],label="L=%f"%allCellLength[i])
plt.xlabel("filling fraction")
plt.ylabel("dataF-saftF")
plt.legend(loc='best')
plt.savefig("deltaF.pdf")
plt.close()

colors=['ro','go','bo','rx','gx','bx']
plt.figure()
size=15.0
for i in range(0,len(allFinfinity)):#ms=size*(0.7)**i
	plt.plot(allFinfinityPercent[i][0],allFinfinityPercent[i][1],colors[i%len(colors)],label="L=%f"%allCellLength[i])
plt.xlabel("filling fraction")
plt.ylabel("(dataF-saftF)/dataF*100.0")
plt.legend(loc='best')
plt.savefig("deltaFpercent.pdf")
plt.close()

print SW.fhs(0,0.1)
print SW.fhs(0,0.2)

コード例 #7

			allFinfinityPercent.append([[], []])
			for i in range(0, len(dataFreeEnergy)):
				allFinfinityPercent[-1][0].append(ffData[i])
				allFinfinityPercent[-1][1].append(dataFreeEnergy[i])
			
###############--------------all ww=1.5  T=infinity  dataF-saftF------##
colors=['ro', 'go', 'bo', 'rx', 'gx', 'bx']
plt.figure()
size=15.0
for i in range(0, len(allFinfinity)):
	plt.plot(allFinfinity[i][0], allFinfinity[i][1], colors[i%len(colors)], label="L=%f"%allCellLength[i])
plt.xlabel("filling fraction")
plt.ylabel("dataF-saftF")
plt.legend(loc='best')
plt.savefig("deltaF.pdf")
plt.close()

colors=['ro', 'go', 'bo', 'rx', 'gx', 'bx']
plt.figure()
size=15.0
for i in range(0, len(allFinfinity)):#ms=size*(0.7)**i
	plt.plot(allFinfinityPercent[i][0], allFinfinityPercent[i][1], colors[i%len(colors)], label="L=%f"%allCellLength[i])
plt.xlabel("filling fraction")
plt.ylabel("(dataF-saftF)/dataF*100.0")
plt.legend(loc='best')
plt.savefig("deltaFpercent.pdf")
plt.close()

print(SW.fhs(0, 0.1))
print(SW.fhs(0, 0.2))