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)
