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)
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)
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)
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)
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)
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)