Ejemplo n.º 1
0
def prop_trace(prop_list, aq, sin=True):
    "Trace of qslash*inv(prop)*(1/12q^2). Supports two forms of q."
    N = len(prop_list)
    prop_inv = inv(sum(prop_list)/N)
    qslash = dw.slash(aq)
    qslash_sin = dw.slash(np.sin(aq))
    qsquared = dw.inner(aq)
    qsquared_sin = dw.inner(np.sin(aq))
    if sin:
        return (1./12)*trace(dot(qslash_sin, prop_inv)).imag/qsquared_sin
    else:
        return (1./12)*trace(dot(qslash, prop_inv)).imag/qsquared
Ejemplo n.º 2
0
 def populate_kinematic_variables(self):
     '''Requires that L, T, and a are defined.'''
     self.ap = dw.ap(self.p, self.tw, self.L, self.T)
     self.ap2 = dw.ap(self.p2, self.tw, self.L, self.T)
     self.aq = dw.aq(self.ap, self.ap2)  # ap - ap2
     self.apSq = dw.inner(self.ap)  # (ap)^2
     self.mu = dw.mu(self.ap, self.a)