Changeset 5015


Ignore:
Timestamp:
08/06/10 09:31:38 (15 years ago)
Author:
seroussi
Message:

changed some nums in tests to put CM at the end

Location:
issm/trunk/test/NightlyRun
Files:
15 added
13 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/test/NightlyRun/test420_nightly.m

    r4998 r5015  
    1 md=tres(md,'dakota');
    2 md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
    3 field_names     ={'importancefactors'};
    4 field_tolerances={1e-12};
     1field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
     2field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
     3};
    54field_values={...
    6          md.results.dakota.importancefactors,...
     5        PatchToVec(md.results.SteadystateSolution.Vx),...
     6        PatchToVec(md.results.SteadystateSolution.Vy),...
     7        PatchToVec(md.results.SteadystateSolution.Vz),...
     8        PatchToVec(md.results.SteadystateSolution.Vel),...
     9        PatchToVec(md.results.SteadystateSolution.Pressure),...
     10        PatchToVec(md.results.SteadystateSolution.Temperature),...
     11        PatchToVec(md.results.SteadystateSolution.MeltingRate),...
    712        };

  • issm/trunk/test/NightlyRun/test423.m

    r4996 r5015  
    22md=geography(md,'../Exp/SquareShelf.exp','');
    33md=parameterize(md,'../Par/SquareSheetShelf.par');
    4 md=extrude(md,3,2);
    54md=setelementstype(md,'macayeal','all');
    65md.cluster='none';
    7 md.dt=0;
    8 md=solve(md,'analysis_type',SteadystateSolutionEnum);
     6
     7%Dakota options
     8md.variables.nuv=normal_uncertain.empty();
     9md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01);
     10md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
     11
     12md.responses=struct();
     13md.responses.rf =response_function.empty();
     14md.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]);
     15
     16md.qmu_method     =dakota_method('nond_l');
     17
     18md.qmu_params.direct=true;
     19md.qmu_params.analysis_driver='diagnostic';
     20md.qmu_params.evaluation_concurrency=1;
     21md.qmu_params.interval_type='forward';
     22md.qmu_params.tabular_graphics_data=true;
     23
     24md.part=0:1:md.numberofgrids-1;
     25md.part=md.part';
     26md.npart=md.numberofgrids;
     27md.qmu_analysis=1;
     28
     29md.eps_rel=10^-10; %tighten for qmu analysese
     30
     31%rm existing results
     32if exist('./qmu')
     33        system('rm -rf qmu/');
     34end
     35
     36md=solve(md,'analysis_type',DiagnosticSolutionEnum);

  • issm/trunk/test/NightlyRun/test423_nightly.m

    r4996 r5015  
    1 field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
    2 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
    3 };
     1md=tres(md,'dakota');
     2md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
     3field_names     ={'importancefactors'};
     4field_tolerances={1e-11};
    45field_values={...
    5         PatchToVec(md.results.SteadystateSolution.Vx),...
    6         PatchToVec(md.results.SteadystateSolution.Vy),...
    7         PatchToVec(md.results.SteadystateSolution.Vz),...
    8         PatchToVec(md.results.SteadystateSolution.Vel),...
    9         PatchToVec(md.results.SteadystateSolution.Pressure),...
    10         PatchToVec(md.results.SteadystateSolution.Temperature),...
    11         PatchToVec(md.results.SteadystateSolution.MeltingRate),...
     6         md.results.dakota.importancefactors,...
    127        };

  • issm/trunk/test/NightlyRun/test424.m

    r4996 r5015  
    22md=geography(md,'../Exp/SquareShelf.exp','');
    33md=parameterize(md,'../Par/SquareSheetShelf.par');
    4 md=extrude(md,3,2);
    54md=setelementstype(md,'macayeal','all');
    65md.cluster=oshostname();
    7 md.dt=0;
    8 md=solve(md,'analysis_type',SteadystateSolutionEnum);
     6
     7%Dakota options
     8md.variables.nuv=normal_uncertain.empty();
     9md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01);
     10md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
     11
     12md.responses=struct();
     13md.responses.rf =response_function.empty();
     14md.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]);
     15
     16md.qmu_method     =dakota_method('nond_l');
     17
     18md.qmu_params.direct=true;
     19md.qmu_params.analysis_driver='diagnostic';
     20md.qmu_params.evaluation_concurrency=1;
     21md.qmu_params.interval_type='forward';
     22md.qmu_params.tabular_graphics_data=true;
     23
     24md.part=0:1:md.numberofgrids-1;
     25md.part=md.part';
     26md.npart=md.numberofgrids;
     27md.qmu_analysis=1;
     28
     29md.eps_rel=10^-10; %tighten for qmu analysese
     30
     31%rm existing results
     32if exist('./qmu')
     33        system('rm -rf qmu/');
     34end
     35
     36md=solve(md,'analysis_type',DiagnosticSolutionEnum);

  • issm/trunk/test/NightlyRun/test424_nightly.m

    r4996 r5015  
    1 field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
    2 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
    3 };
     1md=tres(md,'dakota');
     2md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
     3field_names     ={'importancefactors'};
     4field_tolerances={1e-11};
    45field_values={...
    5         PatchToVec(md.results.SteadystateSolution.Vx),...
    6         PatchToVec(md.results.SteadystateSolution.Vy),...
    7         PatchToVec(md.results.SteadystateSolution.Vz),...
    8         PatchToVec(md.results.SteadystateSolution.Vel),...
    9         PatchToVec(md.results.SteadystateSolution.Pressure),...
    10         PatchToVec(md.results.SteadystateSolution.Temperature),...
    11         PatchToVec(md.results.SteadystateSolution.MeltingRate),...
     6         md.results.dakota.importancefactors,...
    127        };

  • issm/trunk/test/NightlyRun/test425.m

    r4996 r5015  
    22md=geography(md,'../Exp/SquareShelf.exp','');
    33md=parameterize(md,'../Par/SquareSheetShelf.par');
    4 md=extrude(md,3,2);
    5 md=setelementstype(md,'pattyn','all');
     4md=setelementstype(md,'macayeal','all');
    65md.cluster='none';
    7 md.dt=0;
    8 md=solve(md,'analysis_type',SteadystateSolutionEnum);
     6
     7%Dakota options
     8md.variables.nuv=normal_uncertain.empty();
     9md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01);
     10md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
     11
     12md.responses=struct();
     13md.responses.rf =response_function.empty();
     14md.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]);
     15
     16md.qmu_method     =dakota_method('nond_l');
     17
     18md.qmu_params.direct=true;
     19md.qmu_params.analysis_driver='diagnostic';
     20md.qmu_params.evaluation_concurrency=1;
     21md.qmu_params.interval_type='forward';
     22md.qmu_params.tabular_graphics_data=true;
     23
     24%partitioning
     25md.npart=20;
     26md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
     27md.part=md.part-1;
     28md.qmu_analysis=1;
     29
     30md.eps_rel=10^-10; %tighten for qmu analyses
     31
     32%rm existing results
     33if exist('./qmu')
     34        system('rm -rf qmu/');
     35end
     36
     37md=solve(md,'analysis_type',DiagnosticSolutionEnum);

  • issm/trunk/test/NightlyRun/test425_nightly.m

    r5003 r5015  
    1 field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
    2 field_tolerances={1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-06};
     1md=tres(md,'dakota');
     2md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
     3field_names     ={'importancefactors'};
     4field_tolerances={1e-11};
    35field_values={...
    4         PatchToVec(md.results.SteadystateSolution.Vx),...
    5         PatchToVec(md.results.SteadystateSolution.Vy),...
    6         PatchToVec(md.results.SteadystateSolution.Vz),...
    7         PatchToVec(md.results.SteadystateSolution.Vel),...
    8         PatchToVec(md.results.SteadystateSolution.Pressure),...
    9         PatchToVec(md.results.SteadystateSolution.Temperature),...
    10         PatchToVec(md.results.SteadystateSolution.MeltingRate),...
     6         md.results.dakota.importancefactors,...
    117        };

  • issm/trunk/test/NightlyRun/test426.m

    r4996 r5015  
    22md=geography(md,'../Exp/SquareShelf.exp','');
    33md=parameterize(md,'../Par/SquareSheetShelf.par');
    4 md=extrude(md,3,2);
    5 md=setelementstype(md,'pattyn','all');
     4md=setelementstype(md,'macayeal','all');
    65md.cluster=oshostname();
    7 md.dt=0;
    8 md=solve(md,'analysis_type',SteadystateSolutionEnum);
     6
     7%Dakota options
     8md.variables.nuv=normal_uncertain.empty();
     9md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01);
     10md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
     11
     12md.responses=struct();
     13md.responses.rf =response_function.empty();
     14md.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]);
     15
     16md.qmu_method     =dakota_method('nond_l');
     17
     18md.qmu_params.direct=true;
     19md.qmu_params.analysis_driver='diagnostic';
     20md.qmu_params.evaluation_concurrency=1;
     21md.qmu_params.interval_type='forward';
     22md.qmu_params.tabular_graphics_data=true;
     23
     24%partitioning
     25md.npart=20;
     26md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
     27md.part=md.part-1;
     28md.qmu_analysis=1;
     29
     30md.eps_rel=10^-10; %tighten for qmu analyses
     31
     32%rm existing results
     33if exist('./qmu')
     34        system('rm -rf qmu/');
     35end
     36
     37md=solve(md,'analysis_type',DiagnosticSolutionEnum);

  • issm/trunk/test/NightlyRun/test426_nightly.m

    r5008 r5015  
    1 field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
    2 field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09
    3 };
     1md=tres(md,'dakota');
     2md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
     3field_names     ={'importancefactors'};
     4field_tolerances={1e-10};
    45field_values={...
    5         PatchToVec(md.results.SteadystateSolution.Vx),...
    6         PatchToVec(md.results.SteadystateSolution.Vy),...
    7         PatchToVec(md.results.SteadystateSolution.Vz),...
    8         PatchToVec(md.results.SteadystateSolution.Vel),...
    9         PatchToVec(md.results.SteadystateSolution.Pressure),...
    10         PatchToVec(md.results.SteadystateSolution.Temperature),...
    11         PatchToVec(md.results.SteadystateSolution.MeltingRate),...
     6         md.results.dakota.importancefactors,...
    127        };

  • issm/trunk/test/NightlyRun/test427.m

    r4998 r5015  
    55md.cluster='none';
    66
     7%constrain all velocities to 1 m/yr, in the y-direction
     8md.spcvelocity(:,1:3)=1;
     9md.spcvelocity(:,4)=0;
     10md.spcvelocity(:,5)=1*md.yts;
     11md.spcvelocity(:,6)=0;
     12
    713%Dakota options
    814md.variables.nuv=normal_uncertain.empty();
    9 md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01);
    1015md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
    1116
    … …  
    1318md.responses.rf =response_function.empty();
    1419md.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]);
     20md.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]);
     21md.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]);
     22md.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]);
     23md.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]);
     24md.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]);
     25md.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]);
     26
     27%mass flux profiles
     28md.qmu_mass_flux_profiles={'../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp'};
     29md.qmu_mass_flux_profile_directory=pwd;
    1530
    1631md.qmu_method     =dakota_method('nond_l');
    … …  
    2237md.qmu_params.tabular_graphics_data=true;
    2338
    24 md.part=0:1:md.numberofgrids-1;
    25 md.part=md.part';
    26 md.npart=md.numberofgrids;
     39%partitioning
     40md.npart=20;
     41md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
     42md.part=md.part-1;
    2743md.qmu_analysis=1;
    2844
    29 md.eps_rel=10^-10; %tighten for qmu analysese
     45md.eps_rel=10^-10; %tighten for qmu analyses
    3046
    3147%rm existing results

  • issm/trunk/test/NightlyRun/test427_nightly.m

    r4998 r5015  
    11md=tres(md,'dakota');
    2 md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
     2%ok, mass flux of 3 profiles should be -.3, -.3 and the sum, which is -.6
     3%we recover those mass fluxes through the mean of the response.
     4%also, we recover the max velo, which should be 1m/s.
     5%we put all that data in the importancefactors, which we will use to test for success.
     6%also, check that the stddev are 0.
     7md.results.dakota.importancefactors=[];
     8for i=1:7,
     9        md.results.dakota.importancefactors=[md.results.dakota.importancefactors md.results.dakota.dresp_out(i).mean];
     10end
     11for i=1:7,
     12        md.results.dakota.importancefactors=[md.results.dakota.importancefactors md.results.dakota.dresp_out(i).stddev];
     13end
    314field_names     ={'importancefactors'};
    415field_tolerances={1e-11};

  • issm/trunk/test/NightlyRun/test428.m

    r4998 r5015  
    55md.cluster=oshostname();
    66
     7%constrain all velocities to 1 m/yr, in the y-direction
     8md.spcvelocity(:,1:3)=1;
     9md.spcvelocity(:,4)=0;
     10md.spcvelocity(:,5)=1*md.yts;
     11md.spcvelocity(:,6)=0;
     12
    713%Dakota options
    814md.variables.nuv=normal_uncertain.empty();
    9 md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01);
    1015md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
    1116
    … …  
    1318md.responses.rf =response_function.empty();
    1419md.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]);
     20md.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]);
     21md.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]);
     22md.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]);
     23md.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]);
     24md.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]);
     25md.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]);
     26
     27%mass flux profiles
     28md.qmu_mass_flux_profiles={'../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp'};
     29md.qmu_mass_flux_profile_directory=pwd;
    1530
    1631md.qmu_method     =dakota_method('nond_l');
    … …  
    2237md.qmu_params.tabular_graphics_data=true;
    2338
    24 md.part=0:1:md.numberofgrids-1;
    25 md.part=md.part';
    26 md.npart=md.numberofgrids;
     39%partitioning
     40md.npart=20;
     41md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
     42md.part=md.part-1;
    2743md.qmu_analysis=1;
    2844
    29 md.eps_rel=10^-10; %tighten for qmu analysese
     45md.eps_rel=10^-10; %tighten for qmu analyses
    3046
    3147%rm existing results

  • issm/trunk/test/NightlyRun/test428_nightly.m

    r4998 r5015  
    11md=tres(md,'dakota');
    2 md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
     2%ok, mass flux of 3 profiles should be -.3, -.3 and the sum, which is -.6
     3%we recover those mass fluxes through the mean of the response.
     4%also, we recover the max velo, which should be 1m/s.
     5%we put all that data in the importancefactors, which we will use to test for success.
     6%also, check that the stddev are 0.
     7md.results.dakota.importancefactors=[];
     8for i=1:7,
     9        md.results.dakota.importancefactors=[md.results.dakota.importancefactors md.results.dakota.dresp_out(i).mean];
     10end
     11for i=1:7,
     12        md.results.dakota.importancefactors=[md.results.dakota.importancefactors md.results.dakota.dresp_out(i).stddev];
     13end
    314field_names     ={'importancefactors'};
    415field_tolerances={1e-11};
    … …  
    617         md.results.dakota.importancefactors,...
    718        };
     19md.results.dakota.importancefactors
Note: See TracChangeset for help on using the changeset viewer.