def mandelbrot_wrapper(row, col):
x = col/(num_pixels/4.) - 2.
y = row/(num_pixels/4.) - 2.
return ((row, col), P2.mandelbrot(x, y))
mandelbrot_rdd = indices.map(lambda a: mandelbrot_wrapper(*a))
# Now collect the data & plot
mandelbrot_result = mandelbrot_rdd.collect()
plt.grid(False)
# I slightly redefined the draw image function as the original
# implementation annoyed me...I did not want to collect in a draw function!
P2.draw_image(data=mandelbrot_result)
plt.savefig('P2a_mandelbrot.png', dpi=200, bbox_inches='tight')
plt.clf()
# Now create the histogram...I recognize that mandelbrot is computed twice
# but it is for my sanity
summed_rdd = P2.sum_values_for_partitions(mandelbrot_rdd)
summed_result = summed_rdd.collect()
plt.hist(summed_result, bins=np.logspace(3, 8, 20))
sns.rugplot(summed_result, color='red')
plt.gca().set_xscale('log')
plt.xlabel('Total Number of Iterations on Partition')
plt.ylabel('Partition Count')
def mandelbrot_wrapper(row, col):
x = col / (num_pixels / 4.) - 2.
y = row / (num_pixels / 4.) - 2.
return ((row, col), P2.mandelbrot(x, y))
mandelbrot_rdd = indices.map(lambda a: mandelbrot_wrapper(*a))
# Now collect the data & plot
mandelbrot_result = mandelbrot_rdd.collect()
plt.grid(False)
# I slightly redefined the draw image function as the original
# implementation annoyed me...I did not want to collect in a draw function!
P2.draw_image(data=mandelbrot_result)
plt.savefig('P2a_mandelbrot.png', dpi=200, bbox_inches='tight')
plt.clf()
# Now create the histogram...I recognize that mandelbrot is computed twice
# but it is for my sanity
summed_rdd = P2.sum_values_for_partitions(mandelbrot_rdd)
summed_result = summed_rdd.collect()
plt.hist(summed_result, bins=np.logspace(3, 8, 20))
sns.rugplot(summed_result, color='red')
plt.gca().set_xscale('log')
plt.xlabel('Total Number of Iterations on Partition')
plt.ylabel('Partition Count')
