def plot_dUdv(domain, spots, vars, plotter=wireframe): """ Plot the error in the derivative across each row. """ dUdv = center_diff(domain,axis=1)/center_diff(vars) plotter(dUdv, spots, vars) pylab.title("First deriv w.r.t. var") pylab.xlabel("Var") pylab.ylabel("% of strike") pylab.show()
def plot_d2Udv2(domain, spots, vars, plotter=wireframe): """ Plot the error in the second derivative across each row. """ d2Udv2 = center_diff(domain,n=2,axis=1)/center_diff(vars, 2) plotter(d2Udv2, spots, vars) pylab.title("Second deriv w.r.t. var") pylab.xlabel("Var") pylab.ylabel("% of strike") pylab.show()
def plot_dUds_err(domain, analytical, spots, vars, plotter=wireframe): """ Plot the error in the derivative down each columns. """ dUds = center_diff(domain,axis=0)/center_diff(vars) danalyticalds = center_diff(analytical,axis=0)/center_diff(vars) plotter(dUds - danalyticalds, spots, vars) pylab.title("Error in $\Delta$") pylab.xlabel("Var") pylab.ylabel("% of strike")