def evolx(P, tn, PL1, PL2, gold, rk, nj, Qi, ptj, qij, rnz): MU = [ P[0], P[1], P[2], P[3], P[4], P[5], P[6], P[7], P[8], P[9], P[10], P[11], P[12] ] Ntot = P[13] pqr = {} ltn = [] nefft = [] lnefft = [] ntest = [] lntest = [] lth = [] ngt = [] fitn = [] trns = [] flg = 'all-good' gen = {} trs = {} k = 0 for i in gold.keys(): gl = [] if gold[i][0] == 'te': gl = [gold[i][0], gold[i][1]] gen[k] = gl k += 1 if gold[i][0] == 'eff': if nj == 0: gl = [gold[i][0], gold[i][1], gold[i][2]] else: qj = rk.uniform_pos() snew = gold[i][2] + (0.01 * (1.0 - 2.0 * rk.uniform_pos())) if ((qj < Qi) and (snew > 0.0)): gl = [gold[i][0], gold[i][1], snew] else: gl = [gold[i][0], gold[i][1], 0.0] gen[k] = gl k += 1 gold = None for nk in range(tn): trs = Lffit.trates(gen, MU) rij, sr = Lffit.montec(trs, rk) trns.append(rij) if sr != 'TRUE': pqr[nk] = [rij[0], rij[1], gen] ltn.append(Lffit.lent(gen)) lth.append(Lth) wq = Lffit.ft(gen) fitn.append(1.0 - (wq / (wq + wo))) ngt.append(len(gen.keys())) n1x = 0.0 n2x = 0.0 m1x = 0.0 m2x = 0.0 for gk in gen.keys(): if gen[gk][0] == 'eff': n1x += 1 m1x += gen[gk][1] if gen[gk][0] == 'te': n2x += 1 m2x += gen[gk][1] nefft.append(n1x) lnefft.append(m1x) ntest.append(n2x) lntest.append(m2x) gu = {} gu = LffitB.transform(rij, gen, rk, PL1, PL2) gen = {} for i in gu.keys(): gen[i] = gu[i] else: flg = 'sumzero' break if len(gen.keys()) > 2000: flg = 'limit-reached' break #################### import os po = '../RUN' nx = str(rnz) + '/n' + str(ptj) if not os.path.exists(po + nx): os.makedirs(po + nx) efflens, telens = Lffit.splitgen(gen) print len(efflens), len(telens) pthl = po + nx if telens or efflens: pdt.savedata(lth, ltn, ngt, efflens, telens, P, SEED, pthl, fitn, trns, nefft, ntest, lnefft, lntest, tn, JMPS) qij.put(gen)
def evolx(P, tn, PL1, PL2, gold, rk, nj, Qi, ptj, qij, rnz): #print P,tn,PL1,PL2,gold,rk,nj,Qi,ptj MU = [ P[0], P[1], P[2], P[3], P[4], P[5], P[6], P[7], P[8], P[9], P[10], P[11], P[12] ] #MU=[Np,m1,m2,m3,m4,m5,m6,m7,m8,Lth,beta1,beta2,wo] Ntot = P[13] pqr = {} ltn = [] nefft = [] lnefft = [] ntest = [] lntest = [] lth = [] ngt = [] fitn = [] trns = [] flg = 'all-good' gen = {} trs = {} k = 0 for i in gold.keys(): gl = [] if gold[i][0] == 'te': gl = [gold[i][0], gold[i][1]] gen[k] = gl k += 1 if gold[i][0] == 'eff': if nj == 0: gl = [gold[i][0], gold[i][1], gold[i][2]] else: qj = rk.uniform_pos() snew = gold[i][2] + (0.01 * (1.0 - 2.0 * rk.uniform_pos())) if ((qj < Qi) and (snew > 0.0)): gl = [gold[i][0], gold[i][1], snew] else: gl = [gold[i][0], gold[i][1], 0.0] gen[k] = gl k += 1 gold = None for nk in range(tn): trs = Lffit.trates(gen, MU) rij, sr = Lffit.montec(trs, rk) trns.append(rij) if sr != 'TRUE': pqr[nk] = [rij[0], rij[1], gen] ltn.append(Lffit.lent(gen)) lth.append(Lth) wq = Lffit.ft(gen) fitn.append(1.0 - (wq / (wq + wo))) ngt.append(len(gen.keys())) n1x = 0.0 n2x = 0.0 m1x = 0.0 m2x = 0.0 for gk in gen.keys(): if gen[gk][0] == 'eff': n1x += 1 m1x += gen[gk][1] if gen[gk][0] == 'te': n2x += 1 m2x += gen[gk][1] nefft.append(n1x) lnefft.append(m1x) ntest.append(n2x) lntest.append(m2x) gu = {} gu = LffitB.transform(rij, gen, rk, PL1, PL2) gen = {} for i in gu.keys(): gen[i] = gu[i] #gc.collect() else: flg = 'sumzero' break if len(gen.keys()) > 2000: flg = 'limit-reached' break #################### import os po = '../RUN' nx = str(rnz) + '/n' + str(ptj) if not os.path.exists(po + nx): os.makedirs(po + nx) efflens, telens = Lffit.splitgen(gen) print len(efflens), len(telens) pthl = po + nx #pdt.savedata(lth,ltn,ngt,efflens,telens,P,SEED,pthl,wfitn) if telens or efflens: #[lth,ltn,ngt,pqr,gen,flg,fitn,trns] #lth,ltn,ngt,pqr, gen,info,wfitn,trns #print("HELLO") #***********add below tn and number of jumps pdt.savedata(lth, ltn, ngt, efflens, telens, P, SEED, pthl, fitn, trns, nefft, ntest, lnefft, lntest, tn, JMPS) ###*if telens and efflens: ###* pdt.datadisplay(lth,ltn,ngt,efflens,telens,500,50,P,SEED,pthl,fitn,Qi,) #del lth,ltn,ngt,pqr,info,wfitn,trns #gc.collect() ###*else: ###* print("NOT-DATA-TO-SHOW") #################### qij.put(gen)
kj+=1 gl.append(k) if kj==3: gl[2]=0.0 gold[i]=gl #[gen[i][0],gen[i][1]] ############ del gen gen={} lth,ltn,ngt,pqr, gen,info,wfitn,trns = lmainjumps.job(P,SEED,neff, nte,tn, RXCLRPTS, CRKPTS, TES, n, gold,Qi) efflens, telens =Lffit.splitgen(gen) print len(efflens), len(telens) pthl=pthk # pdt.savedata(lth,ltn,ngt,efflens,telens,P,SEED,pthl,wfitn) if telens or efflens: pdt.savedata(lth,ltn,ngt,efflens,telens,P,SEED,pthl,wfitn,trns) if telens and efflens: pdt.datadisplay(lth,ltn,ngt,efflens,telens,500,50,P,SEED,pthl,wfitn) del lth,ltn,ngt,pqr,info,wfitn,trns gc.collect() else: print("NOT-DATA-TO-SHOW") gc.collect() print("CARRY ON") del gen, gold gen={} #gc={} gold={} gc.collect() m3=0.1*m3
def evolx(P,tn,PL1,PL2,gold,rk,nj,Qi,ptj,qij,rnz): #print P,tn,PL1,PL2,gold,rk,nj,Qi,ptj MU=[P[0],P[1],P[2],P[3],P[4],P[5],P[6],P[7],P[8],P[9],P[10],P[11],P[12]] #MU=[Np,m1,m2,m3,m4,m5,m6,m7,m8,Lth,beta1,beta2,wo] Ntot=P[13] pqr={} ltn=[] nefft=[] lnefft=[] ntest=[] lntest=[] lth=[] ngt=[] fitn=[] trns=[] flg='all-good' gen={} trs={} k=0 for i in gold.keys(): gl=[] if gold[i][0]=='te': gl=[gold[i][0],gold[i][1]] gen[k]=gl k+=1 if gold[i][0]=='eff': if nj==0: gl=[gold[i][0],gold[i][1],gold[i][2]] else: qj=rk.uniform_pos() snew=gold[i][2]+(0.01*( 1.0-2.0*rk.uniform_pos() )) if ((qj<Qi) and (snew>0.0)): gl=[gold[i][0],gold[i][1],snew] else: gl=[gold[i][0],gold[i][1],0.0] gen[k]=gl k+=1 gold=None for nk in range(tn): trs=Lffit.trates(gen,MU) rij,sr =Lffit.montec(trs,rk) trns.append(rij) if sr!='TRUE': pqr[nk]=[rij[0],rij[1],gen] ltn.append(Lffit.lent(gen)) lth.append(Lth) wq=Lffit.ft(gen) fitn.append(1.0-(wq/(wq+wo))) ngt.append(len(gen.keys())) n1x=0.0 n2x=0.0 m1x=0.0 m2x=0.0 for gk in gen.keys(): if gen[gk][0]=='eff': n1x+=1 m1x+=gen[gk][1] if gen[gk][0]=='te': n2x+=1 m2x+=gen[gk][1] nefft.append(n1x) lnefft.append(m1x) ntest.append(n2x) lntest.append(m2x) gu={} gu=LffitB.transform(rij,gen,rk, PL1, PL2) gen={} for i in gu.keys(): gen[i]=gu[i] #gc.collect() else: flg='sumzero' break if len(gen.keys())>2000: flg='limit-reached' break #################### import os po='../RUN' nx=str(rnz)+'/n'+str(ptj) if not os.path.exists(po+nx): os.makedirs(po+nx) efflens, telens =Lffit.splitgen(gen) print len(efflens), len(telens) pthl=po+nx #pdt.savedata(lth,ltn,ngt,efflens,telens,P,SEED,pthl,wfitn) if telens or efflens: #[lth,ltn,ngt,pqr,gen,flg,fitn,trns] #lth,ltn,ngt,pqr, gen,info,wfitn,trns #print("HELLO") #***********add below tn and number of jumps pdt.savedata(lth,ltn,ngt,efflens,telens,P,SEED,pthl,fitn,trns,nefft,ntest,lnefft,lntest,tn,JMPS) ###*if telens and efflens: ###* pdt.datadisplay(lth,ltn,ngt,efflens,telens,500,50,P,SEED,pthl,fitn,Qi,) #del lth,ltn,ngt,pqr,info,wfitn,trns #gc.collect() ###*else: ###* print("NOT-DATA-TO-SHOW") #################### qij.put(gen)
gl.append(k) if kj == 3: gl[2] = 0.0 gold[i] = gl #[gen[i][0],gen[i][1]] ############ del gen gen = {} lth, ltn, ngt, pqr, gen, info, wfitn, trns = lmainjumps.job( P, SEED, neff, nte, tn, RXCLRPTS, CRKPTS, TES, n, gold, Qi) efflens, telens = Lffit.splitgen(gen) print len(efflens), len(telens) pthl = pthk # pdt.savedata(lth,ltn,ngt,efflens,telens,P,SEED,pthl,wfitn) if telens or efflens: pdt.savedata(lth, ltn, ngt, efflens, telens, P, SEED, pthl, wfitn, trns) if telens and efflens: pdt.datadisplay(lth, ltn, ngt, efflens, telens, 500, 50, P, SEED, pthl, wfitn) del lth, ltn, ngt, pqr, info, wfitn, trns gc.collect() else: print("NOT-DATA-TO-SHOW") gc.collect() print("CARRY ON") del gen, gold gen = {} #gc={} gold = {} gc.collect()
def evolx(P,tn,PL1,PL2,gold,rk,nj,Qi,ptj,qij,rnz): MU=[P[0],P[1],P[2],P[3],P[4],P[5],P[6],P[7],P[8],P[9],P[10],P[11],P[12]] Ntot=P[13] pqr={} ltn=[] nefft=[] lnefft=[] ntest=[] lntest=[] lth=[] ngt=[] fitn=[] trns=[] flg='all-good' gen={} trs={} k=0 for i in gold.keys(): gl=[] if gold[i][0]=='te': gl=[gold[i][0],gold[i][1]] gen[k]=gl k+=1 if gold[i][0]=='eff': if nj==0: gl=[gold[i][0],gold[i][1],gold[i][2]] else: qj=rk.uniform_pos() snew=gold[i][2]+(0.01*( 1.0-2.0*rk.uniform_pos() )) if ((qj<Qi) and (snew>0.0)): gl=[gold[i][0],gold[i][1],snew] else: gl=[gold[i][0],gold[i][1],0.0] gen[k]=gl k+=1 gold=None for nk in range(tn): trs=Lffit.trates(gen,MU) rij,sr =Lffit.montec(trs,rk) trns.append(rij) if sr!='TRUE': pqr[nk]=[rij[0],rij[1],gen] ltn.append(Lffit.lent(gen)) lth.append(Lth) wq=Lffit.ft(gen) fitn.append(1.0-(wq/(wq+wo))) ngt.append(len(gen.keys())) n1x=0.0 n2x=0.0 m1x=0.0 m2x=0.0 for gk in gen.keys(): if gen[gk][0]=='eff': n1x+=1 m1x+=gen[gk][1] if gen[gk][0]=='te': n2x+=1 m2x+=gen[gk][1] nefft.append(n1x) lnefft.append(m1x) ntest.append(n2x) lntest.append(m2x) gu={} gu=LffitB.transform(rij,gen,rk, PL1, PL2) gen={} for i in gu.keys(): gen[i]=gu[i] else: flg='sumzero' break if len(gen.keys())>2000: flg='limit-reached' break #################### import os po='../RUN' nx=str(rnz)+'/n'+str(ptj) if not os.path.exists(po+nx): os.makedirs(po+nx) efflens, telens =Lffit.splitgen(gen) print len(efflens), len(telens) pthl=po+nx if telens or efflens: pdt.savedata(lth,ltn,ngt,efflens,telens,P,SEED,pthl,fitn,trns,nefft,ntest,lnefft,lntest,tn,JMPS) qij.put(gen)