Python findPhotonGSの例

コード例 #1

ファイル: calcConductivity.py プロジェクト: ce335805/comeChainComeShine

def getPhGS(eta):
    gs = arbOrder.findGS(eta, 3)
    gsJ = eF.J(gs)
    gsT = eF.T(gs)
    gsJ = 0.
    gsT = -2. / np.pi * prms.chainLength

    return phState.findPhotonGS([gsT, gsJ], eta, 3)

コード例 #2

def expectAnnihil(eta):
    gs = arbOrder.findGS(eta, 3)
    gsJ = eF.J(gs)
    gsT = eF.T(gs)

    ptGS = phState.findPhotonGS([gsT, gsJ], eta, 3)
    a = setUpAnnihil(eta)

    return np.dot(np.conj(ptGS), np.dot(a, ptGS))

コード例 #3

def currentGS(eta):
    gs = arbOrder.findGS(eta, 3)
    gsJ = eF.J(gs)
    gsT = eF.T(gs)

    ptGS = phState.findPhotonGS([gsT, gsJ], eta, 3)
    jOp = setupCurrentOperator(gsT, gsJ, eta)

    current = np.dot(np.conj(ptGS), np.dot(jOp, ptGS))
    return current

コード例 #4

def expectSinA(eta):
    gs = arbOrder.findGS(eta, 3)
    gsJ = eF.J(gs)
    gsT = eF.T(gs)

    ptGS = phState.findPhotonGS([gsT, gsJ], eta, 3)
    x = setUpA(eta)
    sinX = sciLin.sinm(x)

    return np.dot(np.conj(ptGS), np.dot(sinX, ptGS))

コード例 #5

ファイル: nonEqGreen.py プロジェクト: ce335805/comeChainComeShine

def getPhGS(eta):
    gs = arbOrder.findGS(eta, 3)
    gsJ = eF.J(gs)
    gsT = eF.T(gs)
    return phState.findPhotonGS([gsT, gsJ], eta, 3)