Exemplo n.º 1
0
    #    ax.set_zlabel(r'$p_{\theta}$', fontsize=16)
    #    plt.title('Standard map')
    #    unit   = 0.5
    #    q_tick = np.arange(0, 2.0+unit, unit)
    #    q_label = [r"$0$",r"$\frac{1}{2}\pi$",r"$\pi$",r"$\frac{3}{2}\pi$",r"$2\pi$"]
    #    ax.set_ylim(0, 2*pi)
    #    ax.set_yticks(q_tick*pi)
    #    ax.set_yticklabels(q_label, fontsize=18)
    #    ax.set_xlim(0, kmax)
    #    if max(ary_p[ntransient:]) > max_p:
    #        max_p = max(ary_p[ntransient:])
    #    if min(ary_p[ntransient:]) < min_p:
    #        min_p = min(ary_p[ntransient:])
    #    margin =  (max_p-min_p)*0.05
    #    ax.set_zlim(min_p-margin, max_p+margin)

    # --- Plot the standard map. ---
    mlab.points3d([
        0.1,
    ] * (nmax - ntransient + 1), ary_q[ntransient:], ary_p[ntransient:])

    plt.show()


if __name__ == '__main__':

    window = SetParameter()
    default_params = [{'theta': 1.5}, {'p': 0.2}, {'kmax': 0.1}]
    commands = [{'OK': static}]
    window.show_setting_window(default_params, commands)