def read(fn,iopt=0,isort=False):
datl = open(fn,'r').readlines()
i=13
psi1 = map(float, datl[i].split()[2:])
sin2psi1 = map(float, datl[i+1].split()[2:])
dspacing1 = map(float, datl[i+2].split()[1:])
ehkl1 = map(float, datl[i+4].split()[1:]) # mu strain * 10^3
i=25
psi2 = map(float, datl[i].split()[2:])
sin2psi2 = map(float, datl[i+1].split()[2:])
dspacing2 = map(float, datl[i+2].split()[1:])
ehkl2 = map(float, datl[i+4].split()[1:]) # mu strain * 10^3
psi = np.array(psi1 + psi2)
# sin2psi = sin2psi1+sin2psi2
sin2psi = np.sign(psi) * np.sin(psi*np.pi/180.)**2
ehkl = ehkl1+ehkl2
ehkl = np.array(ehkl) / 1e3
dspacing = dspacing1 +dspacing2
if iopt==1:
psi=psi1
sin2psi = sin2psi1
ehkl = ehkl1
dspacing=dspacing1
if isort:
sin2psi, dspacing, ehkl, psi = sort(sin2psi, dspacing, ehkl, psi)
return sin2psi, ehkl, dspacing, psi
def read(fn, iopt=0, isort=False):
datl = open(fn, 'r').readlines()
i = 13
psi1 = map(float, datl[i].split()[2:])
sin2psi1 = map(float, datl[i + 1].split()[2:])
dspacing1 = map(float, datl[i + 2].split()[1:])
ehkl1 = map(float, datl[i + 4].split()[1:]) # mu strain * 10^3
i = 25
psi2 = map(float, datl[i].split()[2:])
sin2psi2 = map(float, datl[i + 1].split()[2:])
dspacing2 = map(float, datl[i + 2].split()[1:])
ehkl2 = map(float, datl[i + 4].split()[1:]) # mu strain * 10^3
psi = np.array(psi1 + psi2)
# sin2psi = sin2psi1+sin2psi2
sin2psi = np.sign(psi) * np.sin(psi * np.pi / 180.)**2
ehkl = ehkl1 + ehkl2
ehkl = np.array(ehkl) / 1e3
dspacing = dspacing1 + dspacing2
if iopt == 1:
psi = psi1
sin2psi = sin2psi1
ehkl = ehkl1
dspacing = dspacing1
if isort:
sin2psi, dspacing, ehkl, psi = sort(sin2psi, dspacing, ehkl, psi)
return sin2psi, ehkl, dspacing, psi
def main(fn=None, ref='Bsteel_BB_opt0.txt'):
""" """
markers = ['-o', '--x', '-+', '-^', '-d', '-*']
datl = open(ref, 'r').readlines()
datl = datl[2:]
ind = [0, 2, 11, 16, 18]
fn = []
strain = []
fig = plt.figure(3)
# fig = plt.figure(3,figsize=(12,5))
ax1 = fig.add_subplot(111)
ax1.grid('on')
# ax1 = fig.add_subplot(121); ax1.grid('on')
# ax2 = fig.add_subplot(122); ax2.grid('on')
for i in xrange(len(ind)):
f = datl[ind[i]].split()[0]
epsx = float(datl[ind[i]].split()[1])
epsy = float(datl[ind[i]].split()[2])
eps = epsx + epsy
strain.append(eps)
f = f + 'Data2Phi0.txt'
fn.append(f)
for i in xrange(len(fn)):
sin2psi, ehkl, dspacing, psi = read(fn[i], iopt=1) # d: strain * 10^3
x = np.array(sin2psi).copy()
x, ehkl = sort(x, ehkl)
#x,ind = sort(x,ehkl)
#ehkl = ss(ind, ehkl)
# dspacing = ss(ind, dspacing)
ax1.plot(x,
ehkl,
markers[i],
ms=8,
color='k',
label=r'$\bar{E}^{\mathrm{eff}}=%5.2f$' % strain[i])
# ax2.plot(x, dspacing,markers[i],ms=8,
# label=r'$\bar{E}^{\mathrm{eff}}=%5.2f$'%strain[i])
ax1.legend(loc='best', fancybox=True).get_frame().set_alpha(0.5)
#ax2.legend(loc='best',fancybox=True).get_frame().set_alpha(0.5)
ax1.set_ylabel(r'$\varepsilon(hkl,\phi,\psi) \times 10^3 $',
dict(fontsize=28))
#ax2.set_ylabel(r'$d(hkl,\phi,\psi)$',
#dict(fontsize=28))
ax1.set_xlabel(r'$\sin^2{\psi}$', dict(fontsize=28))
#ax2.set_xlabel(r'$\sin^2{\psi}$',dict(fontsize=28))
fig.tight_layout()
fig.savefig('sin2psi_ex.pdf')
def main(fn=None,ref='Bsteel_BB_opt0.txt'):
""" """
markers = ['-o','--x','-+','-^','-d','-*']
datl = open(ref,'r').readlines()
datl = datl[2:]
ind = [0,2,11,16,18]
fn = []; strain = []
fig = plt.figure(3)
# fig = plt.figure(3,figsize=(12,5))
ax1 = fig.add_subplot(111); ax1.grid('on')
# ax1 = fig.add_subplot(121); ax1.grid('on')
# ax2 = fig.add_subplot(122); ax2.grid('on')
for i in range(len(ind)):
f = datl[ind[i]].split()[0]
epsx = float(datl[ind[i]].split()[1])
epsy = float(datl[ind[i]].split()[2])
eps = epsx+epsy
strain.append(eps)
f = f+'Data2Phi0.txt'
fn.append(f)
for i in range(len(fn)):
sin2psi, ehkl,dspacing,psi = read(fn[i],iopt=1) # d: strain * 10^3
x = np.array(sin2psi).copy()
x,ehkl = sort(x,ehkl)
#x,ind = sort(x,ehkl)
#ehkl = ss(ind, ehkl)
# dspacing = ss(ind, dspacing)
ax1.plot(x, ehkl,markers[i],ms=8,color='k',
label=r'$\bar{E}^{\mathrm{eff}}=%5.2f$'%strain[i])
# ax2.plot(x, dspacing,markers[i],ms=8,
# label=r'$\bar{E}^{\mathrm{eff}}=%5.2f$'%strain[i])
ax1.legend(loc='best',fancybox=True).get_frame().set_alpha(0.5)
#ax2.legend(loc='best',fancybox=True).get_frame().set_alpha(0.5)
ax1.set_ylabel(r'$\varepsilon(hkl,\phi,\psi) \times 10^3 $',
dict(fontsize=28))
#ax2.set_ylabel(r'$d(hkl,\phi,\psi)$',
#dict(fontsize=28))
ax1.set_xlabel(r'$\sin^2{\psi}$',dict(fontsize=28))
#ax2.set_xlabel(r'$\sin^2{\psi}$',dict(fontsize=28))
fig.tight_layout()
fig.savefig('sin2psi_ex.pdf')
