Index: /issm/trunk/test/NightlyRun/qmu/SquareSheetShelf.qmu.in =================================================================== --- /issm/trunk/test/NightlyRun/qmu/SquareSheetShelf.qmu.in (revision 5015) +++ /issm/trunk/test/NightlyRun/qmu/SquareSheetShelf.qmu.in (revision 5015) @@ -0,0 +1,65 @@ +strategy, + single_method + tabular_graphics_data + +method, + nond_local_reliability + distribution cumulative + num_probability_levels = 9 9 9 9 9 9 9 + probability_levels = + 0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 + 0.9999 + 0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 + 0.9999 + 0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 + 0.9999 + 0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 + 0.9999 + 0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 + 0.9999 + 0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 + 0.9999 + 0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 + 0.9999 + +model, + single + +variables, + normal_uncertain = 20 + nuv_means = + 1 1 1 1 1 1 + 1 1 1 1 1 1 + 1 1 1 1 1 1 + 1 1 + nuv_std_deviations = + 0.01 0.01 0.01 0.01 0.01 0.01 + 0.01 0.01 0.01 0.01 0.01 0.01 + 0.01 0.01 0.01 0.01 0.01 0.01 + 0.01 0.01 + nuv_descriptors = + 'DragCoefficient1' 'DragCoefficient2' 'DragCoefficient3' + 'DragCoefficient4' 'DragCoefficient5' 'DragCoefficient6' + 'DragCoefficient7' 'DragCoefficient8' 'DragCoefficient9' + 'DragCoefficient10' 'DragCoefficient11' 'DragCoefficient12' + 'DragCoefficient13' 'DragCoefficient14' 'DragCoefficient15' + 'DragCoefficient16' 'DragCoefficient17' 'DragCoefficient18' + 'DragCoefficient19' 'DragCoefficient20' + +interface, + direct + analysis_driver = 'diagnostic' + failure_capture abort + deactivate evaluation_cache restart_file + +responses, + num_response_functions = 7 + response_descriptors = + 'MaxVel' 'MassFlux1' 'MassFlux2' 'MassFlux3' 'MassFlux4' 'MassFlux5' + 'MassFlux6' + numerical_gradients + method_source dakota + interval_type forward + fd_gradient_step_size = 0.001 + no_hessians + Index: /issm/trunk/test/NightlyRun/qmu/SquareSheetShelf.qmu.out =================================================================== --- /issm/trunk/test/NightlyRun/qmu/SquareSheetShelf.qmu.out (revision 5015) +++ /issm/trunk/test/NightlyRun/qmu/SquareSheetShelf.qmu.out (revision 5015) @@ -0,0 +1,1047 @@ +Writing new restart file dakota.rst +methodName = nond_local_reliability +gradientType = numerical +Numerical gradients using forward differences +to be calculated by the dakota finite difference routine. +hessianType = none + +>>>>> Running nond_local_reliability iterator. + +>>>>> Evaluating response at mean values + +------------------------------------------ +Begin Dakota derivative estimation routine +------------------------------------------ + +>>>>> Initial map for analytic portion of response: + +------------------------------ +Begin Function Evaluation 1 +------------------------------ +Parameters for function evaluation 1: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 1: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[1] + h: + +------------------------------ +Begin Function Evaluation 2 +------------------------------ +Parameters for function evaluation 2: + 1.0010000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 2: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[2] + h: + +------------------------------ +Begin Function Evaluation 3 +------------------------------ +Parameters for function evaluation 3: + 1.0000000000e+00 DragCoefficient1 + 1.0010000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 3: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[3] + h: + +------------------------------ +Begin Function Evaluation 4 +------------------------------ +Parameters for function evaluation 4: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0010000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 4: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[4] + h: + +------------------------------ +Begin Function Evaluation 5 +------------------------------ +Parameters for function evaluation 5: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0010000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 5: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[5] + h: + +------------------------------ +Begin Function Evaluation 6 +------------------------------ +Parameters for function evaluation 6: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0010000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 6: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[6] + h: + +------------------------------ +Begin Function Evaluation 7 +------------------------------ +Parameters for function evaluation 7: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0010000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 7: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[7] + h: + +------------------------------ +Begin Function Evaluation 8 +------------------------------ +Parameters for function evaluation 8: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0010000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 8: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[8] + h: + +------------------------------ +Begin Function Evaluation 9 +------------------------------ +Parameters for function evaluation 9: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0010000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 9: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[9] + h: + +------------------------------ +Begin Function Evaluation 10 +------------------------------ +Parameters for function evaluation 10: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0010000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 10: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[10] + h: + +------------------------------ +Begin Function Evaluation 11 +------------------------------ +Parameters for function evaluation 11: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0010000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 11: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[11] + h: + +------------------------------ +Begin Function Evaluation 12 +------------------------------ +Parameters for function evaluation 12: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0010000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 12: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[12] + h: + +------------------------------ +Begin Function Evaluation 13 +------------------------------ +Parameters for function evaluation 13: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0010000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 13: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[13] + h: + +------------------------------ +Begin Function Evaluation 14 +------------------------------ +Parameters for function evaluation 14: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0010000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 14: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[14] + h: + +------------------------------ +Begin Function Evaluation 15 +------------------------------ +Parameters for function evaluation 15: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0010000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 15: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[15] + h: + +------------------------------ +Begin Function Evaluation 16 +------------------------------ +Parameters for function evaluation 16: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0010000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 16: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[16] + h: + +------------------------------ +Begin Function Evaluation 17 +------------------------------ +Parameters for function evaluation 17: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0010000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 17: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[17] + h: + +------------------------------ +Begin Function Evaluation 18 +------------------------------ +Parameters for function evaluation 18: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0010000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 18: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[18] + h: + +------------------------------ +Begin Function Evaluation 19 +------------------------------ +Parameters for function evaluation 19: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0010000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 19: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[19] + h: + +------------------------------ +Begin Function Evaluation 20 +------------------------------ +Parameters for function evaluation 20: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0010000000e+00 DragCoefficient19 + 1.0000000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 20: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Dakota finite difference gradient evaluation for x[20] + h: + +------------------------------ +Begin Function Evaluation 21 +------------------------------ +Parameters for function evaluation 21: + 1.0000000000e+00 DragCoefficient1 + 1.0000000000e+00 DragCoefficient2 + 1.0000000000e+00 DragCoefficient3 + 1.0000000000e+00 DragCoefficient4 + 1.0000000000e+00 DragCoefficient5 + 1.0000000000e+00 DragCoefficient6 + 1.0000000000e+00 DragCoefficient7 + 1.0000000000e+00 DragCoefficient8 + 1.0000000000e+00 DragCoefficient9 + 1.0000000000e+00 DragCoefficient10 + 1.0000000000e+00 DragCoefficient11 + 1.0000000000e+00 DragCoefficient12 + 1.0000000000e+00 DragCoefficient13 + 1.0000000000e+00 DragCoefficient14 + 1.0000000000e+00 DragCoefficient15 + 1.0000000000e+00 DragCoefficient16 + 1.0000000000e+00 DragCoefficient17 + 1.0000000000e+00 DragCoefficient18 + 1.0000000000e+00 DragCoefficient19 + 1.0010000000e+00 DragCoefficient20 + +Direct function: invoking diagnostic + +Active response data for function evaluation 21: +Active set vector = { 1 1 1 1 1 1 1 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + + +>>>>> Total response returned to iterator: + +Active set vector = { 3 3 3 3 3 3 3 } Deriv vars vector = { 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 } + 3.1536000000e+07 MaxVel + -8.6755499563e+00 MassFlux1 + -8.6755499563e+00 MassFlux2 + -1.7351102342e+01 MassFlux3 + -5.3543305858e-06 MassFlux4 + -2.5421317493e-06 MassFlux5 + -6.3466722620e-07 MassFlux6 + [ 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + ] MaxVel gradient + [ 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + ] MassFlux1 gradient + [ 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + ] MassFlux2 gradient + [ 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + ] MassFlux3 gradient + [ 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + ] MassFlux4 gradient + [ 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + ] MassFlux5 gradient + [ 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00 + ] MassFlux6 gradient + + + +<<<<< Iterator nond_local_reliability completed. +<<<<< Function evaluation summary: 21 total (21 new, 0 duplicate) +----------------------------------------------------------------- +MV Statistics for MaxVel: + Approximate Mean Response = 3.1536000000e+07 + Approximate Standard Deviation of Response = 0.0000000000e+00 + Importance Factors not available. + +Warning: negligible standard deviation renders CDF results suspect. + +Cumulative Distribution Function (CDF) for MaxVel: + Response Level Probability Level Reliability Index General Rel Index + -------------- ----------------- ----------------- ----------------- + 3.1536000000e+07 1.0000000000e-04 3.7190164855e+00 3.7190164855e+00 + 3.1536000000e+07 1.0000000000e-03 3.0902323062e+00 3.0902323062e+00 + 3.1536000000e+07 1.0000000000e-02 2.3263478740e+00 2.3263478740e+00 + 3.1536000000e+07 2.5000000000e-01 6.7448975020e-01 6.7448975020e-01 + 3.1536000000e+07 5.0000000000e-01 -0.0000000000e+00 -0.0000000000e+00 + 3.1536000000e+07 7.5000000000e-01 -6.7448975020e-01 -6.7448975020e-01 + 3.1536000000e+07 9.9000000000e-01 -2.3263478740e+00 -2.3263478740e+00 + 3.1536000000e+07 9.9900000000e-01 -3.0902323062e+00 -3.0902323062e+00 + 3.1536000000e+07 9.9990000000e-01 -3.7190164855e+00 -3.7190164855e+00 +MV Statistics for MassFlux1: + Approximate Mean Response = -8.6755499563e+00 + Approximate Standard Deviation of Response = 0.0000000000e+00 + Importance Factors not available. + +Warning: negligible standard deviation renders CDF results suspect. + +Cumulative Distribution Function (CDF) for MassFlux1: + Response Level Probability Level Reliability Index General Rel Index + -------------- ----------------- ----------------- ----------------- + -8.6755499563e+00 1.0000000000e-04 3.7190164855e+00 3.7190164855e+00 + -8.6755499563e+00 1.0000000000e-03 3.0902323062e+00 3.0902323062e+00 + -8.6755499563e+00 1.0000000000e-02 2.3263478740e+00 2.3263478740e+00 + -8.6755499563e+00 2.5000000000e-01 6.7448975020e-01 6.7448975020e-01 + -8.6755499563e+00 5.0000000000e-01 -0.0000000000e+00 -0.0000000000e+00 + -8.6755499563e+00 7.5000000000e-01 -6.7448975020e-01 -6.7448975020e-01 + -8.6755499563e+00 9.9000000000e-01 -2.3263478740e+00 -2.3263478740e+00 + -8.6755499563e+00 9.9900000000e-01 -3.0902323062e+00 -3.0902323062e+00 + -8.6755499563e+00 9.9990000000e-01 -3.7190164855e+00 -3.7190164855e+00 +MV Statistics for MassFlux2: + Approximate Mean Response = -8.6755499563e+00 + Approximate Standard Deviation of Response = 0.0000000000e+00 + Importance Factors not available. + +Warning: negligible standard deviation renders CDF results suspect. + +Cumulative Distribution Function (CDF) for MassFlux2: + Response Level Probability Level Reliability Index General Rel Index + -------------- ----------------- ----------------- ----------------- + -8.6755499563e+00 1.0000000000e-04 3.7190164855e+00 3.7190164855e+00 + -8.6755499563e+00 1.0000000000e-03 3.0902323062e+00 3.0902323062e+00 + -8.6755499563e+00 1.0000000000e-02 2.3263478740e+00 2.3263478740e+00 + -8.6755499563e+00 2.5000000000e-01 6.7448975020e-01 6.7448975020e-01 + -8.6755499563e+00 5.0000000000e-01 -0.0000000000e+00 -0.0000000000e+00 + -8.6755499563e+00 7.5000000000e-01 -6.7448975020e-01 -6.7448975020e-01 + -8.6755499563e+00 9.9000000000e-01 -2.3263478740e+00 -2.3263478740e+00 + -8.6755499563e+00 9.9900000000e-01 -3.0902323062e+00 -3.0902323062e+00 + -8.6755499563e+00 9.9990000000e-01 -3.7190164855e+00 -3.7190164855e+00 +MV Statistics for MassFlux3: + Approximate Mean Response = -1.7351102342e+01 + Approximate Standard Deviation of Response = 0.0000000000e+00 + Importance Factors not available. + +Warning: negligible standard deviation renders CDF results suspect. + +Cumulative Distribution Function (CDF) for MassFlux3: + Response Level Probability Level Reliability Index General Rel Index + -------------- ----------------- ----------------- ----------------- + -1.7351102342e+01 1.0000000000e-04 3.7190164855e+00 3.7190164855e+00 + -1.7351102342e+01 1.0000000000e-03 3.0902323062e+00 3.0902323062e+00 + -1.7351102342e+01 1.0000000000e-02 2.3263478740e+00 2.3263478740e+00 + -1.7351102342e+01 2.5000000000e-01 6.7448975020e-01 6.7448975020e-01 + -1.7351102342e+01 5.0000000000e-01 -0.0000000000e+00 -0.0000000000e+00 + -1.7351102342e+01 7.5000000000e-01 -6.7448975020e-01 -6.7448975020e-01 + -1.7351102342e+01 9.9000000000e-01 -2.3263478740e+00 -2.3263478740e+00 + -1.7351102342e+01 9.9900000000e-01 -3.0902323062e+00 -3.0902323062e+00 + -1.7351102342e+01 9.9990000000e-01 -3.7190164855e+00 -3.7190164855e+00 +MV Statistics for MassFlux4: + Approximate Mean Response = -5.3543305858e-06 + Approximate Standard Deviation of Response = 0.0000000000e+00 + Importance Factors not available. + +Warning: negligible standard deviation renders CDF results suspect. + +Cumulative Distribution Function (CDF) for MassFlux4: + Response Level Probability Level Reliability Index General Rel Index + -------------- ----------------- ----------------- ----------------- + -5.3543305858e-06 1.0000000000e-04 3.7190164855e+00 3.7190164855e+00 + -5.3543305858e-06 1.0000000000e-03 3.0902323062e+00 3.0902323062e+00 + -5.3543305858e-06 1.0000000000e-02 2.3263478740e+00 2.3263478740e+00 + -5.3543305858e-06 2.5000000000e-01 6.7448975020e-01 6.7448975020e-01 + -5.3543305858e-06 5.0000000000e-01 -0.0000000000e+00 -0.0000000000e+00 + -5.3543305858e-06 7.5000000000e-01 -6.7448975020e-01 -6.7448975020e-01 + -5.3543305858e-06 9.9000000000e-01 -2.3263478740e+00 -2.3263478740e+00 + -5.3543305858e-06 9.9900000000e-01 -3.0902323062e+00 -3.0902323062e+00 + -5.3543305858e-06 9.9990000000e-01 -3.7190164855e+00 -3.7190164855e+00 +MV Statistics for MassFlux5: + Approximate Mean Response = -2.5421317493e-06 + Approximate Standard Deviation of Response = 0.0000000000e+00 + Importance Factors not available. + +Warning: negligible standard deviation renders CDF results suspect. + +Cumulative Distribution Function (CDF) for MassFlux5: + Response Level Probability Level Reliability Index General Rel Index + -------------- ----------------- ----------------- ----------------- + -2.5421317493e-06 1.0000000000e-04 3.7190164855e+00 3.7190164855e+00 + -2.5421317493e-06 1.0000000000e-03 3.0902323062e+00 3.0902323062e+00 + -2.5421317493e-06 1.0000000000e-02 2.3263478740e+00 2.3263478740e+00 + -2.5421317493e-06 2.5000000000e-01 6.7448975020e-01 6.7448975020e-01 + -2.5421317493e-06 5.0000000000e-01 -0.0000000000e+00 -0.0000000000e+00 + -2.5421317493e-06 7.5000000000e-01 -6.7448975020e-01 -6.7448975020e-01 + -2.5421317493e-06 9.9000000000e-01 -2.3263478740e+00 -2.3263478740e+00 + -2.5421317493e-06 9.9900000000e-01 -3.0902323062e+00 -3.0902323062e+00 + -2.5421317493e-06 9.9990000000e-01 -3.7190164855e+00 -3.7190164855e+00 +MV Statistics for MassFlux6: + Approximate Mean Response = -6.3466722620e-07 + Approximate Standard Deviation of Response = 0.0000000000e+00 + Importance Factors not available. + +Warning: negligible standard deviation renders CDF results suspect. + +Cumulative Distribution Function (CDF) for MassFlux6: + Response Level Probability Level Reliability Index General Rel Index + -------------- ----------------- ----------------- ----------------- + -6.3466722620e-07 1.0000000000e-04 3.7190164855e+00 3.7190164855e+00 + -6.3466722620e-07 1.0000000000e-03 3.0902323062e+00 3.0902323062e+00 + -6.3466722620e-07 1.0000000000e-02 2.3263478740e+00 2.3263478740e+00 + -6.3466722620e-07 2.5000000000e-01 6.7448975020e-01 6.7448975020e-01 + -6.3466722620e-07 5.0000000000e-01 -0.0000000000e+00 -0.0000000000e+00 + -6.3466722620e-07 7.5000000000e-01 -6.7448975020e-01 -6.7448975020e-01 + -6.3466722620e-07 9.9000000000e-01 -2.3263478740e+00 -2.3263478740e+00 + -6.3466722620e-07 9.9900000000e-01 -3.0902323062e+00 -3.0902323062e+00 + -6.3466722620e-07 9.9990000000e-01 -3.7190164855e+00 -3.7190164855e+00 +----------------------------------------------------------------- Index: /issm/trunk/test/NightlyRun/test335.m =================================================================== --- /issm/trunk/test/NightlyRun/test335.m (revision 5015) +++ /issm/trunk/test/NightlyRun/test335.m (revision 5015) @@ -0,0 +1,8 @@ +md=mesh(model,'../Exp/Square.exp',180000); +md=geography(md,'',''); +md=parameterize(md,'../Par/SquareSheetConstrained.par'); +md=extrude(md,4,1); +md=setelementstype(md,'hutter','all'); +md.cluster='none'; +md.dt=0; +md=solve(md,'analysis_type',SteadystateSolutionEnum); Index: /issm/trunk/test/NightlyRun/test335_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test335_nightly.m (revision 5015) +++ /issm/trunk/test/NightlyRun/test335_nightly.m (revision 5015) @@ -0,0 +1,12 @@ +field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'}; +field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13 +}; +field_values={... + PatchToVec(md.results.SteadystateSolution.Vx),... + PatchToVec(md.results.SteadystateSolution.Vy),... + PatchToVec(md.results.SteadystateSolution.Vz),... + PatchToVec(md.results.SteadystateSolution.Vel),... + PatchToVec(md.results.SteadystateSolution.Pressure),... + PatchToVec(md.results.SteadystateSolution.Temperature),... + PatchToVec(md.results.SteadystateSolution.MeltingRate),... + }; Index: /issm/trunk/test/NightlyRun/test336.m =================================================================== --- /issm/trunk/test/NightlyRun/test336.m (revision 5015) +++ /issm/trunk/test/NightlyRun/test336.m (revision 5015) @@ -0,0 +1,8 @@ +md=mesh(model,'../Exp/Square.exp',180000); +md=geography(md,'',''); +md=parameterize(md,'../Par/SquareSheetConstrained.par'); +md=extrude(md,4,1); +md=setelementstype(md,'hutter','all'); +md.cluster=oshostname(); +md.dt=0; +md=solve(md,'analysis_type',SteadystateSolutionEnum); Index: /issm/trunk/test/NightlyRun/test336_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test336_nightly.m (revision 5015) +++ /issm/trunk/test/NightlyRun/test336_nightly.m (revision 5015) @@ -0,0 +1,12 @@ +field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'}; +field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13 +}; +field_values={... + PatchToVec(md.results.SteadystateSolution.Vx),... + PatchToVec(md.results.SteadystateSolution.Vy),... + PatchToVec(md.results.SteadystateSolution.Vz),... + PatchToVec(md.results.SteadystateSolution.Vel),... + PatchToVec(md.results.SteadystateSolution.Pressure),... + PatchToVec(md.results.SteadystateSolution.Temperature),... + PatchToVec(md.results.SteadystateSolution.MeltingRate),... + }; Index: /issm/trunk/test/NightlyRun/test419.m =================================================================== --- /issm/trunk/test/NightlyRun/test419.m (revision 5015) +++ /issm/trunk/test/NightlyRun/test419.m (revision 5015) @@ -0,0 +1,8 @@ +md=mesh(model,'../Exp/Square.exp',150000); +md=geography(md,'../Exp/SquareShelf.exp',''); +md=parameterize(md,'../Par/SquareSheetShelf.par'); +md=extrude(md,3,2); +md=setelementstype(md,'macayeal','all'); +md.cluster='none'; +md.dt=0; +md=solve(md,'analysis_type',SteadystateSolutionEnum); Index: /issm/trunk/test/NightlyRun/test419_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test419_nightly.m (revision 5015) +++ /issm/trunk/test/NightlyRun/test419_nightly.m (revision 5015) @@ -0,0 +1,12 @@ +field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'}; +field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13 +}; +field_values={... + PatchToVec(md.results.SteadystateSolution.Vx),... + PatchToVec(md.results.SteadystateSolution.Vy),... + PatchToVec(md.results.SteadystateSolution.Vz),... + PatchToVec(md.results.SteadystateSolution.Vel),... + PatchToVec(md.results.SteadystateSolution.Pressure),... + PatchToVec(md.results.SteadystateSolution.Temperature),... + PatchToVec(md.results.SteadystateSolution.MeltingRate),... + }; Index: /issm/trunk/test/NightlyRun/test420.m =================================================================== --- /issm/trunk/test/NightlyRun/test420.m (revision 5015) +++ /issm/trunk/test/NightlyRun/test420.m (revision 5015) @@ -0,0 +1,8 @@ +md=mesh(model,'../Exp/Square.exp',150000); +md=geography(md,'../Exp/SquareShelf.exp',''); +md=parameterize(md,'../Par/SquareSheetShelf.par'); +md=extrude(md,3,2); +md=setelementstype(md,'macayeal','all'); +md.cluster=oshostname(); +md.dt=0; +md=solve(md,'analysis_type',SteadystateSolutionEnum); Index: /issm/trunk/test/NightlyRun/test420_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test420_nightly.m (revision 5014) +++ /issm/trunk/test/NightlyRun/test420_nightly.m (revision 5015) @@ -1,7 +1,12 @@ -md=tres(md,'dakota'); -md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')'; -field_names ={'importancefactors'}; -field_tolerances={1e-12}; +field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'}; +field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13 +}; field_values={... - md.results.dakota.importancefactors,... + PatchToVec(md.results.SteadystateSolution.Vx),... + PatchToVec(md.results.SteadystateSolution.Vy),... + PatchToVec(md.results.SteadystateSolution.Vz),... + PatchToVec(md.results.SteadystateSolution.Vel),... + PatchToVec(md.results.SteadystateSolution.Pressure),... + PatchToVec(md.results.SteadystateSolution.Temperature),... + PatchToVec(md.results.SteadystateSolution.MeltingRate),... }; Index: /issm/trunk/test/NightlyRun/test421.m =================================================================== --- /issm/trunk/test/NightlyRun/test421.m (revision 5015) +++ /issm/trunk/test/NightlyRun/test421.m (revision 5015) @@ -0,0 +1,8 @@ +md=mesh(model,'../Exp/Square.exp',150000); +md=geography(md,'../Exp/SquareShelf.exp',''); +md=parameterize(md,'../Par/SquareSheetShelf.par'); +md=extrude(md,3,2); +md=setelementstype(md,'pattyn','all'); +md.cluster='none'; +md.dt=0; +md=solve(md,'analysis_type',SteadystateSolutionEnum); Index: /issm/trunk/test/NightlyRun/test421_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test421_nightly.m (revision 5015) +++ /issm/trunk/test/NightlyRun/test421_nightly.m (revision 5015) @@ -0,0 +1,11 @@ +field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'}; +field_tolerances={1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-06}; +field_values={... + PatchToVec(md.results.SteadystateSolution.Vx),... + PatchToVec(md.results.SteadystateSolution.Vy),... + PatchToVec(md.results.SteadystateSolution.Vz),... + PatchToVec(md.results.SteadystateSolution.Vel),... + PatchToVec(md.results.SteadystateSolution.Pressure),... + PatchToVec(md.results.SteadystateSolution.Temperature),... + PatchToVec(md.results.SteadystateSolution.MeltingRate),... + }; Index: /issm/trunk/test/NightlyRun/test422.m =================================================================== --- /issm/trunk/test/NightlyRun/test422.m (revision 5015) +++ /issm/trunk/test/NightlyRun/test422.m (revision 5015) @@ -0,0 +1,8 @@ +md=mesh(model,'../Exp/Square.exp',150000); +md=geography(md,'../Exp/SquareShelf.exp',''); +md=parameterize(md,'../Par/SquareSheetShelf.par'); +md=extrude(md,3,2); +md=setelementstype(md,'pattyn','all'); +md.cluster=oshostname(); +md.dt=0; +md=solve(md,'analysis_type',SteadystateSolutionEnum); Index: /issm/trunk/test/NightlyRun/test422_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test422_nightly.m (revision 5015) +++ /issm/trunk/test/NightlyRun/test422_nightly.m (revision 5015) @@ -0,0 +1,12 @@ +field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'}; +field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09 +}; +field_values={... + PatchToVec(md.results.SteadystateSolution.Vx),... + PatchToVec(md.results.SteadystateSolution.Vy),... + PatchToVec(md.results.SteadystateSolution.Vz),... + PatchToVec(md.results.SteadystateSolution.Vel),... + PatchToVec(md.results.SteadystateSolution.Pressure),... + PatchToVec(md.results.SteadystateSolution.Temperature),... + PatchToVec(md.results.SteadystateSolution.MeltingRate),... + }; Index: /issm/trunk/test/NightlyRun/test423.m =================================================================== --- /issm/trunk/test/NightlyRun/test423.m (revision 5014) +++ /issm/trunk/test/NightlyRun/test423.m (revision 5015) @@ -2,7 +2,35 @@ md=geography(md,'../Exp/SquareShelf.exp',''); md=parameterize(md,'../Par/SquareSheetShelf.par'); -md=extrude(md,3,2); md=setelementstype(md,'macayeal','all'); md.cluster='none'; -md.dt=0; -md=solve(md,'analysis_type',SteadystateSolutionEnum); + +%Dakota options +md.variables.nuv=normal_uncertain.empty(); +md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01); +md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01); + +md.responses=struct(); +md.responses.rf =response_function.empty(); +md.responses.rf (end+1)=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + +md.qmu_method =dakota_method('nond_l'); + +md.qmu_params.direct=true; +md.qmu_params.analysis_driver='diagnostic'; +md.qmu_params.evaluation_concurrency=1; +md.qmu_params.interval_type='forward'; +md.qmu_params.tabular_graphics_data=true; + +md.part=0:1:md.numberofgrids-1; +md.part=md.part'; +md.npart=md.numberofgrids; +md.qmu_analysis=1; + +md.eps_rel=10^-10; %tighten for qmu analysese + +%rm existing results +if exist('./qmu') + system('rm -rf qmu/'); +end + +md=solve(md,'analysis_type',DiagnosticSolutionEnum); Index: /issm/trunk/test/NightlyRun/test423_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test423_nightly.m (revision 5014) +++ /issm/trunk/test/NightlyRun/test423_nightly.m (revision 5015) @@ -1,12 +1,7 @@ -field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'}; -field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13 -}; +md=tres(md,'dakota'); +md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')'; +field_names ={'importancefactors'}; +field_tolerances={1e-11}; field_values={... - PatchToVec(md.results.SteadystateSolution.Vx),... - PatchToVec(md.results.SteadystateSolution.Vy),... - PatchToVec(md.results.SteadystateSolution.Vz),... - PatchToVec(md.results.SteadystateSolution.Vel),... - PatchToVec(md.results.SteadystateSolution.Pressure),... - PatchToVec(md.results.SteadystateSolution.Temperature),... - PatchToVec(md.results.SteadystateSolution.MeltingRate),... + md.results.dakota.importancefactors,... }; Index: /issm/trunk/test/NightlyRun/test424.m =================================================================== --- /issm/trunk/test/NightlyRun/test424.m (revision 5014) +++ /issm/trunk/test/NightlyRun/test424.m (revision 5015) @@ -2,7 +2,35 @@ md=geography(md,'../Exp/SquareShelf.exp',''); md=parameterize(md,'../Par/SquareSheetShelf.par'); -md=extrude(md,3,2); md=setelementstype(md,'macayeal','all'); md.cluster=oshostname(); -md.dt=0; -md=solve(md,'analysis_type',SteadystateSolutionEnum); + +%Dakota options +md.variables.nuv=normal_uncertain.empty(); +md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01); +md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01); + +md.responses=struct(); +md.responses.rf =response_function.empty(); +md.responses.rf (end+1)=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + +md.qmu_method =dakota_method('nond_l'); + +md.qmu_params.direct=true; +md.qmu_params.analysis_driver='diagnostic'; +md.qmu_params.evaluation_concurrency=1; +md.qmu_params.interval_type='forward'; +md.qmu_params.tabular_graphics_data=true; + +md.part=0:1:md.numberofgrids-1; +md.part=md.part'; +md.npart=md.numberofgrids; +md.qmu_analysis=1; + +md.eps_rel=10^-10; %tighten for qmu analysese + +%rm existing results +if exist('./qmu') + system('rm -rf qmu/'); +end + +md=solve(md,'analysis_type',DiagnosticSolutionEnum); Index: /issm/trunk/test/NightlyRun/test424_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test424_nightly.m (revision 5014) +++ /issm/trunk/test/NightlyRun/test424_nightly.m (revision 5015) @@ -1,12 +1,7 @@ -field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'}; -field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13 -}; +md=tres(md,'dakota'); +md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')'; +field_names ={'importancefactors'}; +field_tolerances={1e-11}; field_values={... - PatchToVec(md.results.SteadystateSolution.Vx),... - PatchToVec(md.results.SteadystateSolution.Vy),... - PatchToVec(md.results.SteadystateSolution.Vz),... - PatchToVec(md.results.SteadystateSolution.Vel),... - PatchToVec(md.results.SteadystateSolution.Pressure),... - PatchToVec(md.results.SteadystateSolution.Temperature),... - PatchToVec(md.results.SteadystateSolution.MeltingRate),... + md.results.dakota.importancefactors,... }; Index: /issm/trunk/test/NightlyRun/test425.m =================================================================== --- /issm/trunk/test/NightlyRun/test425.m (revision 5014) +++ /issm/trunk/test/NightlyRun/test425.m (revision 5015) @@ -2,7 +2,36 @@ md=geography(md,'../Exp/SquareShelf.exp',''); md=parameterize(md,'../Par/SquareSheetShelf.par'); -md=extrude(md,3,2); -md=setelementstype(md,'pattyn','all'); +md=setelementstype(md,'macayeal','all'); md.cluster='none'; -md.dt=0; -md=solve(md,'analysis_type',SteadystateSolutionEnum); + +%Dakota options +md.variables.nuv=normal_uncertain.empty(); +md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01); +md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01); + +md.responses=struct(); +md.responses.rf =response_function.empty(); +md.responses.rf (end+1)=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + +md.qmu_method =dakota_method('nond_l'); + +md.qmu_params.direct=true; +md.qmu_params.analysis_driver='diagnostic'; +md.qmu_params.evaluation_concurrency=1; +md.qmu_params.interval_type='forward'; +md.qmu_params.tabular_graphics_data=true; + +%partitioning +md.npart=20; +md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on'); +md.part=md.part-1; +md.qmu_analysis=1; + +md.eps_rel=10^-10; %tighten for qmu analyses + +%rm existing results +if exist('./qmu') + system('rm -rf qmu/'); +end + +md=solve(md,'analysis_type',DiagnosticSolutionEnum); Index: /issm/trunk/test/NightlyRun/test425_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test425_nightly.m (revision 5014) +++ /issm/trunk/test/NightlyRun/test425_nightly.m (revision 5015) @@ -1,11 +1,7 @@ -field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'}; -field_tolerances={1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-06}; +md=tres(md,'dakota'); +md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')'; +field_names ={'importancefactors'}; +field_tolerances={1e-11}; field_values={... - PatchToVec(md.results.SteadystateSolution.Vx),... - PatchToVec(md.results.SteadystateSolution.Vy),... - PatchToVec(md.results.SteadystateSolution.Vz),... - PatchToVec(md.results.SteadystateSolution.Vel),... - PatchToVec(md.results.SteadystateSolution.Pressure),... - PatchToVec(md.results.SteadystateSolution.Temperature),... - PatchToVec(md.results.SteadystateSolution.MeltingRate),... + md.results.dakota.importancefactors,... }; Index: /issm/trunk/test/NightlyRun/test426.m =================================================================== --- /issm/trunk/test/NightlyRun/test426.m (revision 5014) +++ /issm/trunk/test/NightlyRun/test426.m (revision 5015) @@ -2,7 +2,36 @@ md=geography(md,'../Exp/SquareShelf.exp',''); md=parameterize(md,'../Par/SquareSheetShelf.par'); -md=extrude(md,3,2); -md=setelementstype(md,'pattyn','all'); +md=setelementstype(md,'macayeal','all'); md.cluster=oshostname(); -md.dt=0; -md=solve(md,'analysis_type',SteadystateSolutionEnum); + +%Dakota options +md.variables.nuv=normal_uncertain.empty(); +md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01); +md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01); + +md.responses=struct(); +md.responses.rf =response_function.empty(); +md.responses.rf (end+1)=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + +md.qmu_method =dakota_method('nond_l'); + +md.qmu_params.direct=true; +md.qmu_params.analysis_driver='diagnostic'; +md.qmu_params.evaluation_concurrency=1; +md.qmu_params.interval_type='forward'; +md.qmu_params.tabular_graphics_data=true; + +%partitioning +md.npart=20; +md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on'); +md.part=md.part-1; +md.qmu_analysis=1; + +md.eps_rel=10^-10; %tighten for qmu analyses + +%rm existing results +if exist('./qmu') + system('rm -rf qmu/'); +end + +md=solve(md,'analysis_type',DiagnosticSolutionEnum); Index: /issm/trunk/test/NightlyRun/test426_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test426_nightly.m (revision 5014) +++ /issm/trunk/test/NightlyRun/test426_nightly.m (revision 5015) @@ -1,12 +1,7 @@ -field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'}; -field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09 -}; +md=tres(md,'dakota'); +md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')'; +field_names ={'importancefactors'}; +field_tolerances={1e-10}; field_values={... - PatchToVec(md.results.SteadystateSolution.Vx),... - PatchToVec(md.results.SteadystateSolution.Vy),... - PatchToVec(md.results.SteadystateSolution.Vz),... - PatchToVec(md.results.SteadystateSolution.Vel),... - PatchToVec(md.results.SteadystateSolution.Pressure),... - PatchToVec(md.results.SteadystateSolution.Temperature),... - PatchToVec(md.results.SteadystateSolution.MeltingRate),... + md.results.dakota.importancefactors,... }; Index: /issm/trunk/test/NightlyRun/test427.m =================================================================== --- /issm/trunk/test/NightlyRun/test427.m (revision 5014) +++ /issm/trunk/test/NightlyRun/test427.m (revision 5015) @@ -5,7 +5,12 @@ md.cluster='none'; +%constrain all velocities to 1 m/yr, in the y-direction +md.spcvelocity(:,1:3)=1; +md.spcvelocity(:,4)=0; +md.spcvelocity(:,5)=1*md.yts; +md.spcvelocity(:,6)=0; + %Dakota options md.variables.nuv=normal_uncertain.empty(); -md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01); md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01); @@ -13,4 +18,14 @@ md.responses.rf =response_function.empty(); md.responses.rf (end+1)=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); +md.responses.rf (end+1)=response_function('MassFlux1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); +md.responses.rf (end+1)=response_function('MassFlux2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); +md.responses.rf (end+1)=response_function('MassFlux3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); +md.responses.rf (end+1)=response_function('MassFlux4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); +md.responses.rf (end+1)=response_function('MassFlux5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); +md.responses.rf (end+1)=response_function('MassFlux6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + +%mass flux profiles +md.qmu_mass_flux_profiles={'../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp'}; +md.qmu_mass_flux_profile_directory=pwd; md.qmu_method =dakota_method('nond_l'); @@ -22,10 +37,11 @@ md.qmu_params.tabular_graphics_data=true; -md.part=0:1:md.numberofgrids-1; -md.part=md.part'; -md.npart=md.numberofgrids; +%partitioning +md.npart=20; +md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on'); +md.part=md.part-1; md.qmu_analysis=1; -md.eps_rel=10^-10; %tighten for qmu analysese +md.eps_rel=10^-10; %tighten for qmu analyses %rm existing results Index: /issm/trunk/test/NightlyRun/test427_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test427_nightly.m (revision 5014) +++ /issm/trunk/test/NightlyRun/test427_nightly.m (revision 5015) @@ -1,4 +1,15 @@ md=tres(md,'dakota'); -md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')'; +%ok, mass flux of 3 profiles should be -.3, -.3 and the sum, which is -.6 +%we recover those mass fluxes through the mean of the response. +%also, we recover the max velo, which should be 1m/s. +%we put all that data in the importancefactors, which we will use to test for success. +%also, check that the stddev are 0. +md.results.dakota.importancefactors=[]; +for i=1:7, + md.results.dakota.importancefactors=[md.results.dakota.importancefactors md.results.dakota.dresp_out(i).mean]; +end +for i=1:7, + md.results.dakota.importancefactors=[md.results.dakota.importancefactors md.results.dakota.dresp_out(i).stddev]; +end field_names ={'importancefactors'}; field_tolerances={1e-11}; Index: /issm/trunk/test/NightlyRun/test428.m =================================================================== --- /issm/trunk/test/NightlyRun/test428.m (revision 5014) +++ /issm/trunk/test/NightlyRun/test428.m (revision 5015) @@ -5,7 +5,12 @@ md.cluster=oshostname(); +%constrain all velocities to 1 m/yr, in the y-direction +md.spcvelocity(:,1:3)=1; +md.spcvelocity(:,4)=0; +md.spcvelocity(:,5)=1*md.yts; +md.spcvelocity(:,6)=0; + %Dakota options md.variables.nuv=normal_uncertain.empty(); -md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01); md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01); @@ -13,4 +18,14 @@ md.responses.rf =response_function.empty(); md.responses.rf (end+1)=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); +md.responses.rf (end+1)=response_function('MassFlux1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); +md.responses.rf (end+1)=response_function('MassFlux2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); +md.responses.rf (end+1)=response_function('MassFlux3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); +md.responses.rf (end+1)=response_function('MassFlux4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); +md.responses.rf (end+1)=response_function('MassFlux5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); +md.responses.rf (end+1)=response_function('MassFlux6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + +%mass flux profiles +md.qmu_mass_flux_profiles={'../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp'}; +md.qmu_mass_flux_profile_directory=pwd; md.qmu_method =dakota_method('nond_l'); @@ -22,10 +37,11 @@ md.qmu_params.tabular_graphics_data=true; -md.part=0:1:md.numberofgrids-1; -md.part=md.part'; -md.npart=md.numberofgrids; +%partitioning +md.npart=20; +md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on'); +md.part=md.part-1; md.qmu_analysis=1; -md.eps_rel=10^-10; %tighten for qmu analysese +md.eps_rel=10^-10; %tighten for qmu analyses %rm existing results Index: /issm/trunk/test/NightlyRun/test428_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test428_nightly.m (revision 5014) +++ /issm/trunk/test/NightlyRun/test428_nightly.m (revision 5015) @@ -1,4 +1,15 @@ md=tres(md,'dakota'); -md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')'; +%ok, mass flux of 3 profiles should be -.3, -.3 and the sum, which is -.6 +%we recover those mass fluxes through the mean of the response. +%also, we recover the max velo, which should be 1m/s. +%we put all that data in the importancefactors, which we will use to test for success. +%also, check that the stddev are 0. +md.results.dakota.importancefactors=[]; +for i=1:7, + md.results.dakota.importancefactors=[md.results.dakota.importancefactors md.results.dakota.dresp_out(i).mean]; +end +for i=1:7, + md.results.dakota.importancefactors=[md.results.dakota.importancefactors md.results.dakota.dresp_out(i).stddev]; +end field_names ={'importancefactors'}; field_tolerances={1e-11}; @@ -6,2 +17,3 @@ md.results.dakota.importancefactors,... }; +md.results.dakota.importancefactors