Changeset 9652
- Timestamp:
- 09/07/11 10:28:30 (14 years ago)
- Location:
- issm/trunk/test/NightlyRun
- Files:
-
- 12 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/NightlyRun/test235.m
r9641 r9652 56 56 57 57 %partitioning 58 md. npart=md.numberofnodes;59 md =partitioner(md,'package','linear','npart',md.npart);60 md. part=md.part-1;58 md.qmu.numberofpartitions=md.numberofnodes; 59 md.qmu.partitionitioner(md,'package','linear','npart',md.qmu.numberofpartitions); 60 md.qmu.partition=md.qmu.partition-1; 61 61 62 62 %Dakota options 63 63 %variables 64 md. variables.rheology_B=normal_uncertain('scaled_MaterialsRheologyB',1,.5);64 md.qmu.variables.rheology_B=normal_uncertain('scaled_MaterialsRheologyB',1,.5); 65 65 66 66 %responses 67 md. responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);67 md.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 68 68 69 69 %method 70 md.qmu _method =dakota_method('nond_l');70 md.qmu.method =dakota_method('nond_l'); 71 71 72 72 %parameters 73 md.qmu _params.direct=true;74 md.qmu _params.analysis_driver='diagnostic';75 md.qmu _params.evaluation_concurrency=1;76 md.qmu _params.interval_type='forward';73 md.qmu.params.direct=true; 74 md.qmu.params.analysis_driver='diagnostic'; 75 md.qmu.params.evaluation_concurrency=1; 76 md.qmu.params.interval_type='forward'; 77 77 78 78 %imperative! 79 79 md.eps_rel=10^-10; %tighten for qmu analysese 80 md. dakota_analysis=1;80 md.qmu.isdakota=1; 81 81 82 82 %solve -
issm/trunk/test/NightlyRun/test236.m
r9641 r9652 56 56 57 57 %partitioning 58 md. npart=md.numberofnodes;59 md =partitioner(md,'package','linear','npart',md.npart);60 md. part=md.part-1;58 md.qmu.numberofpartitions=md.numberofnodes; 59 md.qmu.partitionitioner(md,'package','linear','npart',md.qmu.numberofpartitions); 60 md.qmu.partition=md.qmu.partition-1; 61 61 62 62 %Dakota options 63 63 %variables 64 md. variables.rheology_B=normal_uncertain('scaled_MaterialsRheologyB',1,.05);64 md.qmu.variables.rheology_B=normal_uncertain('scaled_MaterialsRheologyB',1,.05); 65 65 66 66 %responses 67 md. responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);67 md.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 68 68 69 69 %method 70 md.qmu _method =dakota_method('nond_l');70 md.qmu.method =dakota_method('nond_l'); 71 71 72 72 %parameters 73 md.qmu _params.direct=true;74 md.qmu _params.analysis_driver='diagnostic';75 md.qmu _params.evaluation_concurrency=1;76 md.qmu _params.interval_type='forward';73 md.qmu.params.direct=true; 74 md.qmu.params.analysis_driver='diagnostic'; 75 md.qmu.params.evaluation_concurrency=1; 76 md.qmu.params.interval_type='forward'; 77 77 78 78 %imperative! 79 79 md.eps_rel=10^-10; %tighten for qmu analysese 80 md. dakota_analysis=1;80 md.qmu.isdakota=1; 81 81 82 82 %solve -
issm/trunk/test/NightlyRun/test423.m
r9641 r9652 6 6 7 7 %partitioning 8 md. npart=md.numberofnodes;9 md =partitioner(md,'package','linear','npart',md.npart);10 md. part=md.part-1;11 md. dakota_analysis=1;8 md.qmu.numberofpartitions=md.numberofnodes; 9 md.qmu.partitionitioner(md,'package','linear','npart',md.qmu.numberofpartitions); 10 md.qmu.partition=md.qmu.partition-1; 11 md.qmu.isdakota=1; 12 12 13 13 %Dakota options 14 14 %variables 15 md. variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);16 md. variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,.01);17 %md. variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,.01*ones(md.npart,1));15 md.qmu.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01); 16 md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,.01); 17 %md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,.01*ones(md.qmu.numberofpartitions,1)); 18 18 19 19 %responses 20 md. responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);20 md.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 21 21 22 22 %method 23 md.qmu _method =dakota_method('nond_l');23 md.qmu.method =dakota_method('nond_l'); 24 24 25 25 %parameters 26 md.qmu _params.direct=true;27 md.qmu _params.analysis_driver='diagnostic';28 md.qmu _params.evaluation_concurrency=1;29 md.qmu _params.interval_type='forward';26 md.qmu.params.direct=true; 27 md.qmu.params.analysis_driver='diagnostic'; 28 md.qmu.params.evaluation_concurrency=1; 29 md.qmu.params.interval_type='forward'; 30 30 31 31 -
issm/trunk/test/NightlyRun/test424.m
r9641 r9652 6 6 7 7 %partitioning 8 md. npart=md.numberofnodes;9 md =partitioner(md,'package','linear','npart',md.npart);10 md. part=md.part-1;11 md. dakota_analysis=1;8 md.qmu.numberofpartitions=md.numberofnodes; 9 md.qmu.partitionitioner(md,'package','linear','npart',md.qmu.numberofpartitions); 10 md.qmu.partition=md.qmu.partition-1; 11 md.qmu.isdakota=1; 12 12 13 13 14 14 %Dakota options 15 15 %variables 16 md. variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);17 md. variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01);16 md.qmu.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01); 17 md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01); 18 18 19 19 %responses 20 md. responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);20 md.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 21 21 22 22 %method 23 md.qmu _method =dakota_method('nond_l');23 md.qmu.method =dakota_method('nond_l'); 24 24 25 25 %parameters 26 md.qmu _params.direct=true;27 md.qmu _params.analysis_driver='diagnostic';28 md.qmu _params.evaluation_concurrency=1;29 md.qmu _params.interval_type='forward';26 md.qmu.params.direct=true; 27 md.qmu.params.analysis_driver='diagnostic'; 28 md.qmu.params.evaluation_concurrency=1; 29 md.qmu.params.interval_type='forward'; 30 30 31 31 %imperative! -
issm/trunk/test/NightlyRun/test425.m
r9641 r9652 8 8 9 9 %partitioning 10 md. npart=20;11 md =partitioner(md,'package','chaco','npart',md.npart,'weighting','on');12 md. part=md.part-1;10 md.qmu.numberofpartitions=20; 11 md.qmu.partitionitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on'); 12 md.qmu.partition=md.qmu.partition-1; 13 13 14 14 %variables 15 md. variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);16 md. variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01);15 md.qmu.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01); 16 md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01); 17 17 18 18 %responses 19 md. responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);19 md.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 20 20 21 21 %method 22 md.qmu _method =dakota_method('nond_l');22 md.qmu.method =dakota_method('nond_l'); 23 23 24 24 %parameters 25 md.qmu _params.direct=true;26 md.qmu _params.analysis_driver='diagnostic';27 md.qmu _params.evaluation_concurrency=1;28 md.qmu _params.interval_type='forward';25 md.qmu.params.direct=true; 26 md.qmu.params.analysis_driver='diagnostic'; 27 md.qmu.params.evaluation_concurrency=1; 28 md.qmu.params.interval_type='forward'; 29 29 30 30 31 31 %imperative! 32 32 md.eps_rel=10^-5; %tighten for qmu analyses 33 md. dakota_analysis=1;33 md.qmu.isdakota=1; 34 34 35 35 %solve -
issm/trunk/test/NightlyRun/test426.m
r9641 r9652 8 8 9 9 %partitioning 10 md. npart=20;11 md =partitioner(md,'package','chaco','npart',md.npart,'weighting','on');12 md. part=md.part-1;10 md.qmu.numberofpartitions=20; 11 md.qmu.partitionitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on'); 12 md.qmu.partition=md.qmu.partition-1; 13 13 14 14 %variables 15 md. variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);16 md. variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01);15 md.qmu.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01); 16 md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01); 17 17 18 18 %responses 19 md. responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);19 md.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 20 20 21 21 %method 22 md.qmu _method =dakota_method('nond_l');22 md.qmu.method =dakota_method('nond_l'); 23 23 24 24 %parameters 25 md.qmu _params.direct=true;26 md.qmu _params.analysis_driver='diagnostic';27 md.qmu _params.evaluation_concurrency=1;28 md.qmu _params.interval_type='forward';25 md.qmu.params.direct=true; 26 md.qmu.params.analysis_driver='diagnostic'; 27 md.qmu.params.evaluation_concurrency=1; 28 md.qmu.params.interval_type='forward'; 29 29 30 30 31 31 %imperative! 32 32 md.eps_rel=10^-5; %tighten for qmu analyses 33 md. dakota_analysis=1;33 md.qmu.isdakota=1; 34 34 35 35 %solve -
issm/trunk/test/NightlyRun/test427.m
r9641 r9652 15 15 16 16 %partitioning 17 md. npart=20;18 md =partitioner(md,'package','chaco','npart',md.npart,'weighting','on');19 md. part=md.part-1;17 md.qmu.numberofpartitions=20; 18 md.qmu.partitionitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on'); 19 md.qmu.partition=md.qmu.partition-1; 20 20 21 21 %variables 22 md. variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01);22 md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01); 23 23 24 24 %responses 25 md. responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);26 md. responses.MassFlux1=response_function('indexed_MassFlux_1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);27 md. responses.MassFlux2=response_function('indexed_MassFlux_2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);28 md. responses.MassFlux3=response_function('indexed_MassFlux_3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);29 md. responses.MassFlux4=response_function('indexed_MassFlux_4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);30 md. responses.MassFlux5=response_function('indexed_MassFlux_5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);31 md. responses.massFlux6=response_function('indexed_MassFlux_6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);25 md.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 26 md.qmu.responses.MassFlux1=response_function('indexed_MassFlux_1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 27 md.qmu.responses.MassFlux2=response_function('indexed_MassFlux_2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 28 md.qmu.responses.MassFlux3=response_function('indexed_MassFlux_3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 29 md.qmu.responses.MassFlux4=response_function('indexed_MassFlux_4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 30 md.qmu.responses.MassFlux5=response_function('indexed_MassFlux_5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 31 md.qmu.responses.massFlux6=response_function('indexed_MassFlux_6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 32 32 33 33 %mass flux profiles 34 md.qmu _mass_flux_profiles={'../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp'};35 md.qmu _mass_flux_profile_directory=pwd;34 md.qmu.mass_flux_profiles={'../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp'}; 35 md.qmu.mass_flux_profile_directory=pwd; 36 36 37 37 %method 38 md.qmu _method =dakota_method('nond_l');38 md.qmu.method =dakota_method('nond_l'); 39 39 40 40 %parameters 41 md.qmu _params.direct=true;42 md.qmu _params.analysis_driver='diagnostic';43 md.qmu _params.evaluation_concurrency=1;44 md.qmu _params.interval_type='forward';45 md. dakota_analysis=1;41 md.qmu.params.direct=true; 42 md.qmu.params.analysis_driver='diagnostic'; 43 md.qmu.params.evaluation_concurrency=1; 44 md.qmu.params.interval_type='forward'; 45 md.qmu.isdakota=1; 46 46 md.eps_rel=10^-5; %tighten for qmu analyses 47 47 -
issm/trunk/test/NightlyRun/test428.m
r9641 r9652 15 15 16 16 %partitioning 17 md. npart=20;18 md =partitioner(md,'package','chaco','npart',md.npart,'weighting','on');19 md. part=md.part-1;17 md.qmu.numberofpartitions=20; 18 md.qmu.partitionitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on'); 19 md.qmu.partition=md.qmu.partition-1; 20 20 21 21 %variables 22 md. variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01);22 md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01); 23 23 24 24 %responses 25 md. responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);26 md. responses.MassFlux1=response_function('indexed_MassFlux_1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);27 md. responses.MassFlux2=response_function('indexed_MassFlux_2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);28 md. responses.MassFlux3=response_function('indexed_MassFlux_3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);29 md. responses.MassFlux4=response_function('indexed_MassFlux_4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);30 md. responses.MassFlux5=response_function('indexed_MassFlux_5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);31 md. responses.MassFlux6=response_function('indexed_MassFlux_6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);25 md.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 26 md.qmu.responses.MassFlux1=response_function('indexed_MassFlux_1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 27 md.qmu.responses.MassFlux2=response_function('indexed_MassFlux_2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 28 md.qmu.responses.MassFlux3=response_function('indexed_MassFlux_3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 29 md.qmu.responses.MassFlux4=response_function('indexed_MassFlux_4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 30 md.qmu.responses.MassFlux5=response_function('indexed_MassFlux_5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 31 md.qmu.responses.MassFlux6=response_function('indexed_MassFlux_6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 32 32 33 33 %mass flux profiles 34 md.qmu _mass_flux_profiles={'../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp'};35 md.qmu _mass_flux_profile_directory=pwd;34 md.qmu.mass_flux_profiles={'../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp'}; 35 md.qmu.mass_flux_profile_directory=pwd; 36 36 37 37 %method 38 md.qmu _method =dakota_method('nond_l');38 md.qmu.method =dakota_method('nond_l'); 39 39 40 40 %parameters 41 md.qmu _params.direct=true;42 md.qmu _params.analysis_driver='diagnostic';43 md.qmu _params.evaluation_concurrency=1;44 md.qmu _params.interval_type='forward';45 md. dakota_analysis=1;41 md.qmu.params.direct=true; 42 md.qmu.params.analysis_driver='diagnostic'; 43 md.qmu.params.evaluation_concurrency=1; 44 md.qmu.params.interval_type='forward'; 45 md.qmu.isdakota=1; 46 46 md.eps_rel=10^-5; %tighten for qmu analyses 47 47 -
issm/trunk/test/NightlyRun/test434.m
r9641 r9652 13 13 14 14 %Dakota options 15 md. variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01);15 md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01); 16 16 17 md. responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);18 md. responses.MassFlux1=response_function('indexed_MassFlux_1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);19 md. responses.MassFlux2=response_function('indexed_MassFlux_2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);20 md. responses.MassFlux3=response_function('indexed_MassFlux_3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);21 md. responses.MassFlux4=response_function('indexed_MassFlux_4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);22 md. responses.MassFlux5=response_function('indexed_MassFlux_5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);23 md. responses.MassFlux6=response_function('indexed_MassFlux_6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);17 md.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 18 md.qmu.responses.MassFlux1=response_function('indexed_MassFlux_1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 19 md.qmu.responses.MassFlux2=response_function('indexed_MassFlux_2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 20 md.qmu.responses.MassFlux3=response_function('indexed_MassFlux_3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 21 md.qmu.responses.MassFlux4=response_function('indexed_MassFlux_4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 22 md.qmu.responses.MassFlux5=response_function('indexed_MassFlux_5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 23 md.qmu.responses.MassFlux6=response_function('indexed_MassFlux_6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 24 24 25 25 %mass flux profiles 26 md.qmu _mass_flux_profiles={'../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp'};27 md.qmu _mass_flux_profile_directory=pwd;26 md.qmu.mass_flux_profiles={'../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp'}; 27 md.qmu.mass_flux_profile_directory=pwd; 28 28 29 29 30 30 %% nond_sampling study 31 31 32 md.qmu _method =dakota_method('nond_samp');33 md.qmu _method(end)=dmeth_params_set(md.qmu_method(end),...32 md.qmu.method =dakota_method('nond_samp'); 33 md.qmu.method(end)=dmeth_params_set(md.qmu.method(end),... 34 34 'seed',1234,... 35 35 'samples',20,... … … 37 37 38 38 %% a variety of parameters 39 md.qmu _params.evaluation_concurrency=1;40 md.qmu _params.analysis_driver='';41 md.qmu _params.analysis_components='';39 md.qmu.params.evaluation_concurrency=1; 40 md.qmu.params.analysis_driver=''; 41 md.qmu.params.analysis_components=''; 42 42 43 43 %partitioning 44 md. npart=20;45 md =partitioner(md,'package','chaco','npart',md.npart,'weighting','on');46 md. part=md.part-1;47 md. dakota_analysis=1;44 md.qmu.numberofpartitions=20; 45 md.qmu.partitionitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on'); 46 md.qmu.partition=md.qmu.partition-1; 47 md.qmu.isdakota=1; 48 48 49 49 md.eps_rel=10^-5; %tighten for qmu analyses -
issm/trunk/test/NightlyRun/test435.m
r9641 r9652 7 7 8 8 %partitioning 9 md. npart=100;10 md =partitioner(md,'package','chaco','npart',md.npart);11 md. part=md.part-1;9 md.qmu.numberofpartitions=100; 10 md.qmu.partitionitioner(md,'package','chaco','npart',md.qmu.numberofpartitions); 11 md.qmu.partition=md.qmu.partition-1; 12 12 13 13 vector=(1:1:md.numberofnodes)'; 14 14 vector_on_partition=AreaAverageOntoPartition(md,vector); 15 vector_on_nodes=vector_on_partition(md. part+1);15 vector_on_nodes=vector_on_partition(md.qmu.partition+1); 16 16 17 17 field_names ={'vector_on_nodes'}; -
issm/trunk/test/NightlyRun/test439.m
r9641 r9652 6 6 7 7 %partitioning 8 md. npart=10;9 md =partitioner(md,'package','chaco','npart',md.npart);10 md. part=md.part-1;11 md. dakota_analysis=1;8 md.qmu.numberofpartitions=10; 9 md.qmu.partitionitioner(md,'package','chaco','npart',md.qmu.numberofpartitions); 10 md.qmu.partition=md.qmu.partition-1; 11 md.qmu.isdakota=1; 12 12 13 13 %Dakota options 14 14 %variables 15 md. variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);15 md.qmu.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01); 16 16 17 17 %responses 18 md. responses.MaxVel=response_function('scaled_Thickness',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);18 md.qmu.responses.MaxVel=response_function('scaled_Thickness',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 19 19 20 20 %method 21 md.qmu _method =dakota_method('nond_l');21 md.qmu.method =dakota_method('nond_l'); 22 22 23 23 %parameters 24 md.qmu _params.direct=true;25 md.qmu _params.analysis_driver='diagnostic';26 md.qmu _params.evaluation_concurrency=1;27 md.qmu _params.interval_type='forward';24 md.qmu.params.direct=true; 25 md.qmu.params.analysis_driver='diagnostic'; 26 md.qmu.params.evaluation_concurrency=1; 27 md.qmu.params.interval_type='forward'; 28 28 29 29 … … 36 36 37 37 %test on thickness 38 h=zeros(md. npart,1);39 for i=1:md. npart,38 h=zeros(md.qmu.numberofpartitions,1); 39 for i=1:md.qmu.numberofpartitions, 40 40 h(i)=md.dakotaresults.dresp_out(i).mean; 41 41 end 42 42 43 43 %project onto grid 44 thickness=h(md. part+1);44 thickness=h(md.qmu.partition+1); 45 45 46 46 %Fields and tolerances to track changes -
issm/trunk/test/NightlyRun/test440.m
r9641 r9652 6 6 7 7 %partitioning 8 md. npart=10;9 md =partitioner(md,'package','chaco','npart',md.npart);10 md. part=md.part-1;11 md. dakota_analysis=1;8 md.qmu.numberofpartitions=10; 9 md.qmu.partitionitioner(md,'package','chaco','npart',md.qmu.numberofpartitions); 10 md.qmu.partition=md.qmu.partition-1; 11 md.qmu.isdakota=1; 12 12 13 13 %Dakota options 14 14 %variables 15 md. variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);15 md.qmu.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01); 16 16 17 17 %responses 18 md. responses.MaxVel=response_function('scaled_Thickness',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);18 md.qmu.responses.MaxVel=response_function('scaled_Thickness',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 19 19 20 20 %method 21 md.qmu _method =dakota_method('nond_l');21 md.qmu.method =dakota_method('nond_l'); 22 22 23 23 %parameters 24 md.qmu _params.direct=true;25 md.qmu _params.analysis_driver='diagnostic';26 md.qmu _params.evaluation_concurrency=1;27 md.qmu _params.interval_type='forward';24 md.qmu.params.direct=true; 25 md.qmu.params.analysis_driver='diagnostic'; 26 md.qmu.params.evaluation_concurrency=1; 27 md.qmu.params.interval_type='forward'; 28 28 29 29 … … 36 36 37 37 %test on thickness 38 h=zeros(md. npart,1);39 for i=1:md. npart,38 h=zeros(md.qmu.numberofpartitions,1); 39 for i=1:md.qmu.numberofpartitions, 40 40 h(i)=md.dakotaresults.dresp_out(i).mean; 41 41 end 42 42 43 43 %project onto grid 44 thickness=h(md. part+1);44 thickness=h(md.qmu.partition+1); 45 45 46 46 %Fields and tolerances to track changes
Note:
See TracChangeset
for help on using the changeset viewer.