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
