Python cov 예제들

예제 #1

def set_C(norm=30.):
    ds.clear_cache()
    Cs, iCs = {}, {}
    for k in ks:
        Ndim = ds.x[k].shape[0]
        Cs[k] = oqe.cov(ds.x[k][:, :], ds.w[k][:, ]) + norm * np.identity(Ndim)
        ds.set_C({k: Cs[k]})
        iCs[k] = ds.iC(k)

예제 #2

파일: red_pspec_cov.py 프로젝트: SaulAryehKohn/capo

def set_C(norm=30.):
    ds.clear_cache()
    Cs,iCs = {},{}
    for k in ks:
        Ndim = ds.x[k].shape[0]
        Cs[k] = oqe.cov(ds.x[k][:,:],ds.w[k][:,])+norm*np.identity(Ndim)
        ds.set_C({k:Cs[k]})
        iCs[k] = ds.iC(k)

예제 #3

파일: comb_cov.py 프로젝트: SaulAryehKohn/capo

def set_C(dst,norm=3e-6):
    dst.clear_cache()
    Cs,iCs = {},{}
    #import IPython; IPython.embed()
    for k in dst.x.keys():
        #Cs[k] = sum([oqe.cov(ds.x[k][:,400:],ds.w[k][:,400:])+norm*np.identity(NCHAN) for ki in ks if ki != k])
        #Cs[k] = sum([oqe.cov(ds.x[ki][:,400:],ds.w[ki][:,400:])+norm*np.identity(NCHAN) for ki in ks if ki[2] != k[2]])
        Ndim = dst.x[k].shape[0]
        Cs[k] = oqe.cov(dst.x[k][:,0:],dst.w[k][:,0:])+norm*np.identity(Ndim)
        dst.set_C({k:Cs[k]})
        iCs[k] = dst.iC(k)

예제 #4

파일: corss_pspec_cov_v001.py 프로젝트: nkern/capo

def set_C(dst, norm=3e-6):
    dst.clear_cache()
    Cs, iCs = {}, {}
    #import IPython; IPython.embed()
    for k in dst.x.keys():
        #Cs[k] = sum([oqe.cov(ds.x[k][:,400:],ds.w[k][:,400:])+norm*np.identity(NCHAN) for ki in ks if ki != k])
        #Cs[k] = sum([oqe.cov(ds.x[ki][:,400:],ds.w[ki][:,400:])+norm*np.identity(NCHAN) for ki in ks if ki[2] != k[2]])
        Ndim = dst.x[k].shape[0]
        Cs[k] = oqe.cov(dst.x[k][:, 0:],
                        dst.w[k][:, 0:]) + norm * np.identity(Ndim)
        dst.set_C({k: Cs[k]})
        iCs[k] = dst.iC(k)

예제 #5

파일: omni_pspec_cov_v000.py 프로젝트: AaronParsons/capo

def set_C(norm=3e-6):
    ds.clear_cache()
    Cs,iCs = {},{}
    for k in ks:
        #Cs[k] = sum([capo.oqe.cov(ds.x[k][:,400:],ds.w[k][:,400:])+norm*np.identity(NCHAN) for ki in ks if ki != k])
        #Cs[k] = sum([capo.oqe.cov(ds.x[ki][:,400:],ds.w[ki][:,400:])+norm*np.identity(NCHAN) for ki in ks if ki[2] != k[2]])
        Cs[k] = sum([oqe.cov(ds.x[k][:,400:],ds.w[k][:,400:])+norm*np.identity(NCHAN) for ki in ks if ki[2] != k[2]])
        #w = np.where(ds.w[ki] > 0, 1, 0)
        #Cs[k] = sum([capo.oqe.cov(ds.x[ki][:,400:],w[:,400:])+norm*np.identity(NCHAN) for ki in ks if ki[2] != k[2]])
        #Cs[k] = sum([ds.C(k)+norm*np.identity(NCHAN) for ki in ks if ki != k])
        #Cs[k] = sum([ds.C(k)+norm*np.identity(NCHAN) for ki in data if ki[2] != k[2]])
        ds.set_C({k:Cs[k]})
        iCs[k] = ds.iC(k)

예제 #6

파일: red_pspec_cov.py 프로젝트: karakundert/capo

def set_C(norm=3e-6):
    ds.clear_cache()
    Cs, iCs = {}, {}
    for k in ks:
        #Cs[k] = sum([capo.oqe.cov(ds.x[k][:,400:],ds.w[k][:,400:])+norm*np.identity(NCHAN) for ki in ks if ki != k])
        #Cs[k] = sum([capo.oqe.cov(ds.x[ki][:,400:],ds.w[ki][:,400:])+norm*np.identity(NCHAN) for ki in ks if ki[2] != k[2]])
        Cs[k] = sum([
            oqe.cov(ds.x[k][:, 400:], ds.w[k][:, 400:]) +
            norm * np.identity(NCHAN) for ki in ks if ki[2] != k[2]
        ])
        #w = np.where(ds.w[ki] > 0, 1, 0)
        #Cs[k] = sum([capo.oqe.cov(ds.x[ki][:,400:],w[:,400:])+norm*np.identity(NCHAN) for ki in ks if ki[2] != k[2]])
        #Cs[k] = sum([ds.C(k)+norm*np.identity(NCHAN) for ki in ks if ki != k])
        #Cs[k] = sum([ds.C(k)+norm*np.identity(NCHAN) for ki in data if ki[2] != k[2]])
        ds.set_C({k: Cs[k]})
        iCs[k] = ds.iC(k)