% opening :/home/verlaanm/deltares/src/openda_1_20100401/public/tests/waqua_usedoc_rough_stat9/./stochobserver/obs_rough_20000707.xml
% Starting Algorithm: 
%	className: org.openda.algorithms.Dud
%	dir.: /home/verlaanm/deltares/src/openda_1_20100401/public/tests/waqua_usedoc_rough_stat9/./algorithm
%	config.: dudAlgorithm.xml
%  Algorithm  Dud initialisation (stoch. obs. and stoch. model have been set)
%  Algorithm  configString = dudAlgorithm.xml
% opening :/home/verlaanm/deltares/src/openda_1_20100401/public/tests/waqua_usedoc_rough_stat9/./algorithm/dudAlgorithm.xml
%  Algorithm  Retrieving initial parameters from model
%  Algorithm  Starting optimizer
%  Algorithm  costFunction@class=org.openda.algorithms.SimulationKwadraticCostFunction
%  Algorithm  costFunction@weakParameterConstraint=false
%  Algorithm  costFunction@factor=0.5
%  Algorithm  costFunction@tryParallel=false
%  Algorithm  outerLoop@maxIterations=10
%  Algorithm  outerLoop@absTolerance=0.01
%  Algorithm  outerLoop@relTolerance=0.01
%  Algorithm  outerLoop@relToleranceLinearCost=0.01
%  Algorithm  lineSearch@maxIterations=5
%  Algorithm  lineSearch@maxRelStepSize=10.0
%  Algorithm  lineSearch/backtracking@startIterationNegativeLook=3
%  Algorithm  lineSearch/backtracking@shorteningFactor=0.5
% Application initializing finished
% Initializing application
% Evaluating with parameters 
% Params - containing:
%    A-411= 0
%    A-412= 0
%    A-413= 0
%    A-414= 0
%    A-415= 0
%    A-416= 0
%    A-417= 0
%    A-418= 0
%    A-419= 0
% ========================================================================
% prepare no 1
%  Params: A-411, A-412, A-413, A-414, A-415, A-416, A-417, A-418, A-419
evaluatedParameters{1}	=[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0];
% ========================================================================
% evaluate no. 1
costObserved{1}	=10.475393981479684;
costTotal{1}	=10.475393981479684;
% SimulationKwadraticCostFunction: evaluation 1 : cost = 10.475
% Evaluating with parameters 
% Params - containing:
%    A-411= 0.02
%    A-412= 0
%    A-413= 0
%    A-414= 0
%    A-415= 0
%    A-416= 0
%    A-417= 0
%    A-418= 0
%    A-419= 0
% ========================================================================
% prepare no 2
%  Params: A-411, A-412, A-413, A-414, A-415, A-416, A-417, A-418, A-419
evaluatedParameters{2}	=[0.02, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0];
% ========================================================================
% evaluate no. 2
costObserved{2}	=10.65011372489148;
costTotal{2}	=10.65011372489148;
% SimulationKwadraticCostFunction: evaluation 2 : cost = 10.65
% Evaluating with parameters 
% Params - containing:
%    A-411= 0
%    A-412= 0.02
%    A-413= 0
%    A-414= 0
%    A-415= 0
%    A-416= 0
%    A-417= 0
%    A-418= 0
%    A-419= 0
% ========================================================================
% prepare no 3
%  Params: A-411, A-412, A-413, A-414, A-415, A-416, A-417, A-418, A-419
evaluatedParameters{3}	=[0.0, 0.02, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0];
% ========================================================================
% evaluate no. 3
costObserved{3}	=12.46851057165752;
costTotal{3}	=12.46851057165752;
% SimulationKwadraticCostFunction: evaluation 3 : cost = 12.469
% Evaluating with parameters 
% Params - containing:
%    A-411= 0
%    A-412= 0
%    A-413= 0.02
%    A-414= 0
%    A-415= 0
%    A-416= 0
%    A-417= 0
%    A-418= 0
%    A-419= 0
% ========================================================================
% prepare no 4
%  Params: A-411, A-412, A-413, A-414, A-415, A-416, A-417, A-418, A-419
evaluatedParameters{4}	=[0.0, 0.0, 0.02, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0];
% ========================================================================
% evaluate no. 4
costObserved{4}	=9.812032826321508;
costTotal{4}	=9.812032826321508;
% SimulationKwadraticCostFunction: evaluation 4 : cost = 9.812
% Evaluating with parameters 
% Params - containing:
%    A-411= 0
%    A-412= 0
%    A-413= 0
%    A-414= 0.02
%    A-415= 0
%    A-416= 0
%    A-417= 0
%    A-418= 0
%    A-419= 0
% ========================================================================
% prepare no 5
%  Params: A-411, A-412, A-413, A-414, A-415, A-416, A-417, A-418, A-419
evaluatedParameters{5}	=[0.0, 0.0, 0.0, 0.02, 0.0, 0.0, 0.0, 0.0, 0.0];
% ========================================================================
% evaluate no. 5
costObserved{5}	=7.344127308811886;
costTotal{5}	=7.344127308811886;
% SimulationKwadraticCostFunction: evaluation 5 : cost = 7.344
% Evaluating with parameters 
% Params - containing:
%    A-411= 0
%    A-412= 0
%    A-413= 0
%    A-414= 0
%    A-415= 0.02
%    A-416= 0
%    A-417= 0
%    A-418= 0
%    A-419= 0
% ========================================================================
% prepare no 6
%  Params: A-411, A-412, A-413, A-414, A-415, A-416, A-417, A-418, A-419
evaluatedParameters{6}	=[0.0, 0.0, 0.0, 0.0, 0.02, 0.0, 0.0, 0.0, 0.0];
% ========================================================================
% evaluate no. 6
costObserved{6}	=13.071297847233835;
costTotal{6}	=13.071297847233835;
% SimulationKwadraticCostFunction: evaluation 6 : cost = 13.071
% Evaluating with parameters 
% Params - containing:
%    A-411= 0
%    A-412= 0
%    A-413= 0
%    A-414= 0
%    A-415= 0
%    A-416= 0.02
%    A-417= 0
%    A-418= 0
%    A-419= 0
% ========================================================================
% prepare no 7
%  Params: A-411, A-412, A-413, A-414, A-415, A-416, A-417, A-418, A-419
evaluatedParameters{7}	=[0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.0, 0.0, 0.0];
% ========================================================================
% evaluate no. 7
costObserved{7}	=12.972960255285379;
costTotal{7}	=12.972960255285379;
% SimulationKwadraticCostFunction: evaluation 7 : cost = 12.973
% Evaluating with parameters 
% Params - containing:
%    A-411= 0
%    A-412= 0
%    A-413= 0
%    A-414= 0
%    A-415= 0
%    A-416= 0
%    A-417= 0.02
%    A-418= 0
%    A-419= 0
% ========================================================================
% prepare no 8
%  Params: A-411, A-412, A-413, A-414, A-415, A-416, A-417, A-418, A-419
evaluatedParameters{8}	=[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.0, 0.0];
% ========================================================================
% evaluate no. 8
costObserved{8}	=11.725000827075089;
costTotal{8}	=11.725000827075089;
% SimulationKwadraticCostFunction: evaluation 8 : cost = 11.725
% Evaluating with parameters 
% Params - containing:
%    A-411= 0
%    A-412= 0
%    A-413= 0
%    A-414= 0
%    A-415= 0
%    A-416= 0
%    A-417= 0
%    A-418= 0.02
%    A-419= 0
% ========================================================================
% prepare no 9
%  Params: A-411, A-412, A-413, A-414, A-415, A-416, A-417, A-418, A-419
evaluatedParameters{9}	=[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.0];
% ========================================================================
% evaluate no. 9
costObserved{9}	=10.568471544908093;
costTotal{9}	=10.568471544908093;
% SimulationKwadraticCostFunction: evaluation 9 : cost = 10.568
% Evaluating with parameters 
% Params - containing:
%    A-411= 0
%    A-412= 0
%    A-413= 0
%    A-414= 0
%    A-415= 0
%    A-416= 0
%    A-417= 0
%    A-418= 0
%    A-419= 0.02
% ========================================================================
% prepare no 10
%  Params: A-411, A-412, A-413, A-414, A-415, A-416, A-417, A-418, A-419
evaluatedParameters{10}	=[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02];
% ========================================================================
% evaluate no. 10
costObserved{10}	=9.101302923362338;
costTotal{10}	=9.101302923362338;
% SimulationKwadraticCostFunction: evaluation 10 : cost = 9.101
% Application initialized
% Application starting next step
% ======================================================
% DUD outer iteration no.1
% ======================================================
% -----------------------------------------------------
costs{1}	=[7.344127308811886,9.101302923362338,9.812032826321508,10.475393981479684,10.568471544908093,10.65011372489148,11.725000827075089,12.46851057165752,12.972960255285379,13.071297847233835];
% -----------------------------------------------------
% Start search until improvement,
% Next try p=
%   Params - containing:
%      A-411= 0.03
%      A-412= -0.034
%      A-413= 0.02
%      A-414= 0.038
%      A-415= -0.033
%      A-416= -0.022
%      A-417= -2.083E-5
%      A-418= 0.01
%      A-419= 0.02
%   delta_p= 
%   Params - containing:
%      A-411= 0.03
%      A-412= -0.034
%      A-413= 0.02
%      A-414= 0.018
%      A-415= -0.033
%      A-416= -0.022
%      A-417= -2.083E-5
%      A-418= 0.01
%      A-419= 0.02
% ========================================================================
% prepare no 11
%  Params: A-411, A-412, A-413, A-414, A-415, A-416, A-417, A-418, A-419
evaluatedParameters{11}	=[0.029587231149000637, -0.033620865064649966, 0.020395335706991837, 0.03772841775057813, -0.0334412381419587, -0.021967527190926955, -2.0829786864763332E-5, 0.010296148661485793, 0.01995766371367313];
% ========================================================================
% evaluate no. 11
costObserved{11}	=0.07652926459913593;
costTotal{11}	=0.07652926459913593;
% SimulationKwadraticCostFunction: evaluation 11 : cost = 0.077
linearCost{1}	=1.5172471443244175E-4;
% linearized cost rel error0.010400026390350535 < 0.01
% Error estimate for this outer iteration
%  Params: A-411, A-412, A-413, A-414, A-415, A-416, A-417, A-418, A-419
parameterErrorEstimateStd{1}	=[0.03954303172308691, 0.020227352170501593, 0.01555610458475943, 0.012824326207301463, 0.01870421422210616, 0.017904596563014928, 0.015406020182249836, 0.013618180902773155, 0.011768314048772351];
parameterErrorCorrelations{1}	=[1.0000000000000002,-0.39359911785957047,0.06934356864361108,0.09128048938168472,0.1353331354717972,0.16335348534091787,0.1780303150370984,0.18772681918974105,0.175781246531927;-0.3935991178595701,1.0,-0.19829415078289536,-0.020307722453554673,-0.023564540094464653,-0.024350114827344962,-0.024894565068266223,-0.025645260373723555,-0.023645411702641114;0.06934356864361027,-0.19829415078289475,0.9999999999999999,-0.0599562332725577,0.021656830193513388,0.04012752141061271,0.04985022793238574,0.054984926840662675,0.05232403236211902;0.09128048938168096,-0.020307722453555294,-0.05995623327255952,1.0,-0.06254375772978035,0.03770268205334547,0.051487932247975675,0.05871587710057868,0.056408667259897856;0.13533313547179635,-0.02356454009446432,0.021656830193513655,-0.06254375772977885,1.0,-0.03355448210594744,0.05692881823745838,0.06898153444164698,0.06758970823851675;0.1633534853409183,-0.02435011482734526,0.04012752141061339,0.037702682053346454,-0.03355448210594707,1.0,-0.00989997361843623,0.07139866637700974,0.07274301228165453;0.1780303150370994,-0.024894565068265963,0.0498502279323866,0.05148793224797687,0.05692881823745975,-0.009899973618435547,1.0000000000000002,0.00848809961135801,0.07538300359519264;0.18772681918974068,-0.025645260373723573,0.05498492684066343,0.05871587710058,0.06898153444164687,0.0713986663770093,0.0084880996113569,1.0000000000000002,0.018482839929362747;0.1757812465319267,-0.023645411702641138,0.052324032362120194,0.056408667259899084,0.06758970823851741,0.07274301228165468,0.07538300359519322,0.018482839929362675,0.9999999999999998];
% Application starting next step
% ======================================================
% DUD outer iteration no.2
% ======================================================
% -----------------------------------------------------
costs{2}	=[0.07652926459913593,7.344127308811886,9.101302923362338,9.812032826321508,10.475393981479684,10.568471544908093,10.65011372489148,11.725000827075089,12.46851057165752,12.972960255285379];
% -----------------------------------------------------
% Start search until improvement,
% Next try p=
%   Params - containing:
%      A-411= 0.03
%      A-412= -0.03
%      A-413= 0.02
%      A-414= 0.039
%      A-415= -0.03
%      A-416= -0.02
%      A-417= 3.988E-8
%      A-418= 0.01
%      A-419= 0.02
%   delta_p= 
%   Params - containing:
%      A-411= 2.186E-4
%      A-412= 3.982E-3
%      A-413= -3.783E-4
%      A-414= 1.632E-3
%      A-415= 3.719E-3
%      A-416= 2.041E-3
%      A-417= 2.087E-5
%      A-418= -2.67E-4
%      A-419= 2.492E-5
% ========================================================================
% prepare no 12
%  Params: A-411, A-412, A-413, A-414, A-415, A-416, A-417, A-418, A-419
evaluatedParameters{12}	=[0.02980578784489008, -0.029638825741383296, 0.020017055631134872, 0.039360742854908884, -0.029722375465941937, -0.019926152833719717, 3.9882364867214576E-8, 0.010029196818569529, 0.019982586126616863];
% ========================================================================
% evaluate no. 12
costObserved{12}	=0.001248470372525991;
costTotal{12}	=0.001248470372525991;
% SimulationKwadraticCostFunction: evaluation 12 : cost = 1.248E-3
linearCost{2}	=6.322191167276974E-6;
% linearized cost rel error0.016232363030898873 < 0.01
% Error estimate for this outer iteration
%  Params: A-411, A-412, A-413, A-414, A-415, A-416, A-417, A-418, A-419
parameterErrorEstimateStd{2}	=[0.039526642399144954, 0.020401686404874598, 0.01556212348324743, 0.012913605017515678, 0.016623439287141894, 0.017979150565227076, 0.015405357864103903, 0.013629290595727884, 0.01176737774010323];
parameterErrorCorrelations{2}	=[1.0,-0.404766332490861,0.07107678057222296,0.08152738969653439,0.13229057881135312,0.15444652193436817,0.1778348040161655,0.18888677283450342,0.17550156040424164;-0.4047663324908606,0.9999999999999998,-0.2006884308171553,-0.0038460926094903507,-0.13252546318605893,-0.011978738898909953,-0.03079689295782997,-0.03438078709847298,-0.030666670457466537;0.07107678057222706,-0.20068843081715418,1.0,-0.0628544603248127,0.035242555940183336,0.037253480780821996,0.05058022580817976,0.05624179351525855,0.05318499703634217;0.08152738969653447,-0.003846092609490042,-0.06285446032480532,1.0,-0.13280480773346626,0.04808224404689535,0.04721046514003953,0.051935435261775804,0.05140233101804261;0.13229057881135609,-0.13252546318605218,0.0352425559401884,-0.13280480773344083,0.9999999999999998,-0.09691632623994953,0.056171004483770756,0.07987678323756014,0.0664074771272514;0.15444652193436523,-0.01197873889890718,0.03725348078081503,0.048082244046888614,-0.09691632623996133,1.0000000000000002,-0.01340085501741228,0.06560643823564306,0.06826987763751535;0.17783480401616952,-0.030796892957830935,0.050580225808178926,0.047210465140040625,0.05617100448376887,-0.01340085501741093,1.0,0.00904429958210775,0.07527414373760277;0.18888677283450916,-0.034380787098470825,0.05624179351525855,0.05193543526176848,0.07987678323755518,0.0656064382356541,0.009044299582109824,1.0,0.01911464709467753;0.1755015604042488,-0.030666670457467304,0.0531849970363434,0.05140233101803251,0.06640747712724603,0.06826987763752264,0.07527414373760465,0.019114647094678215,1.0000000000000002];
% Application starting next step
% ======================================================
% DUD outer iteration no.3
% ======================================================
% -----------------------------------------------------
costs{3}	=[0.001248470372525991,0.07652926459913593,7.344127308811886,9.101302923362338,9.812032826321508,10.475393981479684,10.568471544908093,10.65011372489148,11.725000827075089,12.46851057165752];
% -----------------------------------------------------
% Start search until improvement,
% Next try p=
%   Params - containing:
%      A-411= 0.03
%      A-412= -0.03
%      A-413= 0.02
%      A-414= 0.04
%      A-415= -0.03
%      A-416= -0.02
%      A-417= 1.546E-6
%      A-418= 0.01
%      A-419= 0.02
%   delta_p= 
%   Params - containing:
%      A-411= 1.395E-4
%      A-412= -3.012E-4
%      A-413= -8.766E-6
%      A-414= 4.17E-4
%      A-415= -2.239E-4
%      A-416= -5.028E-5
%      A-417= 1.507E-6
%      A-418= -2.351E-5
%      A-419= 1.183E-5
% ========================================================================
% prepare no 13
%  Params: A-411, A-412, A-413, A-414, A-415, A-416, A-417, A-418, A-419
evaluatedParameters{13}	=[0.029945249166477964, -0.02994007059554397, 0.020008289712231, 0.03977775837823853, -0.029946255213741606, -0.019976432845067987, 1.5464407632170055E-6, 0.010005685441680593, 0.0199944185503383];
% ========================================================================
% evaluate no. 13
costObserved{13}	=1.1730707376823815E-4;
costTotal{13}	=1.1730707376823815E-4;
% SimulationKwadraticCostFunction: evaluation 13 : cost = 1.173E-4
linearCost{3}	=1.7251348987142986E-6;
% linearized cost rel error0.09270694234973921 < 0.01
% Error estimate for this outer iteration
%  Params: A-411, A-412, A-413, A-414, A-415, A-416, A-417, A-418, A-419
parameterErrorEstimateStd{3}	=[0.03978718888064361, 0.026760918327161144, 0.015575299217025954, 0.024713415901792563, 0.013813750881964643, 0.02018694798927119, 0.015408379086110799, 0.013705383765312516, 0.011770146411180897];
parameterErrorCorrelations{3}	=[1.0,-0.3929570189668238,0.06419457617363253,0.14940102169824346,0.1759410234815544,-0.025987131284610783,0.18463612292327758,0.17204921751694832,0.17700478120954669;-0.39295701896721363,1.0,-0.12545591898120154,-0.5580258895732944,0.029557928091596325,0.54063721728124,-0.05015154703179197,0.042181500465901293,-0.04481430470750551;0.06419457617358802,-0.12545591898122235,1.0000000000000002,-0.07037907672381734,0.05947255747646065,0.06179953052166225,0.04914865912058722,0.06027893346213852,0.05146958233428059;0.1494010216986684,-0.558025889572993,-0.0703790767240094,1.0,-0.20503939366562413,-0.7868906878584354,0.06415762504758772,-0.06447577920058738,0.05254189902032933;0.17594102348139243,0.029557928091932435,0.059472557476459645,-0.2050393936664,1.0,0.6016659787919341,0.03231942760946087,0.12514345980151056,0.08991564140273;-0.025987131285219983,0.5406372172814402,0.0617995305216282,-0.7868906878589782,0.6016659787915953,1.0,-0.02101091782377194,0.14256430424520908,0.01476794310670653;0.18463612292323822,-0.05015154703146842,0.04914865912060701,0.06415762504716706,0.03231942760957107,-0.021010917823316207,1.0000000000000002,0.00489601128165097,0.07751299982727237;0.17204921751686766,0.04218150046603257,0.060278933462232516,-0.06447577920064387,0.12514345980143043,0.14256430424540778,0.004896011281587756,1.0,0.015336181332141317;0.1770047812095589,-0.04481430470736894,0.05146958233434516,0.0525418990204315,0.08991564140276714,0.01476794310694569,0.07751299982727941,0.015336181332108467,1.0];
% ===================================================================
% SimulationKwadraticCostfunction: optimal results
%     number of evaluations: 13
%     all cost values:
%         [10.475,10.65,12.469,9.812,7.344,13.071,12.973,11.725,10.568,9.101,0.077,1.248E-3,1.173E-4]
%     all parameter values
%         [0,0.02,0,0,0,0,0,0,0,0,0.03,0.03,0.03;0,0,0.02,0,0,0,0,0,0,0,-0.034,-0.03,-0.03;0,0,0,0.02,0,0,0,0,0,0,0.02,0.02,0.02;0,0,0,0,0.02,0,0,0,0,0,0.038,0.039,0.04;0,0,0,0,0,0.02,0,0,0,0,-0.033,-0.03,-0.03;0,0,0,0,0,0,0.02,0,0,0,-0.022,-0.02,-0.02;0,0,0,0,0,0,0,0.02,0,0,-2.083E-5,3.988E-8,1.546E-6;0,0,0,0,0,0,0,0,0.02,0,0.01,0.01,0.01;0,0,0,0,0,0,0,0,0,0.02,0.02,0.02,0.02]
%     number of observations: 808
%     best cost:
%         cost = 1.173E-4
%     best parameters:
%         
%        Params - containing:
%           A-411= 0.03
%           A-412= -0.03
%           A-413= 0.02
%           A-414= 0.04
%           A-415= -0.03
%           A-416= -0.02
%           A-417= 1.546E-6
%           A-418= 0.01
%           A-419= 0.02
% ===================================================================
% Application Done
