Changeset 5193


Ignore:
Timestamp:
08/12/10 12:07:39 (15 years ago)
Author:
Eric.Larour
Message:

Using John's simplified naming conventions for qmu variables and responses

Location:
issm/trunk/test/NightlyRun
Files:
8 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/test/NightlyRun/test231.m

    r5119 r5193  
    77
    88%solver, to handle penalties
    9 md=solversettoasm(md);
    10 %md=solversettomumps(md);
     9%md=solversettoasm(md);
     10md=solversettomumps(md);
    1111
    1212%rift settings
  • issm/trunk/test/NightlyRun/test232.m

    r5119 r5193  
    55md=setelementstype(md,'macayeal','all');
    66md.cluster=oshostname();
    7 md.np=1
     7md.np=2
    88
    99%solver, to handle penalties
     
    1414md.rifts.fill=MelangeEnum();
    1515md.rifts.fraction=0;
    16 md.penalty_lock=2;
     16md.penalty_lock=10;
     17md.min_mechanical_constraints=0;
    1718md.rifts.fractionincrement=.1;
    18 
    1919md=solve(md,'analysis_type',DiagnosticSolutionEnum);
    2020
  • issm/trunk/test/NightlyRun/test423.m

    r5115 r5193  
    66
    77%Dakota options
    8 md.variables.nuv=normal_uncertain.empty();
    9 md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01);
    10 md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
     8%variables
     9md.variables.rho_ice=normal_uncertain('RhoIce',md.rho_ice,0.01);
     10md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01);
    1111
    12 md.responses=struct();
    13 md.responses.rf =response_function.empty();
    14 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]);
     12%responses
     13md.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    1514
     15%method
    1616md.qmu_method     =dakota_method('nond_l');
    1717
     18%parameters
    1819md.qmu_params.direct=true;
    1920md.qmu_params.analysis_driver='diagnostic';
     
    2223md.qmu_params.tabular_graphics_data=true;
    2324
     25%partitioning
    2426md.part=0:1:md.numberofgrids-1;
    2527md.part=md.part';
     
    2729md.qmu_analysis=1;
    2830
     31%imperative!
    2932md.eps_rel=10^-5; %tighten for qmu analysese
    3033
     34%solve
    3135md=solve(md,'analysis_type',DiagnosticSolutionEnum,'overwrite','y');
    3236
  • issm/trunk/test/NightlyRun/test424.m

    r5115 r5193  
    66
    77%Dakota options
    8 md.variables.nuv=normal_uncertain.empty();
    9 md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01);
    10 md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
     8%variables
     9md.variables.rho_ice=normal_uncertain('RhoIce',md.rho_ice,0.01);
     10md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01);
    1111
    12 md.responses=struct();
    13 md.responses.rf =response_function.empty();
    14 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]);
     12%responses
     13md.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    1514
     15%method
    1616md.qmu_method     =dakota_method('nond_l');
    1717
     18%parameters
    1819md.qmu_params.direct=true;
    1920md.qmu_params.analysis_driver='diagnostic';
     
    2223md.qmu_params.tabular_graphics_data=true;
    2324
     25%partitioning
    2426md.part=0:1:md.numberofgrids-1;
    2527md.part=md.part';
     
    2729md.qmu_analysis=1;
    2830
     31%imperative!
    2932md.eps_rel=10^-5; %tighten for qmu analysese
    3033
     34%solve
    3135md=solve(md,'analysis_type',DiagnosticSolutionEnum,'overwrite','y');
    3236
  • issm/trunk/test/NightlyRun/test425.m

    r5115 r5193  
    66
    77%Dakota options
    8 md.variables.nuv=normal_uncertain.empty();
    9 md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01);
    10 md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
    118
    12 md.responses=struct();
    13 md.responses.rf =response_function.empty();
    14 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]);
     9%partitioning
     10md.npart=20;
     11md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
     12md.part=md.part-1;
    1513
     14%variables
     15md.variables.rho_ice=normal_uncertain('RhoIce',md.rho_ice,0.01);
     16md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01);
     17
     18%responses
     19md.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     20
     21%method
    1622md.qmu_method     =dakota_method('nond_l');
    1723
     24%parameters
    1825md.qmu_params.direct=true;
    1926md.qmu_params.analysis_driver='diagnostic';
     
    2229md.qmu_params.tabular_graphics_data=true;
    2330
    24 %partitioning
    25 md.npart=20;
    26 md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
    27 md.part=md.part-1;
     31
     32%imperative!
     33md.eps_rel=10^-5; %tighten for qmu analyses
    2834md.qmu_analysis=1;
    2935
    30 md.eps_rel=10^-5; %tighten for qmu analyses
    31 
     36%solve
    3237md=solve(md,'analysis_type',DiagnosticSolutionEnum,'overwrite','y');
    3338
  • issm/trunk/test/NightlyRun/test426.m

    r5115 r5193  
    66
    77%Dakota options
    8 md.variables.nuv=normal_uncertain.empty();
    9 md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01);
    10 md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
    118
    12 md.responses=struct();
    13 md.responses.rf =response_function.empty();
    14 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]);
     9%partitioning
     10md.npart=20;
     11md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
     12md.part=md.part-1;
    1513
     14%variables
     15md.variables.rho_ice=normal_uncertain('RhoIce',md.rho_ice,0.01);
     16md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01);
     17
     18%responses
     19md.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     20
     21%method
    1622md.qmu_method     =dakota_method('nond_l');
    1723
     24%parameters
    1825md.qmu_params.direct=true;
    1926md.qmu_params.analysis_driver='diagnostic';
     
    2229md.qmu_params.tabular_graphics_data=true;
    2330
    24 %partitioning
    25 md.npart=20;
    26 md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
    27 md.part=md.part-1;
     31
     32%imperative!
     33md.eps_rel=10^-5; %tighten for qmu analyses
    2834md.qmu_analysis=1;
    2935
    30 md.eps_rel=10^-5; %tighten for qmu analyses
    31 
     36%solve
    3237md=solve(md,'analysis_type',DiagnosticSolutionEnum,'overwrite','y');
    3338
  • issm/trunk/test/NightlyRun/test427.m

    r5115 r5193  
    1414
    1515%Dakota options
    16 md.variables.nuv=normal_uncertain.empty();
    17 md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
    1816
    19 md.responses=struct();
    20 md.responses.rf =response_function.empty();
    21 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]);
    22 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]);
    23 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]);
    24 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]);
    25 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]);
    26 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]);
    27 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]);
     17%partitioning
     18md.npart=20;
     19md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
     20md.part=md.part-1;
     21
     22%variables
     23md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01);
     24
     25%responses
     26md.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     27md.responses.MassFlux1=response_function('MassFlux1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     28md.responses.MassFlux2=response_function('MassFlux2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     29md.responses.MassFlux3=response_function('MassFlux3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     30md.responses.MassFlux4=response_function('MassFlux4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     31md.responses.MassFlux5=response_function('MassFlux5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     32md.responses.massFlux6=response_function('MassFlux6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    2833
    2934%mass flux profiles
     
    3136md.qmu_mass_flux_profile_directory=pwd;
    3237
     38%method
    3339md.qmu_method     =dakota_method('nond_l');
    3440
     41%parameters
    3542md.qmu_params.direct=true;
    3643md.qmu_params.analysis_driver='diagnostic';
     
    3845md.qmu_params.interval_type='forward';
    3946md.qmu_params.tabular_graphics_data=true;
    40 
    41 %partitioning
    42 md.npart=20;
    43 md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
    44 md.part=md.part-1;
    4547md.qmu_analysis=1;
    46 
    4748md.eps_rel=10^-5; %tighten for qmu analyses
    4849
     50%solve
    4951md=solve(md,'analysis_type',DiagnosticSolutionEnum,'overwrite','y');
     52md=tres(md,'dakota');
    5053
    5154%Fields and tolerances to track changes
    52 md=tres(md,'dakota');
     55
    5356%ok, mass flux of 3 profiles should be -3 Gt/yr -3 Gt/yr and the sum, which is -6 Gt/yr
    5457%we recover those mass fluxes through the mean of the response.
  • issm/trunk/test/NightlyRun/test428.m

    r5115 r5193  
    1414
    1515%Dakota options
    16 md.variables.nuv=normal_uncertain.empty();
    17 md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
    1816
    19 md.responses=struct();
    20 md.responses.rf =response_function.empty();
    21 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]);
    22 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]);
    23 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]);
    24 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]);
    25 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]);
    26 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]);
    27 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]);
     17%partitioning
     18md.npart=20;
     19md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
     20md.part=md.part-1;
     21
     22%variables
     23md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01);
     24
     25%responses
     26md.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     27md.responses.MassFlux1=response_function('MassFlux1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     28md.responses.MassFlux2=response_function('MassFlux2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     29md.responses.MassFlux3=response_function('MassFlux3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     30md.responses.MassFlux4=response_function('MassFlux4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     31md.responses.MassFlux5=response_function('MassFlux5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     32md.responses.massFlux6=response_function('MassFlux6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    2833
    2934%mass flux profiles
     
    3136md.qmu_mass_flux_profile_directory=pwd;
    3237
     38%method
    3339md.qmu_method     =dakota_method('nond_l');
    3440
     41%parameters
    3542md.qmu_params.direct=true;
    3643md.qmu_params.analysis_driver='diagnostic';
     
    3845md.qmu_params.interval_type='forward';
    3946md.qmu_params.tabular_graphics_data=true;
    40 
    41 %partitioning
    42 md.npart=20;
    43 md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
    44 md.part=md.part-1;
    4547md.qmu_analysis=1;
    46 
    4748md.eps_rel=10^-5; %tighten for qmu analyses
    4849
     50%solve
    4951md=solve(md,'analysis_type',DiagnosticSolutionEnum,'overwrite','y');
     52md=tres(md,'dakota');
    5053
    5154%Fields and tolerances to track changes
    52 md=tres(md,'dakota');
    5355%ok, mass flux of 3 profiles should be -3 Gt/yr -3 Gt/yr and the sum, which is -6 Gt/yr
    5456%we recover those mass fluxes through the mean of the response.
Note: See TracChangeset for help on using the changeset viewer.