raise ValueError("\nusage: python %s dim k qmin qmax pmin pmax " % sys.argv[0]) dim = int(sys.argv[1]) k = float(sys.argv[2]) qmin = float(sys.argv[3]) qmax = float(sys.argv[4]) pmin = float(sys.argv[5]) pmax = float(sys.argv[6]) domain = [[qmin,qmax],[pmin,pmax]] #dim = 10 #domain=[[0.0, 1.0],[-0.5,0.5]] #k=1.0 map = StandardMap(k=k) qmap = Qmap(map, dim, domain) evals, evecs = qmap.eigen() for i, evec in enumerate(evecs): evec.savetxt("eigen_qrep_%i.dat" % i) sample=300 tmax=300 q = np.random.random(sample) p = (np.random.random(sample) -0.5) * (pmax - pmin) * 2 with open("trajectories.dat", "w") as f: for i in range(tmax): pp = p - map.func0(q) qq = q + map.func1(pp) q = qq - np.floor(qq)
if len(sys.argv) != 7: raise ValueError("\nusage: python %s dim k qmin qmax pmin pmax " % sys.argv[0]) dim = int(sys.argv[1]) k = float(sys.argv[2]) qmin = float(sys.argv[3]) qmax = float(sys.argv[4]) pmin = float(sys.argv[5]) pmax = float(sys.argv[6]) tmax = 100 map = StandardMap(k=k) domain = [[qmin,qmax],[pmin,pmax]] qmap = Qmap(map, dim, domain) state = qmap.getState().cs(0.5,0.5) #init[dim/2] = 1.0 + 0.j qmap.setInit(state) tmax=10 for i in range(tmax): state = qmap.getIn() state.savetxt("evolv_qrep_%d.dat" % i) qmap.evol() sample=300 tmax=300 q = np.random.random(sample) p = (np.random.random(sample) -0.5) * (pmax - pmin) * 2
sys.argv[0]) dim = int(sys.argv[1]) k = float(sys.argv[2]) qmin = float(sys.argv[3]) qmax = float(sys.argv[4]) pmin = float(sys.argv[5]) pmax = float(sys.argv[6]) domain = [[qmin, qmax], [pmin, pmax]] #dim = 10 #domain=[[0.0, 1.0],[-0.5,0.5]] #k=1.0 map = StandardMap(k=k) qmap = Qmap(map, dim, domain) evals, evecs = qmap.eigen() for i, evec in enumerate(evecs): evec.savetxt("eigen_qrep_%i.dat" % i) sample = 300 tmax = 300 q = np.random.random(sample) p = (np.random.random(sample) - 0.5) * (pmax - pmin) * 2 with open("trajectories.dat", "w") as f: for i in range(tmax): pp = p - map.func0(q) qq = q + map.func1(pp) q = qq - np.floor(qq) p = pp