def do_plots_c(Ud, Unew):
    ' plot Ud,new and Ud with zoom on the bug '
    pylab.clf()
    pylab.cla()
    f = pylab.figure()
    f.text(0.5, 0.95, '$U_{\\rm d}$ (left) and $U_{\\rm d, new}$ (right) ', horizontalalignment='center')
    pylab.subplot(221)
    pylab.imshow(Ud[0])
    pylab.ylabel('# of cells', size=8)
    pylab.subplot(223)
    pylab.imshow(Ud[1])
    pylab.xlim(1, 32)
    pylab.xlabel('# of cells', size=8)
    pylab.ylabel('# of cells', size=8)
    pylab.subplot(222)
    pylab.imshow(Unew[0])
    pylab.ylabel('# of cells', size=8)
    pylab.subplot(224)
    pylab.imshow(Unew[1])
    pylab.xlim(1, 32)
    pylab.xlabel('# of cells', size=8)
    pylab.ylabel('# of cells', size=8)
    pylab.savefig('plots/item_c_Udnew.png')
plot Ud,new and Ud with zoom on the bug
def do_plots_c(Ud, Unew):
    '  '
    pylab.clf()
    pylab.cla()
    f = pylab.figure()
    f.text(0.5, 0.95, '$U_{\\rm d}$ (left) and $U_{\\rm d, new}$ (right) ', horizontalalignment='center')
    pylab.subplot(221)
    pylab.imshow(Ud[0])
    pylab.ylabel('# of cells', size=8)
    pylab.subplot(223)
    pylab.imshow(Ud[1])
    pylab.xlim(1, 32)
    pylab.xlabel('# of cells', size=8)
    pylab.ylabel('# of cells', size=8)
    pylab.subplot(222)
    pylab.imshow(Unew[0])
    pylab.ylabel('# of cells', size=8)
    pylab.subplot(224)
    pylab.imshow(Unew[1])
    pylab.xlim(1, 32)
    pylab.xlabel('# of cells', size=8)
    pylab.ylabel('# of cells', size=8)
    pylab.savefig('plots/item_c_Udnew.png')