% 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/./model1/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/./model1/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/./model1/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.04313443918374465,-0.6309221702813459,0.1759519292051853,0.4657659569580843,-0.2452863996167396,-0.3179787785265744,0.2676806537476915,0.1946206473664078,-0.25796354297701396,-0.09802109778430443];
residuals{1}	=[0.0,-0.10300849000705582,-0.008407429022448842,0.22805168831317807,-0.018591390128484264,-0.2901175989177405,0.05884531247445479,0.30590623973760833,-0.09851447303582975,-0.29030932920303837,0.12998862887590798];
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}	=17.756295868749987;
costTotal{1}	=17.756295868749987;
% SimulationKwadraticCostFunction: evaluation 1 : cost = 17.756
% Evaluating with parameters  [9,1.571]
% ========================================================================
% no2
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model1/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.04313443918374465,-0.6309221702813459,0.1759519292051853,0.4657659569580843,-0.2452863996167396,-0.3179787785265744,0.2676806537476915,0.1946206473664078,-0.25796354297701396,-0.09802109778430443];
residuals{2}	=[0.0,-0.09671306712973375,0.009290951946216541,0.22335893272753982,-0.046543223229187136,-0.28685161101529827,0.0919540053283418,0.30384779557849934,-0.13337348457273365,-0.2892317218471685,0.16439614154145513];
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}	=18.651382948941716;
costTotal{2}	=18.651382948941716;
% SimulationKwadraticCostFunction: evaluation 2 : cost = 18.651
% Evaluating with parameters  [8,1.696]
% ========================================================================
% no3
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model1/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.04313443918374465,-0.6309221702813459,0.1759519292051853,0.4657659569580843,-0.2452863996167396,-0.3179787785265744,0.2676806537476915,0.1946206473664078,-0.25796354297701396,-0.09802109778430443];
residuals{3}	=[0.0,-0.009631319904631869,-0.015024276675364656,0.018547709503065474,0.019663211925355806,-0.027588885114972345,-0.016967462675909983,0.03396074056172119,0.010022740041215006,-0.03658518314069384,-0.0014548136659924688];
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}	=0.23462856704877713;
costTotal{3}	=0.23462856704877713;
% SimulationKwadraticCostFunction: evaluation 3 : cost = 0.235
%  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}	=[0.23462856704877713,17.756295868749987,18.651382948941716];
% -----------------------------------------------------
% Start search until improvement,
% Next try p=   [8.663,1.71]
%   delta_p=    [0.663,0.014]
% ========================================================================
% no4
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model1/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}	=[8.663227729470337,1.710474267579495];
predicted{4}	=[0.8,-0.04871533545298321,-0.623329637530464,0.18902524808841653,0.4492358459183649,-0.25841736087328293,-0.29429346004168633,0.2761308470105584,0.16701975973254954,-0.25949317150559,-0.06999665346208908];
observed{4}	=[0.8,-0.04313443918374465,-0.6309221702813459,0.1759519292051853,0.4657659569580843,-0.2452863996167396,-0.3179787785265744,0.2676806537476915,0.1946206473664078,-0.25796354297701396,-0.09802109778430443];
residuals{4}	=[0.0,0.005580896269238564,-0.007592532750881897,-0.013073318883231239,0.01653011103971941,0.013130961256543339,-0.023685318484888085,-0.00845019326286689,0.027600887633858262,0.0015296285285760214,-0.028024444322215347];
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}	=0.14436448284808975;
costTotal{4}	=0.14436448284808975;
% SimulationKwadraticCostFunction: evaluation 4 : cost = 0.144
% 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}	=[0.14436448284808975,0.23462856704877713,17.756295868749987];
% -----------------------------------------------------
% Start search until improvement,
% Next try p=   [9.057,1.699]
%   delta_p=    [0.393,-0.011]
% ========================================================================
% no5
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model1/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}	=[9.056627185141876,1.6993903683204152];
predicted{5}	=[0.8,-0.04303478974611259,-0.6318717397479764,0.17559631802055461,0.4674586318105381,-0.24517959039061235,-0.32008675659863206,0.26806876891663833,0.19680560214609158,-0.2588953862726204,-0.1000007750551769];
observed{5}	=[0.8,-0.04313443918374465,-0.6309221702813459,0.1759519292051853,0.4657659569580843,-0.2452863996167396,-0.3179787785265744,0.2676806537476915,0.1946206473664078,-0.25796354297701396,-0.09802109778430443];
residuals{5}	=[0.0,-9.964943763205586E-5,9.495694666304955E-4,3.5561118463067753E-4,-0.0016926748524537794,-1.0680922612724131E-4,0.002107978072057648,-3.881151689468232E-4,-0.0021849547796837787,9.31843295606416E-4,0.0019796772708724764];
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}	=9.035157107104139E-4;
costTotal{5}	=9.035157107104139E-4;
% SimulationKwadraticCostFunction: evaluation 5 : cost = 9.035E-4
% Error estimate for this outer iteration
parameterErrorEstimateStd{2}	=[4.465558725976945,0.04314757949382094];
parameterErrorCorrelations{2}	=[1.0,-0.8446564104810934;-0.8446564104810911,1.0000000000000002];
% Application starting next step
% ======================================================
% DUD outer iteration no.3
% ======================================================
% -----------------------------------------------------
costs{3}	=[9.035157107104139E-4,0.14436448284808975,0.23462856704877713];
% -----------------------------------------------------
% Start search until improvement,
% Next try p=   [9.004,1.7]
%   delta_p=    [-0.052,6.186E-4]
% ========================================================================
% no6
% configstring = OscillatorStochModel.xml
% opening :/home/verlaanm/tudelft/svn_costa/trunk/openda/tests/two_oscillators/./model1/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}	=[9.004484712448102,1.7000089455533829];
predicted{6}	=[0.8,-0.04316935722993659,-0.6309837467730878,0.17601944806912906,0.46584993580789136,-0.24539107920893447,-0.3180553199030431,0.2678144337862502,0.19467083253085882,-0.25811241634984355,-0.0980361973210575];
observed{6}	=[0.8,-0.04313443918374465,-0.6309221702813459,0.1759519292051853,0.4657659569580843,-0.2452863996167396,-0.3179787785265744,0.2676806537476915,0.1946206473664078,-0.25796354297701396,-0.09802109778430443];
residuals{6}	=[0.0,3.4918046191941E-5,6.157649174187707E-5,-6.751886394376627E-5,-8.397884980704262E-5,1.0467959219487288E-4,7.654137646867509E-5,-1.3378003855868403E-4,-5.0185164451016506E-5,1.4887337282959434E-4,1.509953675307707E-5];
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}	=3.8122752208815227E-6;
costTotal{6}	=3.8122752208815227E-6;
% SimulationKwadraticCostFunction: evaluation 6 : cost = 3.812E-6
% Error estimate for this outer iteration
parameterErrorEstimateStd{3}	=[1.5743609516034422,0.02036735626685016];
parameterErrorCorrelations{3}	=[1.0000000000000002,-0.11806660474048626;-0.11806660474048615,0.9999999999999999];
% ===================================================================
% SimulationKwadraticCostfunction: optimal results
%     number of evaluations: 6
%     all cost values:
%         [17.756,18.651,0.235,0.144,9.035E-4,3.812E-6]
%     all parameter values
%         [8,9,8,8.663,9.057,9.004;1.571,1.571,1.696,1.71,1.699,1.7]
%     number of observations: 11
%     best cost:
%         cost = 3.812E-6
%     best parameters:
%                 [9.004,1.7]
% ===================================================================
% Application Done
