def plot_reval2(epsimp1, epsfc1, epsimp2, phi, varphi):
e1data = np.linspace(-3.0, 1.0, 30)
e2data = np.linspace(-3.0, 1.0, 30)
e1data, e2data = np.meshgrid(e1data, e2data)
lambdafun = np.vectorize(
lambda e1, e2: reval2(e1, e2, epsimp1, epsfc1, epsimp2, phi, varphi))
ldata = lambdafun(e1data, e2data)
p3d.plot3d(e1data,
e2data,
ldata,
xlabel=r'$\epsilon_1^i$',
ylabel=r'$\epsilon_2^i$',
zlabel=r'$\lambda^i_2$',
export="lambda2.pdf")
def plot_reval2(epsimp1, epsfc1, epsimp2, phi, varphi):
e1data = np.linspace(-3.0, 1.0, 30)
e2data = np.linspace(-3.0, 1.0, 30)
e1data, e2data = np.meshgrid(e1data, e2data)
lambdafun = np.vectorize(lambda e1, e2:
reval2(e1, e2, epsimp1, epsfc1, epsimp2,
phi, varphi))
ldata = lambdafun(e1data, e2data)
p3d.plot3d(e1data, e2data, ldata,
xlabel=r'$\epsilon_1^i$',
ylabel=r'$\epsilon_2^i$',
zlabel=r'$\lambda^i_2$'
, export="lambda2.pdf"
)
# shrinks, positive enlarges in the corresponding direction.
trim = [-0.1, 0.1, 0.1, 0.0]
# Package the plot options
kwargs = {"xlabel": r'$p_A$',
"ylabel": r'$p_B$',
"elev": elev, "azim": azim, "trim": trim,
"zform": "5.2f", "zlabelshift": 7,
"labeldist": 3}
# ************************************************************************
# Do plots
p3d.plot3d(xdata, ydata, OmegaAdata,
zlabel=r'$\Omega_A$',
plotlabel=r'(a)',
export='Omega-A.pdf',
**kwargs)
p3d.plot3d(xdata, ydata, OmegaBdata,
zlabel=r'$\Omega_B$',
plotlabel=r'(b)',
export='Omega-B.pdf',
**kwargs)
p3d.plot3d(xdata, ydata, OmegaCdata,
zlabel=r'$\Omega_A + \Omega_B$',
plotlabel=r'(c)',
export='Omega-C.pdf',
**kwargs)
"ylabel": r'$p_B$',
"elev": elev,
"azim": azim,
"trim": trim,
"zform": "5.2f",
"zlabelshift": 7,
"labeldist": 3
}
# ************************************************************************
# Do plots
p3d.plot3d(xdata,
ydata,
OmegaAdata,
zlabel=r'$\Omega_A$',
plotlabel=r'(a)',
export='Omega-A.pdf',
**kwargs)
p3d.plot3d(xdata,
ydata,
OmegaBdata,
zlabel=r'$\Omega_B$',
plotlabel=r'(b)',
export='Omega-B.pdf',
**kwargs)
p3d.plot3d(xdata,
ydata,
OmegaCdata,
# shrinks, positive enlarges in the corresponding direction.
trim = [-0.1, 0.1, 0.1, 0.0]
# Package the plot options
kwargs = {"xlabel": r'$p_A$',
"ylabel": r'$p_B$',
"elev": elev, "azim": azim, "trim": trim,
"zform": "5.2f", "zlabelshift": 7,
"labeldist": 3}
# ************************************************************************
# Do plots
p3d.plot3d(xdata, ydata, OmegaAdata,
zlabel=r'$\Omega_A$',
plotlabel=r'(a)',
export='Omega-simple-A.pdf',
**kwargs)
p3d.plot3d(xdata, ydata, OmegaBdata,
zlabel=r'$\Omega_A$',
plotlabel=r'(b)',
export='Omega-simple-B.pdf',
**kwargs)
p3d.plot3d(xdata, ydata, OmegaCdata,
zlabel=r'$\Omega_A$',
plotlabel=r'(c)',
export='Omega-simple-C.pdf',
**kwargs)
"ylabel": r'$p_B$',
"elev": elev,
"azim": azim,
"trim": trim,
"zform": "5.2f",
"zlabelshift": 7,
"labeldist": 3
}
# ************************************************************************
# Do plots
p3d.plot3d(xdata,
ydata,
OmegaAdata,
zlabel=r'$\Omega_A$',
plotlabel=r'(a)',
export='Omega-simple-A.pdf',
**kwargs)
p3d.plot3d(xdata,
ydata,
OmegaBdata,
zlabel=r'$\Omega_A$',
plotlabel=r'(b)',
export='Omega-simple-B.pdf',
**kwargs)
p3d.plot3d(xdata,
ydata,
OmegaCdata,
# Trim output so that the graph is neatly contained in the figure
# area. Vector of left, bottom, right, and top trim values. Negative
# shrinks, positive enlarges in the corresponding direction.
trim = [-0.1, 0.1, 0.1, 0.0]
# Package the plot options
kwargs = {"elev": elev, "azim": azim, "trim": trim,
"zform": "5.0f", "zlabelshift": 3,
"labeldist": 3}
# ************************************************************************
# Do plots
p3d.plot3d(x1data, y1data, c2Apdf,
xlabel=r'$C_2/L_2$',
ylabel=r'$p_A$',
zlabel=r'$C_2/L_2\; \mathrm{pdf.}$',
plotlabel="(a)",
export='caprat2-pdf-A.pdf',
**kwargs)
p3d.plot3d(x2data, y2data, C2Apdf,
xlabel=r'$C_2$',
ylabel=r'$p_A$',
zlabel=r'$C_2\; \mathrm{pdf.}$',
plotlabel="(b)",
export='caplev2-pdf-A.pdf',
**kwargs)
