B }N_$@s`ddlmZddlTddlZddlmZmZddlm Z ddl m Z e j Z Gdd d e j ZdS) )print_function)*N)platformexit)path)thermocsZeZdZdZfddZddd Zd d Zd d ZddZddZ ddZ ddZ Z S)cubicz Interface to cubic cstt|t|j|dd|_t|j|dd|_t|j|dd|_t|j|dd|_ t|j|dd|_ dS) z= Initialize cubic specific function pointers eoslibinitZ init_cubicZthermopack_getkijZthermopack_setkijandjiZthermopack_gethvparamZthermopack_sethvparamN) superr __init__getattrtpget_export_nameeoslibinit_init_cubic s_get_kij s_set_kijs_get_hv_params_set_hv_param)self) __class__../pycThermopack/pyctp/cubic.pyr s zcubic.__init__vdWClassicDefaultFc Cst|d}tt|}t|d} tt|} t|d} tt|} t|d} tt|}t|d}tt|}|rtd}ntd}ttttttttttttg |j_d|j_|| || | |t |||| | | t t| dt| d|_ dS)aInitialize cubic model in thermopack Args: comps (str): Comma separated list of component names eos (str): Equation of state (SRK, PR, ...) mixing (str, optional): Mixture model. Defaults to "vdW". alpha (str, optional): Alpha model. Defaults to "Classic". parameter_reference (str, optional): Which parameters to use?. Defaults to "Default". asciirrN ,) c_char_pencode c_len_typelenc_intPOINTERrargtypesrestypebyrefmaxsplitnc)rcompseosmixingalphaparameter_referenceZ volume_shifteos_ceos_lenmixing_c mixing_lenalpha_c alpha_len comp_string_ccomp_string_len ref_string_cref_string_lenZ vol_shift_crrrinit)sJ        z cubic.initcCsZt|}t|}td}ttttttg|j_d|j_|t|t|t||jS)zGet attractive energy interaction parameter Args: c1 (int): Component one c2 (int): Component two Returns: kij (float): i-j interaction parameter gN)r$c_doubler%rr&r'r(value)rc1c2c1_cc2_ckij_crrrget_kij_s   z cubic.get_kijcCsXt|}t|}t|}ttttttg|j_d|j_|t|t|t|dS)zSet attractive energy interaction parameter Args: c1 (int): Component one c2 (int): Component two kij (float): i-j interaction parameter N)r$r<r%rr&r'r()rr>r?Zkijr@rArBrrrset_kijxs z cubic.set_kijcCsdS)zGet co-volume interaction Args: c1 (int): Component one c2 (int): Component two Returns: lij (float): i-j interaction parameter g?r)rr>r?rrrget_lijs z cubic.get_lijcCs tddS)zSet co-volume interaction Args: c1 (int): Component one c2 (int): Component two lij ([type]): [description] z Setting lijN)print)rr>r?Zlijrrrset_lijsz cubic.set_lijc Cst|}t|}td}td}td}td}td} td} td} td} ttttttttttttttttttttg |j_d|j_|t|t|t|t|t|t|t| t| t| t|  |j|j|j|j| j| j| j| jfS)aGet Huron-Vidal parameters Args: c1 (int): Component one c2 (int): Component two Returns: alpha_ij (float): alpha param i-j alpha_ji (float): alpha param j-i a_ij (float): a param i-j a_ji (float): a param j-i b_ij (float): b param i-j b_ji (float): b param j-i c_ij (float): c param i-j c_ji (float): c param j-i gN)r$r<r%rr&r'r(r=) rr>r?r@rA alpha_ij_c alpha_ji_ca_ij_ca_ji_cb_ij_cb_ji_cc_ij_cc_ji_crrr get_hv_params@  zcubic.get_hv_paramc  Cst|} t|} t|} t|}t|}t|}t|}t|}t| }t| }ttttttttttttttttttttg |j_d|j_|t| t| t| t|t|t|t|t|t|t| dS)aSet Huron-Vidal parameters Args: c1 (int): Component one c2 (int): Component two alpha_ij (float): alpha param i-j alpha_ji (float): alpha param j-i a_ij (float): a param i-j a_ji (float): a param j-i b_ij (float): b param i-j b_ji (float): b param j-i c_ij (float): c param i-j c_ji (float): c param j-i N)r$r<r%rr&r'r()rr>r?Zalpha_ijZalpha_jiZa_ijZa_jiZb_ijZb_jiZc_ijZc_jir@rArHrIrJrKrLrMrNrOrrr set_hv_params> zcubic.set_hv_param)rrrF) __name__ __module__ __qualname____doc__r r;rCrDrErGrPrQ __classcell__rr)rrr s  5  5r ) __future__rctypesnumpynpsysrrosrr rr" thermopackr rrrrs