ee.sampler = function(d,net,k) {
        omcmc = fast.order.sample(d,k,c(0,20000),S=0)
        ee    = NULL
        r     = 6
        tL   = temp.ladder(1,50,6)
        eL   = energy.ladder(min(omcmc@score),min(omcmc@score)+diff(range(omcmc@score))/2,r)
        sL   = .25*(1/r)^(1/tL)
        i    = 1
        while(i <= r) {
                s = r-(i-1)
                ee.n = fast.order.sample(d,k,c(10000,s*10000),H=eL[s],T=tL[s],S=sL[s],X=ee,use.max=FALSE)
#               if(!is.null(ee)) {
#                       first = ee.n@diag@burnin
#                       pop = c(ee@population,ee.n@population[(first+1):length(ee.n@population)])
#                       sco = c(ee@score,ee.n@score[(first+1):length(ee.n@score)])
#                       diag= ee.n@diag
#                       diag@burnin = 0
#                       ee = new("sample",population=pop,temperature=ee.n@temperature,energy=ee.n@energy,diag=diag)
#               } else {
                        first = ee.n@diag@burnin
                        ee.n@population = ee.n@population[(first+1):length(ee.n@population)]
                        ee.n@score      = ee.n@score[(first+1):length(ee.n@population)]
                        ee.n@diag@burnin = 0
                        ee = ee.n
#               }
                ee = landscape(ee,eL[2:length(eL)])
                i =  i + 1
                print(ee)
        }
        os = orderSet(ee.n@population,ee.n@score,ee.n@diag@burnin+1,10)
        ds = fast.dag.wrp(d,k,os,1000,.95)
}