def credentials(cls, args):
"""Credentials handler."""
# Prep
symfony = Symfony(url=args.url, session=args.session)
engine = Engine(size=args.threads, session=args.session)
engine.start()
# Run
cls.log.info(Logs.name)
Logs(symfony, engine).run()
# Save
engine.stop()
if args.output:
print()
EOS.save(symfony, args.output)
def sources(cls, args):
"""Sources handler."""
# Prep
symfony = Symfony(url=args.url, session=args.session)
engine = Engine(size=args.threads, session=args.session)
engine.start()
# Run
cls.log.info(Info.name)
Info(symfony).run()
print()
cls.log.info(Sources.name)
Sources(symfony, engine).run()
# Save
engine.stop()
if args.output:
print()
EOS.save(symfony, args.output)
def test_EOSPuresubstancesInterface_PR1976():
eosname = "Peng and Robinson (1976)"
subs = [methane]
y = [1.0]
k = [[0]]
p = 1e5
t = 150
# test Z
eos_PSI = EOSPureSubstanceInterface(subs, eosname)
eos_old = EOS(MixtureProp(subs, [1.0]), k, eosname)
zsPSI = eos_PSI.getZfromPT(p, t)
zs_old = eos_old.getZfromPT(p, t)
np.testing.assert_allclose(np.min(zs_old), np.min(zsPSI), 1e-5)
np.testing.assert_allclose(np.max(zs_old), np.max(zsPSI), 1e-5)
# test pvp
pvp_aw = subs[0].getPvpAW(t)
pvp_old = eos_old.getPvp(t, pvp_aw)[0]
pvp_PSI = eos_PSI.getPvp(t, pvp_aw)[0]
np.testing.assert_allclose(pvp_PSI, pvp_old, 1e-5)
def scan(cls, args):
"""Scan handler."""
eos = EOS(url=args.url, output=args.output, session=args.session)
eos.run(check_only=args.check_only, threads=args.threads)
plt.plot([1,4,9,16])
fig.canvas.mpl_connect('pick_event', onpick)
"""
plt.close("IME")
UsQ, UsCH = barrios.ix[event.ind][["UsQ", "UsCH"]].values[0]
impedence_matching_exp(q, ch, UsQ, UsCH, showplot=True)
if __name__ == "__main__":
plt.ion()
plt.close('all')
# Import Polystyrene EOS data, from a HYADES EOS file, and change initial
# conditions
ch = EOS("../data/eos_32.dat")
ch._rho0 = 1.05
ch._P0 = 0
# Import Quartz EOS data and use an analytic method bases on Us vs Up
# to calculate the Hugoniot
q = EOS("../data/eos_24.dat",
hugoniot_method="analytic",
upmin=0,
upmax=29,
UsvUp=quartz_UsvUp(model="knudson"))
# Play around with different quarts release models
q.hugoniot.custom_release_model = types.MethodType(quartz_ode_mglr,
q.hugoniot)
q.hugoniot.set_release_model("custom")
def get_context_data(self, **kwargs):
kwargs.update({
'eos': EOS(),
'surveys': self.get_surveys()
})
return super().get_context_data(**kwargs)
from eos import EOS, quartz_UsvUp, quartz_mglr, quartz_ode_mglr
from impediance_match import impedence_matching_exp, monte_carlo_error_prop
from util import MonteCarloVariable
def onpick(event):
plt.close("IME")
UsQ, UsGDP = knudson.ix[event.ind][["UsQ", "UsGDP"]].values[0]
impedence_matching_exp(q, GDP, UsQ, UsGDP, showplot=True)
if __name__ == "__main__":
plt.ion()
plt.close('all')
GDP = EOS("../data/eos_32.dat")
GDP._rho0 = 1.035
GDP._P0 = 0
q = EOS("../data/eos_24.dat",
hugoniot_method="analytic",
upmin=0,
upmax=50,
UsvUp=quartz_UsvUp(model="knudson")
)
# overload the quarts hugoniot with the MGLR model
q.hugoniot.custom_release_model = types.MethodType(quartz_ode_mglr, q.hugoniot)
q.hugoniot.set_release_model("custom")
knudson = pd.read_csv("../data/knudson_test.csv", index_col=False)
UsGDP = knudson.UsGDP.values