def initTriCycle(IGEElement,doPickle=True,samples=10):
fitHist=param.fitHistory()
lumInterval=runLumCycle(maxIter=5)
curParam=param.param()
curParam['lum']=lumInterval['suggestValue']
intervals={'luminterval':lumInterval['interval'],
'vphinterval':initialize.getVphBounds(),
#make sure that no other metal really shoots up
'igeinterval':initialize.getElementBounds(IGEElement,curParam.comp)}
#pdb.set_trace()
lumRange=np.linspace(intervals['luminterval'][0],intervals['luminterval'][1],num=samples)
vphRange=np.linspace(intervals['vphinterval'][0],intervals['vphinterval'][1],num=samples)
IGERange=np.linspace(intervals['igeinterval'][0],intervals['igeinterval'][1],num=samples)
lumMG,vphMG,IGEMG,curParamMG=launcherSteps.launchTriCycle(lumRange,vphRange,IGERange,IGEElement,initParam=curParam)
if doPickle:
pickle.dump(lumMG,file('lumMG0.pkl','w'))
pickle.dump(vphMG,file('vphMG0.pkl','w'))
pickle.dump(IGEMG,file('igeMG0.pkl','w'))
evalTriCycle(lumMG,vphMG,IGEMG,curParamMG,curParam,intervals,fitHist,IGEElement,mode='init')
return fitHist,curParam,intervals
def createRandomParam(randomParamFunc=None):
#creating new paramObject
lumLimits = config.GAConfDict['lumLimits']
vphLimits = config.GAConfDict['vphLimits']
randomParam=param.param()
if config.GAConfDict['lockTiCr']:
randomParam.comp.lockTiCr=True
if config.GAConfDict['lockScTi']:
randomParam.comp.lockScTi=True
if config.GAConfDict['lockVCr']:
randomParam.comp.lockVCr=True
if randomParamFunc==createRandomLogNormalValueBestFit:
fname=glob('*.bf.pkl')[0]
bestFitParam=cPickle.load(file(fname))
#Limiting luminosity and photospheric velocity near best fit value, commented out atm
#randomParam['lum']=random.uniform(bestFitParam['lum']-0.1,bestFitParam['lum']+0.1)
#no limit for lum parameter
#randomParam['lum']=bestFitParam['lum']
randomParam['lum']=random.uniform(*lumLimits)
else:
randomParam['lum']=random.uniform(*lumLimits)
randomParam['vph']=random.uniform(*vphLimits)
rChoice=random.permutation(config.GAConfDict['selElements'])
for element in rChoice:
bounds=initialize.getElementBounds(element,randomParam.comp)
curAbundance=randomParam[element]
#newAbundance=random.uniform(bounds[0],bounds[1])
newAbundance=randomParamFunc(curAbundance,element,bounds)
if any(np.array(bounds)<0): raise Exception('Negative Bounds')
randomParam[element]=newAbundance
randomParam.comp.resetOxygen()
if randomParam['O']<0: pdb.set_trace()
if randomParam[element]<0: pdb.set_trace()
if randomParam.comp.data.has_key('element'): raise Exception()
#randomParam['Ca']=0.01
#randomParam['vph']=11700
return randomParam
def evalTriCycle(lumMG,vphMG,IGEMG,curParamMG,curParam,intervals,fitHist,IGEElement,mode='std'):
#IMPORTANT!!!!!!!! the function DOES NOT need curParamMG, only for adding it to the fithist object
#REMOVE AS SOON AS DEBUGGING for
pp=pprint.PrettyPrinter()
selectSteadyThresh=6
steadyThresh=4
pp.pprint(intervals)
pivots=[0,0,0]
lumBounds=initialize.getLumBounds()
vphBounds=initialize.getVphBounds()
IGEBounds=initialize.getElementBounds(IGEElement,curParam.comp)
lumInterval=dalekSteps.getNextLumInterval(lumMG)
vphInterval=dalekSteps.getNextVphInterval(vphMG)
IGEInterval=dalekSteps.getNextIGEInterval(IGEMG,IGEElement)
lastSuggestLum=np.mean(intervals['luminterval'])
lastSuggestVph=np.mean(intervals['vphinterval'])
lastSuggestIGE=np.mean(intervals['igeinterval'])
print "CurLum=%.4f CurVph=%.2f CurFe0=%.7f"%(curParam['lum'],curParam['vph'],curParam['Fe0'])
print "---------------------------------"
print "LUM:"
pp.pprint(lumInterval)
print
print "VPH:"
pp.pprint(vphInterval)
print
print "IGE:"
pp.pprint(IGEInterval)
#Updating the intervals if the solution is not steady
steadyIncFactor=2
selectSteadyIncFactor=2
curLum=curParam['lum']
closestLumID=np.argmin(lumMG['lum']-curLum)
#intervals=[lumInterval,vphInterval,IGEInterval]
merits=np.argsort([item[2] for item in intervals])
#pivots[merits[2]]+=3
#pivots[merits[1]]+=1
pivots[0]=3
pivots[1]=2
pivots[2]=1
#making sure that we are close in luminosity
print IGEMG['Fe0']
if not (abs((lumMG['lum'][closestLumID]-curLum))<0.05 and abs(lumInterval['merits'][closestLumID])<0.1):
print "lum doesnt allow other parameters atm"
print "closestLum %s closestMerit %s"%(lumMG['lum'][closestLumID]-curLum,lumInterval['merits'][closestLumID])
pivots[1]-3
pivots[2]-3
pivots[0]=10
#Checking if the interval is too small already
if lumInterval['mininterval']: pivots[0]-=5
if vphInterval['mininterval']: pivots[1]-=5
#IGE NOT WORKING properly
#pivots[2]=-10
print "PIVOTS %s"%pivots
#Checking w factors:
print 'vph w %s'%vphMG['w']
print 'ige w %s'%IGEMG['w']
if np.sum(np.abs(vphMG['w']-0.5)>0.1)>4: pivots[1]-=1.5
if np.sum(np.abs(IGEMG['w']-0.5)>0.1)>4: pivots[2]-=1
damping=True
#the factor does not do that much atm BEWARE
dampingFactor=2
if np.argmax(pivots)==0:
#param=copy.copy(lumMG.grid[lumInterval['bestFitID']].param)
if lumInterval['steady']>selectSteadyThresh:
if damping:
curParam['lum']=np.mean([lumInterval['suggestValue'],lastSuggestLum])
intervals.update({'luminterval':[curParam['lum']+i*dampingFactor*lumInterval['dev'] for i in [-1,1]]})
else:
curParam['lum']=lumInterval['suggestValue']
intervals.update({'luminterval':lumInterval['interval']})
else:
print "Not steady solution will increase last interval and try again"
dev=selectSteadyIncFactor*np.abs(intervals['luminterval'][0]-curParam['lum'])
intervals.update({'luminterval':[curParam['lum']-dev,curParam['lum']+dev]})
print "Chose Lum: updated for next run with suggestValue %s and Interval %s"%(curParam['lum'],intervals['luminterval'])
else:
if lumInterval['steady']<steadyThresh:
print "Updating the lum interval as the current solution is not steady"
dev=steadyIncFactor*np.abs(intervals['luminterval'][0]-curParam['lum'])
intervals.update({'luminterval':[curParam['lum']-dev,curParam['lum']+dev]})
if np.argmax(pivots)==1:
#param=copy.copy(vphMG.grid[vphInterval['bestFitID']].param)
if vphInterval['steady']>selectSteadyThresh:
if damping:
curParam['vph']=np.mean([vphInterval['suggestValue'],lastSuggestVph])
intervals.update({'vphinterval':[curParam['vph']+i*dampingFactor*vphInterval['dev'] for i in [-1,1]]})
else:
curParam['vph']=vphInterval['suggestValue']
intervals.update({'vphinterval':vphInterval['interval']})
else:
print "Not steady solution will increase last interval and try again"
dev=selectSteadyIncFactor*np.abs(intervals['vphinterval'][0]-curParam['vph'])
intervals.update({'vphinterval':[curParam['vph']-dev,curParam['vph']+dev]})
print "Chose Vph: updated for next run with suggestValue %s and Interval %s"%(curParam['vph'],intervals['vphinterval'])
else:
if vphInterval['steady']<steadyThresh:
print "Updating the vph interval as the current solution is not steady"
dev=steadyIncFactor*np.abs(intervals['vphinterval'][0]-curParam['vph'])
intervals.update({'vphinterval':[curParam['vph']-dev,curParam['vph']+dev]})
if np.argmax(pivots)==2:
if IGEInterval['steady']>selectSteadyThresh:
if damping:
curParam[IGEElement]=np.mean([IGEInterval['suggestValue'],lastSuggestIGE])
intervals.update({'igeinterval':[curParam[IGEElement]+i*dampingFactor*IGEInterval['dev'] for i in [-1,1]]})
else:
curParam[IGEElement]=IGEInterval['suggestValue']
intervals.update({'igeinterval':IGEInterval['interval']})
else:
print "Not steady solution will increase last interval and try again"
dev=selectSteadyIncFactor*np.abs(intervals['igeinterval'][0]-curParam[IGEElement])
intervals.update({'igeinterval':[curParam[IGEElement]-dev,curParam[IGEElement]+dev]})
print "Chose IGE: updated for next run with suggestValue %s and Interval %s"%(curParam[IGEElement],intervals['igeinterval'])
else:
if IGEInterval['steady']<steadyThresh:
print "Updating the IGE interval as the current solution is not steady"
dev=steadyIncFactor*np.abs(intervals['igeinterval'][0]-curParam[IGEElement])
intervals.update({'igeinterval':[curParam[IGEElement]-dev,curParam[IGEElement]+dev]})
#Checking for crossing Boundaries:
boundCheck=zip(lumBounds,intervals['luminterval'])
lumCheck=[np.max(boundCheck[0]),np.min(boundCheck[1])]
boundCheck=zip(vphBounds,intervals['vphinterval'])
vphCheck=[np.max(boundCheck[0]),np.min(boundCheck[1])]
boundCheck=zip(IGEBounds,intervals['igeinterval'])
igeCheck=[np.max(boundCheck[0]),np.min(boundCheck[1])]
intervals.update({'luminterval':lumCheck,'vphinterval':vphCheck,'igeinterval':igeCheck})
#Readjusting suggest values to be in the middle of the intervals:
curParam['lum']=np.mean(intervals['luminterval'])
curParam['vph']=np.mean(intervals['vphinterval'])
curParam[IGEElement]=np.mean(intervals['igeinterval'])
#Printing the intervals:
pp.pprint(intervals)
#Adding to fitHist
fitHist.addHistItem([lumInterval,vphInterval,IGEInterval],[lumMG,vphMG,IGEMG,curParamMG],pivots,intervals,reallyAdd=False,saveSingle=True)
#pdb.set_trace()
#Break condition
if np.argmax(pivots)>0:
# pdb.set_trace()
return True
else:
return False
