def run_gui_nodisorder(L=24): from pele.gui import run_gui dim=[L,L] system = XYModlelSystem(dim=dim, phi_disorder=0.) system.params.basinhopping.temperature=10. db = system.create_database("xy_%dx%d_nodisorder.sqlite" %(L,L)) run_gui(system, db=db)
def run_gui_nodisorder(L=24): from pele.gui import run_gui dim = [L, L] system = XYModlelSystem(dim=dim, phi_disorder=0.) system.params.basinhopping.temperature = 10. db = system.create_database("xy_%dx%d_nodisorder.sqlite" % (L, L)) run_gui(system, db=db)
def main(): from xy_model_system import XYModlelSystem fname = "gmin.coords" fname = "goldstone_mode.coords" fname = "vortex2.coords" reduced_coords = np.genfromtxt(fname) print(len(reduced_coords)) print np.sqrt(len(reduced_coords)) dim = [24, 24] print "dim", dim system = XYModlelSystem(dim, phi_disorder=0) coords = system.pot.coords_converter.get_full_coords(reduced_coords) x = np.zeros(len(coords)) y = x.copy() dx = x.copy() dy = x.copy() for node in system.pot.G.nodes(): xyz = system.node2xyz(node) i = system.pot.indices[node] x[i] = xyz[0] y[i] = xyz[1] theta = coords[i] dx[i] = np.cos(theta) dy[i] = np.sin(theta) if True: pot = system.get_potential() energies = pot.get_spin_energies(coords) else: energies = None # plt.quiver(x, y, dx, dy, energies, pivot="middle", cmap="YlOrRd_r") if energies is not None: plt.quiver(x, y, dx, dy, energies, pivot="middle", cmap="cool_r") else: plt.quiver(x, y, dx, dy, pivot="middle") d = 0.5 plt.xlim([-d, dim[0] - d]) plt.ylim([-d, dim[1] - d]) plt.tight_layout() plt.xticks([]) plt.yticks([]) plt.gcf().set_facecolor('white') plt.box(on=False) plt.savefig(fname + ".eps", format="eps")