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