Example #1
0
File: core.py Project: liangjj/BdG
 def getTc(self):
     Phi = self.overlaps()
     integrator = Integrator(self.struct)
     T = 0.1
     dims = self.struct.cont_dims
     lims = 2*np.ones(len(dims))
     A = integrator.integrate(lambda x,y,z: self.thermal_weight_LO_k(x, y, z, T),
             lims, dims)
     return A
Example #2
0
File: test.py Project: liangjj/BdG
from parameters import Material, Structure, Parameters, Integrator
import numpy as np
from core import Hamiltonian, System
tin = Material(148.0, 0.000618, 0.666666, 0.1)
Lz = 50.0

struct = Structure(np.array([Lz]))
params = Parameters(tin, struct)
print(params.print_parameters())
ham = Hamiltonian()
integrator = Integrator(struct)
sys = System(ham, params, struct)

print(integrator.integrate(ham.spectrum, [10, 10, 10]))
print(sys.getTc())
print(struct.cont_dims)