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])