ms=ms[sid])
#plt.title('$\gamma$ nodes')
#plt.xlim(-4,4)
#plt.ylim(-4,4)
#plt.legend(loc=3)
parser = argparse.ArgumentParser()
parser.add_argument('method', choices=('gamma', 'pert', 'm', 'k', 'aux',
'int', 'ext'))
io_util.add_arguments(parser, ('boundary', 'N'))
args = parser.parse_args()
# Dummy problem just to get the solver to run
problem = Smooth_Sine()
boundary = problems.boundary.boundaries[args.boundary](problem.R)
problem.boundary = boundary
options = {'problem': problem, 'boundary': boundary, 'N': args.N,
'scheme_order': 4}
solver = ps.ps.PizzaSolver(options)
if args.method == 'gamma':
plot_aux(options={'color': 'black'})
plot_gamma()
plt.axes().set_aspect('equal', 'datalim')
plt.axis('off')
plt.savefig('/Users/sam/Google Drive/research/writeup_v4/union_gamma.pdf',
bbox_inches='tight', transparent=True)
gamma = solver.all_gamma[sid]
x_data, y_data = nodes_to_plottable(gamma)
label_text = '$\gamma_{}$'.format(sid)
plt.plot(x_data, y_data, markers[sid], label=label_text,
mfc='none', mec=colors[sid], mew=1)
#plt.title('$\gamma$ nodes')
plt.xlim(-4,4)
plt.ylim(-4,4)
plt.legend(loc=3)
plt.show()
# Dummy problem just to get the solver to run
problem = Smooth_Sine()
dms = 6
data = [
[
np.pi/6,
{'boundary': 'inner-sine', 'color': 'blue', 'marker': 'o', 'ms': dms},
{'boundary': 'cubic', 'color': 'green', 'marker': '^', 'ms': dms},
],
[
np.pi/2,
{'boundary': 'outer-sine', 'color': 'red', 'marker': 'd', 'ms': dms},
{'boundary': 'sine7', 'color': 'purple', 'marker': '*', 'ms': dms+2.5},
]
]