source: issm/trunk/src/m/qmu/examples/samp_direct.m

Last change on this file was 16137, checked in by Mathieu Morlighem, 12 years ago

merged trunk-jpl and trunk for revision 16135

File size: 2.3 KB
RevLine 
[834]1% set up a sampling study, like might be done in Pig.par
2
3%% a variety of variables
4
5% seems to be a Matlab bug here (on Linux, not WinXP) -- unless
6% the class has been called, "empty" method can not be found
7normal_uncertain;
8continuous_design;
9continuous_state;
10linear_inequality_constraint;
11linear_equality_constraint;
12response_function;
13objective_function;
14least_squares_term;
15nonlinear_inequality_constraint;
16nonlinear_equality_constraint;
17
[9650]18md.qmu.variables=struct();
19md.qmu.variables.nuv=normal_uncertain.empty();
20md.qmu.variables.nuv(end+1)=normal_uncertain('rho_ice',917,45.85);
21md.qmu.variables.nuv(end+1)=normal_uncertain('rho_water',1023,51.15);
22md.qmu.variables.nuv(end+1)=normal_uncertain('heatcapacity',2009,100.45);
23md.qmu.variables.nuv(end+1)=normal_uncertain('thermalconductivity',2.2,0.11);
24md.qmu.variables.nuv(end+1)=normal_uncertain('gravity',9.8,0.49);
25md.qmu.variables.nuv(end+1)=normal_uncertain('thickness',1,0.05);
26md.qmu.variables.nuv(end+1)=normal_uncertain('drag',1,0.05);
[834]27
28%% a variety of responses
29
[9650]30md.qmu.responses=struct();
31md.qmu.responses.rf =response_function.empty();
32md.qmu.responses.rf (end+1)=response_function('min_vx',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
33md.qmu.responses.rf (end+1)=response_function('max_vx',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
34md.qmu.responses.rf (end+1)=response_function('max_abs_vx',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
35md.qmu.responses.rf (end+1)=response_function('min_vy',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
36md.qmu.responses.rf (end+1)=response_function('max_vy',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
37md.qmu.responses.rf (end+1)=response_function('max_abs_vy',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
38md.qmu.responses.rf (end+1)=response_function('min_vel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
39md.qmu.responses.rf (end+1)=response_function('max_vel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
[834]40
41%% nond_sampling study
42
[9650]43md.qmu.method =dakota_method('nond_samp');
44md.qmu.method(end)=dmeth_params_set(md.qmu.method(end),...
[834]45 'seed',1234,...
46 'samples',20,...
47 'sample_type','lhs');
48
49%% a variety of parameters
50
[9650]51md.qmu.params.direct=true;
[16137]52md.qmu.params.analysis_driver='stressbalance';
[9650]53md.qmu.params.evaluation_concurrency=1;
[834]54
[9650]55md.qmu.numberofpartitions=10;
[834]56
57md.qmu
Note: See TracBrowser for help on using the repository browser.