def METHOD_NAME(x, z, Bx, Bz, Bp, Bx2, Bz2, Bp2):
    """
    Plotting utility for supporting comparisons. Not a test.
    """
    plot_defaults()
    f, ax = plt.subplots(3, 4)
    ax[0, 0].contourf(x, z, Bx)
    ax[0, 0].set_title("$B_{x}$ Green's function")
    ax[0, 1].contourf(x, z, Bx2)
    ax[0, 1].set_title("$B_{x}$ Biot-Savart differential function")
    cm = ax[0, 2].contourf(x, z, 100 * (Bx - Bx2) / Bx)
    f.colorbar(cm, ax[0, 3])
    ax[0, 2].set_title("Relative error [%]")
    ax[0, 0].set_aspect("equal")
    ax[0, 1].set_aspect("equal")
    ax[0, 2].set_aspect("equal")
    ax[0, 3].set_aspect(20)
    ax[1, 0].contourf(x, z, Bx)
    ax[1, 0].set_title("$B_{z}$ Green's function")
    ax[1, 1].contourf(x, z, Bx2)
    ax[1, 1].set_title("$B_{z}$ Biot-Savart differential function")
    cm = ax[1, 2].contourf(x, z, 100 * (Bz - Bz2) / Bz)
    f.colorbar(cm, ax[1, 3])
    ax[1, 2].set_title("Relative error [%]")
    ax[1, 0].set_aspect("equal")
    ax[1, 1].set_aspect("equal")
    ax[1, 2].set_aspect("equal")
    ax[1, 3].set_aspect(20)
    ax[2, 0].contourf(x, z, Bp)
    ax[2, 0].set_title("$B_{p}$ Green's function")
    ax[2, 1].contourf(x, z, Bp2)
    ax[2, 1].set_title("$B_{p}$ Biot-Savart differential function")
    cm = ax[2, 2].contourf(x, z, 100 * (Bp - Bp2) / Bp)
    ax[2, 2].set_title("Relative error [%]")
    f.colorbar(cm, ax[2, 3])
    ax[2, 0].set_aspect("equal")
    ax[2, 1].set_aspect("equal")
    ax[2, 2].set_aspect("equal")
    ax[2, 3].set_aspect(20)
    plt.show()
    plt.close(f)