if args.plot == 'effectivity':
title = 'Effectivity'
else:
title = 'Indicator effectivity'
# Plot the axes
s += tkz.get_2d_axes(xmin,
xmax,
ymin,
ymax,
xscale=xscale,
yscale=yscale,
xticks=xticks,
yticks=yticks,
xtick_labels=xtick_labels,
ytick_labels=ytick_labels,
tick_font='normalsize',
tick_frac=0.01,
xlabel_offset=0.085,
label_font='Large',
xlabel='Number of nodes',
ylabel_offset=0.175,
ylabel=title)
s += tkz.get_end_tikz()
fp = open(args.outfile, 'w')
fp.write(s)
fp.close()
color='gray',
line_dim='thick',
xmin=xmin,
xmax=xmax,
ymin=ymin,
ymax=ymax)
# Plot the axes
s += tkz.get_2d_axes(xmin,
xmax,
ymin,
ymax,
xscale=xscale,
yscale=yscale,
xticks=xticks,
yticks=yticks,
xtick_labels=xtick_labels,
ytick_labels=ytick_labels,
tick_font='small',
tick_frac=0.0125,
xlabel_offset=0.1,
xlabel='$\\log_{10}(\\text{Error})$',
ylabel_offset=0.065,
ylabel='Percentage')
for k, d in enumerate(data):
bars = []
ymin = delta * (len(data) - k - 1)
for i in range(d.shape[0]):
bars.append([d[i, 3] + ymin])
s += tkz.get_bar_chart(bars,
yticks = [0, 0.25, 0.5, 0.75, 1.0]
xticks = [ 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5 ]
ylabel_offset = 0.15
length = 0.15
ylegend = 0.025 + len(heuristics)*0.06
xlegend = 3.5
# Get the header info
s = tikz.get_header()
s += tikz.get_begin_tikz(xdim=2, ydim=2, xunit='in', yunit='in')
# Plot the axes
s += tikz.get_2d_axes(xmin, xmax, ymin, ymax,
tick_frac=tick_frac, ylabel_offset=ylabel_offset,
xscale=xscale, yscale=yscale,
xticks=xticks, yticks=yticks,
xlabel='$\\alpha$', ylabel='Fraction of problems')
for k in xrange(len(heuristics)):
tau, rho = get_performance_profile(r[:, k], 1.5*xmax)
s += tikz.get_2d_plot(tau, rho, xscale=xscale, yscale=yscale,
color=colors[k], line_dim='ultra thick',
xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax,
symbol=None)
# Add a label to the legend
xlegend_len = 1.5
ylegend_len = 0.06*(len(heuristics)-1)
s += '\\draw[fill=white] (%f,%f) rectangle (%f,%f);'%(
(xlegend - 0.15)*xscale, (ylegend + 0.04)*yscale,
s += r'\draw[%s, color=%s, fill=%s, fill opacity=0.3]' % (line_dim,
color, color)
s += ' (%f, %f) rectangle (%f, %f);' % (xscale * x1, yscale * y1,
xscale * x2, yscale * y2)
s += r'\draw[font=\%s] (%f, %f) node[right] {%s};' % (
font_size, xscale * (x + xlen), yscale * y, label)
# Plot the axes
s += tkz.get_2d_axes(xmin,
xmax,
ymin,
ymax,
xscale=xscale,
yscale=yscale,
xticks=xticks,
yticks=yticks,
xtick_labels=xtick_labels,
ytick_labels=ytick_labels,
tick_font='small',
tick_frac=0.0125,
xlabel_offset=0.1,
xlabel='Element error',
ylabel_offset=0.075,
ylabel='Percentage')
s += tkz.get_end_tikz()
fp = open(args.outfile, 'w')
fp.write(s)
fp.close()
color='customcolor',
line_dim='ultra thick',
xmin=xmin,
xmax=xmax,
ymin=ymin,
ymax=ymax)
# Plot the axes
s += tikz.get_2d_axes(xmin,
xmax,
ymin,
ymax,
tick_frac=tick_frac,
ylabel_offset=ylabel_offset,
label_font='LARGE',
xscale=xscale,
yscale=yscale,
xticks=xticks,
xtick_labels=xtick_labels,
yticks=yticks,
tick_font='large',
xlabel=r'$\xi$',
ylabel=r'$N(\xi)$')
# Finalize the output string
s += tikz.get_end_tikz()
# Write the whole string to a file
fp = open('basis_functions.tex', 'w')
fp.write(s)
fp.close()
x_sep=x_sep,
xmin=xmin,
xmax=xmax,
ymin=ymin,
ymax=ymax,
line_dim='thick',
xscale=xscale,
yscale=yscale,
bar_width=bar_width,
bar_offset=bar_offset)
# Plot the axes
s += tikz.get_2d_axes(xmin,
xmax,
ymin,
ymax,
tick_frac=tick_frac,
ylabel_offset=ylabel_offset,
xscale=xscale,
yscale=yscale,
xticks=xticks,
yticks=yticks,
xlabel=None,
ylabel=None)
s += tikz.get_end_tikz()
fp = open('bar_plot.tex', 'w')
fp.write(s)
fp.close()
color=colors[k % len(colors)],
symbol=symbols[k % len(symbols)],
symbol_size=0.03,
xscale=xscale,
yscale=yscale)
# Plot the axes
s += tkz.get_2d_axes(xmin,
xmax,
ymin,
ymax,
xscale=xscale,
yscale=yscale,
xticks=xticks,
yticks=yticks,
xtick_labels=xtick_labels,
ytick_labels=ytick_labels,
tick_font='small',
tick_frac=0.01,
xlabel_offset=0.075,
label_font='large',
xlabel='1/$\\rho$',
ylabel_offset=0.165,
ylabel='Aggregation error')
s += tkz.get_end_tikz()
fp = open('results/disk_exact_aggregation_error.tex', 'w')
fp.write(s)
fp.close()
