commented out skopt in optimizers

pull/1/head
titusquah 5 years ago
parent 159a5c8140
commit cab58fda25

@ -2,15 +2,8 @@
<project version="4"> <project version="4">
<component name="ChangeListManager"> <component name="ChangeListManager">
<list default="true" id="f4439dc0-6756-4612-8f7d-596d8949f300" name="Default Changelist" comment=""> <list default="true" id="f4439dc0-6756-4612-8f7d-596d8949f300" name="Default Changelist" comment="">
<change afterPath="$PROJECT_DIR$/data/csvs/may_pitzer_params.csv" afterDir="false" />
<change afterPath="$PROJECT_DIR$/data/csvs/zeroes_removed_PC88A_HCL_NdPrCeLaDySmY.csv" afterDir="false" />
<change afterPath="$PROJECT_DIR$/data/xmls/PC88A_HCL_NdPrCeLaDySmY_w_pitzer.xml" afterDir="false" />
<change afterPath="$PROJECT_DIR$/docs/Examples/iterative_fitter.py" afterDir="false" />
<change afterPath="$PROJECT_DIR$/llepe/objectives.py" afterDir="false" />
<change afterPath="$PROJECT_DIR$/llepe/optimizers.py" afterDir="false" />
<change beforePath="$PROJECT_DIR$/.idea/workspace.xml" beforeDir="false" afterPath="$PROJECT_DIR$/.idea/workspace.xml" afterDir="false" /> <change beforePath="$PROJECT_DIR$/.idea/workspace.xml" beforeDir="false" afterPath="$PROJECT_DIR$/.idea/workspace.xml" afterDir="false" />
<change beforePath="$PROJECT_DIR$/data/xmls/twophase.xml" beforeDir="false" afterPath="$PROJECT_DIR$/data/xmls/twophase.xml" afterDir="false" /> <change beforePath="$PROJECT_DIR$/llepe/optimizers.py" beforeDir="false" afterPath="$PROJECT_DIR$/llepe/optimizers.py" afterDir="false" />
<change beforePath="$PROJECT_DIR$/llepe/__init__.py" beforeDir="false" afterPath="$PROJECT_DIR$/llepe/__init__.py" afterDir="false" />
</list> </list>
<option name="SHOW_DIALOG" value="false" /> <option name="SHOW_DIALOG" value="false" />
<option name="HIGHLIGHT_CONFLICTS" value="true" /> <option name="HIGHLIGHT_CONFLICTS" value="true" />
@ -269,7 +262,14 @@
<option name="project" value="LOCAL" /> <option name="project" value="LOCAL" />
<updated>1594739749122</updated> <updated>1594739749122</updated>
</task> </task>
<option name="localTasksCounter" value="14" /> <task id="LOCAL-00014" summary="added example file to do iterative optimization">
<created>1594767749286</created>
<option name="number" value="00014" />
<option name="presentableId" value="LOCAL-00014" />
<option name="project" value="LOCAL" />
<updated>1594767749286</updated>
</task>
<option name="localTasksCounter" value="15" />
<servers /> <servers />
</component> </component>
<component name="Vcs.Log.Tabs.Properties"> <component name="Vcs.Log.Tabs.Properties">
@ -293,7 +293,8 @@
<MESSAGE value="Edited parity plot to allow color to represent 3rd dimension. Still need to improve colorbar axis name." /> <MESSAGE value="Edited parity plot to allow color to represent 3rd dimension. Still need to improve colorbar axis name." />
<MESSAGE value="updated to LLEPE" /> <MESSAGE value="updated to LLEPE" />
<MESSAGE value="fixed setup file by changing LLEPE to llepe" /> <MESSAGE value="fixed setup file by changing LLEPE to llepe" />
<option name="LAST_COMMIT_MESSAGE" value="fixed setup file by changing LLEPE to llepe" /> <MESSAGE value="added example file to do iterative optimization" />
<option name="LAST_COMMIT_MESSAGE" value="added example file to do iterative optimization" />
</component> </component>
<component name="WindowStateProjectService"> <component name="WindowStateProjectService">
<state x="-1330" y="212" key="#Inspections" timestamp="1590787654691"> <state x="-1330" y="212" key="#Inspections" timestamp="1590787654691">
@ -313,11 +314,11 @@
<screen x="-1920" y="2" width="1920" height="1040" /> <screen x="-1920" y="2" width="1920" height="1040" />
</state> </state>
<state x="-1213" y="379" key="ANALYSIS_DLG_com.intellij.analysis.BaseAnalysisAction$1/0.0.1536.824/-1920.2.1920.1040@-1920.2.1920.1040" timestamp="1590787657711" /> <state x="-1213" y="379" key="ANALYSIS_DLG_com.intellij.analysis.BaseAnalysisAction$1/0.0.1536.824/-1920.2.1920.1040@-1920.2.1920.1040" timestamp="1590787657711" />
<state x="169" y="115" key="CommitChangelistDialog2" timestamp="1594739747678"> <state x="169" y="115" key="CommitChangelistDialog2" timestamp="1594767747695">
<screen x="0" y="0" width="1920" height="1040" /> <screen x="0" y="0" width="1920" height="1040" />
</state> </state>
<state x="-1364" y="117" key="CommitChangelistDialog2/0.0.1536.824/-1920.2.1920.1040@-1920.2.1920.1040" timestamp="1592321237444" /> <state x="-1364" y="117" key="CommitChangelistDialog2/0.0.1536.824/-1920.2.1920.1040@-1920.2.1920.1040" timestamp="1592321237444" />
<state x="169" y="115" key="CommitChangelistDialog2/0.0.1920.1040@0.0.1920.1040" timestamp="1594739747678" /> <state x="169" y="115" key="CommitChangelistDialog2/0.0.1920.1040@0.0.1920.1040" timestamp="1594767747695" />
<state x="-1828" y="94" width="1736" height="856" key="DiffContextDialog" timestamp="1591048879404"> <state x="-1828" y="94" width="1736" height="856" key="DiffContextDialog" timestamp="1591048879404">
<screen x="-1920" y="2" width="1920" height="1040" /> <screen x="-1920" y="2" width="1920" height="1040" />
</state> </state>
@ -367,11 +368,11 @@
<screen x="-1920" y="2" width="1920" height="1040" /> <screen x="-1920" y="2" width="1920" height="1040" />
</state> </state>
<state x="-1245" y="438" width="437" height="168" key="VCS.ChangelistChooser/0.0.1536.824/-1920.2.1920.1040@-1920.2.1920.1040" timestamp="1592321346068" /> <state x="-1245" y="438" width="437" height="168" key="VCS.ChangelistChooser/0.0.1536.824/-1920.2.1920.1040@-1920.2.1920.1040" timestamp="1592321346068" />
<state x="552" y="254" key="Vcs.Push.Dialog.v2" timestamp="1594739751769"> <state x="552" y="254" key="Vcs.Push.Dialog.v2" timestamp="1594767758886">
<screen x="0" y="0" width="1920" height="1040" /> <screen x="0" y="0" width="1920" height="1040" />
</state> </state>
<state x="-1368" y="256" key="Vcs.Push.Dialog.v2/0.0.1536.824/-1920.2.1920.1040@-1920.2.1920.1040" timestamp="1592321242949" /> <state x="-1368" y="256" key="Vcs.Push.Dialog.v2/0.0.1536.824/-1920.2.1920.1040@-1920.2.1920.1040" timestamp="1592321242949" />
<state x="552" y="254" key="Vcs.Push.Dialog.v2/0.0.1920.1040@0.0.1920.1040" timestamp="1594739751769" /> <state x="552" y="254" key="Vcs.Push.Dialog.v2/0.0.1920.1040@0.0.1920.1040" timestamp="1594767758886" />
<state x="263" y="182" width="1283" height="717" key="com.intellij.ide.util.TipDialog" timestamp="1594694030135"> <state x="263" y="182" width="1283" height="717" key="com.intellij.ide.util.TipDialog" timestamp="1594694030135">
<screen x="0" y="0" width="1920" height="1040" /> <screen x="0" y="0" width="1920" height="1040" />
</state> </state>

@ -1,6 +1,6 @@
import scipy.optimize as scipy_opt import scipy.optimize as scipy_opt
from scipy.optimize import minimize from scipy.optimize import minimize
import skopt # import skopt
def dual_anneal_optimizer(objective, x_guess): def dual_anneal_optimizer(objective, x_guess):
bounds = [(1e-1, 1e1)] * len(x_guess) bounds = [(1e-1, 1e1)] * len(x_guess)
@ -21,19 +21,19 @@ def diff_evo_optimizer(objective, x_guess):
return est_parameters, res.fun return est_parameters, res.fun
def forest_lbfgsb_optimizer(objective, x_guess): # def forest_lbfgsb_optimizer(objective, x_guess):
x_guess = list(x_guess) # x_guess = list(x_guess)
bounds = [(1e-1, 1e1)]*len(x_guess) # bounds = [(1e-1, 1e1)]*len(x_guess)
bounds[1] = (1e-1, 2) # bounds[1] = (1e-1, 2)
res = skopt.forest_minimize(objective, # res = skopt.forest_minimize(objective,
bounds, # bounds,
random_state=1, # random_state=1,
acq_func='LCB', # acq_func='LCB',
n_random_starts=30, # n_random_starts=30,
x0=x_guess, # x0=x_guess,
xi=1e-4) # xi=1e-4)
x_guess = res.x # x_guess = res.x
optimizer_kwargs = {"method": 'l-bfgs-b', # optimizer_kwargs = {"method": 'l-bfgs-b',
"bounds": bounds} # "bounds": bounds}
res = minimize(objective, x_guess, **optimizer_kwargs) # res = minimize(objective, x_guess, **optimizer_kwargs)
return res.x, res.fun # return res.x, res.fun

Loading…
Cancel
Save