예제 #1
0
#!/usr/bin/env python3
from flat import Flat
f = Flat(dumping=True,
         m="bigraphene",
         Exc=1,
         Hc=200,
         Ex=0.3,
         omega=5e10,
         T=300,
         n=100,
         dt=1e-13,
         alltime=1e-8)
f.run(True)
예제 #2
0
#!/usr/bin/env python3

import numpy as np
from flat import Flat

freqs = np.linspace(5.25e14, 5.35e14, 3)
waves = lambda freq: [[1, 0, 0, 1e12, 0], [10, 0, 0, freq, 0]]
population = []

for freq in freqs:
    f = Flat(m="bigraphene_0.01_2",
             waves=waves(freq),
             T=300,
             n=100,
             dt=1e-15,
             alltime=1e-9)
    d = f.run()
    population.append(d["population"])

population = np.array(population)
data = np.hstack([freqs.reshape(len(freqs), 1), population[:, :, 0]])
np.savetxt("ps.dat", data)
예제 #3
0
#!/usr/bin/env python3

import sys
import numpy as np
from matplotlib import pyplot as plt
from flat import Flat

Excl = np.linspace(0, 2, 6)
Hcl = np.linspace(20, 5000, 50)
for Exc in Excl:
    line = []
    for Hc in Hcl:
        f = Flat(m="bigraphene_0.01_1",
                 Exc=Exc,
                 Hc=Hc,
                 waves=[[1, 0, 0, 1e12, 0]],
                 T=300,
                 n=100,
                 dt=1e-13,
                 alltime=1e-9)
        line.append(f.run(True)["power"][1][0])
    plt.plot(Hcl, line, label="Exc=%.2f" % Exc)
plt.legend()
plt.gca().set_ylim(bottom=0)
plt.savefig("magnetic.png")
예제 #4
0
파일: cvc.py 프로젝트: citrux/flat
#!/usr/bin/env python3

import sys
import numpy as np
from flat import Flat

result = np.zeros((2, 21))
result[0, :] = np.linspace(-2, 2, 21)
for i, Exc in enumerate(result[0, :]):
    f = Flat(m="graphene", Exc=Exc, T=300, n=100, dt=1e-14, alltime=1e-8)
    d = f.run(True)
    result[1, i] = d["v_y"][0]
np.savetxt("cvc.dat", result.transpose())