%config InlineBackend.figure_format = 'retina'

"""
Explanation: For high dpi displays.
End of explanation
"""


import matplotlib.pyplot as plt
import numpy as np
from uncertainties import unumpy as unp
import pytheos as eos

"""
Explanation: 0. General note
This example compares pressure calculated from pytheos and original publication for the gold scale by Tange 2008.
1. Global setup
End of explanation
"""


eta = np.linspace(1., 0.65, 36)
print(eta)

tange_mgo = eos.periclase.Tange2009()

tange_mgo.print_equations()

tange_mgo.print_equations()

tange_mgo.print_parameters()

v0 = 74.698

tange_mgo.three_r

v = v0 * (eta) 
temp = 3000.

p = tange_mgo.cal_p(v, temp * np.ones_like(v))

"""
Explanation: 3. Compare
End of explanation
"""


print('for T = ', temp)
for eta_i, p_i in zip(eta, p):
    print("{0: .3f} {1: .2f}".format(eta_i, p_i))

"""
Explanation: <img src='./tables/Tange_MgO.png'>
End of explanation
"""


v = tange_mgo.cal_v(p, temp * np.ones_like(p), min_strain=0.6)
print(1.-(v/v0))

"""
Explanation: Make comparison with the Vinet column in the table.
End of explanation
"""


