Beispiel #1
0
import pylab as plt
import TLL as tll
import numpy as np

ls, gs, corr = tll.dvrg(1.5, .5)
plt.figure(1)
plt.subplot(211)
for g in gs:
    plt.plot(ls, np.absolute(g[0]), label=g[1], color=g[2])
plt.legend()
plt.subplot(212)
for cor in corr:
    plt.semilogy(ls, np.absolute(cor[0]), label=cor[1], color=cor[2])
plt.legend()

xpnts = tll.tll(1.5, .5)
plt.figure(2)
plt.xlim([5, 20])
plt.ylim([0, .1])
for tup in xpnts:
    lst = tup[0]
    plt.plot(range(1, 20),
             map(lambda x: (1. / x)**(lst[-1]), range(1, 20)),
             label=tup[1],
             color=tup[2],
             marker=tup[3])
plt.legend()
plt.show()
Beispiel #2
0
import matplotlib.pyplot as plt
import TLL as tll
import numpy as np

plt.rcParams['mathtext.fontset'] = 'stix'

u0 = 1.3
upi = .5
ls, gs, corr = tll.dvrg(u0, upi)
lmax = ls[-1]

g1 = np.array(gs[0][0])
g2 = np.array(gs[1][0])
g3 = np.array(gs[2][0])
g4 = np.array(gs[3][0])
cdws = g2 - g3 - 2.0*g1
cdwb = g2 + g3 - 2.0*g1
sdws = g2 + g3
sdwb = g2 - g3
ss = -g1 - g2
st = g1 - g2

phasegs = [(cdws,'cdws','red'),
           (cdwb,'cdwb','orange'),
           (sdws,'sdws','black'),
           (sdwb,'sdwb','green'),
           (ss,'ss','blue'),
           (st,'st','purple')]
fig = plt.figure()
ax1 = fig.add_subplot(211)
for g in gs:
Beispiel #3
0
import pylab as plt
from numpy import arange
import TLL as tll

# Create a dictionary that maps tuples in phase space to tuples containing
# (label, color, marker) for plotting.
phase = {}
for uzr in arange(0., 1., 1. / 10.):
    for upi in arange(0., 1., 1. / 10.):
        ls, g, corr = tll.dvrg(uzr, upi)
        tup = corr[-1]
        phase[(uzr, upi)] = tup[1:]

for uzr in arange(0., 1., 1. / 10.):
    for upi in arange(1. + 1. / 10., 2. + 1. / 10., 1. / 10.):
        ls, g, corr = tll.dvrg(uzr, upi)
        tup = corr[-1]
        phase[(uzr, upi)] = tup[1:]

for uzr in arange(1., 2. + 1. / 10., 1. / 10.):
    for upi in arange(1. + 1. / 10., 2. + 1. / 10., 1. / 10.):
        ls, g, corr = tll.dvrg(uzr, upi)
        tup = corr[-1]
        phase[(uzr, upi)] = tup[1:]

for uzr in arange(1. + 1. / 10., 2. + 1. / 10., 1. / 10.):
    for upi in arange(0., 1., 1. / 10.):
        xpnt = tll.tll(uzr, upi)
        tup = xpnt[0]
        phase[(uzr, upi)] = tup[1:]
Beispiel #4
0
import pylab as plt
import TLL as tll
import numpy as np

ls, gs, corr = tll.dvrg(1.5, .5)
plt.figure(1)
plt.subplot(211)
for g in gs:
    plt.plot(ls, np.absolute(g[0]), label=g[1], color=g[2])
plt.legend()
plt.subplot(212)
for cor in corr:
    plt.semilogy(ls, np.absolute(cor[0]), label=cor[1], color=cor[2])
plt.legend()

xpnts = tll.tll(1.5, .5)
plt.figure(2)
plt.xlim([5,20])
plt.ylim([0,.1])
for tup in xpnts:
    lst = tup[0]
    plt.plot(range(1,20), map(lambda x: (1./x)**(lst[-1]), range(1,20)),
             label=tup[1], color=tup[2], marker=tup[3])
plt.legend()
plt.show()
Beispiel #5
0
import matplotlib.pyplot as plt
import TLL as tll
import numpy as np

plt.rcParams['mathtext.fontset'] = 'stix'

u0 = 1.3
upi = .5
ls, gs, corr = tll.dvrg(u0, upi)
lmax = ls[-1]

g1 = np.array(gs[0][0])
g2 = np.array(gs[1][0])
g3 = np.array(gs[2][0])
g4 = np.array(gs[3][0])
cdws = g2 - g3 - 2.0 * g1
cdwb = g2 + g3 - 2.0 * g1
sdws = g2 + g3
sdwb = g2 - g3
ss = -g1 - g2
st = g1 - g2

phasegs = [(cdws, 'cdws', 'red'), (cdwb, 'cdwb', 'orange'),
           (sdws, 'sdws', 'black'), (sdwb, 'sdwb', 'green'),
           (ss, 'ss', 'blue'), (st, 'st', 'purple')]
fig = plt.figure()
ax1 = fig.add_subplot(211)
for g in gs:
    ax1.plot(ls,
             g[0],
             label=g[1],
Beispiel #6
0
import pylab as plt
from numpy import arange
import TLL as tll

# Create a dictionary that maps tuples in phase space to tuples containing
# (label, color, marker) for plotting.
phase = {}
for uzr in arange(0., 1., 1./10.):
    for upi in arange(0., 1., 1./10.):
        ls, g, corr = tll.dvrg(uzr,upi)
        tup = corr[-1]
        phase[(uzr,upi)] = tup[1:]

for uzr in arange(0., 1., 1./10.):
    for upi in arange(1.+1./10., 2.+1./10., 1./10.):
        ls, g, corr = tll.dvrg(uzr,upi)
        tup = corr[-1]
        phase[(uzr,upi)] = tup[1:]
    
for uzr in arange(1., 2.+1./10., 1./10.):
    for upi in arange(1.+1./10., 2.+1./10., 1./10.):
        ls, g, corr = tll.dvrg(uzr,upi)
        tup = corr[-1]
        phase[(uzr,upi)] = tup[1:]

for uzr in arange(1.+1./10., 2.+1./10., 1./10.):
    for upi in arange(0., 1., 1./10.):
        xpnt = tll.tll(uzr,upi)
        tup = xpnt[0]
        phase[(uzr,upi)] = tup[1:]