def job(PT, SEED, nte, neff, tn, RXCLRPTS, CRKPTS, TES, nj, gold, Qi):
from math import floor
from pygsl import rng as rn
import sys
import gc
sys.path.append('../LNKMODEL/')
import Lffit
import LffitB
Np = PT[0]
m1 = PT[1] #hgt effs
m2 = PT[2] #hgt te
m3 = PT[3] #eff recomb
m4 = PT[4] #te dup
m5 = PT[5] #eff dup+te
m6 = PT[6] #eff->null
m7 = PT[7] #null ->0
m8 = PT[8] ##eff->0
beta1 = PT[9]
beta2 = PT[10]
wo = PT[11]
#print("parameter vector:"% PT)
av1 = 0.5 * (RXCLRPTS[2] - RXCLRPTS[1])
av2 = 0.5 * (CRKPTS[2] - CRKPTS[1])
av3 = 0.5 * (TES[2] - TES[1])
Lth = floor(
(RXCLRPTS[0] + CRKPTS[0] + TES[0]) * (1.0 / 3.0) * (av1 + av2 + av3))
#print Lth
MU = [Np, m1, m2, m3, m4, m5, m6, m7, m8, Lth, beta1, beta2, wo]
Qj = Qi
#SEED=987654320
rk = rn.rng()
rk.set(SEED)
################
if nj == 0:
gen = {}
trs = {}
gen = Lffit.inigenome(nte, neff, rk, [TES[1], TES[2]],
[CRKPTS[1], RXCLRPTS[2]])
#print gen
#raw_input()
################
if nj > 0:
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':
qj = rk.uniform_pos()
snew = gold[i][2] + (0.01 * (1.0 - 2.0 * rk.uniform_pos()))
if ((qj > Qj) 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
# print gen
# raw_input()
################
Ntot = RXCLRPTS[0] + CRKPTS[0] + TES[0]
pqr = {}
ltn = []
lth = []
ngt = []
fitn = []
trns = []
flg = 'all-good'
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)
#fitn.append(Lffit.ft(gen))
wq = Lffit.ft(gen)
fitn.append(1.0 - (wq / (wq + wo)))
ngt.append(len(gen.keys()))
#del gu
gu = {}
gu = LffitB.transform(rij, gen, rk, [TES[1], TES[2]],
[CRKPTS[1], RXCLRPTS[2]])
del gen
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
#print("JOB COMPLETED")
#print flg
gc.collect()
return [lth, ltn, ngt, pqr, gen, flg, fitn, trns]
av3 = 0.5 * (TES[2] - TES[1])
Lth = floor(
(RXCLRPTS[0] + CRKPTS[0] + TES[0]) * (1.0 / 3.0) * (av1 + av2 + av3))
#print Lth
MU = [Np, m1, m2, m3, m4, m5, m6, m7, m8, Lth, beta1, beta2, wo]
Qj = Qi
#SEED=987654320
rk = rn.rng()
rk.set(SEED)
for nj in range(2):
################
if nj == 0:
gen = {}
trs = {}
gen = Lffit.inigenome(nte, neff, rk, [TES[1], TES[2]],
[CRKPTS[1], RXCLRPTS[2]])
#print gen
#raw_input()
################
if nj > 0:
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':
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)
av2=0.5*(CRKPTS[2]-CRKPTS[1])
av3=0.5*(TES[2]-TES[1])
Lth=floor((RXCLRPTS[0]+CRKPTS[0]+TES[0])*(1.0/3.0)*(av1+av2+av3))
#print Lth
MU=[Np,m1,m2,m3,m4,m5,m6,m7,m8,Lth,beta1,beta2,wo]
Qj=Qi
#SEED=987654320
rk=rn.rng()
rk.set(SEED)
for nj in range(2):
################
if nj==0:
gen={}
trs={}
gen=Lffit.inigenome(nte,neff,rk,[TES[1],TES[2]],[CRKPTS[1], RXCLRPTS[2]])
#print gen
#raw_input()
################
if nj>0:
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':
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)
RUNS = args.r tn = args.DT JMPS = args.j SEED = args.seed m1 = args.probability1 m2 = args.probability2 m3 = args.probability3 m4 = args.probability4 m5 = args.probability5 m6 = args.probability6 m7 = args.probability7 m8 = args.probability8 rxclr = [] crinckler = [] tes = [] rxclr = Lffit.loaddata(pth + 'rxclr.dat') crinckler = Lffit.loaddata(pth + 'crinckler.dat') tes = Lffit.loaddata(pth + 'tes.dat') print len(rxclr), len(crinckler), len(tes) RXCLRPTS = [len(rxclr), min(rxclr), max(rxclr)] CRKPTS = [len(crinckler), min(crinckler), max(crinckler)] TES = [len(tes), min(tes), max(tes)] tes = None rxclr = None crinckler = None print RXCLRPTS, CRKPTS, TES av1 = 0.5 * (RXCLRPTS[2] - RXCLRPTS[1]) av2 = 0.5 * (CRKPTS[2] - CRKPTS[1]) av3 = 0.5 * (TES[2] - TES[1])
from math import floor from pygsl import rng as rn import gc import sys #pth='../' #sys.path.append(pth) import Lffit import os rxclr=[] crinckler=[] tes=[] pth='../../' rxclr=Lffit.loaddata(pth+'rxclr.dat') crinckler=Lffit.loaddata(pth+'crinckler.dat') tes=Lffit.loaddata(pth+'tes.dat') print len(rxclr), len(crinckler), len(tes) #raw_input() RXCLRPTS=[len(rxclr),min(rxclr),max(rxclr)] CRKPTS=[len(crinckler),min(crinckler),max(crinckler)] TES=[len(tes),min(tes),max(tes)] print RXCLRPTS print CRKPTS print TES import lmainjumps import pdt SEED=987654320 nte=5
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)
tn=args.DT ###* JMPS=args.j ###* SEED=args.seed ###* m1=args.probability1 ###* m2=args.probability2 ###* m3=args.probability3 ###* m4=args.probability4 ###* m5=args.probability5 ###* m6=args.probability6 ###* m7=args.probability7 ###* m8=args.probability8 ###* rxclr=[] crinckler=[] tes=[] ###*pth='../../' rxclr=Lffit.loaddata(pth+'rxclr.dat') crinckler=Lffit.loaddata(pth+'crinckler.dat') tes=Lffit.loaddata(pth+'tes.dat') print len(rxclr), len(crinckler), len(tes) RXCLRPTS=[len(rxclr),min(rxclr),max(rxclr)] CRKPTS=[len(crinckler),min(crinckler),max(crinckler)] TES=[len(tes),min(tes),max(tes)] tes= None rxclr=None crinckler=None print RXCLRPTS, CRKPTS, TES av1=0.5*(RXCLRPTS[2]-RXCLRPTS[1]) av2=0.5*(CRKPTS[2]-CRKPTS[1]) av3=0.5*(TES[2]-TES[1])
from math import floor from pygsl import rng as rn import gc import sys #pth='../' #sys.path.append(pth) import Lffit import os rxclr = [] crinckler = [] tes = [] pth = '../../' rxclr = Lffit.loaddata(pth + 'rxclr.dat') crinckler = Lffit.loaddata(pth + 'crinckler.dat') tes = Lffit.loaddata(pth + 'tes.dat') print len(rxclr), len(crinckler), len(tes) #raw_input() RXCLRPTS = [len(rxclr), min(rxclr), max(rxclr)] CRKPTS = [len(crinckler), min(crinckler), max(crinckler)] TES = [len(tes), min(tes), max(tes)] print RXCLRPTS print CRKPTS print TES import lmainjumps import pdt SEED = 987654320 nte = 5
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)
