% Starting Algorithm: 
%	className: org.openda.algorithms.Dud
%	dir.: /home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./algorithm
%	config.: dudAlgorithm.xml
%  Algorithm  configstring = dudAlgorithm.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./algorithm/dudAlgorithm.xml
%  Algorithm  Dud initialized
% Application initializing finished
% Initializing application
%  Algorithm  Dud preparation
%  Algorithm  Retrieving initial parameters from model
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model2/OscillatorStochModel.xml
% simulationTimespan =[0.0,0.05,10.0]
% parameters@names =t_damp,omega parameters=[8.0,1.5708]
% parameterUncertainty@names =t_damp,omega parameterUncertainty=[1.0,0.1257]
% systemNoise = {[0.0,0.0],[0.3,0.3]}
% initialState = [0.8,0.0]
% initialStateUncertainty = [0.8,0.8]
%  Algorithm  Starting optimizer
%  Algorithm  costFunction@class=org.openda.algorithms.SimulationKwadraticCostFunction
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model2/OscillatorStochModel.xml
% simulationTimespan =[0.0,0.05,10.0]
% parameters@names =t_damp,omega parameters=[8.0,1.5708]
% parameterUncertainty@names =t_damp,omega parameterUncertainty=[1.0,0.1257]
% systemNoise = {[0.0,0.0],[0.3,0.3]}
% initialState = [0.8,0.0]
% initialStateUncertainty = [0.8,0.8]
%  Algorithm  costFunction@weakParameterConstraint=false
%  Algorithm  costFunction@factor=0.5
% ========================================================================
% no1
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model2/OscillatorStochModel.xml
% simulationTimespan =[0.0,0.05,10.0]
% parameters@names =t_damp,omega parameters=[8.0,1.5708]
% parameterUncertainty@names =t_damp,omega parameterUncertainty=[1.0,0.1257]
% systemNoise = {[0.0,0.0],[0.3,0.3]}
% initialState = [0.8,0.0]
% initialStateUncertainty = [0.8,0.8]
evaluatedParameters{1}	=[8.0,1.5708];
predicted{1}	=[0.8,0.059874050823311174,-0.622514741258897,-0.05209975910799278,0.4843573470865686,0.04483119930100093,-0.3768240910010292,-0.03822558598991682,0.29313512040223755,0.03234578622602441,-0.2280097266602124];
observed{1}	=[0.8,-0.11351446579670872,-0.5871715953189904,0.32299872426993737,0.3357367607855463,-0.3886669488572505,-0.11092801803286569,0.3512503605862144,-0.05188139665849537,-0.2569942428910992,0.14310793968854146];
residuals{1}	=[0.0,-0.1733885166200199,0.03534314593990662,0.37509848337793017,-0.14862058630102226,-0.43349814815825144,0.26589607296816353,0.3894759465761312,-0.34501651706073294,-0.2893400291171236,0.37111766634875387];
obs_stdev{1}	=[0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1];
costObserved{1}	=47.24474586868815;
costTotal{1}	=47.24474586868815;
% SimulationKwadraticCostFunction: evaluation 1 : cost = 47.245
% Evaluating with parameters  [9,1.571]
% ========================================================================
% no2
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model2/OscillatorStochModel.xml
% simulationTimespan =[0.0,0.05,10.0]
% parameters@names =t_damp,omega parameters=[8.0,1.5708]
% parameterUncertainty@names =t_damp,omega parameterUncertainty=[1.0,0.1257]
% systemNoise = {[0.0,0.0],[0.3,0.3]}
% initialState = [0.8,0.0]
% initialStateUncertainty = [0.8,0.8]
evaluatedParameters{2}	=[9.0,1.5708];
predicted{2}	=[0.8,0.0535786279459891,-0.6402131222275624,-0.04740700352235454,0.5123091801872715,0.041565211398558684,-0.4099327838549162,-0.03616714183080783,0.32799413193914145,0.03126817887015455,-0.26241723932575955];
observed{2}	=[0.8,-0.11351446579670872,-0.5871715953189904,0.32299872426993737,0.3357367607855463,-0.3886669488572505,-0.11092801803286569,0.3512503605862144,-0.05188139665849537,-0.2569942428910992,0.14310793968854146];
residuals{2}	=[0.0,-0.16709309374269782,0.053041526908572,0.3704057277922919,-0.17657241940172513,-0.4302321602558092,0.29900476582205054,0.3874175024170222,-0.37987552859763685,-0.2882624217612537,0.405525179014301];
obs_stdev{2}	=[0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1];
costObserved{2}	=50.77794613756578;
costTotal{2}	=50.77794613756578;
% SimulationKwadraticCostFunction: evaluation 2 : cost = 50.778
% Evaluating with parameters  [8,1.696]
% ========================================================================
% no3
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model2/OscillatorStochModel.xml
% simulationTimespan =[0.0,0.05,10.0]
% parameters@names =t_damp,omega parameters=[8.0,1.5708]
% parameterUncertainty@names =t_damp,omega parameterUncertainty=[1.0,0.1257]
% systemNoise = {[0.0,0.0],[0.3,0.3]}
% initialState = [0.8,0.0]
% initialStateUncertainty = [0.8,0.8]
evaluatedParameters{3}	=[8.0,1.6965];
predicted{3}	=[0.8,-0.03350311927911278,-0.6158978936059812,0.15740421970211982,0.4461027450327285,-0.21769751450176725,-0.30101131585066443,0.23371991318597032,0.1845979073251928,-0.22137835983632012,-0.09656628411831196];
observed{3}	=[0.8,-0.11351446579670872,-0.5871715953189904,0.32299872426993737,0.3357367607855463,-0.3886669488572505,-0.11092801803286569,0.3512503605862144,-0.05188139665849537,-0.2569942428910992,0.14310793968854146];
residuals{3}	=[0.0,-0.08001134651759595,0.028726298286990803,0.16559450456781755,-0.11036598424718219,-0.17096943435548326,0.19008329781779876,0.11753044740024407,-0.23647930398368816,-0.03561588305477906,0.23967422380685344];
obs_stdev{3}	=[0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1];
costObserved{3}	=12.031975288088328;
costTotal{3}	=12.031975288088328;
% SimulationKwadraticCostFunction: evaluation 3 : cost = 12.032
%  Algorithm  outerLoop@maxIterations=10
%  Algorithm  outerLoop@absTolerance=0.01
%  Algorithm  outerLoop@relTolerance=0.01
%  Algorithm  lineSearch@maxIterations=5
%  Algorithm  lineSearch@maxRelStepSize=10.0
%  Algorithm  lineSearch/backtracking@startIterationNegativeLook=3
%  Algorithm  lineSearch/backtracking@shorteningFactor=0.5
% Application initialized
% Application starting next step
% ======================================================
% DUD outer iteration no.1
% ======================================================
% -----------------------------------------------------
costs{1}	=[12.031975288088328,47.24474586868815,50.77794613756578];
% -----------------------------------------------------
% Start search until improvement,
% Next try p=   [3.865,1.768]
%   delta_p=    [-4.135,0.071]
% ========================================================================
% no4
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model2/OscillatorStochModel.xml
% simulationTimespan =[0.0,0.05,10.0]
% parameters@names =t_damp,omega parameters=[8.0,1.5708]
% parameterUncertainty@names =t_damp,omega parameterUncertainty=[1.0,0.1257]
% systemNoise = {[0.0,0.0],[0.3,0.3]}
% initialState = [0.8,0.0]
% initialStateUncertainty = [0.8,0.8]
evaluatedParameters{4}	=[3.8645008782721773,1.7677550591129492];
predicted{4}	=[0.8,-0.019355758465837647,-0.4715034929423244,0.14037763815403573,0.24265249941609046,-0.14997034705642778,-0.10363814732967684,0.11770082218993014,0.029604764281122006,-0.07823831949868218,0.0037350287250344587];
observed{4}	=[0.8,-0.11351446579670872,-0.5871715953189904,0.32299872426993737,0.3357367607855463,-0.3886669488572505,-0.11092801803286569,0.3512503605862144,-0.05188139665849537,-0.2569942428910992,0.14310793968854146];
residuals{4}	=[0.0,-0.09415870733087107,-0.11566810237666603,0.18262108611590164,0.09308426136945586,-0.23869660180082272,-0.007289870703188844,0.23354953839628426,-0.08148616093961739,-0.178755923392417,0.139372910963507];
obs_stdev{4}	=[0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1];
costObserved{4}	=11.692659628836498;
costTotal{4}	=11.692659628836498;
% SimulationKwadraticCostFunction: evaluation 4 : cost = 11.693
% Error estimate for this outer iteration
parameterErrorEstimateStd{1}	=[1.5059320861162961,0.023763314363807318];
parameterErrorCorrelations{1}	=[0.9999999999999998,0.23509589433394898;0.235095894333949,0.9999999999999998];
% Application starting next step
% ======================================================
% DUD outer iteration no.2
% ======================================================
% -----------------------------------------------------
costs{2}	=[11.692659628836498,12.031975288088328,47.24474586868815];
% -----------------------------------------------------
% Start search until improvement,
% Next try p=   [5.176,1.842]
%   delta_p=    [1.311,0.075]
% ========================================================================
% no5
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model2/OscillatorStochModel.xml
% simulationTimespan =[0.0,0.05,10.0]
% parameters@names =t_damp,omega parameters=[8.0,1.5708]
% parameterUncertainty@names =t_damp,omega parameterUncertainty=[1.0,0.1257]
% systemNoise = {[0.0,0.0],[0.3,0.3]}
% initialState = [0.8,0.0]
% initialStateUncertainty = [0.8,0.8]
evaluatedParameters{5}	=[5.17557472191182,1.8424372539680713];
predicted{5}	=[0.8,-0.10329849888247057,-0.49955869208424153,0.2830931215654548,0.218789633563452,-0.2856003932685666,-0.026944508570375912,0.20554267686130054,-0.06929081669312111,-0.11013120070060539,0.09401780029742948];
observed{5}	=[0.8,-0.11351446579670872,-0.5871715953189904,0.32299872426993737,0.3357367607855463,-0.3886669488572505,-0.11092801803286569,0.3512503605862144,-0.05188139665849537,-0.2569942428910992,0.14310793968854146];
residuals{5}	=[0.0,-0.01021596691423815,-0.0876129032347489,0.03990560270448257,0.11694712722209433,-0.10306655558868388,-0.08398350946248978,0.14570768372491386,0.017409420034625735,-0.1468630421904938,0.04909013939111198];
obs_stdev{5}	=[0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1];
costObserved{5}	=4.3118915588890765;
costTotal{5}	=4.3118915588890765;
% SimulationKwadraticCostFunction: evaluation 5 : cost = 4.312
% Error estimate for this outer iteration
parameterErrorEstimateStd{2}	=[1.0038093253933913,0.03127719788351028];
parameterErrorCorrelations{2}	=[0.9999999999999998,-0.6742860472031978;-0.6742860472031978,1.0];
% Application starting next step
% ======================================================
% DUD outer iteration no.3
% ======================================================
% -----------------------------------------------------
costs{3}	=[4.3118915588890765,11.692659628836498,12.031975288088328];
% -----------------------------------------------------
% Start search until improvement,
% Next try p=   [8.306,1.839]
%   delta_p=    [3.13,-3.288E-3]
% ========================================================================
% no6
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model2/OscillatorStochModel.xml
% simulationTimespan =[0.0,0.05,10.0]
% parameters@names =t_damp,omega parameters=[8.0,1.5708]
% parameterUncertainty@names =t_damp,omega parameterUncertainty=[1.0,0.1257]
% systemNoise = {[0.0,0.0],[0.3,0.3]}
% initialState = [0.8,0.0]
% initialStateUncertainty = [0.8,0.8]
evaluatedParameters{6}	=[8.306052605167059,1.8391489653234216];
predicted{6}	=[0.8,-0.14044274535764256,-0.5637274140459242,0.3716132005244781,0.2708937538888926,-0.4176197930862733,-0.01940515678379439,0.33724479147237363,-0.14102226986655234,-0.19972935923991453,0.20339670712609276];
observed{6}	=[0.8,-0.11351446579670872,-0.5871715953189904,0.32299872426993737,0.3357367607855463,-0.3886669488572505,-0.11092801803286569,0.3512503605862144,-0.05188139665849537,-0.2569942428910992,0.14310793968854146];
residuals{6}	=[0.0,0.02692827956093384,-0.023444181273066245,-0.04861447625454074,0.06484300689665373,0.02895284422902278,-0.0915228612490713,0.014005569113840766,0.08914087320805697,-0.05726488365118465,-0.060288767437551294];
obs_stdev{6}	=[0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1];
costObserved{6}	=1.605684983098306;
costTotal{6}	=1.605684983098306;
% SimulationKwadraticCostFunction: evaluation 6 : cost = 1.606
% Error estimate for this outer iteration
parameterErrorEstimateStd{3}	=[1.1453403105440625,0.030063202978196105];
parameterErrorCorrelations{3}	=[1.0,-0.07766907058470954;-0.07766907058470926,1.0];
% Application starting next step
% ======================================================
% DUD outer iteration no.4
% ======================================================
% -----------------------------------------------------
costs{4}	=[1.605684983098306,4.3118915588890765,11.692659628836498];
% -----------------------------------------------------
% Start search until improvement,
% Next try p=   [9.31,1.76]
%   delta_p=    [1.004,-0.079]
% ========================================================================
% no7
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model2/OscillatorStochModel.xml
% simulationTimespan =[0.0,0.05,10.0]
% parameters@names =t_damp,omega parameters=[8.0,1.5708]
% parameterUncertainty@names =t_damp,omega parameterUncertainty=[1.0,0.1257]
% systemNoise = {[0.0,0.0],[0.3,0.3]}
% initialState = [0.8,0.0]
% initialStateUncertainty = [0.8,0.8]
evaluatedParameters{7}	=[9.309779897437522,1.7598917574736115];
predicted{7}	=[0.8,-0.08956415413710442,-0.6156216548542603,0.2765504598026084,0.40483376152019634,-0.3574364177884018,-0.2079523602994104,0.357348302676896,0.04916119027785541,-0.30458063218803605,0.06142104050357018];
observed{7}	=[0.8,-0.11351446579670872,-0.5871715953189904,0.32299872426993737,0.3357367607855463,-0.3886669488572505,-0.11092801803286569,0.3512503605862144,-0.05188139665849537,-0.2569942428910992,0.14310793968854146];
residuals{7}	=[0.0,-0.023950311659604304,0.02845005953526991,0.04644826446732897,-0.06909700073465003,-0.03123053106884871,0.0970243422665447,-0.006097942090681596,-0.10104258693635079,0.04758638929693687,0.08168689918497128];
obs_stdev{7}	=[0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1];
costObserved{7}	=1.8943966189554138;
costTotal{7}	=1.8943966189554138;
% SimulationKwadraticCostFunction: evaluation 7 : cost = 1.894
% Cost is worse! Reducing stepsize.
% Next try p=   [9.31,1.76]
% ========================================================================
% no8
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model2/OscillatorStochModel.xml
% simulationTimespan =[0.0,0.05,10.0]
% parameters@names =t_damp,omega parameters=[8.0,1.5708]
% parameterUncertainty@names =t_damp,omega parameterUncertainty=[1.0,0.1257]
% systemNoise = {[0.0,0.0],[0.3,0.3]}
% initialState = [0.8,0.0]
% initialStateUncertainty = [0.8,0.8]
evaluatedParameters{8}	=[8.807916251302291,1.7995203613985165];
predicted{8}	=[0.8,-0.11541839913490543,-0.5914918260570347,0.3277195720841197,0.3407233243616661,-0.3969483464963261,-0.11330228240009887,0.36147292872637,-0.053782776342122654,-0.26660990206730445,0.14911842698755542];
observed{8}	=[0.8,-0.11351446579670872,-0.5871715953189904,0.32299872426993737,0.3357367607855463,-0.3886669488572505,-0.11092801803286569,0.3512503605862144,-0.05188139665849537,-0.2569942428910992,0.14310793968854146];
residuals{8}	=[0.0,0.0019039333381967033,0.004320230738044217,-0.004720847814182361,-0.004986563576119762,0.008281397639075594,0.0023742643672331892,-0.010222568140155597,0.0019013796836272814,0.009615659176205271,-0.006010487299013961];
obs_stdev{8}	=[0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1];
costObserved{8}	=0.019018162867878868;
costTotal{8}	=0.019018162867878868;
% SimulationKwadraticCostFunction: evaluation 8 : cost = 0.019
% Error estimate for this outer iteration
parameterErrorEstimateStd{4}	=[1.4816870996540394,0.05219621791174397];
parameterErrorCorrelations{4}	=[1.0,-0.6598210574804025;-0.6598210574804031,1.0];
% Application starting next step
% ======================================================
% DUD outer iteration no.5
% ======================================================
% -----------------------------------------------------
costs{5}	=[0.019018162867878868,1.605684983098306,4.3118915588890765];
% -----------------------------------------------------
% Start search until improvement,
% Next try p=   [8.574,1.801]
%   delta_p=    [-0.234,1.512E-3]
% ========================================================================
% no9
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model2/OscillatorStochModel.xml
% simulationTimespan =[0.0,0.05,10.0]
% parameters@names =t_damp,omega parameters=[8.0,1.5708]
% parameterUncertainty@names =t_damp,omega parameterUncertainty=[1.0,0.1257]
% systemNoise = {[0.0,0.0],[0.3,0.3]}
% initialState = [0.8,0.0]
% initialStateUncertainty = [0.8,0.8]
evaluatedParameters{9}	=[8.573734713991831,1.8010324422307218];
predicted{9}	=[0.8,-0.11482288479492929,-0.587667093987248,0.325771759944244,0.33521449294400407,-0.39194738987748234,-0.10884267534260916,0.3538935462030571,-0.055223362127739956,-0.25819440152518536,0.14691121366485543];
observed{9}	=[0.8,-0.11351446579670872,-0.5871715953189904,0.32299872426993737,0.3357367607855463,-0.3886669488572505,-0.11092801803286569,0.3512503605862144,-0.05188139665849537,-0.2569942428910992,0.14310793968854146];
residuals{9}	=[0.0,0.0013084189982205674,4.95498668257599E-4,-0.0027730356743066453,5.222678415422455E-4,0.003280441020231839,-0.0020853426902565225,-0.0026431856168426804,0.003341965469244583,0.0012001586340861792,-0.00380327397631397];
obs_stdev{9}	=[0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1];
costObserved{9}	=0.0029545177133023974;
costTotal{9}	=0.0029545177133023974;
% SimulationKwadraticCostFunction: evaluation 9 : cost = 2.955E-3
% Error estimate for this outer iteration
parameterErrorEstimateStd{5}	=[1.2075743196937008,0.022116512334634512];
parameterErrorCorrelations{5}	=[0.9999999999999998,-0.39069400060531423;-0.39069400060531384,0.9999999999999998];
% Application starting next step
% ======================================================
% DUD outer iteration no.6
% ======================================================
% -----------------------------------------------------
costs{6}	=[0.0029545177133023974,0.019018162867878868,1.605684983098306];
% -----------------------------------------------------
% Start search until improvement,
% Next try p=   [8.482,1.8]
%   delta_p=    [-0.092,-1.068E-3]
% ========================================================================
% no10
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model2/OscillatorStochModel.xml
% simulationTimespan =[0.0,0.05,10.0]
% parameters@names =t_damp,omega parameters=[8.0,1.5708]
% parameterUncertainty@names =t_damp,omega parameterUncertainty=[1.0,0.1257]
% systemNoise = {[0.0,0.0],[0.3,0.3]}
% initialState = [0.8,0.0]
% initialStateUncertainty = [0.8,0.8]
evaluatedParameters{10}	=[8.481574404958938,1.7999642541935335];
predicted{10}	=[0.8,-0.11334852021359293,-0.5869362525910741,0.32261522652636093,0.3355282159604673,-0.3880860677678433,-0.11093322485271499,0.3506152251214009,-0.051602205958077525,-0.25647152562804354,0.14260936233807317];
observed{10}	=[0.8,-0.11351446579670872,-0.5871715953189904,0.32299872426993737,0.3357367607855463,-0.3886669488572505,-0.11092801803286569,0.3512503605862144,-0.05188139665849537,-0.2569942428910992,0.14310793968854146];
residuals{10}	=[0.0,-1.6594558311579244E-4,-2.3534272791636557E-4,3.8349774357643307E-4,2.0854482507903782E-4,-5.808810894071814E-4,5.206819849307154E-6,6.351354648134744E-4,-2.791907004178473E-4,-5.227172630556387E-4,4.985773504682933E-4];
obs_stdev{10}	=[0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1];
costObserved{10}	=8.070463839490167E-5;
costTotal{10}	=8.070463839490167E-5;
% SimulationKwadraticCostFunction: evaluation 10 : cost = 8.07E-5
% Error estimate for this outer iteration
parameterErrorEstimateStd{6}	=[1.6740775365781972,0.02187197093036392];
parameterErrorCorrelations{6}	=[1.0,-0.08149943770759074;-0.08149943770759067,0.9999999999999999];
% ===================================================================
% SimulationKwadraticCostfunction: optimal results
%     number of evaluations: 10
%     all cost values:
%         [47.245,50.778,12.032,11.693,4.312,1.606,1.894,0.019,2.955E-3,8.07E-5]
%     all parameter values
%         [8,9,8,3.865,5.176,8.306,9.31,8.808,8.574,8.482;1.571,1.571,1.696,1.768,1.842,1.839,1.76,1.8,1.801,1.8]
%     number of observations: 11
%     best cost:
%         cost = 8.07E-5
%     best parameters:
%                 [8.482,1.8]
% ===================================================================
% Application Done
