def compare_exp_mod(ntot_psi=21):
import numpy as np
import matplotlib.pyplot as plt
from pepshkl import reader4
from RS import sff_converter
from rs_exp import read_IGSF
from MP.lib import mpl_lib
from MP.mat.mech import FlowCurve as FC
wf = mpl_lib.wide_fig
sff_converter.main(fn='temp.sff',difile=None,itab=True,
ieps0=4, fn_str='STR_STR.OUT')
SF, dum = read_IGSF(fn='temp.sff',fn_str='STR_STR.OUT')
SF.reduce_psi(bounds=[0.,0.5],ntot_psi=ntot_psi)
SF.mask_vol()
SF.plot(ylim=[-2,2])
## biaxial
flow = FC()
exx = open('YJ_Bsteel_BB.sff','r').read().split('\n')[4].split()[1:]
exx = np.array(exx,dtype='float'); eyy = exx[::]; ezz = -(exx+eyy)
flow.epsilon_vm = abs(ezz)
flow.nstp = len(ezz)
SF, dum = read_IGSF(fn='YJ_Bsteel_BB.sff',fc=flow)
SF.plot(mxnphi=3,ylim=[-2,2])
def compare_exp_mod(ntot_psi=21):
import numpy as np
import matplotlib.pyplot as plt
from pepshkl import reader4
from RS import sff_converter
from rs_exp import read_IGSF
from MP.lib import mpl_lib
from MP.mat.mech import FlowCurve as FC
wf = mpl_lib.wide_fig
sff_converter.main(fn='temp.sff',difile=None,itab=True,
ieps0=4, fn_str='STR_STR.OUT')
SF, dum = read_IGSF(fn='temp.sff',fn_str='STR_STR.OUT')
SF.reduce_psi(bounds=[0.,0.5],ntot_psi=ntot_psi)
SF.mask_vol()
SF.plot(ylim=[-2,2])
## biaxial
flow = FC()
exx = open('YJ_Bsteel_BB.sff','r').read().split('\n')[4].split()[1:]
exx = np.array(exx,dtype='float'); eyy = exx[::]; ezz = -(exx+eyy)
flow.epsilon_vm = abs(ezz)
flow.nstp = len(ezz)
SF, dum = read_IGSF(fn='YJ_Bsteel_BB.sff',fc=flow)
SF.plot(mxnphi=3,ylim=[-2,2])
